師兄們的參儀器科學(xué)與技術(shù)

2（P184）解

儀器科學(xué)與技術(shù)冪法的算法5.1的程序如下function[lambda1,x1,k]=mypower(A,v,ep)%A為系數(shù)矩陣，v為初始向量，eplambda=v(iu=v/lambda;%規(guī)格化x1=u;%使用收縮技術(shù)求矩陣的特征值的程序如%function%A為系數(shù)矩陣，lambda1為按模最大的特征值，x1forj=1:n-

H=I-2*w*w'用定義求HouseholderC=(H*B計(jì)算使用原點(diǎn)位移技術(shù)改進(jìn)解的準(zhǔn)確度的程序如下function[lambda1,x1,k]=myyuandian(A,v,p,ep)A為系數(shù)矩陣，v為初始向量，p為原點(diǎn)位移參數(shù)，eplambda=v(iu=v/lambda;%規(guī)格化x1=u;%使用瑞利商法加速收斂的程序如下function[lambda1,x1,k]=myrayleigh(A,v,ep)%A為系數(shù)矩陣，v為初始向量，eplambda=v(iu=v/lambda;%規(guī)格化x1=u;%使用反冪法改進(jìn)解的準(zhǔn)確度的程序如下function[lambdan,xn,k]=myinverse(A,v,ep)%A為系數(shù)矩陣，v為初始向量，eplambda=1/v(iu=lambda*v;%規(guī)格化xn=u;%M=[200;030;00K=[2-10;-12-1;0-1得到結(jié)果為lambda1=x1= - k=lambda= w1=w2=w3=使用原點(diǎn)位移技術(shù)改進(jìn)解的準(zhǔn)確度由于原點(diǎn)位移參數(shù)p不好取，因此在程序中試探得p=0.3時(shí)能加速收斂。M=[200;030;00K=[2-10;-12-1;0-1得到結(jié)果為lambda1=x1= - k=lambda= w1=w2=w3=M=[200;030;00K=[2-10;-12-1;0-1得到結(jié)果為lambda1=x1= - k=w1=w2=w3=M=[200;030;00K=[2-10;-12-1;0-1得到結(jié)果為w1=w2=w3=3（P185）解對(duì)矩陣A 中實(shí)現(xiàn)基本QR算法的程序如下function[lambda,B]=myqr(A,k)%A為系數(shù)矩陣，kfori=1:klambda=diag(A);%B=A最后的擬上三角陣A=[0.50.50.50.5;0.50.5-0.5-0.5;0.5-0.50.5-0.5;0.5-0.5-0.5得到結(jié)果為lambda B 分析：AQR算法后依然保持不變，且特征值均為A對(duì)角