MATLAB的答案

解决系统的常微分方程

55视图(30天)
迭戈Mondaray穆尼奥斯
迭戈Mondaray穆尼奥斯 2023年5月7日15:59
编辑: Torsten 2023年5月8日,今天11。49分在
你好,
我有这个系统的微分方程:
在这两个正在分化尊重adimensional时间定义为:
所有其他变量被定义为:
但是我也有一个第三微分方程是:
有什么方法解决这些方程的数值解算器一样吗?后来我想夫妻这个系统与另一个两个常微分方程,将diferrentiated ' t '和我需要解决所有在一起。
提前谢谢。
4评论
显示3年长的评论
迭戈Mondaray穆尼奥斯
迭戈Mondaray穆尼奥斯 2023年5月7日21:55
我得到我的数据从datos_TH和解决方程:
% % % %简化热工水力模型
清晰的所有
datos_TH;
信谊Tfav (t) Tcav (t) Tm (t)
cond1 = Tfav (0) = = (Tf0-Tm0) / (Tf0-Tm0);
cond2 = Tcav (0) = = (Tc0-Tm0) / (Tf0-Tm0);
气孔导度= [cond1;cond2];
ode1 = diff (Tfav) = = 2 * Bigf * (Tfav-Tcav) / Rf0 + G;
ode2 = diff (Tcav) = = 2 * K * (10 * Rci * (Tfav-Tcav) bi (1) * Tcav) / (1-Rci ^ 2);% Bigc %值药对洛杉矶parametros el莫德罗门好。洛杉矶resultados没有。
常微分方程= [ode1; ode2];
[TcavSol (t) TfavSol (t)) = dsolve(诗赋,气孔导度);
Tf = [];
Tc = [];
n = 1;
t_0 = 50;
time_ad = t_0 * kf / (Rhof * cf * rc0 ^ 2);
为我= 0:0.01:time_ad
双(TfavSol特遣部队(n) =(我))* (Tf0-Tm0) + Tm0;
双(TcavSol Tc (n) =(我))* (Tf0-Tm0) + Tm0;
n = n + 1;
m = n + 1;
结束
%慢化剂温度的颂歌
t1 = linspace (0 t_0 m - 2);
t1 = ' (t1);
Tc = ' (Tc);
Tc_p =适合(t1, Tc,“poly8”);
Tc_p =潜艇(str2sym(公式(Tc_p)), coeffnames (Tc_p) num2cell (coeffvalues (Tc_p)。');
Tc_p =潜艇(Tc_p t);
cond3 = Tm (0) = = Tm0;
ode3 = diff (Tm, t) = = (Sc *嗨(1)* (Tc_p-Tm) 2 * m * cp * (Tm-Tm0)) / (Mc * cp);
TmSol (t) = dsolve (ode3 cond3);
% %的情节功能Adimensional温度
tspan = [0 time_ad];
图
网格在
fplot (TfavSol tspan)
持有在
fplot (TcavSol tspan)
传奇({“Adimesional燃料温度”,“Adimensional熔覆温度”},“位置”,“东”)
包含(“Adimensional时间”)
ylabel (“Adimensional温度”)
标题(“B1 Adimensional温度演化”)
°C % %情节温度
t = linspace (0 t_0 m - 2);
图
情节(t, Tf)
持有在
情节(t, Tc)
传奇({“燃料温度”,包层温度的},“位置”,“东”)
包含(“时间(s)”)
ylabel (“温度°C”)
标题(“维温度演化B1”)
%的阴谋慢化剂温度°C %
tspan = [0 t_0];
图
网格在
fplot (TmSol tspan)
传奇(“慢化剂温度”,“位置”,“东”)
包含(“时间(s)”)
ylabel (“慢化剂温度°C”)
标题(“h1慢化剂温度进化”)
但我这需要相当的时间,我不知道如何´几本系统与其他方程,所以我要解决这个只使用一个数值。

登录置评。

在回答这个问题。

答案(1)

Torsten
Torsten 2023年5月8日在福音11:24
编辑:Torsten 2023年5月8日,今天11。49分在
写下你的三个方程在同一自变量(因此t或τ)。
我认为这三个方程都写在t。
然后你的电话数值
tspan = [tstart往往]
y0 = [Tf0; Tc0 T0];
dTf_dt = @ (t, y)…
dTc_dt = @ (t, y)…
dT_dt = @ (t, y)…
有趣= @ (t、y) [dTf_dt (t、y); dTc_dt (t、y); dT_dt (t, y)];
[时间Y] =数值(乐趣、tspan y0)

类别

找到更多的在常微分方程在帮助中心和文件交换

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!