数值分析实验课秋.pptxVIP

引言本课程是随《数值分析》一起开设的实验课程，为《数值分析》理论课程的实践环节，旨在引导学生利用计算机开展数值试验，掌握数值算法和程序设计的基本原理和技能。学生通过实验训练，逐步掌握数值试验的方法和技巧，培养学生应用计算方法解决工工程计算的能力，以期达到初步的科学计算和研究的目的。 课程介绍教师：肖鸣宇 （mailto:myxiao@myxiao@ 开发工具：MATLAB课时：16课时实验项目：4个上课时间地点：待定信息公布处：计算机学院—师资队伍— 肖鸣宇考核内容和标准出勤：10% （无特殊情况，最少需出勤一次）写报告：50% （按照标准实验报告书写，要求打印程序，若程序注释部分清晰可读则可得高分）电子版程序：10% （按要求发送到老师信箱。格式：*_project#.m，其中*为学号，#为实验项目编号。）现场检测：20% （主动提出来检查的得高分）其它：10% （如完成附加要求等）MATLAB基础创建矩阵：A=[5 2 1; 1,3,4];% 同行元素用“ ”或“,”隔开，行与行间用“;”隔开Zeros()Ones()用小矩阵作元素建立大矩阵 B=[A;[1 2 3]];向量操作等：X=[3:-1:-1];% [3:-1:-1]=[3,2,1,0,-1]X=A(2 ; :);% A的第二列元素(a,b)=size(A)% 矩阵A的大小矩阵运算A*B；%乘法A\b；%左除A/B；%右除 MATLAB基础控制语句If expression StatementsEnd For variable=expression StatementsEnd while expression StatementsEnd Switch x case 1 statement1 case statement2 otherwiseEnd关系运算符==~=|~注释 %续行 …MATLAB基础一个例子g=2;for%20k=1:100%20%20%20%20g=1+1/g;endgg=2;g=2;for%20k=1:100%20%20%20%20g=1+1/g;endgfor k=1:100g=2;for%20k=1:100%20%20%20%20g=1+1/g;endg g=1+1/g;g=2;for%20k=1:100%20%20%20%20g=1+1/g;endgEndg=2;for%20k=1:100%20%20%20%20g=1+1/g;endgg项目1：求解线性代数方程组的直接方法 实验内容：利用高斯消元法求解如下方程组要求输出格式：x1=***, x2=***, x3=***.附加分部分要求编写的程序能解一般的线性方程。输入：方程的系数矩阵。输出：如果没有解，输出：No solution；如果没有唯一解，输出：No unique solution；否则输入x1=***, x2=***, x3=***， x4=***，……现场检查内容：按老师要求输入一组方程的系数矩阵，查看输出结果是否正确。一个古老的实例 《九章算术》今有上禾三秉，中禾二秉，下禾一秉，实三十九斗；上禾二秉，中禾三秉，下禾一秉，实三十四斗；上禾一秉，中禾二秉，下禾三秉，实二十六斗。问上、中、下禾实一秉各几何？答曰：上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。刘辉（/view/2328.htm三国后期 ）为九章算术注释求解步骤做恒等变换：将第一行乘以-2/3加入第二行将第一行乘以-1/3加入第三行将第二行乘以-4/5加入第三行可以解得：一般形式用矩阵表示A和b为程序的输入部分。求解步骤求解步骤求解步骤求解步骤求解步骤程序实现中的细节注意三重循环for k = 1 to n-1… 对第k列,消去对角线以下元素 … (通过每行加上第k行的倍数) for i = k+1 to n… 对第k行以下的每一行i for j = k to n… 第k行的倍数加到第 i 行 A(i,j) = A(i,j) - (A(i,k)/A(k,k)) * A(k,j)第k列第k行以下每行k到n的元素