时间:2008-05-26 11:10:19   来源:论坛整理  作者:  编辑:chinaitzhe

2 x0 M0 Y0
1- x1 2 x1 M1 Y1
………… ……… ………… 乘以 ……… =
1-xn-1 2 xn-1 Mn-1 Yn-1
2 xn Mn Yn


已知： xn Yn
求解： Mn
网友回复: 2 x0 M0 Y0
1- x1 2 x1 M1 Y1
… … … 乘以 … = …
1-xn-1 2 xn-1 Mn-1 Yn-1
2 xn Mn Yn


已知： x0...xn 和y0...Yn
求解： M0... Mn
网友回复:先构造矩阵然后是矩阵的运算
网友回复:写成表达式：

2M0 x0M1 = Y0
(1-x1)M0 2M1 x1M2 = Y1
...
(1-xn-1)Mn-1 2Mn-1 xn-1Mn = Yn-1

xnMn-1 2Mn = Yn
网友回复:直接用lapack
网友回复:2楼：
构造矩阵
矩阵的运算
能否说具体点

网友回复:4楼：
lapack是什么？
网友回复:线性代数问题最好的库lapack.有lapack (CPP的)，CLapack（C语言的），Lapack（F90的)
网友回复:我工程里面一个算法要用到稀疏矩阵，以前没有涉及过，请各位多多指教！
谢谢
网友回复:你的问题需要LUV分解自己做起来很麻烦的，首先你要定义出矩阵，然后是矩阵的运算比如加减乘除nomal转置等等都要自己实现。所以还是推荐使用Lapack.
网友回复:laolaoliu2002 ：Lapack没有使用过，可以现学现用？
网友回复:对你的问题1天就能搞定，呵呵.
网友回复:laolaoliu2002 ：但愿但愿

网友回复:laolaoliu2002 ：Lapack似乎是Linux下的工具包啊
windows下可以使用？
网友回复:请各位高手指点一下，不一定要写出代码
指点一下思路也可以。
谢谢
网友回复:lapackpp-2.5.2
CLAPACK-3.1.1
都有VC的
网友回复:貌似计算方法里的用追赶法解方程组啊
如(设x为未知数)
b[1]*x[1] c[1]x[2]=d[1]
a[2]*x[1] b[2]*x[2] c[2]*x[3]=d[2]
a[3]*x[2] b[3]*x[3] c[3]*x[4]=d[3]
............
..............
a[n-1]*x[n-2] b[n-1]*x[n-1] c[n-1]*x[n]=d[n-1]

a[n]*x[n-1] b[n]*x[n]=d[n]
把他转化成
x[1] q[1]*x[2]=p[1]
x[2] q[2]*x[3]=p[2]
.............
x[n-1] q[n-1]x[n]=p[n-1]
x[n]=p[n]
再迭代回去就可以解出X[n]了

网友回复:http://www.netlib.org/clapack/
网友回复:file CLAPACK-3.1.1-VisualStudio.zip
for CLAPACK, version 3.1.1 for Windows
, (42574025 bytes)
by Julie Langou

网友回复:谢谢大家
网友回复:X*M=Y
X=Y*（M的逆）
且Y,M都为对角矩阵，所以计算起来很简单的
这样就可已用追赶法了

而上面的q，p的计算满足
p[1]=d[1]/b[1]
q[1]=c[1]/b[1]
t[k]=b[k]-a[k]*q[k-1]
p[k]=(d[k]-a[k]*p[k-1])/t[k]
q[k]=c[k]/t[k]
k=2,3,4,........,n
回代得到
x[N]=p[N]
x[k]=p[k]-q[k]*x[k 1]
k=n-1,n-2,n-3.........2,1


貌似方法说的比较复杂哈~~
不过假如知道点线性代数的知识的话相信很快就会看懂的，呵呵
网友回复:高手哪里去乎？
网友回复:哪位高手说清楚点啊！
网友回复:悲哀