| 查看: 464 | 回复: 4 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】4个金币下面程序出错在哪里，请帮忙改正
|
|||
|
程序如下： v=0.45:0.0001:3; G=230000; g=161000; k=111091; f=0.05; Ff=11500; fkp=(0.999*k*g+0.1094*f*k*g)./(0.1094*k+g*v)-G*f; nx=((fkp+g*f-Ff)./(fkp+g*f)).*(1-((0.1094*fkp)./g)-0.001); delta_p=[8000000 15000000]; Cs=0.8e-9; Csm=2.0e-9; u=0.051; n=2200; Cv=0.2e+6; Cvm=0.6e+6; Cf=0.01; Cfm=0.015; Tc=5; Vmax=45e-6; nc=0.9; rd=0.3223; for i=1:1:length(delta_p) for j=1:1:length(v) npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j)); npt(i,j)=1/(1+(Cv*(u*n)/(60*delta_p(i)+Cf+(2*3.141592*Tc)/(delta_p(i)*Vmax))*(1.07882/v(j)))); np(i,j)=npv(i,j)*npt(i,j); nmv(i,j)=1/(1+(Csm*delta_p(i))/(u*(v(j)/(0.377*rd*0.95)))); nmt(i,j)=1-Cfm-(Cvm*(u*(v(j)/(0.377*rd*0.95))/delta_p(i))); nm(i,j)=nmv(i,j)*nmt(i,j); n1(i,j)=np(i,j)*nm(i,j)*nc; end end m=nx*n1(i,j); [Av,H1,H2] = plotyy(v,fkp,v,m,'plot'); set(H1,'LineStyle','-') set(H2,'LineStyle','--') xlabel('行走速度(km/h)') set(get(Av(1),'Ylabel'),'String','牵引力(N)') set(get(Av(2),'Ylabel'),'String','行驶驱动系统总效率') grid on |
回复此楼» 猜你喜欢
确定了，国自然21号放榜
已经有4人回复
-
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
★ ★ ★ ★
xwz118(金币+4,VIP+0):非常感谢 5-4 20:17
xwz118(金币+4,VIP+0):非常感谢 5-4 20:17
|
少用循环，用点乘，运算速度也会大幅提高，错在m=nx*n1(i,j);这句上 v=0.45:0.0001:3; G=230000; g=161000; k=111091; f=0.05; Ff=11500; fkp=(0.999*k*g+0.1094*f*k*g)./(0.1094*k+g*v)-G*f; nx=((fkp+g*f-Ff)./(fkp+g*f)).*(1-((0.1094*fkp)./g)-0.001); delta_p=[8000000 15000000]; Cs=0.8e-9; Csm=2.0e-9; u=0.051; n=2200; Cv=0.2e+6; Cvm=0.6e+6; Cf=0.01; Cfm=0.015; Tc=5; Vmax=45e-6; nc=0.9; rd=0.3223; for i=1:1:length(delta_p) npv(i,  =1-(1.07882.*Cs).*(60.*delta_p(i))./((u*n).*v);npt(i, =1./(1+(Cv.*(u.*n)/(60.*delta_p(i)+Cf+(2.*3.141592.*Tc)./(delta_p(i).*Vmax)).*(1.07882./v)));np(i, =npv(i,.*npt(i,;nmv(i, =1./(1+(Csm.*delta_p(i))./(u.*(v./(0.377.*rd.*0.95))));nmt(i, =1-Cfm-(Cvm.*(u.*(v./(0.377.*rd.*0.95))./delta_p(i)));nm(i, =nmv(i,.*nmt(i,;n1(i, =np(i,.*nm(i,.*nc;m(i, =nx.*n1(i,;end [Av,H1,H2] = plotyy(v,fkp,v,m,'plot'); set(H1,'LineStyle','-') set(H2,'LineStyle','--') xlabel('行走速度(km/h)') set(get(Av(1),'Ylabel'),'String','牵引力(N)') set(get(Av(2),'Ylabel'),'String','行驶驱动系统总效率') grid on [ Last edited by fspdlh on 2009-5-4 at 17:30 ] |
2楼2009-05-04 17:28:20
3楼2009-05-04 19:37:18
4楼2009-05-04 19:39:47
5楼2009-05-04 19:41:51
