原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
MATLAB零维内弹道
实验三 固体火箭发动机零维内弹道计算M文件：function dy=neidandao(t,y);dy=zeros(4,1);rou=y(1);p=y(2);Vc=y(3);e=y(4);d0=0.016;h0=0.08;D0=0.03;rougr=1750;k=1.17;R=300;Tp=3200;b=0.002411;n=0.315;At=pi*16*10^(-6);fai=0.95;ka=0.98;r=b*(p/1.013/10^5)^n;gama=(2/(k+1))^((k+1)/(2*(k-1)))*sqrt(k);c=sqrt(R*Tp)/gama;if e=(D0-d0)/2 Ab=pi*(d0+2*e)*h0;else Ab=0;enddy(1)=(1/Vc)*((rougr-rou)*Ab*r-(fai*p*At)/(c*sqrt(ka)));dy(2)=(1/Vc)*(rougr*Ab*r*k*R*ka*Tp-(fai*p*p*At*k)/(rou*c*sqrt(ka))-p*Ab*r);dy(3)=Ab*r;dy(4)=r;end主程序： [t,y]=ode45(neidandao,[0:0.00001:1.5],[1.29;101300;pi*0.008*0.008*0.08;0]); plot(t,y(:,2))P-t曲线初始段P-t曲线燃烧终了段P-t曲线实验总结这次实验是我们对所学课程即固体火箭发动机零维内弹道计算以及MATLAB软件的一次练习。通过这次实验，我了解了常微分方程组数值解法的一般过程，掌握了用MATLAB软件的具体实现方法，得到了零维内弹道压强曲线，完成了实验的要求。在MATLAB中，实现常微分方程组数值解法的是ode函数（在本实验中用的是ode45），它不需要用户自己编程，使用起来比较简单，总体来说这次实验也完成的比较顺利。
文档评论（0）