空动作业.doc

  1. 1、本文档共4页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
空动作业

syms m1 m2 y p0 p T s r; y=1.4; m1=1:0.1:10; m2=sqrt((1+((y-1)./2).*m1.^2)./(y.*m1.^2-(y-1).*m1.^2)); p=1+2*y/(y+1).*(m1.^2-1); T=(1+2*y/(y+1).*(m1.^2-1)).*(2+(y-1).*m1.^2)./((y+1).*m1.^2); s=8.314*y/(y-1)*log((1+2*y/(y+1).*(m1.^2-1)).*(2+(y-1).*m1.^2)./((y+1).*m1.^2))-8.314*log(1+2*y/(y+1).*(m1.^2-1)); p0=exp(-s/8.314); r=(y+1).*m1.^2./(2+(y-1).*m1.^2); %赋值 plotyy(m1,[p0;m2;m2],m1,[p;T;r],plot); %绘制左边为p0,m2,右边为p,T,ρ的双坐标图 [X,H1,H2]=plotyy(m1,[p0;m2;m2],m1,[p;T;r],plot); %返回双坐标两个纵坐标的句柄 set(get(X(1),Ylabel),String,M2 and p0,2/p0,1); %设置左边纵坐标的句柄属性,改变标注 set(get(X(2),Ylabel),String,T2/T1,p2/p1 and ρ2/ρ1); %设置右边纵坐标的句柄属性,改变标注 set(X(1),ylim,[0,1],ytick,[0:0.1:1]); %设置左边纵坐标的句柄属性,设置区间范围及间隔 set(X(2),ylim,[0,20],ytick,[0:2:20]); %设置右边纵坐标的句柄属性,设置区间范围及间隔 xlabel(M1); %设置横坐标的标注 syms m1 m2 y p0 p T s r;%设值; y=1.4;%假设比热比为1.4; m1=1:0.1:10;%马赫数m1的数值以0.1为步进从1到10变化; m2=sqrt((1+((y-1)./2).*m1.^2)./(y.*m1.^2-(y-1).*m1.^2)); p=1+2*y/(y+1).*(m1.^2-1); T=(1+2*y/(y+1).*(m1.^2-1)).*(2+(y-1).*m1.^2)./((y+1).*m1.^2); s=8.314*y/(y-1)*log((1+2*y/(y+1).*(m1.^2-1)).*(2+(y-1).*m1.^2)./((y+1).*m1.^2))-8.314*log(1+2*y/(y+1).*(m1.^2-1)); p0=exp(-s/8.314); r=(y+1).*m1.^2./(2+(y-1).*m1.^2);%以上为各个变量的赋值 plotyy(m1,p0,m1,p,plot);%绘制左边为p0,右边为p的双坐标图 plotyy(m1,m2,m1,T,plot);%绘制左边为m2,右边为T的双坐标图 plotyy(m1,m2,m1,r,plot);%绘制左边为m2,右边为ρ的双坐标图 hold on; [AX,H3,H4]=plotyy(m1,p0,m1,p,plot);%返回双坐标两个纵坐标的句柄 set(get(AX(1),Ylabel),String,M2 and p0,2/p0,1);%设置左边纵坐标的句柄属性,改变标注 set(get(AX(2),Ylabel),String,T2/T1,p2/p1 and ρ2/ρ1);%设置右边纵坐标的句柄属性,改变标注 set(AX(1),ylim,[0,1],ytick,[0:0.1:1]);%设置左边纵坐标的句柄属性,设置区间范围及间隔 set(AX(2),ylim,[0,20],ytick,[0:2:20]);%设置右边纵坐标的句柄属性,设置区间范围及间隔 [BX,H1,H2]=plotyy(m1,m2,m1,T,plot);%以下属性同上 set(BX(1),ylim,[0,1],ytick,[0:0.1:1]); set(BX(2),ylim,[0,20],ytick,[0:2:20]); [CX,H1,H2]=plotyy(m1,m2,m1,r,plot); set(CX(1),ylim,[0,1],ytick,[0:0.1:1]); set(CX(2),ylim,[0,20],ytick,[0:2:20]); xlabel(M1);%设置横坐标的标注 %%该图形绘制后,会多出一条关于p0的直线,使右边纵坐标产生重叠,在figure中将其删除即可得到和书本上一致的图

文档评论(0)

xy88118 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档