应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 专业学科区 » 数学 » 计算数学 » 中心差分求解二维热传导的matlab程序
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 3382  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

paipai314

金虫 (初入文坛)

[求助] 中心差分求解二维热传导的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)')
运行结果不正确，是方程不收敛么==  求大神指导，求思路
回复此楼

» 收录本帖的淘帖专辑推荐

好帖子

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-09-21 19:56:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cjltxxz

新虫 (初入文坛)
lz，想请教一下您是如何修正的？
赞一下(2人) 回复此楼
9楼2014-09-23 08:47:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

yongcailiu

金虫 (小有名气)
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;
这里面用的公式不对吧？循环体执行的应该是
T(i,j,2)=f*T(i-1,j,1)+f*T(i,j-1,1)+f*T(i,j+1,1)-4*f*T(i,j,1)+f*T(i+1,j,1);

如果结果还不正确，可能对于你给定的方程，中心差分法不满足收敛性或者稳定性了
赞一下 回复此楼
2楼2013-09-21 20:11:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paipai314

金虫 (初入文坛)
引用回帖:
2楼: Originally posted by yongcailiu at 2013-09-21 20:11:05
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;
这里面用的公式不对吧？循环体执行的应该是
T(i,j,2)=f*T(i-1,j,1)+f*T(i,j-1,1)+f*T(i,j+1 ...

也还是不对的，，，
我的思路是，中心差分的时间差有2，所以需要知道前2秒的状态才可以继续算后面的，所以第2秒的状态先用前向差分算。这个思路是不是正确的？有误的话烦请指教一下。或者有没有其他的思路？
赞一下 回复此楼
3楼2013-09-21 20:25:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paipai314

金虫 (初入文坛)
引用回帖:
2楼: Originally posted by yongcailiu at 2013-09-21 20:11:05
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;
这里面用的公式不对吧？循环体执行的应该是
T(i,j,2)=f*T(i-1,j,1)+f*T(i,j-1,1)+f*T(i,j+1 ...

如果是方程不满足收敛性，可以修正么？该如何修正？
赞一下 回复此楼
4楼2013-09-21 20:32:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 086000生物与医药298调剂求助 +6 元元青青 2026-03-31 8/400 2026-04-04 22:15 by lys0704
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 材料专硕 调剂 +13 CXN123456 2026-04-03 13/650 2026-04-04 21:59 by hemengdong
[考研] 085600调剂 +4 1amJJ 2026-04-02 4/200 2026-04-04 21:53 by hemengdong
[考研] 315求调剂 +13 小羊小羊_ 2026-04-02 14/700 2026-04-04 20:30 by 蓝云思雨
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 289-求调剂 +4 这里是_ 2026-04-03 4/200 2026-04-03 14:23 by 1753564080
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 274求调剂 +10 薛定谔的虎。 2026-04-01 10/500 2026-04-03 10:13 by tianyyysss
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 一志愿复旦材料，英一专硕，总分357调剂 +4 1050389037 2026-04-02 5/250 2026-04-02 21:40 by dongzh2009
[考研] 一志愿山东大学，085600，344 +7 魏子per 2026-04-02 8/400 2026-04-02 21:12 by 百灵童888
[考研] 一志愿北京科技材料科学与工程288分，求调剂 +14 是辰啊 2026-04-02 14/700 2026-04-02 21:10 by dongzh2009
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 290求调剂 +5 dfffsar 2026-03-29 5/250 2026-04-01 19:45 by 6781022
[考研] 322求调剂 +8 三水sss 2026-04-01 8/400 2026-04-01 10:19 by 唐沐儿
[考研] 考研材料工程351分调剂 +5 整个好的 2026-03-31 5/250 2026-04-01 09:36 by topgun2009
[考研] 复试调剂 +7 双马尾痞老板2 2026-03-31 7/350 2026-03-31 19:49 by Dyhoer
信息提示
关闭
请填处理意见
关闭