| 查看: 464 | 回复: 3 | |||
| 当前主题已经存档。 | |||
[交流]
求助:MATLB运行时出现Error using ==> mrdivide怎么解决？
|
|||
|
程序如下： v=0.01:0.01:10; delta_p=[8000000 15000000 25000000 35000000]; Cs=0.8e-9; u=0.051; n=2200; Cv=0.2e+6; Cf=0.01; 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+(Cs*delta_p(i))/(u*(v(j)/(0.377*rd)))); nmt(i,j)=1-Cf-(Cv*(u*(v(j)/(0.377*rd))/delta_p(i))); nm(i,j)=nmv(i,j)*nmt(i,j); n(i,j)=np(i,j)*nm(i,j)*nc; end end figure(1) plot(v,n(1,  ,'-b')hold on plot(v,n(2, ,':r')plot(v,n(3, ,'-.b')plot(v,n(4, ,'--k')title('泵—马达系统的理论效率分析') xlabel('钻机行走速度(Km/h)') ylabel('总效率') legend('8Mpa','15Mpa','25Mpa','35Mpa') grid on hold off 运行时出现Error using ==> mrdivide Matrix dimensions must agree. Error in ==> DTH at 14 npv(i,j)=1-(1.07882*Cs)*(60*delta_p(i))/((u*n)*v(j)); |
回复此楼» 猜你喜欢
假如你的研究生提出不合理要求
已经有12人回复
-
实验室接单子
已经有7人回复
-
全日制（定向）博士
已经有5人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
对氯苯硼酸纯化
已经有3人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xwz118(金币+10,VIP+0):谢谢 3-10 19:07
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来交流 3-11 00:00
xwz118(金币+10,VIP+0):谢谢 3-10 19:07
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来交流 3-11 00:00
|
程序开始指定的 n=2200 和你计算的 n 是一个意思吗？ 这两个变更两名重合了！ /////////////////////////////////////////////////////////////// >> v=0.01:0.01:10; delta_p=[8000000 15000000 25000000 35000000]; Cs=0.8e-9; u=0.051; n=2200; Cv=0.2e+6; Cf=0.01; 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+(Cs*delta_p(i))/(u*(v(j)/(0.377*rd)))); nmt(i,j)=1-Cf-(Cv*(u*(v(j)/(0.377*rd))/delta_p(i))); nm(i,j)=nmv(i,j)*nmt(i,j); n1(i,j)=np(i,j)*nm(i,j)*nc; end end figure(1) plot(v,n1(1, ,'-b')hold on plot(v,n1(2, ,':r')plot(v,n1(3, ,'-.b')plot(v,n1(4, ,'--k')title('泵—马达系统的理论效率分析') xlabel('钻机行走速度(Km/h)') ylabel('总效率') legend('8Mpa','15Mpa','25Mpa','35Mpa') grid on hold off >> |
2楼2009-03-10 18:54:34
3楼2009-03-10 19:09:40
4楼2009-03-10 20:12:48