飞行器系统仿真

飞行器系统仿真

ID:39345179

大小:1.31 MB

页数:49页

时间:2019-07-01

飞行器系统仿真_第1页
飞行器系统仿真_第2页
飞行器系统仿真_第3页
飞行器系统仿真_第4页
飞行器系统仿真_第5页
资源描述:

《飞行器系统仿真》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库

1、《飞行器系统仿真与CAD》学习报告第一部分仿真(40)题目1:给定导弹相对于目标的运动学方程组为r(0)=5km,q(0)=60deg,s(0)=30deg,V=1.2Ma,Vm=0.9Ma,1Ma=340m/s,k=2(1)建立系统的方框图模型;(2)用MATLAB语言编写S—函数(3)用窗口菜单对(1),(2)进行仿真,动态显示结果;(4)用命令行对(1),(2)进行仿真,以图形显示结果答:(1)(2)用MATLAB语言编写S函数function[sys,x0,str,ts]=CAD1_sfun(t,x,u,fla

2、g)switchflagcase0[sys,x0,str,ts]=mdlInitializeSizes;case1sys=mdlDerivatives(t,x,u);case3sys=mdlOutputs(t,x,u);case{2,4,9}sys=[];otherwiseerror('unhandledflag=',num2str(flag))endfunction[sys,x0,str,ts]=mdlInitializeSizessizes=simsizes;sizes.NumContStates=3;sizes.

3、NumDiscStates=0;sizes.NumOutputs=3;sizes.NumInputs=0;sizes.DirFeedthrough=1;sizes.NumSampleTimes=1;sys=simsizes(sizes);str=[];x0=[5000,pi/3,pi/6];ts=[00];functionsys=mdlDerivatives(t,x,u)vm=0.9*340;v=1.2*340;k=2;dx(1)=vm*cos(x(2))-v*cos(x(2)-x(3));dx(2)=(v*sin(x

4、(2)-x(3))-vm*sin(x(2)))/x(1);dx(3)=k*dx(2);sys=dx;functionsys=mdlOutputs(t,x,u)sys=x;调用S函数的模型框图(3)框图仿真结果:S函数仿真结果:(4)命令输入clear;clc[tx]=sim('CAD1');hSimulink=figure();subplot(3,1,1);plot(t,x(:,1));grid;ylabel('r');subplot(3,1,2);plot(t,x(:,2));grid;ylabel('q');sub

5、plot(3,1,3);plot(t,x(:,3));grid;ylabel('sigma');[tx]=sim('CAD1_S');hSFun=figure();subplot(3,1,1);plot(t,x(:,1));grid;ylabel('r');subplot(3,1,2);plot(t,x(:,2));grid;ylabel('q');subplot(3,1,3);plot(t,x(:,3));grid;ylabel('sigma');模型仿真结果:S函数仿真结果:题目2:给出动态方程;(1)用MATLA

6、B语言编写S—函数;(2)用命令行gear/adams法对(1)进行仿真,显示曲线x(t=0:100);(3)建立方框图,用RK45仿真50秒,显示曲线答:(1)用MATLAB语言编写S—函数function[sys,x0,str,ts]=CAD2_sfun(t,x,u,flag)switchflagcase0[sys,x0,str,ts]=mdlInitializeSizes;case1sys=mdlDerivatives(t,x,u);case3sys=mdlOutputs(t,x,u);case{2,4,9}sy

7、s=[];otherwiseerror('unhandledflag=',num2str(flag))endfunction[sys,x0,str,ts]=mdlInitializeSizessizes=simsizes;sizes.NumContStates=2;sizes.NumDiscStates=0;sizes.NumOutputs=2;sizes.NumInputs=0;sizes.DirFeedthrough=1;sizes.NumSampleTimes=1;sys=simsizes(sizes);str=

8、[];x0=[1,0];ts=[00];functionsys=mdlDerivatives(t,x,u)dx(1)=x(2);dx(2)=1-t*x(1)-(1-x(1)^2)*x(2);sys=dx;functionsys=mdlOutputs(t,x,u)sys=x;(2)直接调用ode数值积分函数进行仿真,系统微分方程CA

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。