TA的每日心情 | 开心
2022-1-24 15:10 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
1 h/ ^+ J3 r5 P2 s0 i! G2 _
请问大家这个要怎么列呢，其中G、C为4*4的已知矩阵，GL、P为4*1的已知矩阵，T（t-dt）为上一时刻的4*1的温度矩阵，现在想用这个矩阵方程组求t时刻的温度矩阵T（t），其中dt差分步长设定为1e-3，
表示为温度对时间的导数。我的主要程序如下：
. k3 Q* f; T# eT1(:,:,1)=77*ones(4,1);$ d; ~8 p2 I! _: E4 X6 O. K
for i=1：40! [/ m& s& J0 g3 U% y0 U0 c
G*T1(:,:,i+1)==77*GL+P-(C*T1(:,:,i+1)-C*T1(:,:,i))/dt;
) ^3 O( S3 x( e4 ]; U O kend& J3 s, @0 K% x# M
2 X/ D: F7 _; X* ?+ u- Y% C2 A4 _
9 E6 E% V3 P& G- ~1 \# i$ w% _2 r但运行结果是这样
. `# C, Y. z; w- T5 ^
|
|
/1 
发表于 2022-10-13 15:44
收藏
淘帖
支持！
反对！