[求助]低通滤波器如何与信号函数结合
信号与系统的上机题目,第一步是实现一个简单的含噪声的信号,第二步是做一个低通滤波器,参数值已经给定,我用butter和buttord实现的。第三步要求用第二步所实现的低通滤波器对第一步的信号进行滤波,这步我不知道怎么把二者综合,请求帮助,谢谢!第一步的代码:
clc;
clear;
f1=5;
f2=200;
t = 0 : 0.001 : 1;
w = randn(size(t)); %w是生成的噪声信号
y=0.5 + 3*sin(2*pi*f1*t) + 0.5*cos(2*pi*f2*t); % 原始信号是一个低频正弦波和一个高频正弦波的叠加
x=0.5 + 3*sin(2*pi*f1*t) + 0.5*cos(2*pi*f2*t)+w;%最终的信号叠加了噪声
subplot(3,1,2); %以下是输出三个信号的图形
plot(t,w);
subplot(3,1,1);
plot(t,y);
subplot(3,1,3);
plot(t,x);
第二步的代码:
%设计一个Butterworth低通滤波器,通带边界频率为100Hz,阻带边界频率为150Hz...
%通带纹波不大于2dB,阻带衰减不小于30dB
clc;
clear;
Wp=100*2*pi;
Ws=150*2*pi;
ap=2;
as=-30;
[N,Wn]=buttord(Wp,Ws,ap,as,'s'); %Wp通带边缘频率 Ws阻带边缘频率 ap通带最大衰减 as阻带最小衰减
[b,a]=butter(N,Wn,'s');
figure(1);
freqs(b,a); %拉普拉斯变换频率响应
xlabel('freqs');
figure(2);
freqz(b,a);%Z 变换频率响应
xlabel('freqz');
figure(3);
zplane(b,a);%绘制极零图