博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
高斯消元法
阅读量:6311 次
发布时间:2019-06-22

本文共 902 字,大约阅读时间需要 3 分钟。

前言:

利用高斯消元可以求解线性方程组,复杂度 $O(n^3)$

正文:

实现过程有点类似代入消元法

最后将矩阵消成一个倒三角形

最后一行只有一个未知数

倒数第二行有两个,依此类推

所以可以从最后一行解出一个未知数的值

然后往上回带,直至求解出所有未知数

#include
#include
#include
const double eps=1e-10;int n;double ans[111];double mat[111][111];int main(){ scanf("%d",&n); for(int i=1;i<=n;i++) for(int j=1;j<=n+1;j++) scanf("%lf",&mat[i][j]); for(int i=1;i<=n;i++) { int pos=i; for(int j=i+1;j<=n;j++) if(fabs(mat[j][i])>fabs(mat[pos][i])) pos=j;//选出系数绝对值最大的主元,提高精度 if(fabs(mat[pos][i])
=1;i--) { ans[i]=mat[i][n+1]; for(int j=i+1;j<=n;j++) ans[i]-=mat[i][j]*ans[j];//回带 } for(int i=1;i<=n;i++) printf("%.2lf\n",ans[i]); return 0;}

后序:

据说高斯消元好像很少考

或者是我太菜了,总之我现在才会写......

转载于:https://www.cnblogs.com/Vscoder/p/10675589.html

你可能感兴趣的文章
vue-cli脚手架的搭建
查看>>
在网页中加入百度搜索框实例代码
查看>>
在Flex中动态设置icon属性
查看>>
采集音频和摄像头视频并实时H264编码及AAC编码
查看>>
3星|《三联生活周刊》2017年39期:英国皇家助产士学会于2017年5月悄悄修改了政策,不再鼓励孕妇自然分娩了...
查看>>
linux查看命令是由哪个软件包提供的
查看>>
高级Linux工程师常用软件清单
查看>>
堆排序算法
查看>>
folders.cgi占用系统大量资源
查看>>
路由器ospf动态路由配置
查看>>
zabbix监控安装与配置
查看>>
python 异常
查看>>
last_insert_id()获取mysql最后一条记录ID
查看>>
可执行程序找不到lib库地址的处理方法
查看>>
bash数组
查看>>
Richard M. Stallman 给《自由开源软件本地化》写的前言
查看>>
oracle数据库密码过期报错
查看>>
修改mysql数据库的默认编码方式 .
查看>>
zip
查看>>
How to recover from root.sh on 11.2 Grid Infrastructure Failed
查看>>