|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
| | 毕设中遇到的一个问题,急需解决!请多帮忙!!!谢谢!!!!
以下是输入一个微分方程:
%ivpodefun.m
function dxdt=ivpodefun(t,x)
dxdt =zeros(2,1);
dxdt(1)=x(2);
dxdt(2)=-x(1)+x(1)^3-2*x(2);
作图程序如下:
%ode45_example.m
[t,x]=ode45(@ivpodefun,[0 1000],[0.0001;0.0001]);
axis([-50 50 -50 50]);
plot(x(:,1),x(:,2),'-')
我想改变初始点的值作图,即
[t,x]=ode45(@ivpodefun,[0 1000],[-0.0001;-0.0001]);
axis([-50 50 -50 50]);
plot(x(:,1),x(:,2),'-')
[t,x]=ode45(@ivpodefun,[0 1000],[-0.0001;0.0001]);
axis([-50 50 -50 50]);
plot(x(:,1),x(:,2),'-')
[t,x]=ode45(@ivpodefun,[0 1000],[0.0001;-0.0001]);
axis([-50 50 -50 50]);
plot(x(:,1),x(:,2),'-')
然后把这四幅图像拼接起来。
求教!拜谢!!!
|
|
|
|