Matlab解微分方程组

kekek

用Matlab解如下图的微分方程组：我用的是dsolve语句。
这个微分方程组中除了I，都是常数。而I则是一个变量，随着循环会不断改变。
& v( h9 V" p8 B9 }7 W# q9 ?0 [但是我发现dsolve语句解微分方程组，好像不能含有变量。
# l# S& W! _' ^+ g; f7 [& _8 T所以想请问各路大神，我该如何使用Matlab解上图的微分方程组。
. a, s8 u1 R" V) \) e' W# K我编写的程序如下：
" T; V- [% e# u/ o0 fc(1)=0.1;
4 ^: H0 z( j4 \- haa=17.26873*c(1)-3.04292;
8 Y% c4 u( c( m: g  P# R3 Ibb=1.08304*c(1)-0.12098;
/ H; ]8 o/ u' [2 O$ \cc=0.09918*c(1)^(-0.43023);
; o5 B$ {4 H& j6 P% b: A" s" r7 Z  ^z=0.05;
; t4 [7 x4 q/ ^/ `i=aa+bb*exp(cc*z);
[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');
' G& W& u( w+ L3 l8 ~! [

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

来学习一下