[理学]第2章MATLA数值计算功能2-2.ppt

MATLAB 7.X程序设计 第4章 数值计算 主要内容： ①数据分析； ②矩阵分析； ③多项式运算； ④数值插值和拟合； ⑤数值积分和微分 。 4.1 概述 数值计算的特点： 1）数值计算是定义在数值数组的基础上，在计算之前必须定义变量并赋值。 2）数值计算有误差。 3）数值计算的计算速度较快。 4）数值计算除定义了算术运算、关系运算和逻辑运算三种，还提供了许多运算函数。 4.2 数据分析 4.2.1统计分析 1 最大值和最小值max min 2 平均值和中值 mean median 3 求和与求积 sum prod 4 累加和与累乘积cumsum cumprod 5 标准方差 std cov（协方差）corrcoef（相关阵） 6 排序 sort 4.2.2 协方差和相关系数 cov corrcoef 4.2.3 傅立叶变换 fft 4.2 数据分析 数据分析原则： 如果输入是向量，则按整个向量进行运算。 如果输入是矩阵，则按列进行运算。 　通常将需要分析的数据按列进行分类，而用行表示同类数据的不同样本。 4.2 数据分析 4.2.1统计分析 1 最大值和最小值 1）求向量的最大和最小值 2）求矩阵A的最大值和最小值 max(A) 返回一个行向量，向量的第i个元素是矩阵A的第i列上的最大值。 [Y,U]=max(A) 返回行向量Y和U，Y向量记录A的每列的最大值，U向量记录每列最大值的行号。 max(A,[],dim) dim取1或2。dim取1时，该函数和max(A)完全相同；dim取2时，该函数返回一个列向量，其第i个元素是A矩阵的第i行上的最大值。 3）同型的向量或矩阵比较 U=max(A,B) A,B是两个同型的向量或矩阵，结果U是与A,B同型的向量或矩阵，U的每个元素等于A,B对应元素的较大者。 例 4-3 分析下列程序的功能。 x=[4 5 6;1 4 8]; y=[1 7 5;4 5 7]; p=max(x,y) ; P 例4-4 已知x=[1 8 4 2; 9 6 2 5; 3 6 7 1]，从不同维方向求出其平均值和中值。 median(x) median(x,1) %按列方向，求数组的中值 例4-5 已知x=[4 5 6; 1 4 8]，分析矩阵x的每行、每列元素的乘积和全部元素的乘积。 sum(x)或sum(x,1) %求数组各列元素的和 sum(x,2) %求数组各行元素的和 sum(sum(x)) %求数组所有元素的和 例4-6 已知a = [1 2 3; 3 9 6; 4 10 8;  4 0 7]，求矩阵a的每行、每列元 素的乘积和全部元素的累加和。 cumsum(a) %求矩阵的各列元素的累加和 cumsum(a,2) %求矩阵的各行元素的累加和 cumprod(a) cumprod(a,1) %求矩阵的各列元素的累乘积 cumprod(a,2) %求矩阵的各行元素的累乘积 4.2.3 傅立叶变换 函数名 含义 fft 离散傅立叶变换 fft2 二维离散傅立叶变换 fftn n维离散傅立叶变换 fftshift 将零延迟移到频谱中心 ifft 离散傅立叶反变换 ifft2 二维离散傅立叶反变换 ifftn n维离散傅立叶反变换 ifftshift 逆fftshift 例 快速傅立叶变换 给定数学函数 x(t)=12sin(2π×10t+π/4)+5cos(2π×40t) 取N=128，试对t从0~1秒采样，用fft作快速傅立叶变换，绘制相应的振幅-频率图。 分析： 在0-1秒时间范围内采样128点，从而可以确定采样周期和采样频率。由于离散傅立叶变换时的下标应是从0到N-1，故在实际应用时下标应该前移1。又考虑到对离散傅立叶变换来说，其振幅| F(k)|是关于N/2对称的，故只须使k从0到N/2即可。 程序如下： N=128; % 采样点数 T=1; % 采样时间终点 t=linspace(0,T,N); % 给出N个采样时间ti(I=1:N) x=12*sin(2*pi*10*t+pi/4)