clc
c=2.997925e8;
lambda1=[2:0.15:12,13:4:30];
lambda=lambda1*1e-6;
T=1000;
k=1.38054e-23;
h=6.626e-34;
c1=2*pi*h*c^2;
c2=h*c/k;
m1=c1*lambda.^(-5).*exp(-c2./(lambda*T));
figure
plot(lambda1,m1,'LineWidth',2)
xlabel("波长长度(um)")
ylabel("单色辐射度")
%维恩公示
lambda1=7:0.15:30;
lambda=lambda1*1e-6;
m2=2*pi*c*k*T./(lambda.^4);
hold on
plot(lambda1,m2,'LineWidth',2)
%瑞利-金斯公式
lambda1=[0.5:0.15:5,5.5:1:30];
lambda=lambda1*1e-6;
m3=2*pi*h*c^2./lambda.^5./exp(h*c./(lambda*k*T))-1;
hold on
plot(lambda1,m3,'LineWidth',2)
%普朗克公式
legend('维恩公式','瑞利-金斯公式','普朗克公式')
%用
%用频率表示
lambda=0.5:0.1:50;
mu=c./lambda;
%维恩公式
m4&