高斯列主消元法法和LU分解法.doc
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《高斯列主消元法法和LU分解法.doc》由会员分享,可在线阅读,更多相关《高斯列主消元法法和LU分解法.doc(6页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、【精品文档】如有侵权,请联系网站删除,仅供学习与交流高斯列主消元法法和LU分解法.精品文档.高斯列主消元法:#include#include#define N 10int main() int i=0,j=0,r=0,n=0,k=1; float aNN,xN,max=-1,sum1=0,sum2=0;/输入矩阵维数 printf(Please enter n 0n11:n); scanf(%d,&n);/输入矩阵元素,最后一列为b for(i=1;i=n;i+) for(j=1;j=n+1;j+) printf(Please enter a%d%d:n,i,j); scanf(%f,&aij
2、); do for(i=k;imax) max=aik; r=i; if(ark=0) printf(Cannot solve!n); return 0;/换行 for(j=1;j=n+1;j+) sum1=akj; akj=arj; arj=sum1;/将矩阵化为上三角矩阵 for(i=k+1;i=n;i+) for(j=k+1;j=n+1;j+) aij=aij-(aik)*(akj)/(akk); k=k+1; max=-1; while(k=1;k-) for(j=k+1;j=n;j+) sum2=sum2+(akj)*(xj); xk=(akn+1-sum2)/(akk); sum2
3、=0;/输出结果 for(i=1;i=n;i+) printf(x%d=%fn,i,xi); while(1); return 0;LU分解法:#include #include #define N 10 /矩阵大小范围 * 使用已经求出的x,向前计算x(供getx()调用) * float a 矩阵U * float x 方程组解 * int i 解的序号(数组X元素序号) * int n 矩阵大小 * return 公式中需要的和 float getmx(float aNN, float xN, int i, int n) float mx = 0; int r; for(r=i+1; r
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 高斯列主消元法法 LU 解法
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内