《MATLAB 7 X程序设计语言》课件第5章.ppt

图5.17对相角进行纠正其中在l的内部在l的边缘上，采用相同形式的三次外插公式。对于多变量函数U(x，y，z，…)，del2(U)逼近于其中，N为函数U中的变量数。在L=del2(U)中，U为矩形阵列，则L是的离散逼近，矩阵L与U同大小，每个元素为U的元素与相邻四元素均值的差。当U为多维阵列时，L=del2(U)为的逼近。在L=del2(U，h)中，标量h用于指定每个方向的间隔，缺省时h=l；L=del2(U，hx，hy)用于U为矩阵时，hx，hy为指定这两个方向上的间隔；L=del2(U，hx，hy，hz，…)用于U为多维阵列时，hx，hy，hz，…为指定各个方向上的间隔。5.8.3相关1．corrcoef功能：计算相关系数。格式：R=corrcoef(X)R=corrcoef(x,y)[R,P]=corrcoef(…)[R,P,RLO,RUP]=corrcoef(…)[…]=corrcoef(…,param1,val1,param2,val2,…)说明：R=corrcoef(X)可计算出输入矩阵X的相关系数矩阵，其中X的行为观察样本，列为变量。矩阵R=corrcoef(X)与协方差矩阵C=cov(X)的关系为R=corrcoef(X)是协方差函数的零阶滞后，即是xcov(x，coeff)封装成方形阵列的零阶滞后。在R=corrcoef(x，y)中，当x，y为列向量时，则它等同于corrcoef([xy])。[R,P]=corrcoef(…)还可以得到矩阵P，当P(i,j)较小(比如小于0.05)时，则说明该函数得到的R(i,j)是有意义的。[R,P,RLO,RUP]=corrcoef(…)还可以得到矩阵RLO和RUP，它们分别是相关系数95％可信区间的下界和上界值。[…]=corrcoef(…,param1,val1,param2,val2,…)可以指定附加的参量及其数值，这里可以指定的参量包括●?alpha：0～1之间的值，用于指定可信水平：100*(1-alpha)%，缺省值为0.05，即指定95％的可信区间。●?rows：其值可以取all(缺省值)，表示使用所有的行；取complete表示使用没有NaN值的行；取pairwise表示计算R(i,j)时将没有NaN值的行应用于列i或j。2．cov功能：计算协方差阵。格式：C=cov(X)C=cov(x，y)说明：在C=cov(X)中，当X为向量时，C为向量元素的方差；当X为矩阵时，cov(X)为协方差矩阵，其中X的行为观测样本，列为变量，diag(cov(X))为由每一列的方差构成的向量，sqrt(diag(cov(X)))为标准差向量。在C=cov(x，y)中，当x，y为列向量时，则它等效于C=cov([xy])。注意协方差的定义为cov(x1,x2)=E[(x1-1)(x2-2)]其中i=Exi，E表示求数学期望。5.8.4滤波1．filter功能：利用无限冲激响应(IIR)或有限冲激响应(FIR)滤波器对数据进行滤波。格式：y=filter(b,a,X) y=filter(b,a,X,zi,dim)[y,zf]=filter(b,a,X) […]=filter(b,a,X,[],dim)[y,zf]=filter(b,a,X,zi)说明：filter函数可对实的和复的数据序列进行滤波，其滤波器采用标准差分方程的直接Ⅱ型实现。在y=filter(b，a，X)中，由b(s)/a(s)构成滤波器，对输入数据向量X中的数据进行滤波。当X为矩阵时，滤波器对X中的数据按列进行滤波；当X为多维阵列时，滤波器对X中的数据沿第一个非单点维进行滤波。[y,zf]=filter(b，a，X)还可以得到滤波器延迟的最终条件zf，实际上zf为max(size(a)，size(b))。[y,zf]=filter(b,a,x,zi)可输入初始条件zi,向量zi的长度必须是max(1ength(a),length(b))-1。y=filter(b,a,X，zi，dim)和[…]=filter(b，a，X，[]，dim)可指定滤波的维数dim。2．filter2功能：二维数字滤波。格式：Y=filter2(b，X)Y=filter2(b，X，