主题:急求:fortran求laplace求任意n*n矩阵的值的编程,大概过程也行啊
fortranlove1
[专家分:0] 发布于 2013-02-02 01:42:00
急求:fortran语言,用laplace拉普拉斯展开求任意n*n矩阵的值的编程,大概过程也行啊
最后更新于:2013-02-02 18:22:00
回复列表 (共2个回复)
板凳
heizaoniangao [专家分:310] 发布于 2013-02-04 11:09:00
哦,只看到标题,以为帖子内容和标题一样。。。
可以按行或者按列展开,以按行为例
公式是 sum(b_ij*(-1)^(i+j)|M_ij|
有三项 b,-1,M(余子式的行列式)
考虑递推的方法:
n*n矩阵的行列式 = sum(b_ij*(-1)^(i+j)|M_ij|
而余子式的行列式 |M_ij| 又可以用它的余子式的行列式求得,整个递推下去直到余子式是一个数字为止。
总共递推n-1次。
大致方法是这样。
现在都不用这种方法了吧,运算次数太多,都用LU分解或者数值近似什么的
我来回复