理论分析及测试结果.DOC

规范化幕法 廖绮绮计算机04级系2班 实验目的 理解和掌握 规范化幕法 解决矩阵特征值问题计算的问题，了解Matlab，加深和巩固知识。 用户手册 MATLAB 6.5下通过。 由于我是采用用文件读入数据的方法，所以用户在使用之前必须先对文件data_matrix.txt、data_u0.txt进行修改,否则，则以默认值进行运算。（默认值为课本303页的例题3的解答。） 请用户在data_matrix.txt写入线性方程组的系数矩阵A 在data_u0.txt中写入初始迭代值u0，注意，用户必须输入行向量，因为我在程序中已对u0 显式的进行转置。 注意格式应与例子文件一样。 执行（RUN）Mifa.m 首先屏幕会显示用户写入的数据，用户验证一下是否有错。 （采用多位有效数位的显示方法） A = 1.00000000000000 1.00000000000000 0.50000000000000 1.00000000000000 1.00000000000000 0.25000000000000 0.50000000000000 0.25000000000000 2.00000000000000 u0 = 1 1 1 然后用户一直按Enter键即可得出特征值和特征向量： lamda = 2.53652586041718 vector = 0.74822114869438 0.64966114427996 1.00000000000000 实验内容 用 规范化幕法 来演示求解矩阵特征的过程。用Malab软件实现。主要是因为Matlab强大的数值计算能力。 方法 源码： % Created by colent 12/7 % Copyright (C) 2007-2008 XMU 数值分析 Experiment. All rights % reserved.\nAuthors: 廖绮绮(ID. % 幕法 function [ ] = Mifa() format long % 数据显示格式 A = importdata(data_matrix.txt) % 用户输入矩阵 u0 = importdata(data_u0.txt); % 用户输入的初始迭代值u0 ,注意必须是行向量，否则出错！！！！ u0 = u0; % 对用户输入的初始迭代值u0进行转置成列向量 pause; % 等待用户输入 v0 = u0; % v0 u0 = v0/max(v0); % u0 规范化向量 v1 = A * u0; % v1 maxCount = 10000; % 最大容许迭代次数,为了计算出精确值，因此设置得很大，相信Matlab可以处理得很好 count = 0; % 显示迭代次数 for i = 0:maxCount % 最多迭代10000次 tmp = max(v1); % max(v1) u1 = v1 / tmp; % u1 v1 = A * u1; % v1 if ( abs(max(v1) - tmp )eps*abs(max(v1))) % 两次的max(v1)若很相近，则说明已经得出特征值了 break; % 退出 end count = count +1; % 显示迭代次数 v1; % 显示v1 u1; % 显示u1 tmp; % 显示max(v1) end count % 显示迭代次数 lamda = tmp % 特征值 vector = u1 % 特征向量 我按照书上的算法302实现规范化幂法，程序简单，支持文件输入，支持多位数位显示，支持高精度计算。程序完全由个人完成。 有一点必须提一下，就是u0必须是列向量，但用户输入的必须是行向量，因为我在程序中有显式的转置。 我没有设置数组来存储书上算法上的ui (i = 0,1,….) 和 vi (i = 0,1,2,…..)，而是把算出来的数据覆盖掉先前计算过的。 我设置最大迭代次数为10000,主要是为了高精度计算，相信Matl