信号与系统实验报告-实验3--周期信号的频谱分析 下载本文

内容发布更新时间 : 2024/5/20 1:31:07星期一 下面是文章的全部内容请认真阅读。

(6)仿照程序3_2,编写程序Q3_6,计算并绘制出原始信号x1(t) 的波形图,用有限项级数合成的y1(t) 的波形图,以及x1(t) 的幅度频谱和相位频谱的谱线图。 程序如下:

clc,clear,close all T=2;

dt=0.00001; t=-3:dt:3;

x=(t+1).*(ut(t+1)-ut(t))-(t-1).*(ut(t)-ut(t-1)); x1=0;

for m=-2:2

x1=x1+(t+1-m*T).*(ut(t+1-m*T)-ut(t-m*T))-(t-1-m*T).*(ut(t-m*t)-ut(t-1-m*t)); end

w0=2*pi/T; N=10; L=2*N+1;

for k=-N:N;

ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt; end

phi=angle(ak); y=0;

for q=1:L;

y=y+ak(q)*exp(j*(q-1-N)*w0*t); end;

subplot(221)

plot(t,x)%plot x

axis([-3 3 -0.2 1.2]); grid on;

title('The original signal x(t)'); subplot(223)

plot(t,y)%Plot y

axis([-3 3 -0.2 1.2]); grid on;

title('The synthesis signal y(t)'); subplot(222);

xlabel('Time i (sec)'); subplot(222);

实验心得:

通过这次实验,了解了连续时间周期信号的傅里叶级数的物理意义,观察了截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因,了解掌握

了各种典型的连续时间非周期信号的频谱特征。从开始的不了解,到后来通过看书,上网查找资料做出这个实验,我学到了很多东西,虽然花费了不少时间,但是却获得了编程经验。