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

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

%雅可比迭代法解线性方程组；
9 H- a+ @. J, _- v\" d4 N
clc;clear;close;
0 j. c. w. F) X& C1 b& ~
% digits(4);. n8 l2 W, Q: j0 W\" u! M. ^
A=input('请输入系数矩阵：');
$ r' M2 I( r3 n0 T' H, z- a
b=input('请输入b矩阵(列向量)：');7 D5 N8 t/ W$ T9 Q
X0=input('请输入迭代初值(列向量)：');! _\" r. n4 I7 }8 G
k=input('请输入迭代次数：');
7 f/ [: W8 |- l0 c) _
[M1,N1]=size(A);. [4 a# m: h' v C
% X0=[3.1481 4.8 2.037]';
, N, x: \5 O: M& Y\" C2 l# m
% f=X0;$ D! w9 C# g6 T* K
% X=X0;. s. }' l: H: I5 R
% vpa(X0);
5 z2 b5 m5 K1 U* v5 _) Q5 R$ ?7 k
% vpa(X);
; r# c& g0 [# S0 D$ p' a
% vpa(f);2 E% \$ W2 P0 J2 ~ \: x U
OUT=[];& O$ Z4 U* q/ h3 E( ^
% vpa(OUT);+ ~5 j$ y( S4 ^4 {/ w: {! l1 H
%M=[0 -0.2222 0.0370;-0.4 0 -0.1333;-0.0185 -0.0185 0];
9 @\" {1 Z1 k; B# l4 i% @
%求解M矩阵；7 g# O. L% B0 o5 l) f
for m=1:M1
$ a2 {/ M% t c% p4 ]/ K
for n=1:N16 a+ K. y: F. r3 G/ R2 M& }. y8 t
if m==n
1 o6 B1 A$ s* Y& p- w\" H
D(m,n)=A(m,n);' t: V/ j\" U# y2 ^. s$ o# z\" y& i
else
; u# ?, `* f2 p- U* E, \4 U
D(m,n)=0;7 \7 X+ A6 x* c
end
* a% c) o# O' Q6 \
end: C0 Q) R$ `9 p( k' f
end/ k* `8 ?$ f2 G3 g% g: R
LU=A-D6 B$ E0 p: I, S0 m
M=-inv(D)*(LU)3 K: V% {# a0 q3 s- y0 m3 ~ q
f=inv(D)*b
! X# U8 f d( z2 a9 H' h
X=X0;
/ N* j) t( v* Z' x3 e7 ~1 K E
%特征根求解；
, Z. H7 f4 P* J+ @& P$ O
[v,d]=eig(M). P' K. x+ J, z, s$ X7 r
% vpa(M);+ @ t\" _- H; t' U+ v
for i=1:k;4 _8 |0 ]6 u: g) D$ v
X=M*X+f;) G# z4 g* y. a/ C
OUT=[OUT;X'];
: E o) N3 [/ [6 o
K1(i)=i;
. x1 e) z( E L9 r, M
end
7 O' F5 j. c' z# @- R
K1=[0 K1]';
( M0 z; ^8 ?6 T6 G4 |* x
OUT1=[X0';OUT];% S% D* s& N\" N0 t$ ?9 F. F3 }
OUT1=[K1 OUT1]
\" {/ f. g# B5 u5 n8 K1 u1 S
% vpa(OUT1);/ A* m& W7 D6 U8 p- v\" c( H7 @+ Z
xlswrite('C:\Users\LENOVO\Desktop\雅可比迭代法2.xls',OUT1);3 C; h' V# ?5 r9 N0 C; r% d