Question

matlab解微分方程画图？？？

function mat1()
syms y
t=0:0.1:1
y=dsolve('Dy=1-2*y+t','t')
[y]
end

command window：
t =

0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000 0.7000 0.8000 0.9000 1.0000

y =

t/2 + C3/(4*exp(2*t)) + 1/4

y =

t/2 + C3/(4*exp(2*t)) + 1/4

>> plot(t,y)
??? Undefined function or variable 'y'.

Answer 1

存在的问题：

1、解微分方程的代码写在一个函数里，函数有自己独立工作区，运行函数后，变量的值不会自动保存在命令窗口里，所以会提示y未定义（你应该在之前定义过t，否则应该提示t未定义）。

2、微分方程没提供初始条件，因而求解得到的表达式中含有自由系数C3。

3、解微分方程得到的y是t的表达式，需要将t的实际取值代入才可以得到绘图所需数据。

 

参考代码（假定初值为0）：

syms  y
t=0:0.1:1;
y=dsolve('Dy=1-2*y+t','y(0)=0')
y=subs(y);
plot(t,y)

追问

y=subs(y)是什么意思，还有syms不是定义了y么？

追答

y=subs(y)：参见上面说的第3点。你用dsolve求出的是关于t的表达式，subs作用是把当前工作区中存在的变量t的具体值代入表达式。
syms y：参见上面说的第1点。你的y是在函数中定义的，而plot是在命令窗口中，二者工作区不一样。

Answer 2

差初始条件