MATLAB实验报告6要点.pdf

MATLAB 实验报告 1.实验目的：将三种含有频率分别为60Hz，180Hz，240Hz 的信号 （1）、请选择合适的采样频率和采样点数，用FFT 分 析该信号的频谱并正确显示。（2）、请用切比雪夫I 型低通滤 波器原型设计一个IIR 系统，使其输出只含有60Hz 的信号， 并显示滤波后的信号及其频谱。 2.实验代码： wp=0.2*pi; %滤波器的通带截止频率 ws=0.28*pi; %滤波器的阻带截止频率 Rp=1;As=30; %输入滤波器的通阻带衰减指标 ripple=10^(-Rp/20); Attn=10^(-As/20); %转换为模拟原型滤波器指标 Fs=100;T=1/Fs; Omgp=(2/T)*tan(wp/2); Omgs=(2/T)*tan(ws/2); %模拟原型滤波器计算 [n,Omgc]=buttord(Omgp,Omgs,Rp,As,s) %计算阶数n 和截止 频率 [z0,p0,k0]=buttap(n); %归一化切比雪夫1 型原型设计 ba=k0*real(poly(z0)); %求原型滤波器系数b aa=real(poly(p0)); %求原型滤波器系数a [ba1,aa1]=lp2lp(ba,aa,Omgc);%变换为模拟低通滤波器 %注意,以上4 行求滤波器系数ba1、aa1 的程序，可由下一条 程序替代 %[ba1,aa1]=butter(n,Omgc,s); %直接求模拟滤波器系数 %用双线性变换法计算数字滤波器系数 [bd,ad]=bilinear(ba1,aa1,Fs) %双线性变换 [sos,g]=tf2sos(bd,ad) %由直接型变换为级联型 %求数字系统的频率特性 [H,w]=freqz(bd,ad); dbH=20*log10(abs(H)/max(abs(H))); %化为分贝值 % figure(1); plot(w/pi,dbH);title(幅度响应(dB)); ylabel(dB);xlabel(频率(\pi));axis([0,1,-40,5]); T=0.001; n=0:99; t=T*n; x=2*cos(2*pi*60*t)+0.5*cos(2*pi*180*t)+0.6*cos(2*pi*240*t); y=filter(bd,ad,x); figure(2); subplot(2,1,1),plot(t,x); subplot(2,1,2),plot(t,y); y1=fft(y); figure(3); plot(t,abs(y1)); y2=fft(x); figure(4); plot(t,abs(y2)); 3.实验结果：