数字信号处理_上机实验.doc

下载文档

51
0
约4.93千字
约 9页
2017-08-28 发布于河南
举报
版权申诉
保障服务

数字信号处理_上机实验.doc

1、本文档共9页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

8.2.1 FFT频谱分析及应用 clc; N=40;n=0:N-1;t=0.01*n; x=2*cos(4*pi*t)+5*cos(8*pi*t); subplot(2,1,1);stem(n,x,.);title(signal x(n));grid; X=fft(x,N); magx=abs(X(1:1:N/2+1));k=[0:1:N/2];w=2*pi/N*k;y=w/2*pi*100 subplot(2,1,2);plot(w/pi,magx);title(FFT N=128); xlabel(频率(单位:pi));ylabel(|X|);grid N=128;n=0:N-1;t=0.01*n; x=2*cos(4*pi*t)+5*cos(8*pi*t); subplot(2,1,1);stem(t,x,.);title(signal x(n));grid; X=fft(x,N); magx=abs(X(1:1:N/2+1));k=[0:1:N/2];w=2*pi/N*k; subplot(2,1,2);plot(w/pi,magx);title(FFT N=128); xlabel(频率(单位:pi));ylabel(|X|);grid t=0:0.001:0.3; x=sin(2*pi*50*t)+sin(2*pi*120*t); subplot(2,1,1);stem(t,x,.);title(signal x(n));grid; y=x+1.5*randn(1,length(t)); fs=1000; N=512; Y=fft(y,N); P=Y.*conj(Y)/N; p=P(1:1:N/2+1);f=fs*(0:N/2)/N; subplot(2,1,2);plot(f,p);title(FFT N=512);grid; xlabel(频率（单位：Hz）);ylabel(功率谱);grid; clc; fs=32000; s=wavread(s8.wav); N=1024; S=fft(s,N); magx=abs(S(1:1:N/2+1));f=fs*(0:N/2)/N; plot(f,magx);title(FFT); 8.2.2 IIR数字滤波器 clear;clc;close all; Rp=1.2;Rs=20;T=0.001;fp=300;fs=200; wp=2*pi*fp*T;ws=2*pi*fs*T; wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); [n,wn]=cheb1ord(wp1,ws1,Rp,Rs,s);[b,a]=cheby1(n,Rp,wn,high,s); [bz,az]=bilinear(b,a,1/T) [db,mag,pha,grd,w]=freqz_m(bz,az); plot(w/pi,db);axis([0,1,-30,2]);grid; bz =0.0691 -0.2073 0.2073 -0.0691 az =1.0000 1.0327 0.9031 0.3174 clear;clc;close all; %脉冲响应不变法 Rp=1;Rs=25;T=0.001;fp=200;fs=300; wp=2*pi*fp;ws=2*pi*fs; [n,wn]=buttord(wp,ws,Rp,Rs,s); [b,a]=butter(n,wn,s); [bz,az]=impinvar(b,a,1/T) [db,mag,pha,grd,w]=freqz_m(bz,az); plot(w/(2*pi*T),db,Color,r);hold on %双线性变换法 wp=2*pi*fp*T;ws=2*pi*fs*T; wp1=(2/T)*tan(wp/2);ws1=(2/T)*tan(ws/2); [n,wn]=buttord(wp1,ws1,Rp,Rs,s); [b,a]=butter(n,wn,s); [bz,az]=bilinear(b,a,1/T) [db,mag,pha,grd,w]=freqz_m(bz,az); plot(w/(2*pi*T),db);axis([0,500,-80,2]);grid;hold off 脉冲响应不变法 bz =-0.0000 0.0002 0.0153 0.0995 0.1444 0.0611 0.0075 0.0002 0.0000 0 az =1.0000 -1.9199 2.5324 -2.2053 1.3868 -0.6309 0.2045 -0.0450 0.0060 -0.0004 双线性变换法 bz =0.0179