matlab一阶微分方程画图问题

Question

输入：

clear;clc;clf;
close all
syms t y;
V0=0.2;
C0=0.1;
L=0.07;
F=9648485;
k=2.5*10^(-3);
i=4;
y=dsolve('(V0+k*i*t)*Dy+(k*i+i/F*L)*y-i/F=0','y(0)=0.1','t');
t=[0:0.01:50];
plot(t,y);

输出：
错误使用 plot
'Line' 中不支持非数值数据

请问这种错误应该怎么解决？

Answer 1

你好，dsolve出来的结果是关于t的一个符号表达式，不是真正的数值，你可以这样做

clear;clc;clf;
close all
syms t y(t);%定义t以及关于t的函数y
V0=0.2;
C0=0.1;
L=0.07;
F=9648485;
k=2.5*10^(-3);
i=4;
Dy = diff(y,t,1);
y=dsolve((V0+k*i*t)*Dy+(k*i+i/F*L)*y-i/F==0,y(0)==0.1,t);
tval=[0:0.01:50];%定于tval，就是t的数值
yval = eval(subs(y, t, tval));%将y表达式中的t用tval代替，再求出精确值
plot(tval,yval);%画图像

注意：不要把符号t和数值t混在一起写，我这里把数值t用tval表示，数值y用yval表示，只有数值t和数值y才能画图，符号求解的y和t是公式需要，可以用subs以及eval获取数值的y