Matlab实现普朗克函数

最新推荐文章于 2025-04-15 11:19:37 发布

Escape the bug

最新推荐文章于 2025-04-15 11:19:37 发布

阅读量7.4k

点赞数 16

文章标签： matlab

本文链接：https://blog.csdn.net/jnbfknasf113/article/details/108974691

版权

本文分享了一段MATLAB代码，用于绘制不同温度下普朗克函数所表示的光谱辐出度随波长变化的图像。通过设置温度范围和波长计算步长，可以得到清晰直观的图像。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

做论文实验需要普朗克函数，在CSDN发现了这个算法，不错分享了。记录下来~~~~

普朗克函数+图像

在这里插入图片描述
代码：

clear all;
close all;
c1=3.741832;%第一辐射常数
c2=14387.86;%第二辐射常数
for T=270:10:370 % 设置辐射温度(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)
hold on % 在指定位置按给定方式标记对应温度
end
set(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}') % 纵坐标名称及单位
hold off