24小时热门版块排行榜

返回列表

当前主题已经存档。

xwz118

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 787.6
帖子: 358
在线: 2.9小时
虫号: 457256
注册: 2007-11-12
性别: GG
专业: 机构学与机器人

[交流] 【求助】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

回复此楼

» 猜你喜欢

1楼 2009-05-04 08:22:37

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

fspdlh

金虫 (正式写手)

应助: 0 (幼儿园)
金币: 345.2
红花: 3
帖子: 522
在线: 4小时
虫号: 578902
注册: 2008-07-04

★ ★ ★ ★
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 ]