# ode45 solve and taking current state as input

1 Ansicht (letzte 30 Tage)
WANYI ZHANG am 24 Apr. 2019
Kommentiert: Jan am 24 Apr. 2019
Hi
Can anyone help me how to fix this problem?
my questions are blow:
this is my function to be solved:
function dx = CSmacro_derive(x,fluid_params,dQ,dV,h_in,dm_in,dm_out)
dx=zeros(3,1);
v=x(1)/x(2);
h=fluid_params.h;
P=fluid_params.P;
dx(1)=dV;
dx(2)=dm_in-dm_out;
dx(3)=(1/(x(2)*fluid_params.dudT))*...
(dQ+((h_in-h)*dm_in)+(P+fluid_params.dudv)*(v*dx(2)-dx(1)));
end
This my main script:
t=linspace(0,10,10);
dV =0.1*sin(t);
dm_in=0.5;
dm_out=0;
% Initial condition
V=0.3326;
m=1;
T=15 + Kel;
x0=[V;m;T];
% Upstream refrigerant
v_in = 0.3326;
T_in = 15 + Kel;
fluid_in = refprop(v_in,T_in,r410a);
h_in=fluid_in.h;
[~,~,dQ] = convecCool(theta_a,x(3),convec_params);
fluid_params = CSmacro_props(x,r410a);
%this is my question.
%how can I take the current x value as my input to these two functions.
for i=length(t)
dx = @(t,x) CSmacro_derive(x,fluid_params,dQ,dV(i),h_in,dm_in,dm_out);
opts=odeset('RelTol',1e-6);
[t,x] = ode45(dx,t,x0,opts);
end
Thanks!
##### 6 Kommentare4 ältere Kommentare anzeigen4 ältere Kommentare ausblenden
Torsten am 24 Apr. 2019
You changed
dx = @(t,x) CSmacro_derive(x,fluid_params,dQ,dV(i),h_in,dm_in,dm_out);
to
dx = @(t,x) CSmacro_derive(x,dV(i),h_in,dm_in,dm_out);
?
Jan am 24 Apr. 2019
I've formatted the code as code to improve the readability. Please do this by your own in future posts. Thank you.

Melden Sie sich an, um zu kommentieren.

### Kategorien

Mehr zu Programming finden Sie in Help Center und File Exchange

R2019a

### Community Treasure Hunt

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

Start Hunting!

Translated by