| 查看: 550 | 回复: 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)); |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
梦想
已经有4人回复
-
中青基了要发朋友圈吗？
已经有6人回复
-
基金不中，共勉
已经有12人回复
-
能否申诉？
已经有7人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
4楼2009-03-10 20:12:48
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
投票:
10