巴特沃斯低通滤波
%**************************************************************************************** % % 創建兩個信號Mix_Signal_1 和信號 Mix_Signal_2 % %***************************************************************************************Fs = 1000; %采樣率 N = 1000; %采樣點數 n = 0:N-1; t = 0:1/Fs:1-1/Fs; %時間序列 Signal_Original_1 =sin(2*pi*10*t)+sin(2*pi*20*t)+sin(2*pi*30*t); Noise_White_1 = [0.3*randn(1,500), rand(1,500)]; %前500點高斯分部白噪聲,后500點均勻分布白噪聲 Mix_Signal_1 = Signal_Original_1 + Noise_White_1; %構造的混合信號Signal_Original_2= [zeros(1,100),20*ones(1,20),-2*ones(1,30),5*ones(1,80),-5*ones(1,30),9*ones(1,140), -4*ones(1,40), 3*ones(1,220), 12*ones(1,100), 5*ones(1,20), 25*ones(1,30), 7 *ones(1,190)]; Noise_White_2 = 0.5*randn(1,1000); %高斯白噪聲 Mix_Signal_2 = Signal_Original_2 + Noise_White_2; %構造的混合信號%**************************************************************************************** % % 信號Mix_Signal_1 和 Mix_Signal_2 分別作巴特沃斯低通濾波。 % %***************************************************************************************%混合信號 Mix_Signal_1 巴特沃斯低通濾波 figure(1); Wc=2*50/Fs; %截止頻率 50Hz [b,a]=butter(4,Wc); Signal_Filter=filter(b,a,Mix_Signal_1);subplot(4,1,1); %Mix_Signal_1 原始信號 plot(Mix_Signal_1); axis([0,1000,-4,4]); title('原始信號 ');subplot(4,1,2); %Mix_Signal_1 低通濾波濾波后信號 plot(Signal_Filter); axis([0,1000,-4,4]); title('巴特沃斯低通濾波后信號');%混合信號 Mix_Signal_2 巴特沃斯低通濾波 Wc=2*100/Fs; %截止頻率 100Hz [b,a]=butter(4,Wc); Signal_Filter=filter(b,a,Mix_Signal_2);subplot(4,1,3); %Mix_Signal_2 原始信號 plot(Mix_Signal_2); axis([0,1000,-10,30]); title('原始信號 ');subplot(4,1,4); %Mix_Signal_2 低通濾波濾波后信號 plot(Signal_Filter); axis([0,1000,-10,30]); title('巴特沃斯低通濾波后信號');
效果圖:總結
- 上一篇: ADT无法识别飞思卡尔的解决办法
- 下一篇: 无监督学习的简单统一