24小时热门版块排行榜    

查看: 553  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:56 by o3EMxHUAQDtS
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
信息提示
请填处理意见