samp6_8.m

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

%Samp6_8
Fs=1000;    %采样频率
wp=300*2/Fs;ws=200*2/Fs; %根据采样频率将滤波器边界频率进行转换（6-19式）
Rp=1;Rs=20;    %通带波纹和阻带衰减
Nn=128;      %显示滤波器频率特性的数据长度
[N,Wn]=buttord(wp,ws,Rp,Rs);  %求得数字滤波器的最小阶数和截止频率（归一化频率）
[b,a]=butter(N,Wn,'high');  %设计Butterworth高通数字滤波器
figure(1)
[H,f]=freqz(b,a,Nn,Fs);    %用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:127;
dt=1/Fs;t=n*dt;   %时间序列
f1=10;f2=400; %输入信号频率
x=sin(2*pi*f1*t)+0.5*cos(2*pi*f2*t); %输入信号
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')