实验二连续时间信号的频域分析.docx
- 文档编号:15408972
- 上传时间:2023-07-04
- 格式:DOCX
- 页数:20
- 大小:226.96KB
实验二连续时间信号的频域分析.docx
《实验二连续时间信号的频域分析.docx》由会员分享,可在线阅读,更多相关《实验二连续时间信号的频域分析.docx(20页珍藏版)》请在冰点文库上搜索。
实验二连续时间信号的频域分析
一、实验目的
1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法;
2、观察截短傅里叶级数而产生的“Gibbs现象”,了解其特点以及产生的原因;
3、掌握连续时间傅里叶变换的分析方法及其物理意义;
4、学习利用MATLAB语言编写计算CTFS和CTFT的仿真程序。
基本要求:
掌握并深刻理傅里叶变换的物理意义,掌握信号的傅里叶变换的计算方法,掌握利用MATLAB编程完成相关的傅里叶变换的计算。
二、实验原理及方法
1、连续时间周期信号的傅里叶级数CTFS分析
任何一个周期为T1的正弦周期信号,只要满足狄利克利条件,就可以展开成傅里叶级数。
其中三角傅里叶级数为:
9.1
或:
9.2
其中
,称为信号的基本频率,
分别是信号
的直流分量、余弦分量幅度和正弦分量幅度,
为合并同频率项之后各正弦谐波分量的幅度和初相位,它们都是频率
的函数,绘制出它们与
之间的图像,称为信号的频谱图(简称“频谱”),
-
图像为幅度谱,
-
图像为相位谱。
指数形式的傅里叶级数为:
9.3
其中,
为指数形式的傅里叶级数的系数,按如下公式计算:
9.4
假设谐波项数为N,则上面的和成式为:
9.5
显然,N越大,所选项数越多,有限项级数合成的结果越逼近原信号x(t)。
2、连续时间信号傅里叶变换----CTFT
傅里叶变换在信号分析中具有非常重要的意义,它主要是用来进行信号的频谱分析的。
傅里叶变换和其逆变换定义如下:
9.6
9.7
连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。
任意非周期信号,如果满足狄里克利条件,那么,它可以被看作是由无穷多个不同频率(这些频率都是非常的接近)的周期复指数信号ejt的线性组合构成的,每个频率所对应的周期复指数信号ejt称为频率分量其相对幅度为对应频率的|X(j)|之值,其相位为对应频率的X(j)的相位。
给定一个连续时间非周期信号x(t),它的频谱也是连续且非周期的。
对于连续时间周期信号,也可以用傅里变换来表示其频谱,其特点是,连续时间周期信号的傅里叶变换时有冲激序列构成的,是离散的——这是连续时间周期信号的傅里叶变换的基本特征。
3、连续周期信号的傅里叶级数CTFS的MATLAB实现
3.1傅里叶级数的MATLAB计算
设周期信号x(t)的基本周期为T1,且满足狄里克利条件,则其傅里叶级数的系数可由式9.4计算得到。
式9.4重写如下:
基本频率为:
对周期信号进行分析时,我们往往只需对其在一个周期内进行分析即可,通常选择主周期。
假定x1(t)是x(t)中的主周期,则
计算机不能计算无穷多个系数,所以我们假设需要计算的谐波次数为N,则总的系数个数为2N+1个。
在确定了时间范围和时间变化的步长即T1和dt之后,对某一个系数,上述系数的积分公式可以近似为:
对于全部需要的2N+1个系数,上面的计算可以按照矩阵运算实现。
MATLAB实现系数计算的程序如下:
dt=0.01;
T=2;t=-T/2:
dt:
T/2;w0=2*pi/T;
x1=input(‘Typeintheperiodicsignalx(t)overoneperiodx1(t)=’);
N=input(‘TypeinthenumberN=’);
k=-N:
N;L=2*N+1;
ak=x1*exp(-j*k*w0*t’)*dt/T;
例1:
给定一个周期为T1=2s的连续时间周期方波信号,如图所示,其一个周期内的数学表达式为:
%Program2_1
clear,closeall
T=2;dt=0.00001;t=-2:
dt:
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);
end
w0=2*pi/T;
N=10;
L=2*N+1;
fork=-N:
N;
ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t')*dt;
end
phi=anglel(ak);%Evaluatethephaseofak
执行程序Program2_1后,就完成了信号的傅里叶级数的系数的计算,在命令窗口键入
>>ak
命令窗口就可以显示傅里叶级数的21个系数:
在程序Program2_1中添加相应的计算|ak|和绘图语句,就可以绘制出信号的幅度谱和相位谱的谱线图。
3.2周期信号的合成以及Gibbs现象
从傅里叶级数的合成式(Synthesisequation)
可以看出,用无穷多个不同频率和不同振幅的周期复指数信号可以合成一个周期信号。
然而,我们无法用计算机实现对无穷多个周期复指数信号的合成。
但是,用有限项来合成却是可行的,在实际应用中,多半也就是这么做的。
然而,这样做的一个必然结果,就是引入了误差。
如果一个周期信号在一个周期有内断点存在,那么,引入的误差将除了产生纹波之外,还将在断点处产生幅度大约为9%的过冲(Overshot),这种现象被称为吉伯斯现象(Gibbsphenomenon)。
为了能够观察到合成信号与原信号的不同以及Gibbs现象,我们可以利用前面已经计算出的傅里叶级数的系数,计算出截短的傅里叶级数:
这个计算可用L=2N+1次循环来完成:
其中r作为循环次数,x2在循环之前应先清零。
完成这一计算的MATLAB程序为:
x2=0;L=2*N+1;
forr=1:
L;
x2=x2+ak(r)*exp(j*(r-1-N)*w0*t);
end;
完成了所有的计算之后,就可以用绘图函数:
plot()和stem()将计算结果包括x1,x2,abs(ak)和angle(ak)以图形的形式给出,便于我们观察。
例2:
修改程序Program2_1,使之能够用有限项级数合成例1所给的周期方波信号,并绘制出原始周期信号、合成的周期信号、信号的幅度谱和相位谱。
%Program2_2
%ThisprogramisusedtocomputetheFourierseriescoefficientsakofaperiodicsquarewave
clear,closeall
T=2;dt=0.00001;t=-2:
dt:
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
end
w0=2*pi/T;
N=input('TypeinthenumberoftheharmoniccomponentsN=:
');
L=2*N+1;
fork=-N:
1:
N;
ak(N+1+k)=(1/T)*x1*exp(-j*k*w0*t')*dt;
end
phi=angle(ak);
y=0;
forq=1:
L;%Synthesiztheperiodicsignaly(t)fromthefiniteFourierseries
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('Thephasephi(k)ofx(t)'),axis([-N,N,-2,2]),xlabel('Indexk')
4用MATLAB实现CTFT计算
MATLAB进行傅里叶变换有两种方法,一种利用符号运算的方法计算,另一种是数值计算,本实验要求采用数值计算的方法来进行傅里叶变换的计算。
严格来说,用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件,即信号是时限信号(Timelimitedsignal),也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零,这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。
计算机只能处理有限大小和有限数量的数。
采用数值计算算法的理论依据是:
若信号为时限信号,当时间间隔T取得足够小时,上式可演变为:
上式用MATLAB表示为:
X=x*exp(j*t’*w)*T
其中X为信号x(t)的傅里叶变换,w为频率Ω,T为时间步长。
相应的MATLAB程序:
T=0.01;dw=0.1;%时间和频率变化的步长
t=-10:
T:
10;
w=-4*pi:
dw:
4*pi;
X(j)可以按照下面的矩阵运算来进行:
X=x*exp(-j*t’*)*T;%傅里叶变换
X1=abs(X);%计算幅度谱
phai=angle(X);%计算相位谱
为了使计算结果能够直观地表现出来,还需要用绘图函数将时间信号x(t),信号的幅度谱|X(j)|和相位谱X(j)分别以图形的方式表现出来,并对图形加以适当的标注。
三实验内容和要求
1编写程序,绘制下面的信号的波形图:
其中,0=0.5π,要求将一个图形窗口分割成四个子图,分别绘制cos(0t)、cos(30t)、cos(50t)和x(t)的波形图,给图形加title,网格线和x坐标标签,并且程序能够接受从键盘输入的和式中的项数。
clear,closeall,
dt=0.00001;
t=-2:
dt:
4;
w0=0.5*pi;
x1=cos(w0.*t);
x2=cos(3*w0.*t);
x3=cos(5*w0.*t);
N=input('TypeinthenumberoftheharmoniccomponentsN=:
');
x=0
forq=1:
N
x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q;
end
subplot(221)
plot(t,x1)
axis([-2,4,-2,2]);
gridon,
title('signalcos(w0.*t)')
xlabel('thetimet(sec)')
subplot(222)
plot(t,x2)
axis([-2,4,-2,2]);
gridon,
title('signalcos(3*w0.*t)')
xlabel('thetimet(sec)')
subplot(223)
plot(t,x3)
axis([-2,4,-2,2]);
gridon,
title('signalcos(5*w0.*t)')
subplot(224)
plot(t,x)
axis([-2,4,-2,2]);
gridon,
title('signalx(t)')
xlabel('thetimet(sec)')
2给程序Program2_1增加适当的语句,并存盘,使之能够计算例1中的周期方波信号的傅里叶级数的系数,并绘制出信号的幅度谱和相位谱的谱线图。
clear,closeall
T=2;dt=0.00001;t=-2:
dt:
2;
x1=u(t)-u(t-1-dt);x=0;
form=-1:
1%Periodicallyextendx1(t)toformaperiodicsignal
x=x+u(t-m*T)-u(t-1-m*T-dt);
end
w0=2*pi/T;
N=10;%Thenumberoftheharmoniccomponents
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),'k');
axis([-10,10,-2,2]);
gridon;
title('xiangweipu')
xlabel('indexx');
3反复执行程序Program2_2,每次执行该程序时,输入不同的N值,并观察所合成的周期方波信号。
通过观察,描述总结吉伯斯现象的特点。
N=10:
N=20:
N=40:
N=80:
4给定如下周期信号:
4.1仿照程序Program2_1,编写程序,以计算x1(t)的傅里叶级数的系数。
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=-1:
1
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));
end
w0=2*pi/T;
N=input('TypeinthenumberoftheharmoniccomponentsN=:
');
L=2*N+1;
fork=-N:
N
ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt;
end
y=0;
forq=1:
L
y=y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T);
end;
phi=angle(ak);
>>ak
ak=
1至6列
0.0000+0.0000i0.0025-0.0000i0.0000-0.0000i0.0041-0.0000i0.0000+0.0000i0.0081-0.0000i
7至12列
0.0000-0.0000i0.0225-0.0000i0.0000+0.0000i0.2026+0.0000i0.5000+0.0000i0.2026-0.0000i
13至18列
0.0000-0.0000i0.0225+0.0000i0.0000+0.0000i0.0081+0.0000i0.0000-0.0000i0.0041+0.0000i
19至21列
0.0000+0.0000i0.0025+0.0000i0.0000-0.0000i
4.2仿照程序Program2_2,编写程序,计算并绘制出原始信号x1(t)的波形图,用有限项级数合成的x1(t)的波形图,以及x1(t)的幅度频谱和相位频谱的谱线图。
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=-1:
1
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));
end
w0=2*pi/T;
N=input('TypeinthenumberoftheharmoniccomponentsN=:
');
L=2*N+1;
fork=-N:
N
ak(N+1+k)=(1/T)*x*exp(-j*k*w0*t')*dt;
end
y=0;
forq=1:
L
y=y+ak(q)*exp(j*(-(L-1)/2+q-1)*2*pi*t/T);
end;
%phi=angle(ak);
subplot(221);
plot(t,x1);
axis([-4,4,-0.5,1.2]);
gridon;
title('Thesinalx1(t)');
xlabel('Time(sec)');
ylabel('signalx1(t)');
subplot(222);
k=-N:
N;
stem(k,abs(ak),'k');
axis([-10,10,0,2.0]);
gridon;
title('fudupu')
subplot(223);
k=-N:
N;
stem(k,angle(ak),'k');
axis([-10,10,-2,2]);
gridon;
title('xiangweipu');
xlabel('indexx');
subplot(224);
plot(t,y);
axis([-4,4,-0.5,4]);
gridon;
title('Thesignaly(t)');
xlabel('Timet(sec)');
ylabel('signaly(t)');
以下内容,可视个人实验进度决定是否选做
5连续时间非周期信号的傅里叶变换
给定一个时限信号:
(1)给出信号x1(t)的闭式数学表达式为:
x1(t)
x=(t+2).*(u(t+2)-u(t+1))+u(t+1)-u(t-1)+(-t+2).*(u(t-1)-u(t-2))
(2)编写MATLAB程序,计算并绘制出x1(t)信号的时域波形、幅度谱。
clear,closeall
T=0.01;
dt=0.01;t=-2:
dt:
2;
x1=(t+2).*(u(t+2)-u(t+1))+u(t+1)-u(t-1)+(-t+2).*(u(t-1)-u(t-2));
x=fft(x1);
X=fftshift(x);
X1=abs(X);
X2=angle(X);
X3=real(X);
X4=imag(X);
subplot(221)
plot(t,x1)
gridon,
title('原信号');
subplot(222)
plot(t,X1)
gridon,
title('幅度谱');
subplot(223)
plot(t,X2)
gridon,
title('相位谱');
(3)修改上述程序,计算其傅里叶变换,并分别绘制其傅里叶变换的实部、虚部、幅度频谱和相位频谱的图形。
clear,closeall
T=0.01;
dt=0.01;t=-2:
dt:
2;
w=-4*pi:
0.1:
4*pi;
x1=(t+2).*(u(t+2)-u(t+1))+u(t+1)-u(t-1)+(-t+2).*(u(t-1)-u(t-2));
X=x1*exp(-j*t'*w)*T;
X1=abs(X);
X2=angle(X);
X3=real(X);
X4=imag(X);
subplot(221)
plot(w,X3)
gridon,
title('实部');
subplot(222)
plot(w,X4)
gridon,
title('虚部');
subplot(223)
plot(w,X1)
gridon,
title('幅度谱');
subplot(224)
plot(w,X2)
gridon,
title('相位谱');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 连续 时间 信号 分析