数值计算方法上机实验报告
27页1、-数值计算方法上机实验报告实验目的:复习和稳固数值计算方法的根本数学模型,全面掌握运用计算机进展数值计算的具体过程及相关问题。利用计算机语言独立编写、调试数值计算方法程序,培养学生利用计算机和所学理论知识分析解决实际问题的能力。上机练习任务:利用计算机根本C语言编写并调试一系列数值方法计算通用程序,并能正确计算给定题目,掌握调试技能。 掌握文件使用编程技能,如文件的各类操作,数据格式设计、通用程序运行过程中文件输入输出运行方式设计等。一、 各算法的算法原理及计算机程序框图1. 列主元高斯消去法l 算法原理:高斯消去法是利用现行方程组初等变换中的一种变换,即用一个不为零的数乘一个方程后加只另一个方程,使方程组变成同解的上三角方程组,然后再自下而上对上三角方程组求解。列选住院是当高斯消元到第步时,从列的以下包括的各元素中选出绝对值最大的,然后通过行交换将其交换到的位置上。交换系数矩阵中的两行包括常数项,只相当于两个方程的位置交换了,因此,列选主元不影响求解的结果。l 计算机程序框图如上l 源程序:#define N 200#include stdio.h#include math.hFIL
2、E *fp1,*fp2;void LZ() int n,i,j,k=0,l; double d,t,t1; static double *N,aNN;fp1=fopen(a1.t*t,r); fp2=fopen(b1.t*t,w); fscanf(fp1,%d,&n); for(i=0;in;+i) for(j=0;jfabs(d) /*选主元*/ d=aik;l=i; i+;while(in); if(d=0) printf(n输入矩阵有误!n); else /*换行*/ if(l!=k) for(j=k;j=n;j+) t=alj; alj=akj; akj=t; for(j=k+1;j=n;j+) /*正消*/ akj/=akk; for(i=k+1;in;i+) for(j=k+1;j=n;j+) aij-=aik*akj; k+; while(k=0;i-) /*回代*/ t1=0; for(j=i+1;jn;j+) t1+=aij*j; *i=ain-t1; for(i=0;in;i+) fprintf(fp2,n方程组的根为*%d=%lf,i+1,*i); fclose(
3、fp1); fclose(fp2); main() LZ(); l 具体算例及求解结果:用列选主元法求解以下线性方程组输入3 输出结果:方程组的根为*1=6.0000001 2 -3 8 方程组的根为*2=4.000000 2 1 3 22 方程组的根为*3=2.0000003 2 1 28l 输入变量、输出变量说明:输入变量:系数矩阵元素,常向量元素输出变量:解向量元素2. 杜里特尔分解法解线性方程l 算法原理:求解线性方程组时,当对进展杜里特尔分解,则等价于求解,这时可归结为利用递推计算相继求解两个三角形系数矩阵为三角矩阵方程组,用顺代,由求出,再利用回带,由求出。计算机程序框图:源程序:#include stdio.h#include math.hFILE *fp1,*fp2;void main()int i,j,k,N;double s,A200200,B200,*200,y200;static double L200200,U200200;fp1=fopen(a2.t*t,r);fp2=fopen(b2.t*t,w);fscanf(fp1,%d,&N);for(i=0;iN;
《数值计算方法上机实验报告》由会员新**分享,可在线阅读,更多相关《数值计算方法上机实验报告》请在金锄头文库上搜索。
【最新资料】全国高考广东省数学文试卷及答案【精校版】
中学考试数学基础知识归纳
高中语文人教版选修文章写作与修改文笔出彩演练十二 含解析
四年级工作计划
广东省东莞市寮步信义学校七年级语文下册16社戏导学案无答案新人教版
浙江06年10月高等教育自学考试药用植物学试题
多功能爬架施工方案完整(DOC 17页)
人教版八年级上册道德与法治期末复习知识点【最新整理】
小学教学工作总结集合九篇
《吃水不忘挖井人》教学反思
写给消防官兵的感谢信
物业客服一周计划范本(2篇).doc
2011年菏泽市测绘研究院招聘模拟卷(第17期)
小学校长培训班学习心得体会范文五篇
第6章建筑艺术鉴赏
精选面试时简短的自我介绍合集九篇
湿地公园世博模板支撑方案
揭阳芯片成品测试设备项目招商引资方案
2023消防宣传月心得总结(3篇).doc
关于单位档案查阅的介绍信(实用)
2023-07-16 14页
2023-12-04 10页
2023-03-31 24页
2023-02-06 20页
2023-01-23 14页
2023-06-04 23页
2023-11-27 9页
2023-05-12 66页
2023-08-10 7页
2023-01-16 21页