雅可比迭代法

这里写图片描述

这里写图片描述

//节省时间不设置精度
//ps:注意收敛条件
#include <stdio.h>
#include <stdlib.h>
int n = 3;//3阶行列式 
double sum(int i,int a[3][3],double x0[3]){
    double result = 0.0;
    for(int j = 0 ; j < n ; j++){
        if( j != i){
            result += a[i][j] * x0[j];
        }
    }

    return  result;
}
int main(){

    int A[3][3] = {5, 1, 2, 3, 8, -4, -3, 2, 7};
    int b[3] = {21, -23, 15};
    //输出A和b 
    printf("A:\n");
    for(int i = 0; i < n; i++){
        for(int j = 0; j < n; j++){
            printf("%d\t", A[i][j]);
        }
        printf("\n");
    }
    printf("b:\n");
    for(int i = 0; i < n; i++){
        printf("%d\t", b[i]);
    }
    printf("\n");
//  double e;
    double x0[3] = {0,0,0};//初始化向量 
    double t[3] = {0,0,0};//临时变量 
    double x[3] = {0,0,0};//k+1
    int times = 20;//迭代次数 
    for(int  k = 0; k < times; k++){
        for(int i = 0; i < n; i++){

        x[i] =  (b[i] - sum(i,A,x0)) / A[i][i] ;
        t[i] = x[i];    

        }
        printf("\n第%d次迭代:\n",k+1);
        for(int j = 0; j < n; j++){
            x0[j] = t[j];
            printf("x0[%d]=%.6lf\t",j,t[j]);
        }   //更新x0里的值并输出 
    }
    return 0;
}

测试结果:

这里写图片描述

  • 5
    点赞
  • 34
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
Matlab中的可比是一种用于求解线性方程组的代方可比是一种基于分解的方,它通过将方程组表示为对角矩阵和非对角矩阵的乘积形式来进行代计算。 在Matlab中,可以使用以下步骤来实现可比: 1. 将线性方程组表示为系数矩阵A和右侧常数向量b的形式。 2. 初始化一个初始解向量x0。 3. 根据可比代公式$x^{(k+1)} = D^{-1}(b - (L+U)x^{(k)})$,计算下一个近似解$x^{(k+1)}$,其中D是A的对角矩阵,L和U分别是A的严格下三角部分和严格上三角部分。 4. 重复步骤3,直到满足指定的收敛条件(如误差限或代次数)为止。 5. 返回最终的近似解$x^{(k+1)}$。 可比的关键思想是将原始线性方程组转化为一个代过程,通过不断更新解向量来逐步逼近精确解。这种方通常在稀疏矩阵和大规模线性方程组求解中非常有效。 请注意,使用可比求解线性方程组时,收敛性往往取决于系数矩阵A的性质。对于某些特殊的矩阵,可比可能会出现收敛问题,需要采用其他更稳定的方。<span class="em">1</span><span class="em">2</span> #### 引用[.reference_title] - *1* *2* [MATLAB Jacobi 求解线性方程组](https://blog.csdn.net/qq_40264271/article/details/120611627)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值