| 查看: 3715 | 回复: 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)') 运行结果不正确，是方程不收敛么== 求大神指导，求思路 |
回复此楼» 收录本帖的淘帖专辑推荐
 好帖子 |
» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有82人回复
-
感觉是下周放榜了
已经有8人回复
-
filecode=后面第一个是大写字母
已经有10人回复
-
时间戳变了，能看出什么问题？
已经有4人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
10楼2014-09-24 23:36:27
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
