第三章　数值数组及其运算.ppt

* 第三章 数值数组及向量化运算 * 3.1 数值计算的特点和地位 符号运算优点：可以对包含变量字符、参数字符和数字的表达式进行推理、运算，并给出符号结果，与高校教科书中的解析表达式接近。 不足：有很多问题无法解，有些问题求解时间过长 现代主流计算机都只能表示、保存、运算和输出有限精度的数字，处理离散的信息。数值计算正好适应了现代计算机的这种特点，它计算速度快，能处理各种复杂的函数关系。但数值计算是以有限精度数字为基本操作元素，只能用有限长度的数据，以有限的精度表现有限时间和范围内的函数关系。 * 例 3.1-1 已知 求 (1) 符号计算解法 syms t; ft=exp(sin(t)); sx=int(ft, t, 0, 4) Warning: Explicit integral could not be found. sx = int(exp(sin(t)),t = 0 .. 4) (2) 数值计算解法 delt=0.01; x=0:delt:4; y=exp(-sin(x)); sx=delt*cumtrapz(y); plot(x,y, r, LineWidth, 6); hold on; plot(x, sx, .b, MarkerSize, 15); plot(x, ones(size(x)), k); * 3.2 引导 数值数组(Numeric Array)和数组运算(Array Operations)是matlat的核心内容. 数组:一组实数或复数排成的长方阵列. 可以是一维的行或列，二维的矩形, 也可以是三维的若干同维矩形的堆叠,甚至更高维. 数组运算：无论在数组上施加什么运算(加减乘除或函数),总认为那种运算对被运算数组中的每个元素平等地实施同样的操作 matlab 设计数组和数组运算的目的： 1) 使程序简单易读,使程序指令更接近教科书上的数学计算公式. 2) 提高程序的向量化程度, 提高计算效率, 节省计算机开销． 数组运算符: .opr?+, -, .* , ./ , .\ , .^ , fun * 3.3 一维数组的创建和寻访 1. 一维数组的创建: (1) 逐个元素输入法： 如：array=[2, pi/2, sqrt(3), 3+5i] 使用场合及特点 使用场合: 数据元素比较少并都已知.如对少量实验数据的处理可用此种方法. x=[1,2,3,4,5]’ 行向量 列向量 * (2) 冒号生成法： array=a: inc: b 向量 a---数组的第一个元素　 inc---采样点之间的间隔, 即步长. 最后一个元素不一定等于b,其大小为b’=a +inc*[(b-a)/inc]; 步长可以省略, 默认为1; inc可以取正数或负数, 但要注意当取正时,要保证ba, 数组最后一个元素不超过b, 取负时ba, 最后一个元素不小于b. 特点: 等差数列 方便对数据之间的间隔(步长)进行控制.但要注意三个数值之间的关系，可能得到空数组．另外要注意生成的数组的元素的个数．如x=a: (b-a)/n :b (ba)得到n+1个元素的数组. a b inc0 a b inc0 * x=1:5?x=[1,2,3,4,5] y=5:-1:1?y=[5, 4, 3, 2, 1] Z=1:2:7?z=[1,3,5,7] u=1:2:8?u=[1,3,5,7] x=31 , 32 , 33 , 34 , 35 ,36 , 37 , 38, 39 输出一个等比数列： a=1:1:9; x=3.^a * A=1: -1 :9? A=Empty matrix: 1-by-0 B=9: 1 : 1?B =Empty matrix: 1-by-0 B=1: 2 :8?B=[1, 3, 5, 7] C=8: -2 :1?C=[8, 6, 4, 2] * (3) 定数线性采样法 在设定的总点数下, 均匀采样生成一维行数组．格式为： 　　　　　x=linspace (a, b, n) a---数组初始值 b---数组最后值 n---数组总个数 特点：方便对数组的元素的总个数进行控制，最后一个元素肯定是设定之值．三个参数之间不必满足一定条件． 增量的计算---inc=(b-a)/(n-1) 等效于: x=a: (b-a)/(n-1) :b * 2. 一维数组的子数组寻访和赋值 x=[1 2 3 4 5]; x(3) x([1 2 3]), x([1 2 3]) x(1:3)=x([1,2,3]) % 1:3=[1,2,3] x(3:end) % end表示最后一个下标 x(3:-1