| 查看: 297 | 回复: 5 | |||
| 当前主题已经存档。 | |||
[交流]
【求助】请大家帮我分析一下下面这个程序哪里出错了
|
|||
|
请高手给我看看这个程序出错在哪里？如果改正，谢谢！ v=0.45:0.0001:3; G=230000; g=207000; k=[66655 83318 111091]; f=0.08; for x=1:1:length(k) fkp(x)=(0.999*k(x)*g+0.1094*f*k(x)*g)./(0.1094*k(x)+g*v)-G*f; end delta_p=[8000000 15000000 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.91)))); 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*0.8; end end [Av,H1,H2] = plotyy(v,fkp,v,n1,'plot'); set(H1,'LineStyle','--') set(H2,'LineStyle',':') title('潜孔钻机牵引力特性曲线') xlabel('潜孔钻机行走速度(km/h)') set(get(Av(1),'Ylabel'),'String','牵引力(N)') set(get(Av(2),'Ylabel'),'String','行走系统总效率') legend('8MPa,15MPa') grid on |
回复此楼» 猜你喜欢
过年走亲戚时感受到了所开私家车的鄙视链
已经有10人回复
-
今年春晚有几个节目很不错，点赞！
已经有12人回复
-
情人节自我反思：在爱情中有过遗憾吗？
已经有13人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有12人回复
-
基金正文30页指的是报告正文还是整个申请书
已经有5人回复
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
2楼2009-04-12 14:37:06
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
★ ★ ★ ★
kuhailangyu(金币+2,VIP+0):替楼主感谢！ 4-12 15:15
xwz118(金币+2,VIP+0):我想把这些曲线都在同个图里一起显示 4-12 15:48
kuhailangyu(金币+2,VIP+0):替楼主感谢！ 4-12 15:15
xwz118(金币+2,VIP+0):我想把这些曲线都在同个图里一起显示 4-12 15:48
|
正好我在弄类似的东东，帮你改了改 v=0.45:0.0001:3; G=230000; g=207000; k=[66655 83318 111091]; f=0.08; j=length(v) for x=1:1:length(k) for y=1:j fkp(x,y)=(0.999*k(x)*g+0.1094*f*k(x)*g)./(0.1094*k(x)+g*v(j))-G*f; end end delta_p=[8000000 15000000 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.91)))); 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*0.8; end end [Av,H1,H2] = plotyy(v,fkp(1,  ,v,n1,'plot');set(H1,'LineStyle','--') set(H2,'LineStyle',':') title('潜孔钻机牵引力特性曲线') xlabel('潜孔钻机行走速度(km/h)') set(get(Av(1),'Ylabel'),'String','牵引力(N)') set(get(Av(2),'Ylabel'),'String','行走系统总效率') legend('8MPa,15MPa') grid on 注意这里啊[Av,H1,H2] = plotyy(v,fkp(1, ,v,n1,'plot');fkp现在是3行x列的数组了 所以打印的时候，要分开打印，由于不知道你干什么，只好简单这样写了，你看看，有啥问题，大家再探讨！ |
3楼2009-04-12 14:41:50
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
4楼2009-04-12 15:53:23
5楼2009-04-12 16:02:24
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
6楼2009-04-12 20:16:07