Scribd
Upload
7 views4 pages

Aryan

The document outlines a MATLAB script designed to model blackbody radiation and verify Wien's Displacement Law. It calculates the energy density of blackbody radiation across a range of wavelengths for temperatures between 3500K and 5000K, plotting the results and checking the accuracy against Wien's constant. The script provides visual outputs and indicates whether the verification is successful based on a specified error threshold.

Uploaded by

Lokesh Sharma
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
7 views4 pages

Aryan

The document outlines a MATLAB script designed to model blackbody radiation and verify Wien's Displacement Law. It calculates the energy density of blackbody radiation across a range of wavelengths for temperatures between 3500K and 5000K, plotting the results and checking the accuracy against Wien's constant. The script provides visual outputs and indicates whether the verification is successful based on a specified error threshold.

Uploaded by

Lokesh Sharma
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 4

Aim-:Write a Matlab for a very famous Blackbody

Radiation and verify Wien,s Displacement rule.


Theory-:

1.Blackbody Radiation Editor Snapshot-:


clc
clear all
format long
h = 6.626 * 10^-34;
c = 3 * 10^8;
k = 1.38 * 10^-23;
b = 2.898 * 10^-3;
lamb = 1e-9:50e-9:2000e-9;

for T = 3500:500:5000
E = ((8*pi*h*c)./(lamb.^5)) .* (1./(exp((h*c)./(lamb.*k.*T) -
1 )));
[mmax , i ] = max (E);
lambmax = lamb(i);
s = lambmax*T;
error = abs(s-b);
if error<=0.001
disp("verified")
else
disp("notverified")
end
plot(lamb,E,'linewidth',3);
hold on
xlabel( 'Wavelength(m)','fontsize',10)
ylabel( 'Energy Density','fontsize',10)
title ( 'Blackbody Radiation Energy Density v/s
Wavelength','fontsize',10)
plot (lambmax , mmax , 'o' , 'Markersize',8)
text(lambmax , mmax , ['T=' , num2str(T), 'K'] , 'fontsize' ,
10)
hold on
end

Graph-:
Result-:

You might also like