Blackbody Radiation Spectrum from Wien's Law and Planck's law
version 1.0.0.0 (1.49 KB) by
Sathyanarayan Rao
This code calculates blackbody radiation spectrum at different temperatures
Wien's law and Plancks law for black-body radiation is compared and plotted.
Cite As
Sathyanarayan Rao (2021). Blackbody Radiation Spectrum from Wien's Law and Planck's law (https://www.mathworks.com/matlabcentral/fileexchange/48253-blackbody-radiation-spectrum-from-wien-s-law-and-planck-s-law), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (7)
MATLAB Release Compatibility
Created with
R2010a
Compatible with any release
Platform Compatibility
Windows macOS LinuxCategories
Tags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Just wanted to re-affirm that this script did exactly what I expected and Stephan Suckow's suggestions worked as well.
A couple minor suggestions:
1.) Before entering the for loop, initialize the size of I1 and I2 using the following:
I1 = zeros(length(Lam),length(T));
I2 = zeros(length(Lam),length(T));
2.) Change for loop to iteration variable and range to:
for ii = 1:length(T)
The first suggestion may help speed things up. The second one does two things 1) change the for loop iteration variable from i to ii (MATLAB uses i for imaginary numbers) and 2) allows one to change the length of T and not have to modify the for loop range.
Thank you.
Does exatly what it claims. Minor improvement suggestion:
Line 33: plot(Lam,I1(:,i),'b','linewidth',2)
Line 42: line([1.2e-5 1.4e-5], [5e8 5e8],'color','b','linewidth',2)
This way the legend matches the line color for all temperatures in newer Matlab versions.