采集数据,低通滤波,得到数据低频部分
![数据分析2 数据滤波_第1张图片](http://img.e-com-net.com/image/info8/b1ab41f2c19f4a319d568819638a879f.jpg)
clear;
Fs=1000;
f1=2;f3=5;
Wp = 2*pi*f1/Fs; %通带截止频率,f1,f3的值都要小于 Fs/2
Ws = 2*pi*f3/Fs;%阻带截止频率,f1,f3的值都要小于 Fs/2
Rp = 1; %通带内的衰减不超过Rp,这个自定大致定义
Rs = 30;%阻带内的衰减不小于Rs,这个自定大致定义
[n,Wn] = buttord(Wp/pi,Ws/pi,Rp,Rs);%巴特沃斯数字滤波器最小阶数选择函数
[b,a] = butter(n,Wn);%巴特沃斯数字滤波器
[h,w]=freqz(b,a,512,Fs); %计算滤波器的频率响应
figure(1);subplot(311);plot(w,abs(h))%,'LineWidth',1绘制滤波器的幅频响应图
grid on,xlabel('频率(Hz)');title('滤波器的频率响应');
axis([0,20,0,2]);
t=0.001:0.001:2;
load('x.mat');
%x=x+0.005*randn(2,2000); %加随机噪声
x2=0.0002*sin(2*pi*16*t)+0.001*sin(2*pi*10*t);
x(1,:)=x(1,:)+x2;%加随机周期信号
y=filtfilt(b,a,x');
hold on;
subplot(312);hua_fft(x(1,:),1000);
subplot(313);hua_fft(y(:,1),1000);title('滤波后功率谱密度');
figure(2);
plot(t,x,t,y);