雅克比法求对称矩阵的特征值（MATLAB程序设计与应用）课设

雅克比法求对称矩阵的特征值(MATLAB程序设计与应用)课设，里面有2份不同的课设，都已近调试过了。要在MATLAB7.0上运行

%下面是m文件Jacobi1.m的代码，以经典Jacobi的的方法计算特征值

function D=Jacobi1（Aepsilon）
%function keshe（h1h2h3）
if ~isempty（A）
   %        -本函数使用经典Jacobi方法来计算一个对称矩阵的特征值
%        -A 是一个 n*n 矩阵的形参，由用户输入
%        -epsilon 是允许误差，由用户输入，用来控制第k次的Givens矩阵的F范数
%        -D是经过一些列的迭代后满足允许误差的矩阵的对角元素
i=1;
D=A;      
[nn]=size（A）;
flag=1;                     %设定标志位
E=abs（D-diag（diag（D）））;     %E是第n次迭代以后矩阵除去对角元的矩阵
[m1p]=max（E）;              %寻找矩阵E中最大元
[m2q]=max（m1）;
p=p（q）;                     
if p==q                     %处理矩阵本身就是一个对角阵的特殊情况
    D=diag（D）;
else
    while（flag==1）
      y=abs（D（qq）-D（pp））;
      if y==0               %处理D（pp）=D（qq）的特殊情况
         x=2*D（pq）;
         c=sin（pi/4）;
         s=sign（x）*c;
      else
         x=sign（D（pp）-D（qq））*2*D（pq）;
         c=sqrt（（1+y/（sqrt（x^2+y^2）））/2）;
         s=x/（2*c*（sqrt（x^2+y^2）））;
      end
      
     G=eye（n）;              
     G（pp）=c;               %利用cos和sin的值来构造Givens变换阵
     G（qq）=c;
     G（pq）=s;
     G（qp）=-s;
     D=G*D;
     D=D*G‘;
     E=abs（D-diag（diag（D）））;
     [m1 p]=max（abs（D-diag（diag（D））））;
     [m2 q]=max（m1）;
     p=p（q）;
     k（i）=i;
     temp=sqrt（sum（sum（E.^2）））;
     r（i）=temp;
     i=i+1;
     if（temp        flag=0;
     end
    end
  D=diag（D）;
end
%plot（kr）;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件      13824  2012-06-13 21:03  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\1.xls

     文件      13824  2012-06-13 21:08  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\2.xls

     文件      13824  2012-06-13 02:21  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\3.xls

     文件       1436  2012-06-13 21:17  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\Jacobi1.m

     文件       2677  2012-06-13 03:27  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\untitled1.fig

     文件       8019  2012-06-13 21:18  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\untitled1.m

     文件      26920  2012-06-13 03:28  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪\图像 2.png

     文件      13824  2012-06-18 16:42  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\1.xls

     文件      13824  2012-06-18 16:42  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\2.xls

     文件       1436  2012-06-18 16:42  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\Jacobi1.m

     文件       3103  2012-06-18 17:16  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\untitled.fig

     文件       8059  2012-06-19 18:20  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\untitled.m

     文件     167424  2012-06-18 13:38  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\任务书10.doc

     文件      34660  2012-06-18 17:22  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心\图像 1.png

     目录          0  2012-06-13 21:08  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值——尚少陪

     目录          0  2012-06-28 10:15  雅克比法求对称矩阵的特征值\雅克比法求对称矩阵的特征值—舒心

     目录          0  2012-06-26 09:20  雅克比法求对称矩阵的特征值

----------- ---------  ---------- -----  ----

               322854                    17