參與評論

評論

登錄后參與評論

全部評論(1條)

• 

  love_xue_58 2017-05-16 18:34:06

  樓主你好 濾波要先設(shè)計濾波器，b，a代表濾波器設(shè)計的參數(shù)。一般可以使用巴特沃斯濾波器。 例如低通濾波器的設(shè)計代碼為: data = importdata('t2.txt'); Time = data(:，1); SA = data(:，2); dtt=diff(Time); dt=sum(dtt)/length(dtt); Fs=round(1/dt);%采樣頻率 Wp = 5/(Fs/2); %通帶截止頻率，這個自定大致定義 Ws = 10/(Fs/2);%阻帶截止頻率，這個自定大致定義 Rp = 2; %通帶內(nèi)的衰減不超過Rp，這個自定大致定義 Rs = 40;%阻帶內(nèi)的衰減不小于Rs，這個自定大致定義 [n，Wn] = buttord(Wp，Ws，Rp，Rs);%巴特沃斯數(shù)字濾波器Z小階數(shù)選擇函數(shù) [b，a] = butter(n，Wn);%巴特沃斯數(shù)字濾波器 [h，w]=freqz(b，a，512，F(xiàn)s); %計算濾波器的頻率響應(yīng) plot(w，abs(h))%，'LineWidth'，1繪制濾波器的幅頻響應(yīng)圖 %************************************************************************** %對輸入的信號進行濾波 RollAf=filtfilt(b，a，RollA);%filtfilt這個函數(shù)是0相位濾波，沒有偏移。filter有偏移。 ************************************************************************** %% 濾波結(jié)果繪圖 figure subplot(2，2，1) H=plot(Time，RollA，Time，RollAf，'r--');%，'linewidth'，3 set(H(2)，'linewidth'，2) 其中;Wp，Ws是一元向量時，則設(shè)計的是低通或高通濾波器，若Wp，Ws是二元向量，則設(shè)計帶通或帶阻濾波器。 可以使用fft函數(shù)。 希望對你有幫助，

  贊(15)

  回復(fù)(0)

  評論

  評論

  登錄后參與評論