SIRT,LSQR,SVD算法程序中,M,N分别系数矩阵A的行数和列数(对于不同的方程,需自己改动此参数),反演得X分别放在文件SIRT_X.TXT,LSQR_X.TXT,SVD_X.TXT中(由程序自动生成)。
另:运行svd程序时,找到svd.c,打开编译运行即可。
附带一个A*X=B的测试文件,其中a.txt是测试方程矩阵的系数矩阵,x.txt为已知解,b.txt为测试方程右边的常数项。分别用上面提供的三个反演算法程序,计算一下,将a.txt ,b.txt作为已知输入程序,算出解X,再与x.txt 中的已知真实解比较,即可看出哪种算法的精度高。
1