Matlab滤波信号处理函数
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Matlab滤波信号处理函数
2009-12-04 19:32:22| 分类:matlab方法| 标签:|字号大中小订阅
1 conv
功能:求卷积。
格式:c = conv(a,b)
说明:c = conv(a,b)返回向量a、b的卷积c。
举例:a = [1 2 3]
b = [4 5 6]
c = conv(a,b)
c=
4 13 28 27 18
2 impz
功能:数字滤波器的冲激响应。
格式:[h,t] = impz(b,a)
[h,t] = impz(b,a,n)
[h,t] = impz(b,a,n,Fs)
impz(b,a)
impz(...)
说明:[h,t] = impz(b,a)返回系统(b,a)的冲激响应h和相应的时间轴向量t,b、a分别为系统传递函数的分子和分母系数向量。
[h,t] = impz(b,a,n)返回指定的n点冲激响应
[h,t] = impz(b,a,n,Fs)指定了冲激响应采样点的频率间隔1/Fs。Fs 为相对频率,
缺省值为1。
impz(b,a)和impz(...)绘制冲激响应的图形。
举例:计算线性系统(b,a)的冲激响应,结果见图1.4.1。
b =[0.2 0.1 0.3 0.1 0.2];
a =[1 ?.1 1.55 ?.7 0.3];
impz(b,a,50)
3 zplane
功能:离散系统的零极点图。
格式:zplane(z,p)
zplane(b,a)
说明:zplane(z,p)和zplane(b,a)绘制系统的零极点图,用“o”表示零点,“x”表示
极点。z、p分别为零点和极点向量,b、a分别为系统传递函数的分子和分母
系数向量。
举例:计算线性系统(b,a)的零点和极点,结果见图1.4.2。
b =[0.2 0.1 0.3 0.1 0.2];
a =[1.0 -1.1 1.5 -0.7 0.3];
zplane(b,a)
4 abs
功能:求幅值。
格式:y = abs(x)
说明:y = abs(x)返回复数向量x的幅值向量y。
举例计算一个谐波信号(由几个单一频率信号组成的信号)傅立叶变换的幅值,
结果见图1.4.3。
t = (0:99)/100;
x = sin(2*pi*15*t) + sin(2*pi*40*t);
y = fft(x);
m = abs(y);
f = (0:length(y)-1)'/length(y)*100;
plot(f,m),grid
5 angle
功能:求复数的相角。
格式:p = angle(h)
说明:p = angle(h)返回复数向量h的相位向量p。
举例:计算谐波信号的相位,结果见图1.4.4。
t = (0:99)/100;
x = sin(2*pi*15*t) + sin(2*pi*40*t);
y = fft(x);
p = unwrap(angle(y));
f = (0:length(y)-1)?length(y)*100;
plot(f,p)
6 freqs
功能:模拟滤波器的频率响应(幅频和相频)。
格式:h = freqs(b,a,w)
[h,w] = freqs(b,a)
[h,w] = freqs(b,a,n)
freqs(b,a)
说明:h = freqs(b,a,w)计算系统(b,a)在指定频率点向量w上的频率响应h,b、a分
别为系统传递函数的分子和分母系数向量。
[h,w] = freqs(b,a)自动选取200个频率点计算频率响应h,w 为频率点向量。
[h,w] = freqs(b,a,n)计算n个点的频率响应。
freqs(b,a)绘制频率响应曲线。
magnitude————幅度
phase(degrees)————相位
frequency---------------------------频率
举例:计算模拟滤波器(b,a)的频率响应,结果见图1.4.5。
a = [1 0.4 1];
b = [0.2 0.3 1];
w = logspace(-1,1);
freqs(b,a,w)
7 freqz
功能:数字滤波器的频率响应(幅频和相频)。
格式:[h,w] = freqz(b,a,n)
[h,f] = freqz(b,a,n,Fs)
[h,w] = freqz(b,a,n,'whole')
[h,f] = freqz(b,a,n,'whole',Fs)
h = freqz(b,a,w)
h = freqz(b,a,f,Fs)
freqz(b,a)
说明:[h,w] = freqz(b,a,n)计算数字滤波器(b,a)的n点频率响应h,w为相应的频率点向量。
[h,f] = freqz(b,a,n,Fs)指定采样频率Fs,根据采样定理,h在[0,Fs/2]之间均匀
选取n个频率点,f为频率点向量。
[h,w] = freqz(b,a,n,'whole')在[0,2π]之间选取n个频率点。
[h,f] = freqz(b,a,n,'whole',Fs) 在[0,Fs]之间均匀选取n个频率点。
h = freqz(b,a,w)和h = freqz(b,a,f,Fs)计算指定频率点向量的频率响应。
freqz(b,a)绘制频率响应曲线。(计算并且赋值时不绘曲线)magnitude————幅度