24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

百里道

木虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1797
散金: 5
帖子: 46
在线: 11.1小时
虫号: 1070800
注册: 2010-08-07
性别: GG
专业: 工程热物理相关交叉领域

[求助] matlab函数调用与数组的问题

for i=1:I
for j=1:J
      for n=0:N-1
         f(i,j,n)=feq(u,n,rho);
      end
end
end

function result=feq(u,n,rho)
%u(i,j,1)是坐标位置为（i,j）处的水平速度，u(i,j,2)是竖直方向的速度
w=[4/9,1/9,1/9,1/9,1/9,1/36,1/36,1/36,1/36];
e=[0,0;1,0;0,1;-1,0;0,-1;1,1;-1,1;-1,-1;1,-1];
for i=1:I-1
for j=1:J-1
      for n=1:N
      eu=e(n,1).*u(:,:,1)+e(n,2).*u(:,:,2);
      uv=u(:,:,1).*u(:,:,1)+u(:,:,2).*u(:,:,2);
%uv是水平方向与竖直方向速度的平方和
      result=w(n)*rho*(1+3*eu+4.5*eu.*eu-1.5*uv);
      end
end
end
feq这么调用对不？eu和uv的写法对不？
本人初学，函数调用数组这方面不明白，请指点

回复此楼

» 猜你喜欢

论文终于录用啦！满足毕业条件了已经有13人回复
2025年遐想已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复
求个博导看看已经有18人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2012-10-22 16:25:16

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

按你的意思我理解是对应一个(i,j)，有两个速度u(i,j,1),u(i,j,2)，对应8个方向，分别是w的8行，因此一个点应该计算得到8个feq
所以子函数写成function result=feq(u,rho)
其中一句改成feq(i,j,n)=w(n)*rho*(1+3*eu+4.5*eu.*eu-1.5*uv);
这样你在主函数里就不用再循环了。此外，其实用矩阵计算可以少用循环的，会快些。