求采用ode23解下列微分方程的程序!
来源:学生作业帮 编辑:作业帮 分类:数学作业 时间:2024/05/23 14:08:10
求采用ode23解下列微分方程的程序!
160*D2y=5120-4*y-0.05*(Dy)^2-0.2*Dy
y(0)=-160,Dy(0)=0,
其中y是关于时间t的函数,尽量画出x的图像.
160*D2y=5120-4*y-0.05*(Dy)^2-0.2*Dy
y(0)=-160,Dy(0)=0,
其中y是关于时间t的函数,尽量画出x的图像.
该方程为2阶微分方程,可以化为两个1阶方程构成的方程组来求解.
方法为:设x1=y,x2=y'
先建立myode.m文件用来描述该系统
%%% 以下为myode.m的内容
function dx=myode(t,x)
dx=zeros(2,1);
dx=[x(2);
(5120-4*x(1)-0.05*(x(2)^2-0.2*x(2)))/160;];
%%%%%%%%%%%%%%%%%%%%%%%%
%%以下为画出y关于t的图像的测试程序test.m的内容
Tspan = [0 5]; % Solve from t=1 to t=5
IC = [-160,0]; %
[T X] = ode23('myode',Tspan,IC); % Solve ODE
plot(T,X(:,1))
xlabel('t')
ylabel('y(t)')
以下为运行结果:
方法为:设x1=y,x2=y'
先建立myode.m文件用来描述该系统
%%% 以下为myode.m的内容
function dx=myode(t,x)
dx=zeros(2,1);
dx=[x(2);
(5120-4*x(1)-0.05*(x(2)^2-0.2*x(2)))/160;];
%%%%%%%%%%%%%%%%%%%%%%%%
%%以下为画出y关于t的图像的测试程序test.m的内容
Tspan = [0 5]; % Solve from t=1 to t=5
IC = [-160,0]; %
[T X] = ode23('myode',Tspan,IC); % Solve ODE
plot(T,X(:,1))
xlabel('t')
ylabel('y(t)')
以下为运行结果: