实验三快速傅里叶变化及其应用.doc

1．高斯序列 clear x=0:15 p=8 q=2或者4或者8 y=exp(-(x-p).^2/q) subplot(1,2,1) stem(x,y) xk=fft(y,16) subplot(1,2,2) stem(x,abs(xk))  分析： P不变，随着q值的增大，时域信号幅值变化缓慢时域幅度对应变大，低频分量变多，频域信号频谱泄露程度减小，幅度最大值不变。 x=0:15 p=8或者13或者14 q=8 y=exp(-(x-p).^2/q) subplot(1,2,1) stem(x,y) xk=fft(y,16) subplot(1,2,2) stem(x,abs(xk)) 分析： 可见，当q不变，随着p的增大，会在时间轴移位对应右移，时域信号幅值不变。当实验中p=13、14时，明显出现泄漏。 结论： X（n）中的参数p为高斯序列的峰值位置，q则表示高斯序列峰的尖锐度，（即峰值边沿的陡峭度）。q值越大，时域图中图象越平缓，序列变化越慢；其幅频特性图中高频分量越少，频谱越窄，越不容易产生混叠。p值越大，序列右移，在规定的窗口内有效值被截断的越多。因为窗口截断会造成窗口泄露，所以我们可以在幅频特性图中看到，随着p值的变大，高频分量会增加。易出现泄露，当p=13时，特别是p=14时，产生了明显的泄露与混叠。 2. 衰减的正弦序列 clear x=0:15 a=0.1 f=0.0625或者0.4375或者0.5625 y=exp(-(a.*x)).*sin(2.*pi.*f.*x) subplot(1,2,1) stem(x,y) xk=fft(y,16) subplot(1,2,2) stem(x,abs(xk))  结论： 图中的幅频特性图：当f=0.0625时，没有产生明显的混叠和泄露，当f=0.4375和f=0.5625时，产生了混叠，是因为不满足奈奎斯特采样定理的缘故。图中后两个序列的时域图：因为0.4375+0.5625=1，关于0.5对称，其幅频特性是完全相同的。 3、三角波和反三角波 三角波序列 x=0:7; y1=zeros(size(x)); y2=zeros(size(x)); y3=zeros(size(x)); N=length(x); for k=1:N if x(k)4x(k)=0; y1(k)=x(k); elseif x(k)=4x(k)8 ; y2(k)=8-x(k); else y3(k)=0; end end y=y1+y2+y3; subplot(1,2,1) stem(x,y) xk=fft(y,N) subplot(1,2,2) stem(x,abs(xk)) 反三角波序列 x=0:7; y1=zeros(size(x)); y2=zeros(size(x)); y3=zeros(size(x)); N=length(x); for k=1:N if x(k)4x(k)=0; y1(k)=4-x(k); elseif x(k)=4x(k)8 ; y2(k)=x(k)-4; else y3(k)=0; end end y=y1+y2+y3; subplot(1,2,1) stem(x,y) xk=fft(y,N) subplot(1,2,2) stem(x,abs(xk)) N=32 三角波序列 clear x=0:31; y1=zeros(size(x)); y2=zeros(size(x)); y3=zeros(size(x)); N=length(x); for k=1:N if x(k)4x(k)=0; y1(k)=x(k); elseif x(k)=4x(k)8 ; y2(k)=8-x(k); else y3(k)=0; end end y=y1+y2+y3; xk=fft(y,N) stem(x,abs(xk)) 反三角波序列 clear x=0:31; y1=zeros(size(x)); y2=zeros(size(x)); y3=zeros(size(x)); N=length(x); for k=1:N if x(k)4x(k)=0; y1(k)=4-x(k); elseif x(k)=4x(k)8 ; y2(k)=x(k)-4; else y3(k)=0; end end y=y1+y2+y3; xk=fft(y,N) stem(x,abs(xk)) 分析： 原来未补零时候，正反三角波的幅频特性相同，因为两者的时域只差一个相位。 补零后，反三角波的高频分量增多，对信号末尾补零加长整数个周期可以对原信号达到细化频谱的作用。反三角波的边沿比较陡峭，因此它的幅频特性曲线中高频分量比较多。这是因为栅栏效应，当N=8时，一些谱线被挡住。 4. 一个连续信号含两个