网站建设费用 多少,天眼免费查询个人,支持wordpress的主机,网页制作总结心得1、用脉冲响应不变法设计一个Butterworth低通数字滤波器#xff0c;通带截止频率为0.4π #xff0c;通带波纹Rp小于3dB#xff0c;阻带边界频率为0.6π#xff0c;阻带衰减大于15dB#xff0c;采样频率Fs10000Hz。假设一个信号
其中f11000Hz,f24000Hz。试将原信号与通…
1、用脉冲响应不变法设计一个Butterworth低通数字滤波器通带截止频率为0.4π 通带波纹Rp小于3dB阻带边界频率为0.6π阻带衰减大于15dB采样频率Fs10000Hz。假设一个信号
其中f11000Hz,f24000Hz。试将原信号与通过该滤波器的输出信号进行比较。 用经典法设计wp0.4*pi;ws0.6*pi;Rp3;Rs15;
T0.0001;Nn128; %采样间隔
Wpwp/T;Wsws/T; %得到模拟滤波器的频率—采用脉冲响应不变法的频率转换形式
[N,Wn]buttord(Wp,Ws,Rp,Rs,s); %计算模拟滤波器的最小阶数
[z,p,k]buttap(N); %设计低通原型数字滤波器
[Bap,Aap]zp2tf(z,p,k); %零点极点增益形式转换为传递函数形式
[b,a]lp2lp(Bap,Aap,Wn); %低通滤波器频率转换
%[b,a]butter(N,Wn,’s’);
[bz,az]impinvar(b,a,1/T); %脉冲响应不变法设计数字滤波器传递函数
figure(1)
[H,f]freqz(bz,az,Nn,1/T); %输出幅频响应和相频响应
subplot(2,1,1),plot(f,20*log10(abs(H)));
xlabel(频率/Hz);ylabel(振幅/dB);grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel(频率/Hz);ylabel(相位/^o);grid on;
figure(2)
f15;f230; %输入信号含有的频率
N100; %数据点数
n0:N-1;tn*T; %时间序列
xsin(2*pi*f1*t)0.5*cos(2*pi*f2*t); %输入信号
subplot(2,1,1),plot(t,x),title(输入信号)
yfiltfilt(bz,az,x); %对信号进行滤波
y1filter(bz,az,x);
subplot(2,1,2),plot(t,y,t,y1,:),
title(输出信号),xlabel(时间/s)
legend( filtfilt , filter)2、设计一个Butterworth高通数字滤波器通带边界频率为300Hz,阻带边界频率为200Hz通带波纹小于1dB阻带衰减大于20dB,采样频率为1000Hz。试绘出其频率特性图并编制一测试该滤波器满足性能的例子。用完全设计法Fs1000;%采样频率
Wp300*2/Fs;
Ws200*2/Fs;
Rp1;
Rs20;Nn128;
[N,wn]buttord(Wp,Ws,Rp,Rs);
[b,a]butter(N,wn,high) %N为滤波器的阶数
figure(1)
[H,f]freqz(b,a,Nn,Fs); %求得滤波器的频率特性
subplot(2,1,1),plot(f,20*log10(abs(H)));
xlabel(频率/Hz);
ylabel(振幅/dB);
grid on;
subplot(2,1,2),plot(f,180/pi*unwrap(angle(H)))
xlabel(频率/Hz);
ylabel(相位/^o);
