数值计算实验报告1(2011.10.08).docVIP

实 验 报 告 实验课程名称 数值计算方法 实验项目名称 Guass消元法 年 级 2009 专 业 数学与应用数学 学生姓名 学 号 理 学 院 实验时间： 20 年 月 日 学生实验室守则 一、按教学安排准时到实验室上实验课，不得迟到、早退和旷课。 二、进入实验室必须遵守实验室的各项规章制度，保持室内安静、整洁，不准在室内打闹、喧哗、吸烟、吃食物、随地吐痰、乱扔杂物，不准做与实验内容无关的事，非实验用品一律不准带进实验室。 三、实验前必须做好预习（或按要求写好预习报告），未做预习者不准参加实验。 四、实验必须服从教师的安排和指导，认真按规程操作，未经教师允许不得擅自动用仪器设备，特别是与本实验无关的仪器设备和设施，如擅自动用或违反操作规程造成损坏，应按规定赔偿，严重者给予纪律处分。 五、实验中要节约水、电、气及其它消耗材料。 六、细心观察、如实记录实验现象和结果，不得抄袭或随意更改原始记录和数据，不得擅离操作岗位和干扰他人实验。 七、使用易燃、易爆、腐蚀性、有毒有害物品或接触带电设备进行实验，应特别注意规范操作，注意防护；若发生意外，要保持冷静，并及时向指导教师和管理人员报告，不得自行处理。仪器设备发生故障和损坏，应立即停止实验，并主动向指导教师报告，不得自行拆卸查看和拼装。 八、实验完毕，应清理好实验仪器设备并放回原位，清扫好实验现场，经指导教师检查认可并将实验记录交指导教师检查签字后方可离去。 九、无故不参加实验者，应写出检查，提出申请并缴纳相应的实验费及材料消耗费，经批准后，方可补做。 十、自选实验，应事先预约，拟订出实验方案，经实验室主任同意后，在指导教师或实验技术人员的指导下进行。 十一、实验室内一切物品未经允许严禁带出室外，确需带出，必须经过批准并办理手续。 学生所在学院：理学院 专业：数学与应用数学 班级：091 姓 名 学 号 实验组 一 实验时间 指导教师 成 绩 实验项目名称 Guass列主元消元法 实验目的及要求： 在掌握gauss消去法的基础上，要求学生如何增加选主元的技术。 要求在计算机上实现整个算法。 实验（或算法）原理： 根据线性代数知识，当Det（A）！=0时，方程组存在唯一解。对增广矩阵（A，b）施行行初等变换，化A为上三角矩阵A，同时b=b，这时与增广矩阵（A b）相应的线性方程组为上三角形方程组Ax=b。其解为： ， 这是原方程组的解。实现上述求解过程的方法称为Guass消元法。 实验硬件及软件平台：计算机 VC++6.0软件 实验步骤： 输入代码，检查无误后，运行程序 实验内容（包括实验具体内容、算法分析、源代码等等）： #includestdio.h void main() { int i,j,f,n; float a[15][15],b[15],x[15],l,k; printf(Input n:); scanf(%d,n); printf(Input a[][]:); for(i=0;in;i++) for(f=0;fn;f++) scanf(%f,a[i][f]); printf(Input b[]:); for(i=0;in;i++) scanf(%f,b[i]); for(i=0;in-1;i++) { for(j=i+1;jn;j++) { l=-a[j][i]/a[i][i]; b[j]=b[j]+b[i]*l; for(f=i;fn;f++) a[j][f]=a[j][f]+a[i][f]*l; } } x[n-1]=b[n-1]/a[n-1][n-1]; for(i=n-2;i=0;i--) { //x[i]=b[i]; k=0; for(f=i+1;fn;f++) k=k-a[i][f]*x[f]; x[i]=(k+b[i])/a[i][i]; } printf(The answer is ); for(i=0;in;i++) printf(%.4f ,x[i]); } 实验结果与讨论： 指导教师意见： 签名：