数值分析实验讲义.doc

《数值分析》实验讲义 卜兵 数值分析实验一 Matlab基本操作 1. 学习目的 1) 熟悉Matlab的运行环境及各种窗口 2) 掌握Matlab的矩阵变量类型,矩阵输入和矩阵的基本运算 3) 掌握命令及函数文件的作用及区别,并编写简单的M文件 4) 能熟练的向查寻目录中添加新目录,掌握常用的Matlab系统命令 5) 掌握Matlab的控制语句 6) 熟悉数组运算 7) Matlab图形处理功能 8) Matlab程序初步设计 2. 学习内容 (1) Matlab启动与环境设置 (2) Matlab基本运算操作 (3) Matlab的文件 例1: 编写命令文件 demo1 完成以下操作 建数组 a=[1,2,3,...,20],b=[1,3,5,...,39],并求 a, b内积 操作 1) 主窗口点击新建按钮 2) 在弹出的文本编辑窗口添加 a=1:20 b=1:2:39 sum=a*b 3) 单击保存按钮将文件命名为 demo1保存在例1新建文件夹中 4) 在Command Window中输入 demo1并回车 (4) 数组运算(相同类型的运算) 1) ’:’引用 *A (:,n) 矩阵A 的n列所有元素 A=rand(4,5); A (:,3)=(1:4)’ %引用的为一列向量 *A (m,:) 矩阵A 的m行所有元素 A (4,:)=2:6 *A (:) 矩阵A 所有元素 A (:) 2) 维 *reshape(X,M,N,P,..)将已知矩阵X为 M*N*P..矩阵 a=1:12; b=reshape(a,2,6) *用’:’引用 a=zeros(3,4); a(:)=1:12 %Matlab矩阵元素按列存储 a(4) a(1,2) 3) ‘.’运算 同类型矩阵元素对应元素运算 * “.*” ,“./”与”.\’运算 a=[1 2 3;2 3 4;3 4 5]; b=[1 1 1;2 2 2;3 3 3]; a.*b %a,b对应元素相乘 a*b %a,b矩阵相乘 a.\b %a对应元素做分母 a./b %b对应元素做分母 * “.^”与^ b=[1 1 1;2 2 2;3 3 3]; b^3 b.^3 b*b*b %等于b^3 例 :编写函数文件demo3实现 sgn函数功能 操作:1)新建 M 文件,并编辑如下 function val=demo3(x) if x0 val=1; elseif x0 val=-1; else