Matlab 讲义ch3-1.docVIP

第三章 数值数组及向量化运算 数值数组（Numeric Array）和数组运算（Array Operations）始终是MATLAB的核心内容。本书从第3章起，全部注意力将集中于数值数组及其运算。 本章系统阐述：数组浮点算法的特点；一、二维数值数组的创建和寻访；数组运算和向量化编程；实现数组运算的基本函数；常用标准数组生成函数和数组构作技法；非数NaN、“空”数组概念和应用；关系和逻辑操作。 3.1 数值计算的特点和地位 【例3.1-1】已知，求。 （1）符号计算解法 syms t x ft=t^2*cos(t) %定义函数 sx=int(ft,t,0,x) %符号积分 ft = t^2*cos(t) sx = x^2*sin(x)-2*sin(x)+2*x*cos(x) 说明： 符号计算可以得到和手工数学推导相同的结果； 可以得到精确的解析结果 （2）数值计算解法 dt=0.05; %定义采样间隔 t=0:dt:5; %定义数值积分的计算区域，数值计算只在有限的区间的采样点进行 Ft=t.^2.*cos(t); %定义计算的函数 Sx=dt*cumtrapz(Ft); %数值积分，计算从0开始到每个采样点位置的区间内，Ft曲线 %下的面积，此面积由宽度为t的小梯形面积累加而成； %下面分步看一下计算结果： whos bar(t,Ft) %Ft函数的数值 hold on plot(t,Sx,.r,MarkerSize,12) %数值积分结果 axis(0,6,-25,10) hold off Name Size Bytes Class Ft 1x101 808 double array Sx 1x101 808 double array dt 1x1 8 double array t 1x101 808 double array xqxqxq1234 1x1 1 logical array Grand total is 305 elements using 2433 bytes 说明： 进行数值计算，必须首先确定一组自变量采样点，即确定一个自变量的取值区间，并设定取样间隔； 数值计算的表达式都是在意志的数值采样点进行计算； 数值计算的结果也是离散的； 再将数值计算的结果扩展到更大的计算区域时不一定成立； 直接观察计算结果的离散数值难以看出函数关系，但图形有利于加深对函数的认识。 【例3.1-2】已知，求。 （1）符号计算解法 syms t x ft=exp(-sin(t)) sx=int(ft,t,0,4) ft = exp(-sin(t)) Warning: Explicit integral could not be found. In sym.int at 58 sx = int(exp(-sin(t)),t = 0 .. 4) 说明：符号计算的能力有限，不是在任何情况下都可以得到计算结果； （2）数值计算解法 dt=0.05; t=0:dt:4; Ft=exp(-sin(t)); Sx=dt*cumtrapz(Ft); Sx(end) bar(t,Ft) hold on plot(t,Sx,.r,MarkerSize,15) hold off xlabel(x) legend(Ft,Sx) hold off ans = 3.0632 说明：数值计算的应用非常广泛，只要给定计算区间、采用点和计算表达式，总是可以得到数值计算的结果。 3.2 数值数组的创建和寻访 3.2.1 一维数组的创建 1递增/减型一维数组的创建 主要特征： （1）数组元素值的大小按照递增或者递减的顺序排列，元素之间的“差”是固定的或者“等步长”的。 （2）主要用作函数的自变量，for循环中的自变量等。 创建方法： （1）“冒号”生成法 格式： x = a : inc :b 说明：x----变量名称； a/ b-----起点/终点 inc -----步长（正数或者负数） 例子： (1)i = 0:1:10