作业帮 > 综合 > 作业

matlab解微分方程组并画图问题

来源:学生作业帮 编辑:百度作业网作业帮 分类:综合作业 时间:2024/06/16 17:48:28
matlab解微分方程组并画图问题
方程组如下 diff(x,t,2)-a*diff(y,t)+b*x=0,diff(y,t,2)-a*diff(x,t)+b*y=0
a=0.05 b=0.15
刚学matlab的小白QAAAQ 完全不知道用什么命令解
跪求指导QAAAAAAQ
x0=4 y0=4 Dx0=0 Dy0=0 绘图 t 由0到100
matlab解微分方程组并画图问题
有求微分方程的解,还要有 初值才行
再问: 补充了抱歉
再答: 解析解解法s=dsolve('D2x-0.05*Dy+0.15*x=0','D2y-0.05*Dx+0.15*y=0',...
    'x(0)=4','y(0)=4','Dx(0)=0','Dy(0)=0');
x=simple(s.x);
y=simple(s.y);
subplot(211),ezplot(x,[0 100]);ylabel('x');
subplot(212),ezplot(y,[0 100]);ylabel('y');

数值解法ts=[0 100];
z0=[4 4 0 0];
a=0.05;b=0.15;
fun=@(t,z) [z(3);z(4);a*z(4)-b*z(1);a*z(3)-b*z(2)];
[t,z]=ode45(fun,ts,z0);
x=z(:,1);
y=z(:,2);
subplot(211),plot(t,x);xlabel('t');ylabel('x');
subplot(212),plot(t,y);xlabel('t');ylabel('y');
再问: 好像我给的初始条件有点问题...请教下作以x-y为坐标的图的方法就好...感谢
再答: 你给的方程,初值x,y完全对称,结果x,y是一样的
在数值算法最后plot(x,y)就可以画x-y图了