?? appl7_5.m
字號:
%Appl7_5
clear all;
wp=2/25;N=400; %通帶邊界頻率(歸一化頻率)和濾波器階數
dt=0.02; %中國數字地震臺網的采樣間隔為0.02s,采樣頻率為50Hz
b=fir1(N,wp,hamming(N+1)); %設計FIR低通濾波器
figure(1)
[H,f]=freqz(b,1,512,1/dt); %求出頻率特性
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;
load cz02.txt %加載數據
x=cz02'; %輸入信號
y=filtfilt(b,1,x); %產生零相位輸出,采用filtfilt 對輸入信號濾波
t=[0:length(x)-1]*dt;
figure(2)
subplot (2,1,1), plot (t,x),title ('輸入信號')%繪出輸入信號波形
subplot (2,1,2),plot (t,y) %繪出輸出信號波形
title('輸出信號'),xlabel('時間/s')
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -