Matlab解微分方程组

kekek

用Matlab解如下图的微分方程组：我用的是dsolve语句。
8 V5 d1 R! e- f' j4 f; u( X这个微分方程组中除了I，都是常数。而I则是一个变量，随着循环会不断改变。
5 R7 Z6 @( r3 @  ~( v; s) l但是我发现dsolve语句解微分方程组，好像不能含有变量。
0 V9 d* X( o0 c2 }4 e; r. D所以想请问各路大神，我该如何使用Matlab解上图的微分方程组。
# W" d  A+ K) ^. E& A我编写的程序如下：
c(1)=0.1;
/ Z8 p9 P; p2 p! G) [7 Baa=17.26873*c(1)-3.04292;
bb=1.08304*c(1)-0.12098;
cc=0.09918*c(1)^(-0.43023);
z=0.05;
! ~7 Z. f$ S( ci=aa+bb*exp(cc*z);
8 |0 J* c  L, J- J[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');
: J  m0 O7 x) j; n. ]

cichishia

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

最好尽量避免使用i(单位虚数)这样的MATLAB已经内置了其他用法的变量名。

cichishia

来学习一下