| 查看: 242 | 回复: 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 ] |
回复此楼» 猜你喜欢
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有6人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
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