| 查看: 645 | 回复: 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 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有31人回复
-
澳大利亚 Murdoch University 全奖博士招生（3个名额）地质化工冶金领域
已经有72人回复
-
FileCode能看出啥？
已经有29人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有15人回复
-
关于Filecode分析方法
已经有12人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有11人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有10人回复
-
阿姨
已经有4人回复
-
一作与独作在应聘高校教师时区别大吗
已经有4人回复
-
静等基金结果
已经有21人回复
蓝谙
至尊木虫 (职业作家)
- 应助: 3 (幼儿园)
- 金币: 10510.4
- 散金: 1340
- 红花: 20
- 沙发: 7
- 帖子: 3250
- 在线: 280.3小时
- 虫号: 3885761
- 注册: 2015-05-23
- 性别: MM
- 专业: 计算机应用技术
2楼2015-11-28 22:51:13
