作业帮 > 数学 > 作业

matlab中常微分方程符号解和数值解区别?

来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/08 19:32:49
matlab中常微分方程符号解和数值解区别?
这两个道题有什么区别 我只知道符号解用dsolve数值解用ode45等.图中的方程用数值的方法怎么算啊 结果和用dsolve计算的结果有什么区别?
dsolve('x^2*D2y+4*x*Dy+2*y=0','y(1)=2','Dy(1)=-3','x')
 
ans =
 
(x + 1)/x^2
数值的方法:结果是一列向量
令y1=y,y2=y1'=y',把原方程转化成微分方程组:
y1'=y2
y2'=-2y-4*y*y(1),
(x1(0),x2(0))=(0,0).
根据上述一阶微分方程组编写M函数文件myfun.m
function dy=myfun(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-2*y(1)/x^2-4*y(2)/x;
[x,y]=ode45('myfun',[1 10],[2 -3]);plot(x,y);