数字信号处理-扩展题试题.doc

【研讨题目】 扩展题 7．本题研究连续周期信号频谱的近似计算问题。 周期为T0的连续时间周期信号x(t)可用Fourier级数表示为 其中 X(n?0)称为连续时间周期信号x(t)的频谱函数。称为信号的基频（基波），称为信号的谐波。 如果信号x(t)函数表达式已知，则可由积分得出信号的频谱。如果信号x(t)函数表达式未知，或者x(t)函数表达式非常复杂，则很难由积分得信号的频谱。本题的目的就是研究如何利用DFT近似计算连续时间周期信号的频谱。 (1)若在信号x(t)的一个周期T0内抽样N个点，即， T为抽样周期(间隔)，可获得序列x[k] 试分析序列x[k]的DFT与连续时间周期信号x(t)的频谱X(n?0)的关系； (2)由(1)的结论，给出由DFT近似计算周期信号频谱X(n?0)的方案； (3)周期信号x(t)的周期T0=1，x(t)在区间[0，1]的表达式为 x(t)=20t2(1?t)4cos(12?t) (a)试画出信号x(t)在区间[0，1]的波形； (b)若要用10次以内的谐波近似表示x(t)，试给出计算方案，并计算出近似表示的误差。讨论出现误差的原因及减小误差的方法。 （1） 【题目分析】 由“时域抽样，频域周期化”可知，序列x[k]的DFT结果X[m]是连续时间周期信号x(t)的频谱X(n?0)周期化后的结果 【理论推导】 由DFT计算所得结果X[m]与连续周期信号频谱X(n?0)的关系。 因为 又因为 所以 因为 , 所以 所以 令n=m+rN;m=0,1,2,…Ｎ－１，ｒ为整数，则上式为 化简可得 由IDFT定义 所以x[k]的DFT与连续时间信号x(t)的频谱X(nω0)的关系如下： （2）【计算方案】 对连续周期信号以抽样间隔 T= T0 /N 进行抽样，得到 x[k] ；利用 FFT 函数对 x[k]作N点FFT 运算，得到 X[m]；最后求得连续周期信号的频谱为 X(nw0)=X[m]/N。 T0=1; t=0:0.01:T0; xt=20*t.^2.*(1-t).^4.*cos(12*pi*t); plot(t,xt); title(原信号); 【仿真结果】 （3）（b）nω0)发生周期化，再乘以N得到的，由于对于一般的连续周期信号的频谱，当n→±∞时，┃X（nw0）┃→0，所以把X[m]近似为较低的谐波分量，产生的误差很主要在于低次谐波分量与高次谐波分量发生混叠。若要减小误差，则应如果在一个周期内取较多的抽样点N，提高与低频分量发生混叠的高频分量的频率值，高频分量如果足够高，也就是幅度值可以趋近于零，就可以进行有效的近似，其混叠误差也会有效地减少。 同时，在进行信号的重建时，谐波分量选取的次数的多少对误差也有影响。对于一般的周期信号，通常时域是无限的，而我们只能选择次数较低的谐波分量（占有原信号的主要能量）进行信号重建，这是不可避免的会发生泄漏，产生误差。 【仿真程序】 T0=1; t=0:0.01:T0; xt=20*t.^2.*(1-t).^4.*cos(12*pi*t); xm=fft(xt)/length(t); N=10; w0=2*pi; xN=xm(1); for n=1:N; xN=xN+2*real(xm(n+1)*exp(j*n*w0*t)); end; subplot(2,1,1); plot(t,xN); error=(xN-xt); subplot(2,1,2); plot(t,error); title(误差); 【仿真结果】 N=4;谐波为4 N=6谐波为6 N=8谐波为8 N=10谐波为10 N=12谐波为12 N=14谐波为14 N=16谐波为16 （4）【题目分析】 如果周期信号x(t)是带限信号，即信号的最高频率分量为M?0(是正整数)，试确定在一个周期内的最少抽样点N，使得在频谱的计算过程当中不存在混叠误差。与抽样定理给出的结论比较，发表你的看法。 由推出的理论关系式，频谱发生的 是周期化，原信号x(t)带限，最高频率分量为Mω0，考虑到x(t)是实信号，则也会有-Mwo分量，当N=2M时，频谱会发生混叠，所以选取N2M或N≥2M+1，这时再将频谱周期化则不会发生混叠，结果与抽样定理类似。 【结果分析】 取谐波次数的多少对近似计算的影响，由以上仿真结果图像可看出，在10次谐波以内随着谐波次数增加，误差逐渐减小。但是当谐波分量超过14 后，误差该变量很小.误差分析要给出定量的结果，如平均误差，最大误差等。 【自主学习内容】