普朗克公式-matlab程序

普朗克公式matlab程序clearall;c1=3.741832;%第一辐射常数c2=14387.86;%第二辐射常数forT=1000:100:1800%设置辐射温度(K)l=0.0001:0.02:15;%设置波长范围及计算步长M=1e+4.*(c1./(l.^5)./(exp(c2./(l.*T))-1));%计算指定温度光谱辐出度plot(l,M,'-b','LineWidth',1.4)%绘制光谱辐出度曲线maxM=max(M)%找出指定温度最大光谱辐出度i=find(maxM==M);%找峰值波长点text(l(i+20),M(i+20),[num2str(T),'K'],'VerticalAlignment',...'baseline','HorizontalAlignment','left','fontsize',7)holdon%在指定位置按给定方式标记对应温度endset(gca,'XTick',[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15])%设置横坐标点xlabel('\lambda/\mum')%横坐标名称及单位ylabel('M_{b\lambda}/W\cdotcm^{-2}\cdot\mum^{-1}')%纵坐标名称及单位holdoff012345678910111213141505101520251000K1100K1200K1300K1400K1500K1600K1700K1800K/mMb/Wcm-2m-1