作业帮 > 数学 > 作业

求采用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的图像.
该方程为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)')
 
以下为运行结果: