| 查看: 308 | 回复: 1 | |||
| 当前主题已经存档。 | |||
ylp835铁虫 (初入文坛)
|
[交流]
【求助】matlab 求助
|
||
|
我用MATLAB编了一个小程序，运行到us（k）一行时出现错误“Error using ==> mtimes Inner matrix dimensions must agree.” 本人是新手，希望高手指点一下（急等，谢谢）。以下是我的程序： v1=0.00000112;g=9.8;H=0.4; dp=0.00056; s=2.6; cl=0.178; cd=0.45; v=0.24:0.04:0.4; d=0.07:0.02:0.11; k=1:length(v)*length(d); for i=1:length(v) for j=1:length(d) r(k)=v(i)*d(j)/v1; r2=log10(r(k)); a=0.69*r2.^2-6.4*r2+15; m=2.269*r2.^2-20.39*r2+44.42; e=2.7*d(j)/H-0.22; f=v(i)/(sqrt(g*H)); b=0.2212; ya=3.0*(-15.4*f+4.16)*exp(-0.42*b); us(k)=100*sqrt((4*dp*(2*d(j)*g*(s-1)-e*ya*(1-m)*(v(i))^2)*sin(a))/(3*d(j)*(cl*tan(a)+cd))) k=k+1; end end plot(r,us) hold on r=[ 15000 18750 19285 23571 24107 25000 29464 32142 39285]; ys=[10.71 8.92 8.31 9.95 8.59 7.47 8.82 7.71 7.6]; plot(r,ys) h=legend('calculated values','test values',2); xlabel('Reynolds number (Re)'); ylabel('seepage welocity (cm/s)'); grid on [ Last edited by nono2009 on 2010-1-28 at 17:19 ] |
回复此楼» 猜你喜欢
5个%2F，非固定段22个字母
已经有10人回复
-
小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子
已经有16人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有3人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
国社科通讯评审
已经有6人回复
-
面上的水很深
已经有21人回复
mjq240301
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 4276.6
- 红花: 2
- 帖子: 166
- 在线: 32.6小时
- 虫号: 540976
- 注册: 2008-04-07
- 专业: 工程地质学
nono2009(金币+2):谢谢应助！ 2010-01-28 19:36
ylp835(金币+10):谢谢帮助，在线吗？能否在请教一下？ 2010-01-28 19:37
ylp835(金币+10):谢谢帮助，在线吗？能否在请教一下？ 2010-01-28 19:37
|
我觉得你在处理r（k）的时候没处理好，我不知道你最终想要得到一个什么结果，我就稍微做了一下改变，不知道可不可以，你先看看吧！ v1=0.00000112;g=9.8;H=0.4; dp=0.00056; s=2.6; cl=0.178; cd=0.45; v=0.24:0.04:0.4; d=0.07:0.02:0.11; r=1:length(v)*length(d); k=1; for i=1:length(v) for j=1:length(d) r(k)=v(i)*d(j)/v1; r2=log10(r(k)); a=0.69*r2.^2-6.4*r2+15; m=2.269*r2.^2-20.39*r2+44.42; e=2.7*d(j)/H-0.22; f=v(i)/(sqrt(g*H)); b=0.2212; ya=3.0*(-15.4*f+4.16)*exp(-0.42*b); us(k)=100*sqrt((4*dp*(2*d(j)*g*(s-1)-e*ya*(1-m)*(v(i))^2)*sin(a))/(3*d(j)*(cl*tan(a)+cd))); k=k+1; end end |
2楼2010-01-28 17:54:58
10