clear%清除变量k=1.38054e-23;%玻尔兹曼常数h=6.626e-34;%普朗克常数c=2.997925e8;%光速sigma=5.6688e-008;%斯特潘常数b=0.0029;%维恩常数t=98.6;%热力学温度向量n=length(t);%向量长...
clear %清除变量
k=1.38054e-23; %玻尔兹曼常数
h=6.626e-34; %普朗克常数
c=2.997925e8; %光速
sigma=5.6688e-008; %斯特潘常数
b=0.0029; %维恩常数
t=98.6; %热力学温度向量
n=length(t); %向量长度
lambda=[0:0.01:35]*1e-5; %波长向量
lambda(1)=eps; %给零加一小量使分母不为零
[T,L]=meshgrid(t,lambda); %波长和温度矩阵
M=2*pi*h*c^2./(exp(h*c./(k*T.*L))-1)./L.^5;%单色辐射能力
figure %创建图形窗口
plot(lambda*1e5,M) %画曲线
hl=legend([repmat('\itT\rm=',n,1),num2str(t'),repmat('K',n,1)]);%标记图例
fs=16; %字体大小
set(hl,'fontsize',fs) %设置图例大小
grid on %加网格
title(' ','fontsize',fs)%标题
xlabel('波长\it\lambda\rm/\