| 查看: 675 | 回复: 1 | ||
蓝谙至尊木虫 (职业作家)
|
[求助]
自己写了一个二维热传导方程的向前差分法，提示有两个错误，但不知怎么改。帮帮忙啦~
|
|
function u=uu(x,y)%边界函数 u=(x.^2)*sin(y)-(y.^2)*sin(x); end function [x,y,u]=heat(x0,xn,y0,yn,t0,tn,h1,h2,a)%二维热传导方程的向前差分 n1=(xn-x0)/h1;%x y的点的个数 n2=(tn-t0)/h2;%t 的点的个数 for i=1:n1 for j=1:n1 u1(i,j,1)=0;%第一层初始化 end end x(1)=x0;y(1)=y0; for i=1:n1 x(i+1)=x(i)+h1; y(i+1)=y(i)+h1;%求出x y的坐标 end for i=1:n1 for j=1:n2 u1(1,i,j)=feval(uu,x(1),y(i));%x=x0的边值 end end for i=1:n1 for j=1:n2 u1(n1,i,j)=feval(uu,x(n1),y(i));%x=xn的边值 end end for i=1:n1 for j=1:n2 u1(i,1,j)=feval(uu,x(i),y(1));%y=y0的边值 end end for i=1:n1 for j=1:n2 u1(i,n1,j)=feval(uu,x(i),y(xn));%y=yn的边值 end end r=a*h2/(h1^2); for k=2  n2-1)for i=2 n1-1)for j=2 n1-1)u(i,j,k+1)=u(i,j,k)+r*(u(i+1,j,k)-4*u(i,j,k)+u(i-1,j,k)+u(i,j+1,k)+u(i,j-1,k)) end end end for i=1:n1 for j=1:n1 u(i,j)=u1(i,j,n2); end end |
回复此楼» 猜你喜欢
各位基础老师们，看完中央《从严治党》的决定后，对学术界更加有信心了吧
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有7人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有10人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
蓝谙
至尊木虫 (职业作家)
- 应助: 3 (幼儿园)
- 金币: 10511.2
- 散金: 1340
- 红花: 20
- 沙发: 7
- 帖子: 3250
- 在线: 280.3小时
- 虫号: 3885761
- 注册: 2015-05-23
- 性别: MM
- 专业: 计算机应用技术
2楼2015-11-28 22:51:13
