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

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

%雅可比迭代法解线性方程组；* P\" l' Y$ W7 g
clc;clear;close;
% v( e% ~8 J8 s8 Z$ G5 l+ Z3 ? `
% digits(4);
7 V' @. b, r i. _8 S% r
A=input('请输入系数矩阵：');# X3 {1 q6 W) m* l
b=input('请输入b矩阵(列向量)：');
+ y. j7 B# x( D- k
X0=input('请输入迭代初值(列向量)：');
$ _4 _2 Q, `: e& L
k=input('请输入迭代次数：');+ {2 i) D# h( K0 [% _$ f
[M1,N1]=size(A);
; {8 h1 X- W8 s G2 [6 y/ F5 b; F9 R
% X0=[3.1481 4.8 2.037]';& a8 z) P* k( s* Z. I\" k6 c5 T
% f=X0;
% o* S1 x' o6 L
% X=X0;/ F+ ? ^* P: j* P7 _ Y
% vpa(X0);* T) ]0 K, n: p% ~& j) d
% vpa(X);( R3 `- E1 b, @ s
% vpa(f);1 w# X) F; h \0 A& g2 ~
OUT=[];
7 g L% K; P- [9 F
% vpa(OUT);
) M- v\" H- z- [) g, a
%M=[0 -0.2222 0.0370;-0.4 0 -0.1333;-0.0185 -0.0185 0];
5 R: C! Q4 d, M; F- d7 o
%求解M矩阵；5 `3 i) u+ E+ o2 [
for m=1:M1/ x. \\" z9 Q\" L\" H
for n=1:N1* _: _! x% O9 b
if m==n! P' w8 M V1 a5 N) L1 s
D(m,n)=A(m,n);
( g5 N6 _4 o' e! D4 ]( b$ V1 I
else1 w) o. _8 Q( ^
D(m,n)=0;2 p- H- {) _/ S+ F' E
end
) w$ L% G5 e\" m5 y/ ?
end; y+ n: x\" `9 w& Z
end2 d+ i$ W/ A7 Z
LU=A-D1 {$ w% D. i5 m9 ?* j, n
M=-inv(D)*(LU)
, ]' y# w* v\" \\" e# o4 }
f=inv(D)*b
# V1 _ v4 }2 X* I O
X=X0;
# C& @) a\" E; D* [' K
%特征根求解；
# M4 V4 Y8 O \
[v,d]=eig(M)\" G/ U9 }2 `) z4 ?2 m+ l3 L
% vpa(M);
. f. v0 i1 }* B# J# B8 I- P9 T
for i=1:k;
& d! h* P0 K! U8 v2 p/ y3 y
X=M*X+f;
8 @8 x4 |; \/ E! k c( q
OUT=[OUT;X'];/ W6 ]- l9 s; r% ~\" ~
K1(i)=i;1 ~) g9 |& K% Y6 z4 [
end1 [7 f* B6 [4 O5 a) I0 [& K& C3 x
K1=[0 K1]';: w8 i2 ^( A& d! ?+ M
OUT1=[X0';OUT];) W) K! p! m6 q- N( N
OUT1=[K1 OUT1]
0 P+ c% J! U* K$ d
% vpa(OUT1);+ I j% s\" D# b6 w5 Z3 P
xlswrite('C:\Users\LENOVO\Desktop\雅可比迭代法2.xls',OUT1);
& d0 y2 i/ E$ {: x$ d