数值分析实验课（2008年秋）.pptVIP

Tutorial 1 CSE, CUHK 数值分析实验课（2008年秋） 电子课件一 ——课程介绍及预备知识 引言 本课程是随《数值分析》一起开设的实验课程，为《数值分析》理论课程的实践环节，旨在引导学生利用计算机开展数值试验，掌握数值算法和程序设计的基本原理和技能。学生通过实验训练，逐步掌握数值试验的方法和技巧，培养学生应用计算方法解决工工程计算的能力，以期达到初步的科学计算和研究的目的。  课程介绍 教师：肖鸣宇 （myxiao@uestc.edu.cn 开发工具：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基础 MATLAB基础 项目1：求解线性代数方程组的直接方法 实验内容：利用高斯消元法求解如下方程组要求输出格式：x1=***, x2=***, x3=***. 附加分部分 要求编写的程序能解一般的线性方程。输入：方程的系数矩阵。输出：如果没有解，输出：No solution；如果没有唯一解，输出：No unique solution；否则输入x1=***, x2=***, x3=***， x4=***，…… 现场检查内容：按老师要求输入一组方程的系数矩阵，查看输出结果是否正确。 一个古老的实例 《九章算术》今有上禾三秉，中禾二秉，下禾一秉，实三十九斗；上禾二秉，中禾三秉，下禾一秉，实三十四斗；上禾一秉，中禾二秉，下禾三秉，实二十六斗。问上、中、下禾实一秉各几何？答曰：上禾一秉九斗四分斗之一。中禾一秉四斗四分斗之一。下禾一秉二斗四分斗之三。 求解步骤 一般形式 用矩阵表示 求解步骤 求解步骤 求解步骤 求解步骤 求解步骤 程序实现中的细节 * * 计算机科学与工程学院 示范性软件学院 电子科技大学 关系运算符 == ~= | ~ 注释 % 续行 … 控制语句 If expression Statements End For variable=expression Statements End while expression Statements End Switch x case 1 statement1 case statement2 otherwise End 一个例子 g=2; for k=1:100 g=1+1/g; End g 刘辉（三国后期 ）为九章算术注释 做恒等变换：将第一行乘以-2/3加入第二行将第一行乘以-1/3加入第三行 将第二行乘以-4/5加入第三行 可以解得： A和b为程序的输入部分。