| 查看: 482 | 回复: 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 |
回复此楼» 猜你喜欢
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
★ ★ ★ ★
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
30