LU分解法(C语言)_第1页
LU分解法(C语言)_第2页
LU分解法(C语言)_第3页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

1、LU分解法求解线性方程。#includeVoid solve (float l 100，float u 100，float b ，float x ，int n)int i，j；Float t、s1、S2；float y100；for(I=1)；I=n；I) /*启动初始回滚流程*/ S1=0；for(j=1)；j=1；I-)/*启动第二代流程*/S2=0；for(j=n；地；地。J-)t=-uIj；S2=S2 t * xj；xI=(yIS2)/uII；Void main()float a 100 100，l 100 100，u 100 100，x 100，b100；Int i、j、n、r、k；

2、Float s1、S2；for(I=1)；I=99I )/*将所有阵列设定为0，并将l矩阵的对角值设定为1*/for(j=1)；J=99j)Lij=0，uIj=0；if(j=I)lIj=1；printf( input n : n )；/*输入方程式数目*/scanf(“% d”，n)；printf( input array a : n )；/*读取原始矩阵A*/for(I=1)；I=n；I)for(j=1)；j=n；j)Scanf(%f ，aIj)；printf( input array b : n )；/*读取列矩阵B*/for(I=1)；I=n；I)scanf(“% f”，bI)；for(

3、r=1)；r=n；R )/*解决方案矩阵l和U*/for(I=r)；I=n；I)S1=0；for(k=1；k=r-1；k)S1=S1 lrk* ukI；urI=arI-S1；for(I=r 1)；I=n；I) S2=0；for(k=1；k=r-1；k)S2=S2 lIk* ukr；lIr=(aIr-S2)/urr；printf( array l : n )；/*输出矩阵L*/for(I=1)；I=n；I)for(j=1)；j=n；j)Printf(%7.3f ，lIj)；printf(“ n”)；printf( array u : n )；/*输出矩阵U*/for(I=1)；I=n；I)for(j=1)；j=n；j)Printf(%7.3f ，uIj)；printf(“ n”)；Solve(l、u、b、x、n)；Printf(使用 :n解决 n )；for(I=1)；I=n；I)Printf(x%d=%fn ，I，xI)；运行结果：Input n:3Input array A:2 2 34 7 7-2 4 5Input array B:3 1 -7阵列l :1.000 0.000 0.0002.000 1.000 0.000-1.000 2.000 1.000阵列u:2.000 2.000 3.0000.000 3.000 1.0000.

人人文库> 全部分类> 应用文书 > 事务文书

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论