查看: 2393|回复: 0

[代码资源] 雅可比迭代法解线性方程组

[复制链接]

字体大小: 正常放大

seven3069

18 主题	15 听众	945 积分

升级 86.25%

TA的每日心情

	奋斗 2018-11-27 09:23

签到天数: 395 天

[LV.9]以坛为家II

群组: 2016国赛备战群组

电梯直达

1^#

发表于 2016-10-12 22:38 |只看该作者 |正序浏览

|招呼Ta 关注Ta

雅可比迭代法解线性方程组：

%雅可比迭代法解线性方程组；
5 r6 l, ^3 ]6 v- ?8 V
clc;clear;close;\" r% I' p+ a2 x2 h! Y9 W
% digits(4);. p: y7 `5 t2 r9 [\" k
A=input('请输入系数矩阵：');
, ^+ R' ~) F7 b+ x& C! c* I$ }3 M
b=input('请输入b矩阵(列向量)：');9 X/ O\" E, \- j3 O% l
X0=input('请输入迭代初值(列向量)：');' g }1 b% E( u7 ]2 _& G
k=input('请输入迭代次数：');
* Z& ?7 Z' C& X+ l( \, C. a
[M1,N1]=size(A);1 H7 L4 i% u$ Y, W& Y; K% o, L
% X0=[3.1481 4.8 2.037]';. D& ~# i$ a# R5 J. I7 T5 h; S. d
% f=X0;9 C5 R! k: H' T
% X=X0;
4 k# \$ _. ]0 F
% vpa(X0);
9 B3 Z. P; b8 [+ c5 f
% vpa(X);
8 `9 y4 C% P1 C
% vpa(f);
- K0 H/ M7 D7 k. N' j3 t\" |! W
OUT=[];
8 E1 y! M- h7 d) v. s
% vpa(OUT);
! P% W! _9 G+ W) s3 E5 L$ M$ d
%M=[0 -0.2222 0.0370;-0.4 0 -0.1333;-0.0185 -0.0185 0];
2 i* A5 `$ [1 ^! w7 ?( Z: @
%求解M矩阵；, w) L: ?# j1 q- [) I9 E, n
for m=1:M1
: S; S7 o. J: }* G
for n=1:N1& ~+ \3 C3 X1 z& \! r! s; o3 b6 `
if m==n; e7 ]& E. I# J7 C5 z\" t7 y
D(m,n)=A(m,n);( A' ?7 W% i5 S; Q- N
else6 T: z7 \ o\" L6 S8 K
D(m,n)=0;
3 Q& S E& U, H5 e0 E; x
end& i7 n j: A5 n0 [& ]
end0 w) I5 S+ ~- c' h( o
end* c; ?# V6 g V% O
LU=A-D
! ?: y; o\" @8 r$ g/ Y\" h
M=-inv(D)*(LU)8 k* u+ h( o4 y6 o- Z% ]
f=inv(D)*b, M2 p\" S7 v% o8 ^
X=X0;
3 F3 e( v' a- Q6 q
%特征根求解；
, b2 A. m\" |+ T# A+ I* j* C
[v,d]=eig(M)
5 E6 l2 p6 D3 ?) m: F8 m
% vpa(M);
for i=1:k;
# W; o$ h* r9 v7 ]0 ?
X=M*X+f;
' m4 J3 b5 ]$ q# U
OUT=[OUT;X'];7 ~7 W5 x8 k, X x9 o
K1(i)=i;
1 C6 l+ W+ w9 ~) z0 U% l
end. ~7 |, ^/ x8 e# B' W+ d: Z
K1=[0 K1]';
\" U/ K9 w/ p# }8 w/ E: U
OUT1=[X0';OUT];; j5 n- }/ _/ `5 g4 P, K
OUT1=[K1 OUT1]
$ g$ b+ K4 W2 M) ` _6 C
% vpa(OUT1);( P, E- _\" c\" Z
xlswrite('C:\Users\LENOVO\Desktop\雅可比迭代法2.xls',OUT1);
% I\" u3 O1 E% D% e% B