飞行器控制实验概述.docx

Harbin Institute of Technology 飞行器制导与控制 实验报告 专 业： 自动化 班 级： 学 号： 1120410333 姓 名： 设计时间： 2015/12/12 上机实验1： 使用四阶龙格库塔法求解微分方程  MACROBUTTON MTPlaceRef \* MERGEFORMAT  SEQ MTEqn \h \* MERGEFORMAT ( SEQ MTEqn \c \* Arabic \* MERGEFORMAT 1) 先定义参数，初值条件可以自己任取。 源程序： function [x,y] = M1(fun,x0,xt,y0,PointNum) if nargin4 | PointNum=0 PointNum=100; end if nargin3 y0=0; end y(1,:)=y0(:); h=(xt-x0)/(PointNum-1); x=x0+[0:(PointNum)]*h; for k=1:(PointNum) f1=h*feval(fun,x(k),y(k,:)); f1=f1(:); f2=h*feval(fun,x(k)+h/2,y(k,:)); f2=f2(:); f3=h*feval(fun,x(k)+h/2,y(k,:)); f3=f3(:); f4=h*feval(fun,x(k)+h,y(k,:)); f4=f4(:); y(k+1,:)=y(k,:)+(f1+2*(f2+f3)+f4)/6; end 2、运行文件： x0=0; xt=2; Num=100; h=(xt-x0)/(Num-1); x=x0+[0:Num]*h; a=1; yt=1-exp(-a*x); fun=inline(-y+1,x,y); y0=0; PointNum=100; [xr,yr]=M1(fun,x0,xt,y0,Num); M1_x=xr M1_y=yr plot(x,yt,k,xr,yr,r-) legend(jiexi,Runge-Kutta,2) 3、实验结果： 0 0 0.0202 0.0202 0.0404 0.0400 0.0606 0.0594 0.0808 0.0784 0.1010 0.0970 0.1212 0.1153 0.1414 0.1331 0.1616 0.1506 0.1818 0.1678 0.2020 0.1846 0.2222 0.2011 0.2424 0.2172 0.2626 0.2330 0.2828 0.2485 0.3030 0.2637 0.3232 0.2786 0.3434 0.2932 0.3636 0.3074 0.3838 0.3214 0.4040 0.3351 0.4242 0.3486 0.4444 0.3617 0.4646 0.3746 0.4848 0.3873 0.5051 0.3996 0.5253 0.4118 0.5455 0.4236 0.5657 0.4353 0.5859 0.4467 0.6061 0.4579 0.6263 0.4688 0.6465 0.4796 0.6667 0.4901 0.6869 0.5004 0.7071 0.5105 0.7273 0.5204 0.7475 0.5300 0.7677 0.5395 0.7879 0.5488 0.8081 0.5580 0.8283 0.5669 0.8485 0.5756 0.8687 0.5842 0.8889 0.5926 0.9091 0.6008 0.9293 0.6089 0.9495 0.6168 0.9697 0.6245 0.9899 0.6321 1.0101 0.63961.0303 0.6468 1.0505 0.6540 1.0707 0.6610 1.0909 0.6678 1.1111 0.6745 1.1