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

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

%雅可比迭代法解线性方程组；' a$ F4 t _4 D( U
clc;clear;close;
8 B6 E! t3 G2 F5 p+ t# | O' l0 r2 ?
% digits(4);+ x' A5 M5 x$ V$ B: W
A=input('请输入系数矩阵：');2 C\" K% ]) q- u* ]
b=input('请输入b矩阵(列向量)：');
- v' r9 m, G) b8 D' J; i& X
X0=input('请输入迭代初值(列向量)：');( L i% G5 o9 L z# F6 n
k=input('请输入迭代次数：');
! W& a J. O1 T% a& G X. _
[M1,N1]=size(A);% j& x* O/ g4 M3 ]4 F2 \! x
% X0=[3.1481 4.8 2.037]';; ?, M% i. M9 M7 k/ M$ g
% f=X0;9 M! W7 I ~2 A. ^( f, f
% X=X0;9 q o8 A2 J$ A' y
% vpa(X0);
$ f1 r% V; | G% ~4 V& K# o' C
% vpa(X);
6 P! Q- H( q5 [/ y' ?
% vpa(f);
) A9 E2 S$ h* ~
OUT=[];# m( O. i6 @% h# [6 B; |
% vpa(OUT);9 j2 s7 D- d; {
%M=[0 -0.2222 0.0370;-0.4 0 -0.1333;-0.0185 -0.0185 0];
4 b$ f$ G6 o9 O) O- n
%求解M矩阵；
% E+ k! {! t+ w6 D: s
for m=1:M1& {+ _4 M( U& B {( w8 O- I\" B
for n=1:N1
- R6 ]0 y. ]% m; H+ X
if m==n9 b\" ~1 n# { q. s
D(m,n)=A(m,n);( S7 D& T\" q7 P- P% v& I
else+ m! m9 m- }6 v: Z# z0 I, G
D(m,n)=0;6 |8 ~( G+ K4 l, i% P9 z- P) s$ J8 v
end
+ J! i* q' u& e: c
end
& \9 i1 j6 b- g3 O- J
end\" h; O+ u1 C4 N& W, U9 a
LU=A-D
4 B! X, V6 H$ ~2 \) g
M=-inv(D)*(LU)
6 s9 f$ q) I( Z, f
f=inv(D)*b( x& D. j+ [( X3 }4 `5 a
X=X0;3 t6 n' b0 S* E
%特征根求解；
0 C& ]! r% m8 D5 ^! q
[v,d]=eig(M)5 I: \4 A8 n3 e
% vpa(M);
' h6 D\" o, x5 }3 i7 B# \ n3 {
for i=1:k;
$ E6 j# A5 a7 j7 t# f
X=M*X+f;
# C/ O# F4 u- H2 l! U
OUT=[OUT;X'];
' h3 _2 y5 q8 K. K- Z
K1(i)=i;
. ~8 M+ d- X7 V3 r! S
end: T0 x8 N5 h9 \3 ~7 {. q& w3 s% |; D
K1=[0 K1]';
6 _: A8 b) d1 M% Q2 j9 E6 ~6 A
OUT1=[X0';OUT];, l$ ^2 g4 o9 J& |\" u2 @% T
OUT1=[K1 OUT1]: P, r3 h! V\" Q5 n& |. _8 v8 Y: V
% vpa(OUT1);7 w3 |& y( c' | o
xlswrite('C:\Users\LENOVO\Desktop\雅可比迭代法2.xls',OUT1);
- Y: {! K/ d; X' x5 e1 h9 p& F* X