雅克比(Jacobi)迭代法解线性方程组(Matlab程序)

原创 2012年03月30日 18:36:30

%---雅克比迭代法-----
%---Jacobi iteration method
clear;clc;
% A=[10,-1,-2;-1,10,-2;-1,-1,5];
% b=[72,83,42]';
A=[ 28,-3,0,0,0;
    -3,38,-10,0,-5;
    -10,0,25,-15,0;
    0,0,-15,45,0;
    0,-5,0,0,30];
b=[10,0,0,0,0]';
N=length(b);
fprintf('库函数计算结果:');
x=inv(A)*b   %库函数计算结果
B=zeros(N,N);g=zeros(N,1);
x=zeros(N,1);%迭代初始值
eps=0.001;%相邻解的距离小于该数时,结束迭代
% for i=1:N
%     B(i,:)=-A(i,:)/A(i,i);
%     B(i,i)=0;       %迭代矩阵
%     g(i)=b(i)/A(i,i);
% end
D=diag(diag(A));
E=-tril(A,-1);%下三角
F=-triu(A,1);%上三角
B=inv(D)*(E+F);g=inv(D)*b;

for k=1:100 %最大迭代次数为100
    fprintf('第%d次迭代:',k);
    y=B*x+g;
    fprintf('\n与上次计算结果的距离(2范数):%f \n',norm(x-y)^2);
    if norm(x-y)<eps
        break;
    end
    x=y
end
x

 

数值分析 jacobi迭代法求解线性方程组 MATLAB程序实现

数值分析 jacobi迭代法MATLAB程序实现   函数内容(保存为jacobi.m文件,在MATLAB中运行): %function [G,d,x,N]=jacobi(A,b) %Jacob...
  • jingmiaa
  • jingmiaa
  • 2015年11月06日 19:01
  • 5772

matlab 雅可比矩阵求取

syms x1 x2 ;%定义变量 a = [x1 x2]; f = [x1*x1+t*x1;x2*cos(x2)]; %定义函数,以矩阵的形式 x = jacobian(f,a); % 求取雅可比...
  • u012784288
  • u012784288
  • 2016年10月08日 23:00
  • 3575

【数值分析】迭代法解方程:牛顿迭代法、Jacobi迭代法

本科课程参见:《软件学院那些课》 牛顿迭代公式 设已知方程f(x)=0的近似根x0 ,则在x0附近f(x)可用一阶泰勒多项式近似代替.因此, 方程f(x)=0可近似地表示为p(x)=0。用x1...
  • xiaowei_cqu
  • xiaowei_cqu
  • 2013年02月18日 14:31
  • 56542

MATLAB样例之雅克比迭代法

  • 2013年11月07日 21:08
  • 15KB
  • 下载

雅克比迭代算法

#include #include #include #include using namespace std; const int n=3;          //设置方程组的维 float a[n...
  • andyhou
  • andyhou
  • 2007年01月25日 21:50
  • 1242

Jacobi迭代与SOR迭代求解希尔伯特矩阵

给出线性方程组 Hn*x = b,其中系数矩阵Hn为希尔伯特矩阵:      假设 x ∗ =(1, 1, . . . , 1)T,b = Hnx ∗。若取 n = 6,8, 10,分别用 Jaco...
  • qqlu_did
  • qqlu_did
  • 2015年04月27日 01:10
  • 2112

Jacobi迭代法解线性方程matlAB源码

  • 2015年12月07日 22:10
  • 1KB
  • 下载

Jacobi迭代与SOR迭代求解希尔伯特矩阵

给出线性方程组 Hn*x = b,其中系数矩阵Hn为希尔伯特矩阵:      假设 x ∗ =(1, 1, . . . , 1)T,b = Hnx ∗。若取 n = 6,8, 10,分别用 Jaco...
  • qqlu_did
  • qqlu_did
  • 2015年04月27日 01:10
  • 2112

【数值分析】迭代法解方程:牛顿迭代法、Jacobi迭代法

本科课程参见:《软件学院那些课》 牛顿迭代公式 设已知方程f(x)=0的近似根x0 ,则在x0附近f(x)可用一阶泰勒多项式近似代替.因此, 方程f(x)=0可近似地表示为p(x)=0。用x1...
  • xiaowei_cqu
  • xiaowei_cqu
  • 2013年02月18日 14:31
  • 56542

线性方程组——Jacobi迭代和G_S迭代

1)方程组Ax=b 用简单迭代法(Jacobi迭代法)和G-S迭代法分别解方程组。精确到小数点后5位,最大迭代次数N=100,说明是否收敛 Jacobi迭代代码如下 Jacobi.m f...
  • john_bian
  • john_bian
  • 2016年11月09日 15:02
  • 1735
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:雅克比(Jacobi)迭代法解线性方程组(Matlab程序)
举报原因:
原因补充:

(最多只允许输入30个字)