| 查看: 3691 | 回复: 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)') 运行结果不正确，是方程不收敛么== 求大神指导，求思路 |
回复此楼» 收录本帖的淘帖专辑推荐
 好帖子 |
» 猜你喜欢
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有4人回复
-
帮忙看看fileCode
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab数值求解二阶常微分方程
已经有9人回复
- 求助matlab编程：一维热传导问题
已经有5人回复
- 请教MATLAB中的LMI求解问题，急！！！
已经有7人回复
- 求高手解答Matlab求解微分方程的初始值和边界值问题
已经有5人回复
- 求解，热传递方程
已经有6人回复
- 哪位高手能用matlab编程求解一维瞬态热传导方程，先谢谢啦
已经有7人回复
- 分数阶微分方程的求解（用差分法，matlab编程求解）
已经有9人回复
- 有限差分离散二阶偏微分方程后如何编译matlab程序解方程组？
已经有3人回复
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- MATLAB程序求解
已经有12人回复
- matlab中有关微分方程求解的几个小问题
已经有4人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- Matlab求解偏微分方程组
已经有22人回复
- matlaB求解多项式反函数遇阻
已经有7人回复
- 新手求助Matlab求解系数
已经有10人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- 求高人指点 MATLAB 编程求解 洛伦兹的最大李雅普诺夫指数程序！！！！
已经有5人回复
- 【求助】matlab求解一个简单的积分（急！！！）
已经有7人回复
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已经有22人回复
- 【求助】积分微分方程matlab求解
已经有6人回复
- 【求助】用MATLAB写程序求解简单的正弦函数
已经有4人回复
- 【求助】请教虫友关于matlab差分方程的求解和画图程序
已经有4人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
- 【求助】偏微分方程的基本解
已经有5人回复
7楼2013-09-22 22:25:02
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
