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的写法对不？
本人初学，函数调用数组这方面不明白，请指点

回复此楼

» 猜你喜欢

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

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

百里道

木虫 (初入文坛)

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

引用回帖:

10楼: Originally posted by csgt0 at 2012-10-25 09:42:05
按你的意思我理解是对应一个(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*e ...

这是改过后的子函数
function result=feq(u,rho)
% double eu,uv,feq;
I=256;J=256;
N=9;
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(i,j,1)+e(n,2).*u(i,j,2);
uv=u(i,j,1).*u(i,j,1)+u(i,j,2).*u(i,j,2);
  result=w(n)*rho*(1+3*eu+4.5*eu.*eu-1.5*uv);
      end
end
end
现在我在主函数用 f(i,j,