马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
如何把M文件发布成独立的运行程序?
我的M文件:生成一个界面,上面一个文本框用于用户输入数字,另一个按纽,用户按下按纽后,画出曲线.
clear
FigWin=figure('Position',[80 80 800 600],...
'Name','传动比',...
'NumberTitle','off' );
AxesHandle=axes('Position',[0.4 0.2 0.5 0.7],...
'Box','on');
EditBox=uicontrol(FigWin,...
'Style','edit',...
'Position',[50 200 100 30],...
'String','3.54',...
'CallBack',...
['NumStr=get(EditBox,''String'');'...
'il=dtr2num(NumStr);'
]);
push1=uicontrol(FigWin,...
'Style','pushbutton',...
'Position',[50 150 100 30],...
'String','观察一下',...
'CallBack',...
['subplot(AxesHandle);'...
'V=0.123;'...
'Pe=6.1;'... % 额定功率,kW,(7500 rpm)
'Ttq=8;'... % 最大扭矩,Nm,(6000rpm)
'nmin=1500;'... % 最低空载稳定转速(怠速),rpm
'ig=[2.643,1.778,1.316,1.045,0.875];'...% 1档至5档传动比
'i0=3; '... % 初级传动比 'il=3.6; '... % 链条传动比即末级传动比
'NumStr=get(EditBox,''String'');'...
'il=str2num(NumStr);'...
'me=117; '... % 总质量,kg
'm0=75; '... % 整备质量,kg
'K=0.4111; '... %空气阻力系数与迎风面积的乘积即:Cd*A,由实验与计算得到
'r=0.315; '... % 轮胎半径,m
'L=1.4; '... % 摩托车轴距 m
'aa=0.5; '... % 摩托车质心至前轴距离 m
'bb=0.9; '... % 摩托车质心至后轴距离 m
'nmin=2000; '... % 发动机最低转速 r/min
'nmax=9000; '... % 发动机最高转速 r/min
'nt=0.95; '... % 传动系机械效率
'f0=0.018; '... % 滚动阻力系数
'a0=0; '... % 道路负荷增加系数
'g=9.8; '... % 重力加速度 m/s^2
'gama=7; '... % 燃油重度 N/L
'fii=0.75; '... % 附着系数
'hg=0.5; '... % 摩托车质心高度 m---------------------------------------------------------------------
'dt1=0.024;'...
'dt2=0.0015;'...
'dt=1+dt1+dt2*ig.^2; '...
'nm=[8512,7980,7486,7020,6490,6000,5499,4985,4485];'...
'Ttqm=[7.58,8.7,9.13,9.45,9.44,9.73,9.37,9.14,8.84];'...
'bem=[2.71,2.55,2.38,2.43,2.26,2.11,1.97,1.72,1.52]; '... %发动机实验测试数据;
'P=polyfit(nm/1000,Ttqm,4);'... % 拟合发动机使用外特性曲线
'nn=nmin:nmax;'...
'Ttq_nh=polyval(P,nn/1000);'...
'Ft=(ig''*Ttq_nh*i0*nt*il)./r; '... % 驱动力 N
'G=g*me; '... %摩托车重力 N
'Ffi=(G*fii*(bb+hg*f0))/(L+fii*hg); '... % 摩托车附着力 N
'Ftmax=Ffi; '... % 有效最大驱动力 N
'Fwx=find(Ft>Ftmax);'...
'Ft(Fwx)=Ftmax;'...
'ua=(1./ig)''*(0.377*r*nn)/(i0*il); '... % 车速 km/h
'Fw=(K*ua.^2)/21.15; '... % 空气阻力 N
'f=f0*(1+(ua.^2/19440)); '... % 摩托车滚动阻力系数
'Ff=G*f; '... % 滚动阻力 N
'Fw_5=(K*(ua(5,:)).^2)/21.15; '... % 第5档空气阻力 N
'uaa=[ua;ua(5,:)];'...
'ftt=[Ft;Ff(5,:)+Fw_5];'... %
'Pe=Ttq_nh.*nn/9549; '... %摩托车发动机功率 kw
'Pf_Pw=(1/nt)*((Ff(5,:)+Fw(5,:)).*ua(5,:)/3600); '... %摩托车阻力功率 kw
'plot(ua(5,:),Pe,ua(5,:),Pf_Pw);'...
'xlabel(''ua ( km/h )'');'...
'ylabel('' P ( kW ) '');'...
'text(90,7.5,''5档 Pe'');'...
'text(100,9,''(Pf+Pw)/nt'');'...
'hold on;'...
'grid on;'
]); |