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

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

%雅可比迭代法解线性方程组；
) _6 j$ S! a# K: N: T6 o
clc;clear;close;
. Q- M, |! V& L
% digits(4);/ A' d+ z# y+ Y- ], }2 m5 h. C1 G
A=input('请输入系数矩阵：');
/ ?7 k/ q* O8 X, d: y
b=input('请输入b矩阵(列向量)：');0 K! T\" D# J. \8 t
X0=input('请输入迭代初值(列向量)：');
6 v T. ]; H- R/ e: v! L
k=input('请输入迭代次数：');6 A0 [# C' p: R% x
[M1,N1]=size(A); d3 N% D' y# I+ R0 O- H
% X0=[3.1481 4.8 2.037]';
9 z$ ?; D8 y% h J8 [
% f=X0;
/ P: v1 m2 ]\" ]! Y- i- @
% X=X0;3 V- H) a5 ]$ o( U# o- p
% vpa(X0);
2 J' Z6 F9 h7 s; m$ B7 I4 ^\" _
% vpa(X);* G% M# K1 G\" B6 w
% vpa(f);
0 j$ N! N5 W2 n. C( W( \& x! t
OUT=[];4 \6 Y% y\" Z$ v) Q2 e9 \
% vpa(OUT);
% a- B' R0 V& T\" v! J9 Q
%M=[0 -0.2222 0.0370;-0.4 0 -0.1333;-0.0185 -0.0185 0];9 b4 O: t. O4 Z- m K; F. a
%求解M矩阵；# Y' h0 _5 ^0 R _, M6 s& e\" y9 [1 d
for m=1:M1! `. ]3 g. R* O' z) a
for n=1:N1$ l4 A% Y, |# d( X$ @
if m==n. J6 m* K9 t$ {: d; G$ y
D(m,n)=A(m,n);
# n- D- _1 v4 Q; c' a/ m( w9 ^- x
else
4 [2 [4 f1 {1 F' k& V
D(m,n)=0;
; }- ^' a# H* ^
end/ \ m, [. J$ X\" e! _& R3 M
end
1 g/ z: A# r6 Z* D, a9 r* `
end
, L5 I' e/ _; u- _/ y
LU=A-D) R7 j/ b! ]7 l
M=-inv(D)*(LU)/ j5 B% x0 O) S- R
f=inv(D)*b
) o0 A( J8 D$ Q
X=X0;* ^6 ~/ g! j; K8 w* @' ]8 Y2 \% v
%特征根求解；6 `7 B: t4 {2 X: I
[v,d]=eig(M)
* h- W\" p+ e/ |2 i6 g0 m
% vpa(M);\" u7 X! g! x2 c0 K3 M' [
for i=1:k; v) e) d7 {& a8 o
X=M*X+f;- |! j( I' j4 X: O* e z( U
OUT=[OUT;X'];
7 _6 P2 c\" p# p A9 _3 W
K1(i)=i;4 c, I2 c) b: C! s0 @6 L9 C
end
2 X c+ q% e' U6 K8 \0 p
K1=[0 K1]';( c8 [\" V* q. b
OUT1=[X0';OUT];
% E' I4 h% z1 |7 v* }
OUT1=[K1 OUT1]; q+ ~1 B( g3 F: O: b' n! U
% vpa(OUT1);' C, y\" ?/ j; E% V) r
xlswrite('C:\Users\LENOVO\Desktop\雅可比迭代法2.xls',OUT1);
/ S/ f9 ~5 I\" D* P& Z5 K