查看: 2551|回复: 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

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

%雅可比迭代法解线性方程组；2 Q# ~% D. v( W9 b
clc;clear;close;5 E, R. r4 V* T& O7 a; O+ V
% digits(4);, r/ i' s+ P2 W# n
A=input('请输入系数矩阵：');
7 l' q0 g1 a& d+ Q& v% B
b=input('请输入b矩阵(列向量)：');
$ P( S9 ?/ i9 x- o& u
X0=input('请输入迭代初值(列向量)：');
! r6 A! g4 g9 I5 g2 H9 ^1 l2 ]
k=input('请输入迭代次数：');
( r; g9 Y( e) Z: t4 R* y( h
[M1,N1]=size(A);9 G3 i& g9 e* q+ G+ ~6 D
% X0=[3.1481 4.8 2.037]';
3 ^& m/ r7 a% d# V8 e, g7 M* F2 U& X; Q
% f=X0;
- t* t6 b: _. k) J8 M r7 s
% X=X0;
( }/ f( [4 `, Y/ z2 k\" S5 O
% vpa(X0);
$ t# K& D0 Y& p$ R
% vpa(X);0 |. A, M- {; D \2 D! X) e& R+ H
% vpa(f);7 y- h/ p' ]0 s5 a
OUT=[];8 C1 @9 `; |/ Q% \0 `: I; p, g# j
% vpa(OUT);
, q0 n+ Q( Y- |+ j
%M=[0 -0.2222 0.0370;-0.4 0 -0.1333;-0.0185 -0.0185 0];9 l$ ?) |5 D5 ~8 x& ]1 c! f
%求解M矩阵；
6 f9 K8 d\" C( A& w0 Y
for m=1:M12 s: v( V9 o2 a
for n=1:N1
6 r( u+ w: P; i( C9 z
if m==n+ m' M2 W! U. y& X# t1 a' Y) n
D(m,n)=A(m,n);. u1 Q9 M: |. X: c
else1 h+ f* H2 P+ Z
D(m,n)=0;
3 Z4 |6 E\" t/ Z( I
end$ F, T5 w- g2 I
end
5 n) k8 u3 V2 M$ S
end
& @# j# ~+ c3 a# `. G# ?5 J+ R
LU=A-D; ]+ b1 a* n8 B4 t0 r
M=-inv(D)*(LU)' O+ |3 v& a8 r9 l8 C
f=inv(D)*b
2 T+ K. n6 y- d& j* o
X=X0;9 { [\" Q4 R' ^# K. ]\" H( E1 C
%特征根求解；
7 _& D5 {3 m2 t; h8 m* a
[v,d]=eig(M)
1 s- Y. M4 c% o: o, e! ~
% vpa(M);
8 a% e) @7 Z, k* X( ]* B% K
for i=1:k;
6 J, Z. B$ j! N: B' b; |
X=M*X+f;
& o- A8 W+ M. N+ C4 w
OUT=[OUT;X'];0 r1 S u3 J. F; N, {
K1(i)=i;
/ |) \, b+ L9 m. X1 A, w- E( j. j
end
/ [. _7 L7 F' n/ ~7 N
K1=[0 K1]';% g: I# I1 y& W\" t2 m. h\" ]
OUT1=[X0';OUT];7 b) P- I3 x; ^/ O3 j0 P5 \
OUT1=[K1 OUT1]- l9 e+ b( |' Q. J1 f
% vpa(OUT1);* n$ o8 }3 i9 }: |3 q
xlswrite('C:\Users\LENOVO\Desktop\雅可比迭代法2.xls',OUT1);
8 h0 G1 |; }3 q$ g5 K, c3 V6 P& g