数值分析 2.2 平方根法

最新推荐文章于 2022-11-05 13:12:46 发布

mathfinder

最新推荐文章于 2022-11-05 13:12:46 发布

阅读量1.9k

点赞数

分类专栏：数值分析算法实现文章标签：数值分析

本文链接：https://blog.csdn.net/leader880022/article/details/44887965

版权

本文介绍了一种针对系数矩阵为正定矩阵的方程组求解方法——平方根法。通过乔列斯基分解，有效地计算G矩阵，并在回带过程中求得解向量x。代码示例展示了如何实现这一算法。

摘要由CSDN通过智能技术生成

///对于系数矩阵是正定矩阵的方程组可以用平方根法
#include <iostream>
#include <cstdio>
#include <math.h>
using namespace std;
const int MAXN = 1000;
double a[MAXN][MAXN];
double b[MAXN];
double x[MAXN];
double y[MAXN];
int main()
{
    int n;
    //freopen("in.txt","r",stdin);
    while(cin >> n){
        for(int i = 0;i<n;i++)cin >> b[i];
        for(int i = 0;i<n;i++){
            for(int j = 0;j<n;j++){
                cin >> a[i][j];
            }
        }
        //memset(l,0,sizeof(l));
        //memset(u,0,sizeof(u));
        ///乔列斯基分解
        ///为了省空间把G矩阵存放在a数组当中
        for(int k = 0;k<n;k++){///G转置的第k列开始
            for(int j = 0;j<k;j++){
                a[k][k]-=a[k][j]*a[k][j];///先求主元素
            }
            a[k][k] = sqrt(a[k][k]);
            ///把G矩阵的第k列元素依次求出来
            for(int i = k+1;i<n;i++){