MALTAB滤波相关函数解析

最新推荐文章于 2022-02-24 23:00:49 发布

yundanfengqing_nuc

最新推荐文章于 2022-02-24 23:00:49 发布

阅读量991

点赞数

分类专栏： MATLAB

本文链接：https://blog.csdn.net/yundanfengqing_nuc/article/details/45307791

版权

MATLAB 专栏收录该内容

54 篇文章 16 订阅

订阅专栏

1、filter函数

1）filter：一维数字滤波器

2）使用方法：

y = filter(b,a,X) [y,zf] = filter(b,a,X) [y,zf] = filter(b,a,X,zi) y = filter(b,a,X,zi,dim) [...]= filter(b,a,X,[],dim)

3）使用方法解析：

3.1）y = filter(b,a,X)，输入X为滤波前序列，y为滤波结果序列，b/a提供滤波器系数，b为分子，a为

分母，整个滤波过程通过下面差分方程实现：

a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb) - a(2)*y(n-1) - ... - a(na+1)*y(n-na)

3.2）[y,zf] = filter(b,a,X)，输入X为滤波前序列，y为滤波结果序列，b/a提供滤波器系数，b为分子，

a为分母。

《Simulink与信号处理》并输入Zi指定X的初始状态，Zf为最终状态矢量。

《Simulink与信号处理》filter(B,A,X,[],DIM) 或 filter(B,A,X,Zi,DIM)指定X的维数DIM进行操作。

参考文献：http://wenda.haosou.com/q/1364969317060621

2、psd函数

PSD 是做功率谱密度的函数。

x 是信号;Nfft快速傅里叶变换点数;

fs是采样频率;window是加的窗函数;

noverlap是指没有重叠率(根据采样定理可以算出最小采样频率);

dflag好像是判断前边这个noverlap是否有重叠，有的话就假(不继续)，如果没有重叠为真(继续做)。

找本信号处理的书看看应该就知道了。个人浅见，只为抛砖引玉。

参考文献：http://wenda.haosou.com/q/1378519180076336

3、spectrum函数

spectrum：谱估计函数

使用方法：Hs = spectrum.estmethod(input1,...)，estmethod是求谱的各种方法，包括周期图法等，

此行命令创建一个对象Hs，可以用set或get函数设置或取得该对象的属性，然后可以用psd来估计功率谱密度并作图。

例子：
Fs = 1000;
t = 0:1/Fs:.3;
x=cos(2*pi*t*200)+randn(size(t));
Hs=spectrum.periodogram;
psd(Hs,x,'Fs',Fs)

4、fir1函数

功能:设计标准频率响应的基于窗函数的FIR滤波器。

语法:b=fir1(n，Wn);

b=fir1(n，Wn，‘ftytpe’);

b=fir1(n，Wn，Window);

b=fir1(n，Wn，‘ftype’，Window);

说明:fir1函数可以实现加窗线形相位FIR数字滤波器设计，它可以设计出标准的低通、高通、带通和带阻滤波器。b=fir1(n，Wn)可得到n阶低通，截至频率为Wn的汉明加窗线形相位FIR滤波器，0≤Wn≤1，Wn=1相当于0.5fs。滤波器系数包含在b中，可表示为

当Wn=[W1 W2]时，fir1函数可得到带通滤波器，其通带为W1＜w＜W2。

当ftype=high时，设计高通FIR滤波器;

当ftype=stop时，设计带阻滤波器。

在设计高通和带阻滤波器时，由于对奇次阶的滤波器，其在Nyquist频率处的频率响应为零，不适合构成高通和带阻滤波器。因此fir1函数总是使用阶数为偶数的滤波器，当输入的阶数为奇数时，fir1函数会自动将阶数加1。b=fir1(n，Wn，Window)利用参数Window来指定滤波器采用的窗函数类型。其默认值为汉明窗。

b=fir1(n，Wn，‘ftype’，Window)可利用ftype和Window参数，设计各种滤波器。

5、fir2函数

功能:设计任意频率响应的基于窗函数的FIR滤波器。

语法:b=fir2(n，f，m);

b=fir2(n，f，m，Window);

b=fir2(n，f，m，npt);

b=fir2(n，f，m，npt，window);

b=fir2(n，f，m，npt，lap);

b=fir2(n，f，m，npt，lap，Window);

说明:fir2函数可以用于设计有任意频率响应的加窗FIR滤波器，对标准的低通、带通、高通和带阻滤波器的设计可使用fir1函数。

b=fir2(n，f，m)可设计出一个n阶的FIR滤波器，其滤波器的频率特性由参数f和m决定。参数f为频率点矢量，且f∈[0，1]，f=1对应于0.5fs。矢量f按升序排列，且第一个元素必须是0，最后一个必须为1，并可以包含重复的频率点。矢量m中包含了与f相对应的期望得到的滤波器的幅度。

b=fir2(n，f，m，Window)中用参数Window来指定使用的窗函数类型，默认值为汉明窗。