Info

Diese Frage ist geschlossen. Öffnen Sie sie erneut, um sie zu bearbeiten oder zu beantworten.

hi guys, i am having a problem with coupled differential equation i dont know where i am going wrong. pls have a look at my code and post some suggestions..

1 Ansicht (letzte 30 Tage)
function d = fun1(t,T)
f = 0.0001; % area of sensor
rosoil = 0.001;
lamdasoil = 0.55;
alphaCon = (lamdasoil / (rosoil /2));
Ks = ( alphaCon * f);
p = 1300;
Cp = 600;
Csoil = (f * rosoil * Cp * p);
Em_ast = 0.9;
sigma = 5.670*10^-8;
Em_sensor = 0.9;
alphas = 0.2;
Ta = 328;
qs_cosfi_sd = 0;
Fs = (0.25 * 0.29);
b = ((1/((1/Em_ast)+(((1/Em_sensor)-1)*(f/Fs)))));
d(1) = (((qs_cosfi_sd)*(alphas*f))-((b)*(f*sigma*(((T(1))^4)-(Ta^4))))+ (Ks*(T(1)-T(2))))/(Csoil);
d(2) = ((Ks*(T(1)-T(2)))+(Ks*(T(2)- T(3))))/(Csoil);
d(3) = ((Ks*(T(2)-T(3)))+(Ks*(T(3)-T(4))))/(Csoil);
d(4) = (Ks*(T(3)-T(4)))/(Csoil);
d = [ d(1); d(2);d(3);d(4)];
%%%[t d]=ode45('asteroid',[0 3600],173);

Antworten (1)

Matt Kindig
Matt Kindig am 19 Nov. 2012
So what exactly is wrong with the function? What is the problem?

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by