| 查看: 312 | 回复: 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 ] |
回复此楼» 猜你喜欢
这样的filecode谁见过
已经有20人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有8人回复
-
filecode与中标关系的预测
已经有5人回复
-
关于代码变化问题，想知道的进来
已经有17人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
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
