| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 407 人关注过本帖
标题:各位大侠帮帮忙啊,高斯列主元回代最后一个解为什么会出错啊
只看楼主 加入收藏
xfirrl
Rank: 1
等 级:新手上路
帖 子:1
专家分:0
注 册:2009-11-4
结帖率:0
收藏
已结贴  问题点数:20 回复次数:1 
各位大侠帮帮忙啊,高斯列主元回代最后一个解为什么会出错啊
数值分析里面用高斯列主元方法解方程的,后8个解求得正确,最后一个解出错(也就是回代的最后一步出错了),我实在找不出原因,大家帮帮忙啊

#include<iostream.h>
#include<math.h>
#include<iomanip.h>
void main(void)
{
    double a[9][10]={{31,-13,0,0,0,-10,0,0,0,-15},{-13,35,-9,0,-11,0,0,0,0,27},{0,-9,
        31,-10,0,0,0,0,0,-23},{0,0,-10,79,-30,0,0,0,-9,0},{0,0,0,-30,57,-7,0,-5,0,-20},
    {0,0,0,0,-7,47,-30,0,0,12},{0,0,0,0,0,-30,41,0,0,-7},{0,0,0,0,-5,0,0,27,-2,7},{
        0,0,0,-9,0,0,0,-2,29,10}};
    double x[8]={0,0,0};
    double z=0;//交换两行时用的中间值//
    int k=0;//列主元时记录最大主元的所在行//
    for(int i=0;i<=7;i++)//找最大主元//
    {
        for(int j=1;j+i<=8;j++)
        {
            if(fabs(a[i+j][i])>fabs(a[i][i]))
                k=i+j;
            else
                k=i;
        }
        for(int v=i;v<=9;v++)//交换两行//
        {
            z=a[i][v];
            a[i][v]=a[k][v];
            a[k][v]=z;
        }
        for(int w=i+1;w<=8;w++)//列主元后消元,w为计算时的行空竹量//
        {
            double H=-a[w][i]/a[i][i];
            for(int x=0;i+x<=9;x++)//i+x为计算的列控制量//
            {
                a[w][i+x]=a[w][i+x]+a[i][i+x]*H;
            }
        }
    }
        for(int m=0;m<=8;m++)//输出处理后的矩阵//
        {
            for(int n=0;n<=9;n++)
            {
                if(n!=0&&n%9==0)
                cout<<setw(14)<<a[m][n]<<endl;
            else
                cout<<setw(14)<<a[m][n];
            
            }
        }
        cout<<a[0][0]<<endl;//检测a[0][0]//
        x[8]=a[8][9]/a[8][8];
        for(int n=7;n>=0;--n)//回代过程开始//
        {
            double g=0;
            for(int m=n+1;m<=8;m++)
            {
                g+=a[n][m]*x[m];
            }
            x[n]=(a[n][9]-g)/a[n][n];
            cout<<setw(18)<<a[n][n]<<setw(18)<<g<<setw(18)<<x[n]<<endl;
        }//回代结束//
        cout<<a[0][0]<<endl;//检测a[0]//
        for(int t=0;t<=8;t++)
            cout<<setw(14)<<x[t];
        cout<<endl;
}
运行结果那个-30.8937是错的(应该是-0.28923),后面几个是对的。
搜索更多相关主题的帖子: 高斯 
2009-11-04 16:42
flyingcloude
Rank: 10Rank: 10Rank: 10
等 级:青峰侠
威 望:6
帖 子:598
专家分:1512
注 册:2008-1-13
收藏
得分:20 
回复 楼主 xfirrl
没发现你说的那个数

你能学会你想学会的任何东西,这不是你能不能学会的问题,而是你想不想学的问题
2009-11-05 21:11
快速回复:各位大侠帮帮忙啊,高斯列主元回代最后一个解为什么会出错啊
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.025229 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved