| 查看: 1107 | 回复: 2 | ||
liu_xiao_jin金虫 (正式写手)
|
[求助]
一次奇异的线性代数方程组怎么求解已有2人参与
|
|
在求解一个具有第二类边界条件的泊松方程时，离散后得到的线性代数方程组是一次奇异的，因为原微分方程组的解就是可以相差一个任意常数的。 简单的去掉1行1列以后，确实可以求解了，但解的误差非常大，而且这个误差跟去掉的行有关。从理论上分析如果是个齐次线性代数方程组，去掉任意的一行一列都可以，但是如果是非齐次的话就不是随便去掉1行就可以的了。 请教下，有没有方法可以直接求解奇异的线性代数方程组，我记得以前看过一本，其中说过有一种迭代法可以做到，但这次要用了，找了好几天都没有找到。 |
回复此楼» 猜你喜欢
心脉受损
已经有5人回复
-
博士读完未来一定会好吗
已经有15人回复
-
Springer期刊投稿求助
已经有4人回复
-
读博
已经有3人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
我想了解一下一些立体图形的矩阵方程表示？比如球的矩阵方程表示形式？立体椭圆球？等
已经有7人回复- 求助：含多个未知系数的线性/矩阵方程组的解法
已经有4人回复
- 求问这个方程能解吗AX+|BX|=c
已经有3人回复
- 有虫友知道哪些sci杂志能刊发基础数学代数学方向的文章的投稿，
已经有4人回复
- 解线性代数方程组
已经有10人回复
- 【课件】矩阵论（南航）
已经有268人回复
匿名
用户注销 (职业作家)
Matlab专家
- 应助: 223 (大学生)
- 金币: 2891.3
- 散金: 5302
- 红花: 59
- 沙发: 10
- 帖子: 3308
- 在线: 945小时
- 虫号: 0
- 注册: 2012-05-01
- 性别: GG
- 专业: 力学
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liu_xiao_jin: 金币+50, ★★★很有帮助 2015-01-10 15:10:49
感谢参与，应助指数 +1
liu_xiao_jin: 金币+50, ★★★很有帮助 2015-01-10 15:10:49
 本帖仅楼主可见 |
2楼2015-01-09 23:00:54
pippi6
铁杆木虫 (著名写手)
工程和科学数值计算咨询
- 数学EPI: 6
- 应助: 413 (硕士)
- 贵宾: 0.002
- 金币: 7116.5
- 散金: 15
- 红花: 63
- 帖子: 1639
- 在线: 798.9小时
- 虫号: 2469437
- 注册: 2013-05-14
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
liu_xiao_jin: 金币+50, ★★★很有帮助, 谢谢！ 2015-01-10 15:11:14
感谢参与，应助指数 +1
liu_xiao_jin: 金币+50, ★★★很有帮助, 谢谢！ 2015-01-10 15:11:14
|
首先，你讨论的应该是一个Laplacian方程，而不是Poisson方程，否则就不是齐次的了。 其实方程不是奇异的。问题是你少一个条件。使用迭代法时很好处理。解决的办法就是要每次迭代时做一个平均，然后把每点函数都减去这个平均值，保证平均电位为0。在使用SOR时这么处理最好。其他迭代应该也行。但直接法可能不行，矩阵需要处理。但格点数目大时也不可能直接法。 |
3楼2015-01-09 23:26:33