资源描述：
《数字信号处理-matlab函数》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、这些是我在学习《数字信号处理》时整理的一些函数，这样在用时可以很方便的查找，做matlab时还可以直接将函数copy过去。希望大家都能够把自己的一些好的东西放到网上和大家一起分享，不要总是做“拿来主义者”Aabs(Xk32)：求模值angle(Xk32)：求相位axis([a,b,c,d]):(a,b)定义域和(c,d)值域范围B[Z,P,K]=buttap(N):该格式用于计算N阶巴特沃斯归一化（3dB截止频率）模拟低通原型滤波器系统函数的零极点和增益因子。返回长度为N的列向量Z和P，分别给出N个零点和极点的位置，K表示滤波器增益。得到的系统函数为如下形式：式子中，Z(k)和P(k)分别为Z
2、和P的第k个元素。如果要从计算得到的零、极点得到系统函数的分子和分母多项式系数向量B和A，可以调用结构转换函数[B,A]=zp2tf(Z,P,K).[N,wc]=buttord(wp,ws,Rp,As):该格式用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc。调用参数wp和ws分别为数字滤波器的通带边界频率和阻带边界频率的归一化值，要求0<=wp<=1,0<=ws<=1，1表示数字频率（对应模拟频率，为采样频率）。Rp和As分别为通带最大衰减和阻带最小衰减（dB）.当ws<=wp时为高通滤波器；当wp和ws为二元矢量时为带通或带阻滤波器，这时wc也是二元向量。N和wc作为butter函数
3、的调用参量。[N,wc]=buttord(wp,ws,Rp,As,’s’):该格式用于计算巴特沃斯模拟滤波器的阶数N和3dB截止频率wc、wp、ws和wc是实际模拟角频率（rad/s）。其他参数与[N,wc]=buttord(wp,ws,Rp,As)相同。（同样的设计指标用巴特沃斯数字的和模拟的调用公式计算的阶数、wc可能会不一样，但都是正确的）[B,A]=butter(N,wc,’ftype’):计算N阶巴特沃斯数字滤波器系统函数分子和分母多项式的系数向量B和A。调用参数N和wc分别为巴特沃斯数字滤波器的阶数和3dB截止频率的归一化值（关于归一化），一般按格式[N,wc]=buttord(w
4、p,ws,Rp,As)调用函数buttord计算N和wc。由系数向量B和A可以写出数字滤波器系统函数：式子中，B(k)和A(k)分别为向量B和A的第k个元素[B,A]=butter(N,wc,’ftype’,’s’):计算N阶巴特沃斯模拟滤波器系统函数分子和分母多项式的系数向量B和A。调用参数N和wc分别为巴特沃斯模拟滤波器的阶数和3dB截止频率（实际频率）。由系数向量B和A可以写出模拟滤波器系统函数：Ftype=high时，设计3dB截止频率为wc的高通滤波器，缺省时默认为低通滤波器。Ftype=stop时，设计3dB截止频率为wc的带阻滤波器，缺省时默认为带通滤波器，当然此时的wc和ws为
5、二元向量。Cceil(A)：朝正无穷方向舍入[z,p,k]=cheb1ap(N,Rp):切比雪夫1型滤波器，其中Rp,As,N,z,p,k,B,A均和巴特沃斯相同，不过参数wpo是切比雪夫1型滤波器的通带截止功率，而不是3dB的截止频率。[N,wpo]=cheb1ord(wp,ws,Rp,As):参见巴特沃斯[N,wpo]=cheb1ord(wp,ws,Rp,As,’s’):参见巴特沃斯[B,A]=cheby1(N,Rp,wpo,’ftype’):参见巴特沃斯[B,A]=cheby1(N,Rp,wpo,’ftype’,’s’):参见巴特沃斯[z,p,G]=cheb2ap(N,Rs):切比雪夫2
6、型滤波器,Rs为阻带最小衰减，wso为阻带截止频率[N,wso]=cheb2ord(wp,ws,Rp,As):参见巴特沃斯[N,wso]=cheb2ord(wp,ws,Rp,As,’s’):参见巴特沃斯[B,A]=cheby2(N,Rp,wso,’ftype’):参见巴特沃斯[B,A]=cheby2(N,Rp,wso,’ftype’,’s’):参见巴特沃斯conv(xn,yn):卷积Ddisp:显示函数,例如disp(‘welcome’);Eexp(x).:e的多少次方，注意还有一个点不要少[z,p,k]=ellipap(N,Rp,As):椭圆滤波器,wpo为通带边界频率[N,wpo]=ell
7、ipord(wp,ws,Rp,As):参见巴特沃斯[N,wpo]=ellipord(wp,ws,Rp,As,’s’):参见巴特沃斯[B,A]=ellip(N,Rp,As,wpo,’ftype’):参见巴特沃斯[B,A]=ellip(N,Rp,As,wpo,’ftype’,’s’):参见巴特沃斯ezplot:ezmesh:Ffft(xn)：默认长度为length（xn）的离散傅里叶变换fft(xn,
