appl6_1.m

MATLAB7.x数字信号处理 光盘内容

%Appl6_1
load ml031025ns.txt   %加载地震波形记录
dt=0.02;   %中国数字地震台网的采样间隔为0.02s,采样频率为50Hz
x=ml031025ns';
wp=0.5*2*dt;  ws=0.2*2*dt; %根据采样频率将滤波器边界频率进行转换
Rp=1;Rs=30;    %通带波纹和阻带衰减
Nn=128;      %显示滤波器频率特性的数据长度
[N,Wn]=buttord(wp,ws,Rp,Rs);  %求得滤波器的最小阶数和截止频率（归一化频率）
[b,a]=butter(N,Wn,'high');  %设计Butterworth高通数字滤波器
figure(1)
[H,f]=freqz(b,a,Nn,1/dt);    %用Nn点求得频率特性
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;
n=0:length(x)-1;
t=n*dt;   %时间序列
figure(2)
subplot(2,1,1),plot(t,x); title('输入信号')  %绘制输入信号
y=filter(b,a,x);   %对输入信号进行滤波
subplot(2,1,2),plot(t,y),title('输出信号') %绘制输出信号
xlabel('时间/s')