内容发布更新时间 : 2024/11/6 7:10:33星期一 下面是文章的全部内容请认真阅读。
(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现象”,了解其特点以及产生的原因,了解掌握
了各种典型的连续时间非周期信号的频谱特征。从开始的不了解,到后来通过看书,上网查找资料做出这个实验,我学到了很多东西,虽然花费了不少时间,但是却获得了编程经验。