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

回复此楼

» 猜你喜欢

求助:我三月中下旬出站，青基依托单位怎么办？已经有9人回复
Cas 72-43-5需要30g，定制合成，能接单的留言已经有8人回复
北京211副教授，35岁，想重新出发，去国外做博后，怎么样？已经有8人回复
磺酰氟产物，毕不了业了！已经有5人回复
论文终于录用啦！满足毕业条件了已经有25人回复
2026年机械制造与材料应用国际会议（ICMMMA 2026）已经有3人回复
自荐读博已经有3人回复
不自信的我已经有5人回复
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球已经有8人回复

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

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

百里道

木虫 (初入文坛)

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

引用回帖:

17楼: Originally posted by csgt0 at 2012-10-30 17:17:09
这是你的三维矩阵维度的问题，
result(p,rho1)算出来的是一个行向量，放到三维矩阵中应该是用A=(i,:,j)来表示的
所以你跟他加减的项应该也是这样的，为了方便一般都用列向量放第一个最方便。
另外你这一句的问题 ...

非常感谢，不知道怎么才是按列排放，我只好定义了f=zeros(N,I+1,J+1),再没有出现问题
另外，我想请教函数调用的问题
对于下式
f(I+1,j,n)=feq(I+1,j,n)+f(I,j,n)-feq(I,j,n)
（feq只是说明坐标）一个式子里出现两次要调用的函数result，按照你前面所说的
主程序里面对于第（i,j）点就用
p=[u(i,j,1),u(i,j,2)];
rho1=rho(i,j);
feq(i,j,：）=result(p,rho1);

所以我写成
   p2=[u(1,I+1,j),u(2,I+1,j)];
   rho2=rho(I+1,j);
   r2=result(p2,rho2);

   p3=[u(1,I,j),u(2,I,j)];
   rho3=rho(I,j);
   r3=result(p3,rho3);

   f(n,I+1,j)=r2(n)+f(n,I,j)-r3(n);

可以运行程序，不过要求的feq的点比较多，所以这个过程很长，有没有简单的办法调用？