附录1
Axb
1、利用高斯消元法
clearformatratAra
d
5bra
d
51m
sizeAfori1m1
umbi
t2striforji1mAjAjAiAjiAiie
dAe
dxmAm
Ammforim111xiAi
Aii1mxi1mAiie
dx
2、利用高斯列主元消元法
clearformatratAra
d
5bra
d
51m
sizeAfori1m1
umbi
t2stridisp第
umb次选列主元后的增广矩阵tempmaxabsAimiabfi
dabsAimitemptempoAa1i1Aa1i1AiAitempodisp第
umb次消元后的增广矩阵forji1mAjAjAiAjiAiie
dAe
d
f回代过程
disp回代求解xmAm
Ammforim111xiAi
Aii1mxi1mAiie
dx3、LU分解Ara
d
5LUPluALPL
f附录2
5Ara
d
fori1
AiisumabsAie
dbra
d
1DdiagdiagALtrilA1UtriuA1MDLUfDbRhomaxabseigMR1e07switchsig
1Rhocase1dispTheJacobia
Methodis
otapplicableotherwisex1
ormr
d09
1k1whilek50
xk1Mxkfif
ormxk1xkRkk1elseXxk1IterNkbreake
de
de
dYAbER_1
ormYX1ER_2
ormYX2ER_i
f
ormYXi
f
改变谱半径的大小看收敛速度的变化
5Ara
d
fori1

fAiisumabsAie
dbra
d
1DdiagdiagALtrilA1UtriuA1MDLUfDbRhomaxabseigM2R1e07switchsig
1Rhocase1dispTheJacobia
Methodis
otapplicableotherwisex1
ormr
d09
1k1whilek50
xk1Mxkfif
ormxk1xkRkk1elseXxk1IterNkbreake
de
de
dYAbER_1
ormYX1ER_2
ormYX2ER_i
f
ormYXi
f
5Ara
d
fori1
AiisumabsAie
dbra
d
1DdiagdiagALtrilA1UtriuA1MDLUfDbRhomaxabseigM3
fR1e07switchsig
1Rhocase1dispTheJacobia
Methodis
otapplicableotherwisex1
ormr
d09
1k1whilek50
xk1Mxkfif
ormxk1xkRkk1elseXxk1IterNkbreake
de
de
dYAbER_1
ormYX1ER_2
ormYX2ER_i
f
ormYXi
f

5C
ormr
d14
ARqrCAA3eye
Co
dAco
dA2Xo
es
1bAXYAbDdiagdiagALtrilA1UtriuA1fork1100Wk02kGskDWkL1WkDWkURhokmaxabseigGske
dMi
Rhokmi
RhoIterNkwWK
fGSDwL1wDwUFSDwLwbr1e06
附录3
fclearallx3013Fu
x2x610x4fu
Fu
xplotxfu
x0xLi
eWidth3holdo
titleSolvetoalgebraicequatr