数值分析第一次大作业 .docx

数值分析第一次大作业王昕飏 一、 算法的设计方案A中多处为0，所以可将A501501以压缩存储为C5501来节省空间。由题中120，则501为a，1=b，否则相反。由于s是A中按模最小的特征值，所以可由反幂法得到lamdas。采用doolittle分解法将A分解成L、U。A的条件数cond(A)2=|A的按模最大特征值/A的按模最小特征值|=|a/lamdas|。在1中将A进行LU分解时，得到的矩阵U的对角线乘积就是A的行列式。求与每个k最接近的特征值，通过原点平移的方法，B=A-k*I。则此时与其最接近的特征值ik变成B按模最小的特征值，于是通过反幂法求解再加上k于是就得到ik。二、 源代码#include#include#includest