Plot the Planck's Radiation curves Below
조회 수: 22 (최근 30일)
이전 댓글 표시
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
답변 (1개)
Torsten
2022년 9월 18일
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
참고 항목
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!