| 查看: 5836 | 回复: 6 | |||
guoleistones至尊木虫 (著名写手)
|
[求助]
MATLAB中第二类边界条件的一维热传导问题。。
|
| 第二类边界条件的一维热传导问题。。。请问在matlab中除了PDETOOL这种方法外,还有别的什么方法吗?因为我的求解是在嵌套循环语句里,用pdetool好像不太适合。有对偏微分方程比较熟悉的大神吗?求指教。 |
» 猜你喜欢
自荐读博
已经有8人回复
投稿Elsevier的杂志(返修),总是在选择OA和subscription界面被踢皮球
已经有8人回复
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐,对您同样有帮助:
一维非稳态热传导, 温度随时间变化模型
已经有12人回复
偏微分方程的MATLAB解法
已经有415人回复
Fluent热传导和对流传热边界条件设置
已经有24人回复
有限差分法 第二类边界条件
已经有12人回复
关于二维热传导方程!
已经有4人回复
请教下可以用UDF编写第二类边界条件么?
已经有4人回复
跟大家请教一个一维热传导方程建立的问题(高分求助)
已经有3人回复
第二类边界条件的追赶法
已经有5人回复
热传导第二类边界条件(傅里叶定理)
已经有6人回复
求助这个热传导方程怎么解啊
已经有6人回复
MATLAB 解偏微分方程时的边界条件传递问题
已经有4人回复
热传导方程中各个变量对结果的影响哪个更大?
已经有4人回复
问个热传导问题的计算
已经有15人回复
求助matlab编程:一维热传导问题
已经有5人回复
求解,热传递方程
已经有6人回复
哪位高手能用matlab编程求解一维瞬态热传导方程,先谢谢啦
已经有7人回复
有限元第二类边界条件 matlab仿真怎么编
已经有3人回复
有限元方法解抛物方程时边界条件该怎么加?
已经有4人回复
matlab数值求解边界条件微分方程组
已经有7人回复
[交流]分享Matlab写的PML吸收边界条件。代码很精简
已经有11人回复
想用comsol解泊松方程,却不知道怎么设边界条件
已经有7人回复

feixiaolin
荣誉版主 (文坛精英)
-

专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
2楼2014-12-05 06:50:54
guoleistones
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 15203.9
- 帖子: 1150
- 在线: 159.7小时
- 虫号: 1323310
- 注册: 2011-06-15
- 性别: GG
- 专业: 爆炸与冲击动力学

3楼2014-12-05 18:39:20
feixiaolin
荣誉版主 (文坛精英)
-

专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
4楼2014-12-05 21:28:42
guoleistones
至尊木虫 (著名写手)
- 应助: 4 (幼儿园)
- 金币: 15203.9
- 帖子: 1150
- 在线: 159.7小时
- 虫号: 1323310
- 注册: 2011-06-15
- 性别: GG
- 专业: 爆炸与冲击动力学

5楼2014-12-08 06:22:46
feixiaolin
荣誉版主 (文坛精英)
-

专家经验: +518 - 应助: 942 (博后)
- 贵宾: 1.275
- 金币: 3430
- 散金: 58785
- 红花: 532
- 沙发: 11
- 帖子: 24215
- 在线: 2601.8小时
- 虫号: 2139575
- 注册: 2012-11-21
- 专业: 光学信息获取与处理
- 管辖: 数学
6楼2014-12-08 13:00:45
Math露珠
木虫 (知名作家)
- 应助: 10 (幼儿园)
- 金币: 7841.6
- 散金: 2548
- 红花: 123
- 沙发: 33
- 帖子: 8870
- 在线: 885.7小时
- 虫号: 4010226
- 注册: 2015-08-07
- 专业: 常微分方程与动力系统
|
楼主这个问题解决了吗?我也有同样的问题。 matlab写ux(0.t)=vx(0,t)=0,ux(l*pi.t)=vx(l*pi,t)=0,Neumann条件。 不是调用pde函数句柄,是用for end 循环语句写。我写的不对 %boundary value for i=1:n for j=1:p u(i+1,j)=u(i,j)+[0+0.1*cos((i-1)*h)]*h; v(i+1,j)=v(i,j)+[0+0.1*cos((i-1)*h)]*h; u(i+1+l*pi,j)=u(l*pi,j)+[0+0.1*cos((i-1)*h)]*h; v(i+1+l*pi,j)=v(l*pi,j)+[0+0.1*cos((i-1)*h)]*h; end end提示??? Attempted to access u(9.42478,1); index must be a positive integer or logical. Error in ==> u(i+1+l*pi,j)=u(l*pi,j)+[0+0.1*cos((i-1)*h)]*h; |
7楼2017-05-09 15:31:53







回复此楼