怎么用MATLAB解微分方程组?方程如下:初值只有一个,t=0时,dvx/dt=0;dvy/dt=7.98.麻烦写出详细的程序步骤dx/dt-7.76y^2-109.36x^2=0dy/dt-7.76x^2-109.36y^2=-7.98
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/25 22:27:48
怎么用MATLAB解微分方程组?方程如下:初值只有一个,t=0时,dvx/dt=0;dvy/dt=7.98.麻烦写出详细的程序步骤dx/dt-7.76y^2-109.36x^2=0dy/dt-7.76x^2-109.36y^2=-7.98
怎么用MATLAB解微分方程组?方程如下:
初值只有一个,t=0时,dvx/dt=0;dvy/dt=7.98.麻烦写出详细的程序步骤
dx/dt-7.76y^2-109.36x^2=0
dy/dt-7.76x^2-109.36y^2=-7.98
怎么用MATLAB解微分方程组?方程如下:初值只有一个,t=0时,dvx/dt=0;dvy/dt=7.98.麻烦写出详细的程序步骤dx/dt-7.76y^2-109.36x^2=0dy/dt-7.76x^2-109.36y^2=-7.98
我直接怀疑你给的方程有问题
function func
y0=[0 7.98];
[t,y]=ode45(@fun,[0,.0011],y0);
plot(t,y)
function dy=fun(t,y)
dy=zeros(2,1);
dy(1)=7.76*y(2)^2+109.36*y(1)^2;
dy(2)=7.76*y(1)^2+109.36*y(2)^2-7.98;
ts=0:0.1:10;%自己定时间范围
y0=[0 7.98];
[t,y]=ode45('fun',ts,y0);
plot(t,y(:,1),ts,y(:,2))%画vx,vy图像
%----------以下函数文件--fun.m---
function dy=fun(t,y)
dy=zeros(2,1);
dy=[7.76*y(2)^2+109.36*y(1)^2;...
7.76*y(1)^2+109.36*y(2)^2+7.98];