48【大学】求解平方根【北理工乐学C语言作业】

小明上来大学,学习了用求解平方根的迭代公式。

其中:p 为 x 的近似平方根,e 是结果允许的误差。

输入:x, e

输出:p  保留小数点后 8 位。

#include <stdio.h>  
double Sqrt(double x, double p, double e);  
double Abs(double n);  
int main() {  
    double x, p, e;  
    scanf("%lf %lf", &x, &e);  
    p = x + 100;  
    printf("%.8f\n", Sqrt(x,p,e));  
}  
double Sqrt(double x, double p, double e)  
{  
    if (Abs(p * p - x) < e) return p;  
    else return Sqrt(x, (p + x / p) * 0.5, e);  
}  
double Abs(double n)  
{  
    return n >= 0 ? n : -n;  
}  

  • 9
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
平方根法(Cholesky分解法)是求解对称正定矩阵的线性方程组的一种有效方法。对于一个对称正定矩阵A,可以将其分解为LL^T的形式,其中L是一个下三角矩阵,L^T是其转置矩阵,然后就可以通过前代和后代求解方程组。 以下是C语言实现改进的平方根求解方程组的代码: ```c #include <stdio.h> #include <stdlib.h> #include <math.h> #define N 3 //方程组的阶数 int cholesky(double a[][N], double b[], double x[]) { int i, j, k; double sum; double l[N][N]; //下三角矩阵 //Cholesky分解 for (i = 0; i < N; i++) { sum = 0.0; for (k = 0; k < i; k++) sum += l[i][k] * l[i][k]; l[i][i] = sqrt(a[i][i] - sum); for (j = i + 1; j < N; j++) { sum = 0.0; for (k = 0; k < i; k++) sum += l[j][k] * l[i][k]; l[j][i] = (a[j][i] - sum) / l[i][i]; } } //前代求解Ly=b for (i = 0; i < N; i++) { sum = b[i]; for (j = 0; j < i; j++) sum -= l[i][j] * x[j]; x[i] = sum / l[i][i]; } //后代求解L^Tx=y for (i = N - 1; i >= 0; i--) { sum = x[i]; for (j = i + 1; j < N; j++) sum -= l[j][i] * x[j]; x[i] = sum / l[i][i]; } return 0; } int main() { double a[N][N] = {{4, -2, 2}, {-2, 3, -1}, {2, -1, 3}}; //系数矩阵 double b[N] = {7, 4, 10}; //常数向量 double x[N]; //解向量 int i; cholesky(a, b, x); printf("The solution is:\n"); for (i = 0; i < N; i++) printf("x[%d] = %f\n", i, x[i]); return 0; } ``` 这里以一个3阶方程组为例进行演示,通过调用cholesky函数求解方程组。函数中的a为系数矩阵,b为常数向量,x为解向量。程序输出的结果为方程组的解。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值