| 查看: 540 | 回复: 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)); |
回复此楼» 猜你喜欢
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
-
好奇怪的filecode
已经有5人回复
hitzhang
木虫 (正式写手)
- 仿真EPI: 1
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1356.7
- 散金: 969
- 红花: 8
- 帖子: 863
- 在线: 226.5小时
- 虫号: 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
100