北航设计研究生数值分析编程大作业1.docVIP

........................................................................ 专业技术资料 数值分析大作业 算法设计方案 矩阵初始化 矩阵的下半带宽r=2,上半带宽s=2，设置矩阵，在矩阵C中检索矩阵A中的带内元素的方法是：。这样所需要的存储单元数大大减少，从而极大提高了运算效率。 利用幂法求出 幂法迭代格式： 当时，迭代终止。 首先对于矩阵A利用幂法迭代求出一个，然后求出矩阵B，其中(为单位矩阵），对矩阵B进行幂法迭代，求出，之后令，比较，大者为，小者为。 利用反幂法求出 反幂法迭代格式： 当时，迭代终止，。 每迭代一次都要求解一次线性方程组，求解过程为： 作分解 对于执行 求解（数组b先是存放原方程组右端向量，后来存放中间向量y) 使用反幂法，直接可以求得矩阵按模最小的特征值。 求与数最接近的特征值，对矩阵实行反幂法，即可求出对应的。 求出A的条件数和行列式 根据，其中分子分母分别对应按模最大和最小的特征值。 的计算：由于,其中为下三角矩阵，且对角线元素为1，故，所以有，又为上三角矩阵，故为对其对角线上各元素的乘积，最后可得。 程序源代码 （1）定义所需要的函数： #