一. 选用fir1
选用了窗函数规划办法
语法如下:
各项参数含义:
b:回来的FIR滤波器单位脉冲呼应,脉冲呼应为偶对称,长度为n+1
n:滤波器的阶数,需求留意的是,规划出的滤波器长度为n+1
Wn:滤波器的截止频率(归一化)。需求留意的是,wn的取值规模为0《wn《1,1对应为信号采样频率的1/2;
假如wn是单个数值,且ftype参数为low,则表明规划截止频率为wn的低通滤波器,如ftype参数为high,则表明规划截止频率为wn的高通滤波器;
假如wn是由两个数组成的向量[wn1,wn2],ftype为stop,则表明规划带阻滤波器,ftype为bandpass,则表明规划带通滤波器;
假如wn是由多个数组成的向量,则表明依据ftype的值规划多个通带或阻带规模的滤波器,ftype为DC-1,表明规划的第一个频带为通带,ftype为DC-0,表明规划的第一个频带为阻带。
4.window:指定运用的窗函数向量,默以为海明窗(Hamming),最常用的窗函数有汉宁窗(Hanning),海明窗(Hamming),布拉克曼窗(Blackman)和凯塞窗(Kaiser),可通过输入help window查找(当凯塞窗函数的贝塔值取7.856时,凯塞窗具有最好的功能)
5.noscale:指定是否归一化滤波器的起伏
留意:fir1函数智能挑选滤波器的截止频率和阶数,不能挑选滤波器通带,阻带衰减,过渡带宽等参数。
二. 选用kaiserord函数规划
语法如下:[n,wn,beta,filtype]=kaiserord(f,a,dev,fs)
参数含义:
f及fs:f是一个向量,其间的元素为待规划滤波器的过渡带的起始点和完毕点。假如没有fs参数,f中元素的取值规模为0~1,即相关于采样频率一半的归一化频率;假如有fs参数,则fs为信号采样频率,f中元素即为实践的截止频率。如,规划滤波器的过渡带宽为1000~1200hz,2000~2100hz,信号采样频率为8000hz;如没有设置fs参数,则f=[0.25 0.3 0.5 0.525];如设置fs为8000,则f=[1000 1200 2000 2100]
a:指定这些频率段的抱负起伏值。通带为1,阻带为0.a与f相对应。以上比如,应该设置为a=[1 0 1],则表明规划带阻滤波器。
dev:指定通带或阻带内的容许差错。同样是上面的比如,要求通带容许差错为0.01,阻带容许差错为0.02,则dev=[0.01 0.02 0.01]
n:n为滤波器的最小阶数
wn:得到滤波器的截止频率点
beta:核算得到的beta值
ftype:取得滤波器的类型参数
三. 选用fir2函数规划
运用fir2函数,能够完结恣意呼应滤波器的规划,即滤波器的起伏频率呼应在指定的频段规模内有不同的幅值。,如在0~0.1的抱负幅值为1,0.2~0.4的幅值为0.5等。
语法如下:
参数含义:
n及b:滤波器的阶数。留意,依据FIR滤波器的结构特色,当规划的滤波器在归一化频率为1处的起伏值部位0时,n不能为奇数。
f及m:f的取值为0~1,对应于滤波器的归一化频率,m用于设置对应频段规模内的抱负幅值。如要求的滤波器在0~0.125内幅值为1,0.125~0.25为0.5,0.25~0.5为0.25,0.5~1为0.125,则f=[0 0.125 0.125 0.25 0.25 0.5 0.5 1],m=[1 1 0.5 0.5 0.25 0.25 0.125 0.125]
windom:指定窗函数的品种,其长度为滤波器长度n+1,默以为海明窗。
npt:指定在对起伏呼应进行插值的插值点个数,其默以为512
lap:指定对起伏呼应进行插值时,关于不接连点转变成接连时的点数,其默认值是25
差异:
1.stem和plot
stem是离散函数,plot是接连函数
常用的信号处理函数:
1. 随机信号函数
随机信号函数:Rand(1,N)和randn(1,N)
2. 方波信号函数
Square(T)和Square(T,DUTY),其实DUTY为正值区域在一个周期的份额,如DUTY=50时,发生方波。
3. 锯齿波信号函数
Sawtooth(T)和sawtooth(T,WIDTH)
4. 正弦波信号函数
Sin(2*pi*f*t)
常用信号剖析级处理函数
1. 滤波函数
Filter(b,a,x):分子项系数顺次从大到小排列成一个行矩阵b,分母项系数顺次从大到小排列成一个行矩阵a,x为输入信号序列。
2. 单位抽样呼应函数
Impz(b,a,p),其间a,b同上,p为核算的数据点数
3. 频率抽样呼应函数
Freqz(b,a,n,Fs): 其间a,b同上,n为在[0 Fs/2]规模内核算的频率点数量,Fs为采样频率。
4. 零极点增益函数
Root函数核算体系的零极点,zplane函数画出体系的零极点。