连续时间信号的谱实验案例.docVIP

实验二 连续时间信号的谱分析 实验目的： 1. 观察周期信号的合成过程，进一步理解信号的傅里叶级数分解特性。 2. 观察和分析典型周期连续时间信号频谱的特性。 3. 观察和分析典型非周期连续时间信号频谱的特性。 实验内容： 1. 已知如图2.7所示周期矩形脉冲信号 图2.7 (a) 求的傅里叶系数及傅里叶级数表达式（笔算）。 (b) 绘出由前N次谐波合成的信号近似波形，回答思考题（p2.1）。 (c) 绘出的频谱，观察参数A、T和τ变化时对频谱波形的影响，回答思考题（p2.2）。 傅里叶级数表达式： 其中N为奇数。 程序代码： 1(b) t=-1.5:0.01:1.5; N=input(N); N20 x=zeros(size(t)); for n=1:2:N x=x+(2/(pi*n))*(-1).^((n+3)/2).*cos(2*pi*n*t); end x=x+1/2; plot(t,x); 1(c) N=10; n1=-N:-1; c1=sin(n1*pi/2)/pi./n1; c0=1/2; n2=1:N; c2=sin(n2*pi/2)/pi./n2; cn=[c1 c0 c2]; n=-N:N; subplot(211); stem(n,abs(cn),filled); subplot(212); stem(n,angle(cn),filled); 实验结果： 1(b) 结果分析： ①从理论上说，将信号的时间函数x（t）用傅里叶级数表示时，理论上需要无限多项才能逼近原函数波形。但是在一定误差下，只需保留若干项就可以近似的表示原函数波形。 本实验中，N取20，即取了10项，这样可以再试验结果的图中看出，其已经相当逼近于方波信号。 ②从图中可见，在不连续点附近，部分和有起伏，其峰值几乎与N值无关。峰值的最大值是不连续点处高度的1.09倍，即超量9%。在不连续点上，级数收敛于x（t）两个相位角，说明周期方波信号的傅里叶系数是一个偶函数，这与其傅里叶系数公式 是一致的。 还有就是，可以看出，其振幅频谱的主峰峰值是很高的，而其两侧的值，是不断减小，并且在以后，振幅值几乎为0，这也就为我们用有限项傅里叶级数来逼近提供了有力的依据。 最后就是根据数值分析，可得，振幅频谱的主峰值为,主峰两侧第一个零点为，主峰宽度为，谱线间隔,从0到间的谱线数为。 思考题： P2.1吉伯斯现象的含义是什么？产生这种现象的原因是什么？ 答：用傅里叶级数去逼近方波脉冲时，在不连续点附近，部分和有起伏，其峰值几乎与N值无关。峰值的最大值是不连续点处高度的1.09倍，即超量9%。在不连续点上，级数收敛于x（t）的变化，其频谱结构将如何变化？ 答：特点： （1）周期性信号的频谱是离散的。 （2）理论上，周期信号的谐波含量是无限多，其频谱包括无限多条谱线。但是高次谐波虽然有时起伏，但总的趋势是逐渐减小的。 随着周期矩形脉冲信号增大，频谱的主峰高度增大，反之，则减小。 若T不变，随着增大，有效带宽将减小，反之，将增大。 若不变，随着T的增大，有效带宽不变，但里面所包含的谱线数将减小。 2. 已知如图2.8所示矩形脉冲信号 图2.8 (a) 绘出的频谱，观察参数A和τ变化时对频谱波形的影响，回答思考题（p2.3）。 (b) 绘出的频谱，观察其频谱波形，回答思考题（p2.3）。 程序代码： 2(a) clear syms t x=sym(heaviside(t+1)-heaviside(t-1)); X=fourier(x) X = exp(i*w)*(pi*dirac(w) - i/w) - (pi*dirac(w) - i/w)/exp(i*w) w=-20:0.01:20; X=exp(i*w).*(pi.*dirac(w) - i./w) - (pi.*dirac(w) - i./w)./exp(i*w); plot(w,X); 2(b) clear syms t x=sym(heaviside(t)-heaviside(t-2)); X=fourier(x) X = pi*dirac(w) - (pi*dirac(w) - i/w)/exp(2*i*w) - i/w w=-20:0.01:20; X=pi.*dirac(w) - (pi.*dirac(w) - i./w)./exp(2*i*w) - i./w; subplot(211); plot(w,abs(X)); title(Magnitude of X); subpl