| 查看: 636 | 回复: 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 |
回复此楼» 猜你喜欢
确定了，国自然21号放榜
已经有4人回复
-
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
蓝谙
至尊木虫 (职业作家)
- 应助: 3 (幼儿园)
- 金币: 10510.4
- 散金: 1340
- 红花: 20
- 沙发: 7
- 帖子: 3250
- 在线: 280.3小时
- 虫号: 3885761
- 注册: 2015-05-23
- 性别: MM
- 专业: 计算机应用技术
2楼2015-11-28 22:51:13
100