MATLAB实验四.docVIP

实验名称：常微分方程的应用 姓名： 学号：105012009 实验日期：2011年5月10日 实验目的： 1、复习常微分方程的相关知识点 2、学习MATLAB有关求解常微分方程（组）的指令 3、学会建模实验研究产品销售量等实际问题的微分方程建模及其求解方法 实验项目： 求解微分方程组，并绘制Appolo卫星轨迹图 2、广告效应 实验背景： 在商业活动中，常常会需要制定销售策略来预测公司产品的销售量。这就需要建立数学模型来预测销量。 在军事方面,常常会要求改进现有的舰对舰导弹系统。这也需要建立数学模型来改进导弹系统。 实验具体过程： 1、题目：已知Appol卫星的运动轨迹满足下面的方程 其中 试在初值下求解，并绘制Appolo卫星轨迹。 理论推导或编程说明： 首先引入辅助变量 化为一阶方程组 程序（或命令）： 先写M函数eg6_4.m function f=eg6_4fun(t,x） f=[x(2);2*x(3)+x(1)-((1-1/82.45)*(x(1)+1/82.45))/(sqrt((x(1)+1/82.45)^2+x(3)^2))^3-(1/82.45*(x(1)-1+1/82.45))/(sqrt((x(1)+1-1/82.45)^2+x(3)^2))^3; x(4);-2*x(2)+x(3)-((1-1/82.45)*x(3))/(sqrt((x(1)+1/82.45)^2+x(3)^2))^3-(1/82.45*x(3))/(sqrt((x(1)+1-1/82.45)^2+x(3)^2))^3]; 再在指令窗口执行 x0=[1.2;0;0;-1; [t,x]=ode45(@ eg6_4fun,[0 24],x0); ??? Error: File: D:\matlab\work\eg6_4fun.m Line: 1 Column: 24 Missing variable or function. Error in == funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in == ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... plot(x(:,1),x(:,3)); 对实验题目的解答： 解： 首先引入辅助变量 化为一阶方程组 先写M函数eg6_4.m function f=eg6_4fun(t,x） f=[x(2);2*x(3)+x(1)-((1-1/82.45)*(x(1)+1/82.45))/(sqrt((x(1)+1/82.45)^2+x(3)^2))^3-(1/82.45*(x(1)-1+1/82.45))/(sqrt((x(1)+1-1/82.45)^2+x(3)^2))^3; x(4);-2*x(2)+x(3)-((1-1/82.45)*x(3))/(sqrt((x(1)+1/82.45)^2+x(3)^2))^3-(1/82.45*x(3))/(sqrt((x(1)+1-1/82.45)^2+x(3)^2))^3]; 再在指令窗口执行 x0=[1.2;0;0;-1; [t,x]=ode45(@ eg6_4fun,[0 24],x0); ??? Error: File: D:\matlab\work\eg6_4fun.m Line: 1 Column: 24 Missing variable or function. Error in == funfun\private\odearguments at 110 f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0. Error in == ode45 at 173 [neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ... plot(x(:,1),x(:,3)); 最后可得图： 改进程序或思考： 由以上程序可知其存在一定问题，故可将其进行改进。如下： clear;dfun=inline([x(2);2*x(4)+x(1)-((1-1/82.45)*(x(1)+1/82.45))/(sqrt((x(1)+1/82.45)^2+x(3)^2))^3-(1/82.45*(x(1)-1+1/82.45))/(sqrt(