24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1108  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liu_xiao_jin

金虫 (正式写手)

[求助] 一次奇异的线性代数方程组怎么求解已有2人参与

在求解一个具有第二类边界条件的泊松方程时,离散后得到的线性代数方程组是一次奇异的,因为原微分方程组的解就是可以相差一个任意常数的。
简单的去掉1行1列以后,确实可以求解了,但解的误差非常大,而且这个误差跟去掉的行有关。从理论上分析如果是个齐次线性代数方程组,去掉任意的一行一列都可以,但是如果是非齐次的话就不是随便去掉1行就可以的了。

请教下,有没有方法可以直接求解奇异的线性代数方程组,我记得以前看过一本,其中说过有一种迭代法可以做到,但这次要用了,找了好几天都没有找到。
回复此楼
得之泰然,失之淡然!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pippi6

铁杆木虫 (著名写手)

工程和科学数值计算咨询

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liu_xiao_jin: 金币+50, ★★★很有帮助, 谢谢! 2015-01-10 15:11:14
首先,你讨论的应该是一个Laplacian方程,而不是Poisson方程,否则就不是齐次的了。

其实方程不是奇异的。问题是你少一个条件。使用迭代法时很好处理。解决的办法就是要每次迭代时做一个平均,然后把每点函数都减去这个平均值,保证平均电位为0。在使用SOR时这么处理最好。其他迭代应该也行。但直接法可能不行,矩阵需要处理。但格点数目大时也不可能直接法。
3楼2015-01-09 23:26:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

匿名

用户注销 (职业作家)

Matlab专家

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
liu_xiao_jin: 金币+50, ★★★很有帮助 2015-01-10 15:10:49
本帖仅楼主可见
2楼2015-01-09 23:00:54
已阅   申请数学EPI   回复此楼   编辑   查看我的主页
信息提示
请填处理意见