MATLAB ode45求解微分方​程组，最后plot图​画不出来。

syms y1(t) y2(t) y3(t) y4(t) y5(t) y6(t)

kappa=20;

eqns=[diff(y1,t)==-kappa*(3*y1-y2-y6-1/6*(4200*sin(t+pi)+4200)),...

diff(y2,t)==-kappa*(2*y2-y1-y3), diff(y3,t)==-kappa*(2*y3-y2-y4),...

diff(y4,t)==-kappa*(2*y4-y3-y5), diff(y5,t)==-kappa*(2*y5-y4-y6),...

diff(y6,t)==-kappa*(2*y6-y5-y1)];

cond=[y1(t)==0 y2(t)==0 y3(t)==0 y4(t)==0 y5(t)==0 y6(t)==0];

[y1sol(t),y2sol(t),y3sol(t),y4sol(t),y5sol(t),y6sol(t)]=dsolve(eqns)

t=0:12;

plot(t,y1sol(t))

错误原因：

plot(t,y1sol(t))

错误使用 plot

数据必须为可转换为双精度值的数值、日期时间、持续时间或数组。

Answer 1

[金沙遇到注单异常风控部门审核不给提现怎么办？【微81923899】]

仅供参考，没有看到ode45函数

syms y1(t) y2(t) y3(t) y4(t) y5(t) y6(t)
kappa=20;
eqns=[diff(y1,t)==-kappa*(3*y1-y2-y6-1/6*(4200*sin(t+pi)+4200)),...
    diff(y2,t)==-kappa*(2*y2-y1-y3), diff(y3,t)==-kappa*(2*y3-y2-y4),...
    diff(y4,t)==-kappa*(2*y4-y3-y5), diff(y5,t)==-kappa*(2*y5-y4-y6),...
    diff(y6,t)==-kappa*(2*y6-y5-y1)];
cond=[y1(0)==0 y2(0)==0 y3(0)==0 y4(0)==0 y5(0)==0 y6(0)==0];
[y1sol(t),y2sol(t),y3sol(t),y4sol(t),y5sol(t),y6sol(t)]=dsolve(eqns,cond);
t=0:12;
plot(t,y1sol(t))