1、普朗克公式 matlab 程序clear all;c1=3.741832;%第一辐射常数c2=14387.86;%第二辐射常数for T=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)
2、,K,VerticalAlignment,.baseline,HorizontalAlignment,left,fontsize,7)hold on % 在指定位置按给定方式标记对应温度endset(gca,XTick,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15) %设置横坐标点xlabel(lambda / mum) % 横坐标名称及单位ylabel(M_blambda / Wcdotcm-2cdotmum-1) % 纵坐标名称及单位hold off0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 1505101520251000K1100K1200K1300K1400K1500K1600K1700K1800K / mMb/ Wcm-2m-1
