Integrate planck function between two intervals

c1 = 27.742*10.^8;
c2 = 1.43878*10.^4;
sigma = 5.67*10.^-8;
T = 5780;
lamda_1 = 8;
lamda_2 = 13;
fun = @(lambda) c1./lambda.^5.*(exp(c2./(lambda * T))-1);
F = integral(fun,lamda_1,lamda_2)./(sigma.*T^4)

4 个评论
显示 2更早的评论隐藏 2更早的评论

Sarvjeet Singh 2020-9-27

Why F vs T graph are not ploted , is shows blank figure and how i can get all the values of F and T in the workspace as it shows only last value.

c1 = 27.742*10.^8;

c2 = 1.43878*10.^4;

sigma = 5.67*10.^-8;

T = 5780;

lamda_1 = 8;

lamda_2 = 13;

for T=-10:5:30

fun = @(lambda) c1./lambda.^5.*(exp(c2./(lambda * T+273))-1);

F = integral(fun,lamda_1,lamda_2)./(sigma.*(T+273)^4)

plot(T,F)

hold on

grid

end

xlabel('Temperature')

ylabel('Fraction in atmospheric window)')

title('Fraction of the radiation emitted in the atmospheric window as a function of temperature')

legend('T=-10','T=-5','T=0','T=5','T=10','T=15','T=20','T=25','T=30')

hold off

Ameer Hamza 2020-9-27

在 MATLAB Online 中打开

You are not specifying the marker type in plot() and also you are not saving vaues of F at each iteration

c1 = 27.742*10.^8;
c2 = 1.43878*10.^4;
sigma = 5.67*10.^-8;
lamda_1 = 8;
lamda_2 = 13;
Tv = -10:5:30;
F = zeros(size(Tv));
count = 1;
for T=Tv
fun = @(lambda) c1./lambda.^5.*(exp(c2./(lambda * T+273))-1);
F(count) = integral(fun,lamda_1,lamda_2)./(sigma.*(T+273)^4);
plot(T, F(count), '+', 'MarkerSize', 10, 'LineWidth', 2)
hold on
grid
count = count + 1;
end
xlabel('Temperature')
ylabel('Fraction in atmospheric window)')
title('Fraction of the radiation emitted in the atmospheric window as a function of temperature')
legend('T=-10','T=-5','T=0','T=5','T=10','T=15','T=20','T=25','T=30')
hold off