%下面是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