Question

matlab ode45 求方程（1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行，正确解答

我编的函数文件：function dy=sdre(x,y)
dy=[y(2);2*x/(1+x^2)*y(2)];

程序 [x,y]=ode45('sdre',[0 10],[1 3]);
plot(x,y)

Answer 1

1．建立M文件vdp1001a.m如下：function dy=vdp1001a(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=2*x*y(2)/(1+x^2); 2．取x0=0，xf=10，输入命令： [X,Y]=ode45('vdp1001a',[0 10],[1 3]);plot(X,Y(:,1),'-')3．结果如图

Answer 2

你写的程序只有一处错误
ode45中要取函数句柄，
[x,y]=ode(@sdre,[0 10],[1 3]);
我已仿真成功。有问题再问我

Answer 3

．建立M文件vdp1001a.m如下：
function dy=vdp1001a(x,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=2*x*y(2)/(1+x^2);

2．取x0=0，xf=10，输入命令：

[X,Y]=ode45('vdp1001a',[0 10],[1 3]);
plot(X,Y(:,1),'-'