Plot the Planck's Radiation curves Below
17 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I tried to plot both radiation curves attached using Planck's and wien's equation as shown below.
Can you help correct this?
%---------Black Body Radiation Using Planck's Equation
%---------Black Body Radiation Using Wien's Equation
clc;
close all;
clear all;
%-----------------------------------%
c=3*10^8; % speed of light in vaccum
h=6.625*10.^-34; % Planck constant
k=1.38*10.^-23; % Boltzmann constant
T=(300); % Temperatures in Kelvin
V= (0:5:30).*1e2; %The wavenumber(Hundred)
Lam=1./V; %The wavelength
figure(1)
%------------Wiens equation
I300= ((2*h*c*c*(V.^3))./(exp(h*c.*V)./(k*T))-1);
%------------Plancks equation
I2(:,T)=(2*h*c*c)./((Lam.^5).*(exp((h.*c)./(k*T.*Lam))-1));
%-----------------------------------------------------
plot(V,I300,'k','linewidth',1)
hold on
plot(Lam,I2(:,T),'r','linewidth',1)
xlabel('Wavenumber','fontsize',14)
ylabel('Spectral Irradiance [W m^{-2} sr^{-1} nmm^{-1}]','fontsize',14)
title('Plancks Radiance Law','fontsize',14)
fh = figure(1);
set(fh, 'color', 'white');
0 commentaires
Réponses (1)
Torsten
le 18 Sep 2022
x-axis and y-axis for the two plots have very different scales. So there must be an error in your formulae if we compare it to the .png file you attached.
%---------Black Body Radiation Using Planck's Equation
%---------Black Body Radiation Using Wien's Equation
clc;
close all;
clear all;
%-----------------------------------%
c=3*10^8; % speed of light in vaccum
h=6.625*10.^-34; % Planck constant
k=1.38*10.^-23; % Boltzmann constant
T=(300); % Temperatures in Kelvin
V= (0:5:30).*1e2; %The wavenumber(Hundred)
Lam=1./V; %The wavelength
figure(1)
%------------Wiens equation
I300= ((2*h*c*c*(V.^3))./(exp(h*c.*V)./(k*T))-1);
%------------Plancks equation
I2=(2*h*c*c)./((Lam.^5).*(exp((h.*c)./(k*T.*Lam))-1));
%-----------------------------------------------------
figure(1)
plot(V,I300,'k','linewidth',1)
xlabel('Wavenumber','fontsize',14)
ylabel('Spectral Irradiance [W m^{-2} sr^{-1} nmm^{-1}]','fontsize',14)
title('Plancks Radiance Law','fontsize',14)
figure(2)
plot(Lam,I2,'r','linewidth',1)
xlabel('Wavenumber','fontsize',14)
ylabel('Spectral Irradiance [W m^{-2} sr^{-1} nmm^{-1}]','fontsize',14)
title('Plancks Radiance Law','fontsize',14)
%fh = figure(1);
%set(fh, 'color', 'white');
2 commentaires
Voir également
Catégories
En savoir plus sur General Physics dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!