西安电子(丁美玉三版)数字信号处理8.pptVIP

第8章　上 机 实 验 2011-9-5 　　title(′(e) h1(n)与R8(n)的卷积y21(n)′)　　subplot(2， 2， 3)； y=′h2(n)′； tstem(h2n， y)； 　　　%调用函数tstem绘图　　title(′(f) 系统单位脉冲响应h2(n)′)　　subplot(2， 2， 4)； y=′y22(n)′； tstem(y22n， y)； 　　title(′(g) h2(n)与R8(n)的卷积y22(n)′)　　%=====================================　　%内容3： 谐振器分析　　un=ones(1， 256)； %产生信号un　　n=0： 255； 　　xsin=sin(0.014*n)+sin(0.4*n) ； %产生正弦信号 8.1.3 实验结果与波形　　实验结果与波形如图8.1.1所示。 　　在数字信号处理的应用中， 只要涉及时域采样或者频域采样， 都必须服从这两个采样理论的要点。 　 对比上面叙述的时域采样原理和频域采样原理， 得到一个有用的结论， 即两个采样理论具有对偶性： “时域采样频谱周期延拓， 频域采样时域信号周期延拓”。 因此把这两部分内容放在一起进行实验。 8.2.2 实验程序清单　　1. 时域采样理论的验证程序清单　　% 时域采样理论验证程序exp2a.m　　Tp=64/1000； %观察时间Tp=64微秒　　%产生M长采样序列x(n)　　% Fs=1000； T=1/Fs； 　　Fs=1000； T=1/Fs； 　　M=Tp*Fs； n=0： M－1； 　　A=444.128； alph=pi*50*2∧0.5； omega=pi*50*2∧0.5； 　　xnt=A*exp(－alph*n*T).*sin(omega*n*T)； 　　Xk=T*fft(xnt， M)； %M点FFT［xnt)］ 　　2. 频域采样理论的验证程序清单　　%频域采样理论验证程序exp2b.m　　M=27； N=32； n=0： M； 　　%产生M长三角波序列x(n)　　xa=0： floor(M/2)； xb= ceil(M/2)-1： -1： 0； 　　xn=［xa， xb］； 　　Xk=fft(xn， 1024)； 　　　%1024点FFT［x(n)］， 用于近似序列x(n)的TF　　X32k=fft(xn， 32)；%32点FFT［x(n)］　　x32n=ifft(X32k)；%32点IFFT［X32(k)］得到x32(n)　　X16k=X32k(1： 2： N)；%隔点抽取X32k得到X16(K) 　　n1=0： N/2－1； 　　subplot(3， 2， 4)； stem(n1， x16n， ′.′)； box on　　title(′(d) 16点IDFT［X_1_6(k)］′)； xlabel(′n′)； 　　　　ylabel(′x_1_6(n)′)； axis(［0， 32， 0， 20］)　　k=0： N－1； 　　subplot(3， 2， 5)； stem(k， abs(X32k)， ′.′)； box on　title(′(e) 32点频域采样′)； xlabel(′k′)； 　　ylabel(′|X_3_2(k)|′)； 　　　　axis(［0， 16， 0， 200］)　　n1=0： N－1； 　　subplot(3， 2， 6)； stem(n1， x32n， ′.′)； box on　　title(′(f) 32点IDFT［X_3_2(k)］′)； xlabel(′n′)；　　ylabel(′x_3_2(n)′)； axis(［0， 32， 0， 20］) 　　对模拟信号进行谱分析时， 首先要按照采样定理将其变成时域离散信号。 如果是模拟周期信号， 也应该选取整数倍周期的长度， 经过采样后形成周期序列， 按照周期序列的谱分析进行。 　　3． 实验内容及步骤　　（1） 对以下序列进行谱分析: 　　　　　　　x1(n)=R4(n)　　　　　　　　n+1 0≤n≤3　　　　　　　　8－n 4≤n≤7　　　　　　　　0 　 其它n 　　　　　　　4－n 0≤n≤3　 　　　　　　n－3 4≤n≤7　　 　　　　　0 　 其它n 　　 选择FFT的变换