| 查看: 237 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】3个金币，下面这个程序出错在哪里？谢谢
|
|||
|
我想用下面的程序在一张图里画3条曲线，分开画都能运行，但合在一起时就不行了？均以V为x轴，分别画出plot（v，fkp）；plot（v，m）；plot（v，nx） 程序中我用了分部的HOLD ON ,但没有成功。 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 0 0 0]; 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(1,j); [Av,H1,H2] = plotyy(v,fkp,v,m,'plot'); hold on plot(v,nx,'-b') hold off set(H1,'LineStyle','-') set(H2,'LineStyle','--') xlabel('行走速度(km/h)') set(get(Av(1),'Ylabel'),'String','牵引力(N)') set(get(Av(2),'Ylabel'),'String','行驶驱动系统总效率') grid on |
回复此楼» 猜你喜欢
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
2楼2009-05-03 20:29:03
cenwanglai
荣誉版主 (知名作家)
老和山猥琐派九段
- 应助: 46 (小学生)
- 贵宾: 8.842
- 金币: 7440.4
- 散金: 9654
- 红花: 118
- 沙发: 23
- 帖子: 5306
- 在线: 1961.4小时
- 虫号: 537452
- 注册: 2008-04-01
- 性别: GG
- 专业: 理论和计算化学
- 管辖: 仿真模拟
3楼2009-05-03 20:42:48
4楼2009-05-04 08:29:58