24小时热门版块排行榜    

查看: 541  |  回复: 3
当前主题已经存档。

xwz118

金虫 (正式写手)

[交流] 求助: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));
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hitzhang

木虫 (正式写手)

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xwz118

金虫 (正式写手)

谢谢高手呀
3楼2009-03-10 19:09:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fspdlh

金虫 (正式写手)


sunxiao(金币+1,VIP+0):谢谢参与,欢迎常来 3-11 00:00
一看你的程序就知运行起来一定挺慢,多用点乘点除,你这程序用不着循环
4楼2009-03-10 20:12:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xwz118 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] FileCode能看出啥? +6 要乐观耀哥 2026-08-10 12/600 2026-08-10 18:37 by wmfsnow
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] filecode +7 documentary 2026-08-10 7/350 2026-08-10 18:32 by GZUZM
[基金申请] 为什么网上很多人说本周 12号出结果 +3 瞬息宇宙 2026-08-10 3/150 2026-08-10 18:27 by Tide man
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见