实验三周期信号的频谱分析实验报告Word格式文档下载.docx
- 文档编号:16803443
- 上传时间:2022-11-26
- 格式:DOCX
- 页数:12
- 大小:96.24KB
实验三周期信号的频谱分析实验报告Word格式文档下载.docx
《实验三周期信号的频谱分析实验报告Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《实验三周期信号的频谱分析实验报告Word格式文档下载.docx(12页珍藏版)》请在冰豆网上搜索。
dt:
4;
%Specify theintervaloftime
w0=0、5*pi;
x1=cos(w0、*t);
x2=cos(3*w0、*t);
x3=cos(5*w0、*t);
N=input('Typeinthe numberoftheharmonic ponentsN=');
x=0;
forq=1:
N;
x=x+(sin(q*(pi/2))、*cos(q*w0*t))/q;
end
subplot(221)
plot(t,x1)%Plotx1
axis([-24-22]);
gridon,
title('
signalcos(w0、*t)'
)
subplot(222)
plot(t,x2)%Plotx2
axis([-24-2 2]);
gridon,
title('
signalcos(3*w0、*t))'
subplot(223)
plot(t,x3)%Plot x3
axis([-24-2 2])
gridon,
title('
signalcos(5*w0、*t))'
执行程序Q3_1所得到得图形如下:
Q3-2给程序Program3_1增加适当得语句,并以Q3_2存盘,使之能够计算例题1中得周期方波信号得傅里叶级数得系数,并绘制出信号得幅度谱与相位谱得谱线图。
通过增加适当得语句修改Program3_1而成得程序Q3_2抄写如下:
%Program3_1
clear, closeall
T =2;
dt=0、00001;
t=-2:
2;
x1=u(t)- u(t-1-dt);
x=0;
form= -1:
1 %Periodicallyextendx1(t)to form a periodicsignal
x=x + u(t-m*T)-u(t-1-m*T-dt);
end
w0=2*pi/T;
N= 10;
% The numberofthe harmonic ponents
L= 2*N+1;
fork= -N:
N;
%EvaluatetheFourierseriescoefficientsak
ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t'
)*dt;
end
phi=angle(ak);
%Evaluatethephaseofak
subplot(211)'
k=-10:
10;
stem(k,abs(ak),'
k'
);
axis([-10,10,0,0、6]);
gridon;
title('
fudupu');
subplot(212);
k = -10:
10
stem(k,angle(ak),'
);
axis([-10,10,-2,2]);
gridon;
titie('
xiangweipu'
xlabel('Frequencyindexx'
执行程序Q3_2得到得图形
Q3-3反复执行程序Program3_2,每次执行该程序时,输入不同得N值,并观察所合成得周期方波信号。
通过观察,您了解得吉伯斯现象得特点就是:
%Program3_3
%This programisused toputethe Fourier seriescoefficientsakof aperiodicsquarewave
clear,closeall
T= 2;
dt=0、00001;
t=-2:
2;
x1=u(t)-u(t-1-dt);
x= 0;
form=-1:
1
x =x+u(t-m*T)-u(t-1-m*T-dt);
% Periodicallyextendx1(t)toformaperiodicsignal
w0 =2*pi/T;
N =input('
Typeinthe number oftheharmonicponentsN=:
'
L =2*N+1;
fork=-N:
1:
ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t')*dt;
end
phi = angle(ak);
y=0;
forq= 1:
L;
%Synthesiztheperiodic signaly(t) fromthefiniteFourier series
y = y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T);
end;
subplot(221),
plot(t,x),title('
Theoriginalsignalx(t)'
), axis([-2,2,-0、2,1、2]),
subplot(223),
plot(t,y),title('Thesynthesissignaly(t)'),axis([-2,2,-0、2,1、2]), xlabel('
Timet'
),
subplot(222)
k=-N:
N;
stem(k,abs(ak),'
k、'
), title('
Theamplitude|ak|ofx(t)'
),axis([-N,N,-0、1,0、6])
subplot(224)
stem(k,phi,'
r、'
), title('The phasephi(k)ofx(t)'
), axis([-N,N,-2,2]), xlabel('
Indexk'
N=1
N=2
通过观察我们了解到:
如果一个周期信号在一个周期有内断点存在,那么,引入得误差将除了产生纹波之外,还将在断点处产生幅度大约为9%得过冲(Overshot),这种现象被称为吉伯斯现象(Gibbsphenomenon)。
即信号在不连续点附近存在一个幅度大约为9%得过冲,且所选谐波次数越多,过冲点越向不连续点靠近。
4、周期信号得傅里叶级数与GIBBS现象
给定如下两个周期信号:
Q3-4仿照程序Program3_1,编写程序Q3_4,以计算x1(t)得傅里叶级数得系数。
程序Q3_4如下:
clc,clear,closeall
T=2;
dt=0、00001;
t=-3:
dt:
3;
x=(t+1)、*(u(t+1)-u(t))-(t-1)、*(u(t)-u(t-1));
x1=0;
form=-2:
2
x1=x1+(t+1-m*T)、*(u(t+1-m*T)-u(t-m*T))-(t-1-m*T)、*(u(t-m*T)-u(t-1-m*T));
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);
plot(t,x1);
axis([-4 401、2]);
grid on;
Thesignal x1(t)'
xlabel('
Timet(sec)'
ylabel('
signal x1(t)'
执行程序Q3_4所得到得x1(t)得傅里叶级数得ak从-10到10共21个系数如下:
Q3-5仿照程序Program3_1,编写程序Q3_5,以计算x2(t)得傅里叶级数得系数(不绘图)。
程序Q3_5如下:
clc,clear,closeall
T=2;
dt=0、00001;
t=-3:
x=u(t+0、2)-u(t-0、2-dt);
x2=0;
form=-1:
x2=x2+u(t+0、2-m*T)-u(t-0、2-m*T)-u(t-0、2-m*t-dt);
w0=2*pi/T;
N=10;
L=2*N+1;
for k=-N:
ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t'
)*dt;
phi=angle(ak);
plot(t,x2);
axis([-2、5 2、50 1、2]);
gridon;
title('Thesignalx2(t)'
xlabel('
Time t(sec)');
ylabel('
signalx2(t)'
执行程序Q3_5所得到得x2(t)得傅里叶级数得ak从-10到10共21个系数如下:
与您手工计算得ak相比较,就是否相同,如有不同,就是何原因造成得?
Q3-6 仿照程序Program3_2,编写程序Q3_6,计算并绘制出原始信号x1(t)得波形图,用有限项级数合成得y1(t) 得波形图,以及x1(t)得幅度频谱与相位频谱得谱线图。
编写程序Q3_6如下:
%Program Q3_6
%ThisprogramisusedtoevaluatetheFourierseriercoefficientsakof aperiodicsquare
clc,clear,close all
dt=0、00001;
t=-3:
dt:
3;
x=(t+1)、*(u(t+1)-u(t))-(t-1)、*(u(t)-u(t-1));
for m=-2:
2%Periodicallyextendx1(t) toformqperiodic signal
x1=x1+(t+1-m*T)、*(u(t+1-m*T)-u(t-m*T))-(t-1-m*T)、*(u(t-m*t)-u(t-1-m*t));
w0=2*pi/T;
N=10;
%thenumberof theharmonicponents
L=2*N+1;
fork=-N:
ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t'
phi=angle(ak);
%Evaluatethephaseofsk
y=0;
forq=1:
L;
%Synthesizthe periodicsignaly(t) fromthefinite Fourierseries
y=y+ak(q)*exp(j*(q-1-N)*w0*t);
end;
subplot(221)
plot(t,x)%plot x
axis([-3 3-0、21、2]);
gridon;
title('
Theoriginalsignalx(t)'
subplot(223)
plot(t,y)%Ploty
axis([-3 3-0、21、2]);
gridon;
title('Thesynthesissignaly(t)'
subplot(222);
xlabel('
Timei(sec)');
subplot(222);
k=-N:
stem(k,abs(ak),'
k');
axis([-N N-0、10、6]);
title('Theamplitudespectrumofx(t)'
subplot(224);
k=-N:
stem(k,phi,'
axis([-NN-22]);
title('Thephasespectrumofx(t)'
xlabel('
Frequencyindexk'
执行程序Q3_6,输入N= 10所得到得图形如下:
反复执行程序Q3_6,输入不同得N值,观察合成得信号波形中,就是否会产生Gibbs现象?
为什么?
假定输入N=10,得到图形如下:
所以不会产生Gibbs现象,即与N值无关。
给定两个时限信号:
实验体会与心得:
在实验得过程中,掌握连续时间周期信号得傅里叶级数得物理意义与分析方法,观察截短傅里叶级数而产生得“Gibbs现象”,了解其特点以及产生得原因,掌握各种典型得连续时间非周期信号得频谱特征。
发现自己在上课时候完全就是一窍不通,可能就是因为自己练得不够所以在下来得学习中,我认为实练永远就是自己要去做得功课,即使自己现在还不会,但我坚信孰能生巧,自己一定能够学好这门科目。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 周期 信号 频谱 分析 报告