EDA365电子论坛网
标题:
Matlab解微分方程组
[打印本页]
作者:
kekek
时间:
2020-12-8 10:23
标题:
Matlab解微分方程组
用Matlab解如下图的微分方程组：
我用的是dsolve语句。
. W% J! w \4 V5 k) L
这个微分方程组中除了I，都是常数。而I则是一个变量，随着循环会不断改变。
4 w- |7 g* g- i0 Q) v4 F
但是我发现dsolve语句解微分方程组，好像不能含有变量。
1 j1 A L5 A$ u2 p
所以想请问各路大神，我该如何使用Matlab解上图的微分方程组。
/ W/ k& ~- V$ [( K/ g6 w
我编写的程序如下：
' j! K1 A/ n7 P" K
c(1)=0.1;
% K, [6 Z$ O8 X- _+ B1 p: ^
aa=17.26873*c(1)-3.04292;
/ s4 T k& w4 ]- A! D' \1 X$ q
bb=1.08304*c(1)-0.12098;
/ P; T2 J3 w. `$ `4 b' |" k
cc=0.09918*c(1)^(-0.43023);
9 _- @- A \ O; m5 y
z=0.05;
$ Z, q; R$ g0 a+ V" M4 T8 \
i=aa+bb*exp(cc*z);
# G7 H; a* C1 K) Y3 E" [) r9 }
[x y]=dsolve('Dx=-(7.253*10^(-4)*i+9.968*10^(-6))*x+(0.009485-9.968*10^(-6))*y+9.968*10^(-6),Dy=7.253*10^(-4)*i*x-(0.009485+1.703*10^(-8)*i)*y','x(0)=1,y(0)=0','t');
2 a) Z3 R- f v( t9 f$ A
作者:
cichishia
时间:
2020-12-8 11:06
Support for character vector or string inputs will be removed in a future release. Instead, use syms to declare variables and replace inputs such as dsolve('Dy = y') with syms y(t); dsolve(diff(y,t) == y).
作者:
Uifhjvv
时间:
2020-12-8 13:44
最好尽量避免使用i(单位虚数)这样的MATLAB已经内置了其他用法的变量名。
作者:
cichishia
时间:
2020-12-8 14:54
来学习一下
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2