| 查看: 3850 | 回复: 12 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
中心差分求解二维热传导的matlab程序
|
||
|
clear;close all; L=100;S=120;M=20;N=24;timeMax=50000; Currtime=100;DL=[0:L/M:L];DS=[0:S/N:S]; T=ones(M,N,timeMax)*30; T(1:M,1,1:1:timeMax)=50; T(1:M,N,1:1:timeMax)=50; T(1,1:N,1:1:timeMax)=80; T(M,1:N,1:1:timeMax)=100; f=0.1; for i=2:M-1 for j=2:N-1 T(i,j,2)=f*T(i-1,j,1)+f*T(i,j-1,1)+(1-4*f)*T(i,j,1)+f*T(i+1,j,1); end; end; for t=3:timeMax for i=2:M-1 for j=2:N-1 T(i,j,t)=f*T(i-1,j,t-1)+f*T(i,j-1,t-1)-4*f*T(i,j,t-1)+f*T(i+1,j,t-1)+f*T(i,j+1,t-1)+T(i,j,t-2); end; end; end; [x,y]=meshgrid(1:M,1:N); mesh(T(1:M,1:N,Currtime)); title(['2D Temperature Field, t=',num2str(Currtime-1),'s']); xlabel('x length(L/M)'); ylabel('y length(L/N)'); zlabel('Temperature(^oC)') 运行结果不正确，是方程不收敛么== 求大神指导，求思路 |
回复此楼» 猜你喜欢
40+中年人失业后，还能找到好工作吗？
已经有15人回复
-
二审4个多月了还是Awaiting Reviewer Scores
已经有5人回复
-
申请科研项目纯靠实力能成功吗？
已经有12人回复
-
用NBS溴化时溶剂除了CCl4外还有什么
已经有6人回复
-
极限拉扯，如何抉择
已经有15人回复
-
OpenAI又攻克数学难题了，离AI做有机合成还有多远？
已经有13人回复
-
课题组会帮你通过考核么？
已经有7人回复
-
招收博士生一名
已经有3人回复
12楼2015-11-07 22:01:44
yongcailiu
金虫 (小有名气)
- 应助: 63 (初中生)
- 金币: 2283.9
- 红花: 10
- 帖子: 181
- 在线: 250.4小时
- 虫号: 1406563
- 注册: 2011-09-18
- 性别: GG
- 专业: 计算数学与科学工程计算
2楼2013-09-21 20:11:05
3楼2013-09-21 20:25:40
4楼2013-09-21 20:32:37