| 查看: 1436 | 回复: 7 | |||
| 当前主题已经存档。 | |||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 小木虫总管 的 5 个金币 | |||
小木虫总管金虫 (正式写手)
|
[交流]
【求助】怎么用matlab编追赶法的程序
|
||
|
自己在解决二维的差分问题的时候遇到了点问题，请教版上的同学怎么编程？ 问题是一个二维的区域划分成8*22个网格，已知函数关系为 a*x(i-1,j)+b*x(i,j)+c*x(i+1,j)=x(i,j-1) i,j是网格节点，i=1到8,j=1到22 abc是已知的数 边界条件是 i=0时x(0,j)=常数 j=1到22 x(i,22)-x(i,21)=常数 i=1到8 请教大家怎么用追赶法编程？谢谢了 |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有293人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
zeoliters
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 0.902
- 金币: 2088.4
- 散金: 16
- 帖子: 445
- 在线: 60.6小时
- 虫号: 772181
- 注册: 2009-05-16
- 专业: 计算化学，量子化学
2楼2009-09-16 19:00:33
redcrag
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 1778.2
- 红花: 1
- 帖子: 1163
- 在线: 45小时
- 虫号: 254694
- 注册: 2006-05-27
- 性别: GG
- 专业: 光学
3楼2009-09-17 08:55:08
4楼2009-09-17 21:21:42
redcrag
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 1778.2
- 红花: 1
- 帖子: 1163
- 在线: 45小时
- 虫号: 254694
- 注册: 2006-05-27
- 性别: GG
- 专业: 光学
5楼2009-09-18 07:58:38
★ ★
zeoliters(金币+2,VIP+0):谢谢回帖交流！ 9-19 16:04
zeoliters(金币+2,VIP+0):谢谢回帖交流！ 9-19 16:04
|
你的右端项是已知的吗？若是，那就是Au=b，A为你的三对角矩阵，u是待求，b是已知。这样的程序在百度或是google里输入“追赶法”会找到一些，你可以试试。或者你给我你的信箱我发一个给你。若右端项也是未知（我认为是未知的），则将右端项移到左边，写成Au=b的形式，考虑到边界值则b不为零向量，此时的A也就不是三对角矩阵了，不能用追赶法，但可以直接求逆(u=inv(A)*b)，或者用一些迭代法（比如雅克比迭代、高斯-赛德尔迭代、GMRES等等）。建议你先看一些数值解的书，比如我建议的那本。 [ Last edited by 边际浪子 on 2009-9-19 at 19:28 ] |
6楼2009-09-19 19:22:39
小木虫总管
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1363.4
- 散金: 887
- 红花: 3
- 帖子: 450
- 在线: 103.1小时
- 虫号: 845704
- 注册: 2009-09-11
- 专业: 金属结构材料
7楼2009-09-21 16:53:03
8楼2009-10-30 14:16:42
