虫虫首页|资源下载|资源专辑|精品软件|心愿单
m 登录|注册

您现在的位置是:虫虫下载站 > 资源下载 > 源码 > C语言编写雅可比迭代

C语言编写雅可比迭代

  • 资源大小:13 K
  • 上传时间: 2019-10-13
  • 下载次数:0次
  • 浏览次数:5次
  • 上传用户:大萌萌撒
  • 资源积分:2 下载积分
  • 标      签: C语言 编写 迭代

资 源 简 介

# include<stdio.h>
# include<math.h>
# define N 3

main(){
    float NF2(float *x,float *y);
    float A[N][N]={{10,-1,-2},{-1,10,-2},{-1,-1,5}};
    float b[N]={7.2,8.3,4.2},sum=0;
    float x[N]= {0,0,0},y[N]={0},x0[N]={};
    int i,j,n=0;
    for(i=0;i<N;i++)
    {
        x[i]=x0[i];
    }
    for(n=0;;n++){
                //计算下一个值
    for(i=0;i<N;i++){
        sum=0;
        for(j=0;j<N;j++){
            if(j!=i){
                sum=sum+A[i][j]*x[j];
            }
        }
        y[i]=(1/A[i][i])*(b[i]-sum);
        //sum=0;
    }
    //判断误差大小
        if(NF2(x,y)>0.01){
                for(i=0;i<N;i++){
        x[i]=y[i];
    }
    }
            else
            break;
    }
    printf("经过%d次雅可比迭代解出方程组的解:\n",n+1);
    for(i=0;i<N;i++){
        printf("%f      ",y[i]);
    }
}
//求两个向量差的二范数函数
float NF2(float *x,float *y){
int i;
float z,sum1=0;
for(i=0;i<N;i++){
    sum1=sum1+pow(y[i]-x[i],2);
}
z=sqrt(sum1);
return z;
}

您 可 能 感 兴 趣 的

相 关 代 码

相 关 资 源

该 用 户 还 上 传

VIP VIP
  • 865394974 5小时前 开通了:月度VIP会员

  • 2789 15小时前 开通了:月度VIP会员

  • a2533505 1天前 开通了:月度VIP会员

  • ccacy 1天前 开通了:月度VIP会员

  • Lyq000624 2天前 开通了:月度VIP会员

  • yuhao129 2天前 开通了:月度VIP会员

  • zhuchen66 2天前 开通了:月度VIP会员

  • aaooaa 2天前 开通了:月度VIP会员

  • kelecn 5天前 开通了:月度VIP会员

  • pxp1234 6天前 开通了:月度VIP会员