连续时间信号的卷积及信号的频域分析报告
- 格式:doc
- 大小:771.00 KB
- 文档页数:8
一、实验目的1. 理解卷积的概念及其物理意义。
2. 掌握卷积运算的原理和方法。
3. 通过实验加深对卷积运算在实际应用中的理解。
二、实验原理1. 卷积的定义:卷积是一种线性运算,它描述了两个信号在时域上的相互作用。
对于两个连续时间信号f(t)和g(t),它们的卷积定义为:F(t) = ∫f(τ)g(t-τ)dτ其中,F(t)是卷积结果,f(τ)是信号f(t)的任意时刻的值,g(t-τ)是信号g(t)在时刻t-τ的值。
2. 卷积的性质:卷积具有交换律、结合律和分配律等性质。
其中,交换律是指f(t)和g(t)的卷积与g(t)和f(t)的卷积相等;结合律是指三个信号f(t)、g(t)和h(t)的卷积可以分别进行两两卷积后再进行一次卷积;分配律是指一个信号与两个信号的卷积等于该信号分别与两个信号卷积后的和。
三、实验内容1. 实验一:连续时间信号卷积实验(1)选用信号:选取两个连续时间信号f(t)和g(t),其中f(t)为矩形脉冲信号,g(t)为指数衰减信号。
(2)卷积计算:根据卷积的定义,计算f(t)和g(t)的卷积F(t)。
(3)结果分析:观察F(t)的波形,分析卷积结果的物理意义。
2. 实验二:离散时间信号卷积实验(1)选用信号:选取两个离散时间信号f[n]和g[n],其中f[n]为单位阶跃信号,g[n]为矩形脉冲信号。
(2)卷积计算:根据离散时间信号卷积的定义,计算f[n]和g[n]的卷积F[n]。
(3)结果分析:观察F[n]的波形,分析卷积结果的物理意义。
3. 实验三:MATLAB仿真实验(1)选用信号:选取两个连续时间信号f(t)和g(t),其中f(t)为正弦信号,g(t)为余弦信号。
(2)MATLAB编程:利用MATLAB的信号处理工具箱,编写程序实现f(t)和g(t)的卷积运算。
(3)结果分析:观察MATLAB仿真得到的卷积结果,分析其物理意义。
四、实验结果与分析1. 实验一:连续时间信号卷积实验(1)实验结果:通过计算得到f(t)和g(t)的卷积F(t)的波形。
实验一 连续时间信号分析一、实验目的(一)掌握使用Matlab 表示连续时间信号1、学会运用Matlab 表示常用连续时间信号的方法2、观察并熟悉常用信号的波形和特性(二)掌握使用Matlab 进行连续时间信号的相关运算1、学会运用Matlab 进行连续时间信号的时移、反褶和尺度变换2、学会运用Matlab 进行连续时间信号微分、积分运算3、学会运用Matlab 进行连续时间信号相加、相乘运算4、学会运用Matlab 进行连续时间信号卷积运算二、实验条件装用Matlab R2015a 的电脑。
三、实验内容1、利用Matlab 命令画出下列连续信号的波形图。
(1))4/3t (2cos π+ 程序:t=-3:0.01:3; ft=2*cos(3*t+pi/4); plot(t,ft)图像:(2))t (u )e 2(t--程序:t=-6:0.01:6; ut=(t>=0);ft=(2-1*exp(-t)).*ut; plot(t,ft)图像:(3))]2()(u )][t (cos 1[--+t u t π 程序:t=-6:0.01:6; ut=(t>=0); ut2=(t>=2);ft=(1+cos(pi*t)).*(ut-ut2); plot(t,ft)图像:2、利用Matlab 命令画出复信号)4/t (j 2e )t (f π+=的实部、虚部、模和辐角。
程序:t=0:0.01:20;ft=2*exp(1j*(t+pi/4));subplot(2,2,1);plot(t,real(ft));title('ʵ²¿');axis([-0.5,20,-2.5,2.5]); subplot(2,2,2);plot(t,imag(ft));title('Ð鲿');axis([-0.5,20,-2.5,2.5]); subplot(2,2,3);plot(t,abs(ft));title('Ä£');axis([-0.5,20,-0.5,2.5]); subplot(2,2,4);plot(t,angle(ft));title('·ø½Ç');axis([-0.5,20,-3.5,3.5]);图像:3、已知信号的波形如下图所示:试用Matlab 命令画出()()()()2332----t f t f t f t f ,,,的波形图。
课程实验报告题目:连续时间信号的卷积及信号的频域分析学院通信与信息工程学院学生姓名班级学号指导教师开课学院通信与信息工程学院日期实验内容:(一)连续时间信号的卷积问题1:用计算机算卷积是把连续信号进行采样,得到一个个离散数值,然后用数值计算代替连续信号的卷积,请推导数值计算与连续信号的卷积之间的关系。
答:x(t)和y(t)为两个连续信号,进行采样后,得到离散值x(ne)与y(ne)(e为取样时间间隔),他们可以近似表示成x和y的函数值。
那么两函数卷积结果 h (t )=⎰∞∞--drr t y r x )()( ,即 h (ke )=∑∞-∞=-n e ne ke h ne x )()(,用k 替代ke h (k )=∑∞-∞=-n en k y n x )()( ,因此,若令e 无限趋于零,那么h (k )的极限值即是两函数卷积函数值。
上机题1.已知两个信号)2()1()(1---=t t t x εε和)1()()(2--=t t t x εε,试分别画出)(),(21t x t x 和卷积)()()(21t x t x t y *=的波形。
>> T=0.01; >> t1=1;t2=2; >> t3=0;t4=1; >> t=0:T:t2+t4;>> x1=ones(size(t)).*((t>t1)-(t>t2)); >> x2=ones(size(t)).*((t>t3)-(t>t4)); >> y=conv(x1,x2)*T;>> subplot(3,1,1),plot(t,x1); >> ylabel('x1(t)');>> subplot(3,1,2),plot(t,x2); >> ylabel('x2(t)');>> subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1)); >> ylabel('y(t)=x1*x2'); >> xlabel('----->t/s');上机题2.已知两个信号)()(t e t x t ε-=和)()(2/t te t h t ε-=,试用数值计算法求卷积,并分别画出)(),(t h t x 和卷积)()()(t h t x t y *=的波形。
实验一连续时间信号的时域和频域分析一. 实验目的:1. 熟悉MATLAB 软件平台。
2. 掌握MATLAB 编程方法、常用语句和可视化绘图技术。
3. 编程实现常用信号及其运算MATLAB 实现方法。
4. 编程实现常用信号的频域分析。
二. 实验原理:1、连续时间信号的描述:(1)向量表示法连续信号是指自变量的取值范围是连续的,且对于一切自变量的取值,除了有若干个不连续点之外,信号都有确定的值与之对应。
严格来说,MATLAB 并不能处理连续信号,而是用等时间间隔点的样值来近似表示连续信号。
当取样时间间隔足够小时,这些离散的样值就能较好地近似连续信号。
矩阵是MATLAB 进行数据处理的基本单元,矩阵运算是MATLAB 最重要的运算。
通常意义上的数量(也称为标量)在MATLAB 系统中是作为1×1 的矩阵来处理的,而向量实际上是仅有一行或者一列的矩阵。
通常用向量表示信号的时间取值范围,如t = -5:5,但信号x(t)、向量t 本身的下标都是从1 开始的,因此必须用一个与向量x 等长的定位时间变量t,以及向量x,才能完整地表示序列x(t)。
在MATLAB 可视化绘图中,对于以t 为自变量的连续信号,在绘图时统一用plot 函数;而对n 为自变量的离散序列,在绘图时统一用stem 函数。
(2)符号运算表示法符号对象(Symbolic Objects 不同于普通的数值计算)是Matlab 中的一种特殊数据类型,它可以用来表示符号变量、表达式以及矩阵,利用符号对象能够在不考虑符号所对应的具体数值的情况下能够进行代数分析和符号计算(symbolic math operations),例如解代数方程、微分方程、进行矩阵运算等。
符号对象需要通过sym 或syms 函数来指定, 普通的数字转换成符号类型后也可以被作为符号对象来处理.我们可以用一个简单的例子来表明数值计算和符号计算的区别: 2/5+1/3 的结果为0.7333(double 类型数值运算), 而sym(2)/sym(5)+sym(1)/sym(3)的结果为11/15, 且这里11/15 仍然是属于sym 类型, 是符号数。
一、实验目的1. 理解信号卷积的概念及其物理意义。
2. 掌握信号卷积的计算方法,包括连续卷积和离散卷积。
3. 分析卷积运算在信号处理中的应用,如信号滤波、信号重构等。
二、实验原理1. 信号卷积的概念信号卷积是指两个信号x(t)和h(t)的乘积在时间域上的积分。
卷积运算可以描述信号之间的相互作用和影响,对于信号处理、通信系统、控制系统等领域具有重要的应用。
2. 卷积的数学表示(1)连续卷积设x(t)和h(t)为两个连续信号,它们的卷积y(t)可以表示为:y(t) = ∫[x(τ)h(t-τ)]dτ(2)离散卷积设x[n]和h[n]为两个离散信号,它们的卷积y[n]可以表示为:y[n] = ∑[x[k]h[n-k]]3. 卷积的性质(1)交换律:x(t) h(t) = h(t) x(t)(2)结合律:(x(t) h(t)) g(t) = x(t) (h(t) g(t))(3)分配律:x(t) (h(t) + g(t)) = x(t) h(t) + x(t) g(t)(4)卷积的导数:d/dt(x(t) h(t)) = x(t) d/dt(h(t))三、实验仪器与设备1. 双踪示波器2. 信号源3. 信号处理模块4. 计算机5. MATLAB软件四、实验内容与步骤1. 连续信号卷积实验(1)选择两个连续信号,如方波信号和三角波信号。
(2)利用示波器观察两个信号的波形。
(3)通过计算机计算两个信号的卷积,并观察卷积结果的波形。
2. 离散信号卷积实验(1)选择两个离散信号,如单位阶跃信号和单位冲激信号。
(2)利用示波器观察两个信号的波形。
(3)通过计算机计算两个信号的卷积,并观察卷积结果的波形。
3. 卷积运算在信号处理中的应用实验(1)信号滤波:选择一个信号,如含噪声的信号,通过卷积运算实现滤波操作,去除噪声。
(2)信号重构:选择一个信号,如被压缩的信号,通过卷积运算实现信号重构,恢复原始信号。
五、实验结果与分析1. 连续信号卷积实验结果通过实验,我们可以观察到连续信号卷积的结果。
信号与系统实验报告目录1. 内容概要 (2)1.1 研究背景 (3)1.2 研究目的 (4)1.3 研究意义 (4)2. 实验原理 (5)2.1 信号与系统基本概念 (7)2.2 信号的分类与表示 (8)2.3 系统的分类与表示 (9)2.4 信号与系统的运算法则 (11)3. 实验内容及步骤 (12)3.1 实验一 (13)3.1.1 实验目的 (14)3.1.2 实验仪器和设备 (15)3.1.4 实验数据记录与分析 (16)3.2 实验二 (16)3.2.1 实验目的 (17)3.2.2 实验仪器和设备 (18)3.2.3 实验步骤 (19)3.2.4 实验数据记录与分析 (19)3.3 实验三 (20)3.3.1 实验目的 (21)3.3.2 实验仪器和设备 (22)3.3.3 实验步骤 (23)3.3.4 实验数据记录与分析 (24)3.4 实验四 (26)3.4.1 实验目的 (27)3.4.2 实验仪器和设备 (27)3.4.4 实验数据记录与分析 (29)4. 结果与讨论 (29)4.1 实验结果汇总 (31)4.2 结果分析与讨论 (32)4.3 结果与理论知识的对比与验证 (33)1. 内容概要本实验报告旨在总结和回顾在信号与系统课程中所进行的实验内容,通过实践操作加深对理论知识的理解和应用能力。
实验涵盖了信号分析、信号处理方法以及系统响应等多个方面。
实验一:信号的基本特性与运算。
学生掌握了信号的表示方法,包括连续时间信号和离散时间信号,以及信号的基本运算规则,如加法、减法、乘法和除法。
实验二:信号的时间域分析。
在本实验中,学生学习了信号的波形变换、信号的卷积以及信号的频谱分析等基本概念和方法,利用MATLAB工具进行了实际的信号处理。
实验三:系统的时域分析。
学生了解了线性时不变系统的动态响应特性,包括零状态响应、阶跃响应以及脉冲响应,并学会了利用MATLAB进行系统响应的计算和分析。
信号与系统公式总结在信号与系统的学习过程中,公式总结是非常重要的,它可以帮助我们更好地理解和掌握知识。
下面将对信号与系统中常见的公式进行总结,希望能够对大家的学习有所帮助。
一、基本概念公式总结。
1. 信号的分类:连续时间信号,x(t)。
离散时间信号,x[n]2. 基本信号:单位冲激函数,δ(t)或δ[n]阶跃函数,u(t)或u[n]3. 基本性质:奇偶性,x(t) = x(-t),x[n] = x[-n]周期性,x(t) = x(t+T),x[n] = x[n+N]二、时域分析公式总结。
1. 基本运算:时移性质,x(t-t0)或x[n-n0]反褶性质,x(-t)或x[-n]放大缩小,Ax(t)或Ax[n]2. 基本运算公式:加法,x1(t) + x2(t)或x1[n] + x2[n]乘法,x1(t)x2(t)或x1[n]x2[n]三、频域分析公式总结。
1. 傅里叶变换:连续时间信号,X(ω) = ∫x(t)e^(-jωt)dt。
离散时间信号,X(e^jω) = Σx[n]e^(-jωn)。
2. 傅里叶变换性质:线性性质,aX1(ω) + bX2(ω)。
时移性质,x(t-t0)对应X(ω)e^(-jωt0)。
频移性质,x(t)e^(jω0t)对应X(ω-ω0)。
四、系统分析公式总结。
1. 系统性质:线性性,y(t) = ax1(t) + bx2(t)。
时不变性,y(t) = x(t-t0)对应h(t-t0)。
2. 系统时域分析:离散卷积,y[n] = Σx[k]h[n-k]连续卷积,y(t) = ∫x(τ)h(t-τ)dτ。
3. 系统频域分析:系统函数,H(ω) = Y(ω)/X(ω)。
五、采样定理公式总结。
1. 采样定理:连续信号采样,x(t)对应x[n],x[n] = x(nT)。
重建滤波器,h(t) = Tsinc(πt/T)。
六、傅里叶级数公式总结。
1. 傅里叶级数:周期信号的傅里叶级数展开。
实验三 连续时间LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义;2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用;3、学习和掌握幅度特性、相位特性以及群延时的物理意义;4、掌握用MA TLAB 语言进行系统频响特性分析的方法。
基本要求:掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法,深刻理解LTI 系统的频率响应特性的物理意义,理解滤波和滤波器的概念,掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。
二、实验原理及方法1 连续时间LTI 系统的频率响应所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response ),是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。
上图中x(t)、y(t)分别为系统的时域激励信号和响应信号,h(t)是系统的单位冲激响应,它们三者之间的关系为:)(*)()(t h t x t y =,由傅里叶变换的时域卷积定理可得到:)()()(ωωωj H j X j Y =3.1或者: )()()(ωωωj X j Y j H =3.2)(ωj H 为系统的频域数学模型,它实际上就是系统的单位冲激响应h(t)的傅里叶变换。
即⎰∞∞--=dt e t h j H tj ωω)()(3.3由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换,如果h(t)是收敛的,或者说是绝对可积(Absolutly integrabel )的话,那么H(j ω)一定存在,而且H(j ω)通常是复数,因此,也可以表示成复数的不同表达形式。
在研究系统的频率响应时,更多的是把它表示成极坐标形式:)()()(ωϕωωj ej H j H = 3.4上式中,)j (ωH 称为幅度频率相应(Magnitude response ),反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,)(ωϕ称为相位特性(Phase response ),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。
电子信息工程系实验报告课程名称:数字信号处理Array实验项目名称:用DFT(FFT)对连续信号进行频谱分析实验时间:班级:通信姓名: xxp 学号:一、实验目的:1.掌握用DFT(FFT)对模拟信号进行谱分析的方法,理解可能出现的分析误差及其原因,以便在实际中正确应用FFT。
2.熟悉应用FFT实现两个序列的线性卷积的方法。
二、实验原理:1.用DFT(FFT)对连续信号进行频谱分析用DFT(FFT)对模拟信号做谱分析是一种近似的谱分析。
首先一般的模拟信号(周期信号除外)的频谱是连续谱,而用FFT做谱分析得到的是数字谱,因此应该取FFT的点数多一些,用它的包络作为模拟信号的近似谱。
另外,如果模拟信号不是严格的带限信号,会因为频谱混叠现象引起谱分析的误差,这种情况下可以预先将模拟信号进行预滤,或者尽量将采样频率取高一些。
最后要注意一般的模拟信号是无限长的,分析时要截断,截断的长度与对模拟信号进行频谱分析的分辨率有关。
如果要进行谱分析的模拟信号是周期信号,最好选择观测时间是信号周期的整数倍,如果不知道信号的周期,要尽量选择观测时间长一些,以减少截断效应的影响。
在运用DFT(FFT)对模拟信号进行谱分析的过程中主要可能产生以下三种误差:(1) 混叠现象对模拟信号进行谱分析时首先要对其采样,当采样速率不满足Nyquist定理时,就会发生频谱混叠,使得采样后的信号序列频谱不能真实的反映原模拟信号的频谱。
避免混叠现象的唯一方法是保证采样速率足够高,使频谱混叠现象不致出现,即在确定采样频率之前,必须对频谱的性质有所了解,在一般情况下,为了保证高于折叠频率的分量不会出现,在采样前,先用低通模拟滤波器对信号进行滤波。
(2) 截断效应实际中我们往往用截短的序列来近似很长的甚至是无限长的序列,这样可以使用较短的DFT来对信号进行频谱分析,这种截短等价于给原信号序列乘以一个矩形窗函数,也相当于在频域将信号的频谱和矩形窗函数的频谱卷积,所得的频谱是原序列频谱的扩展。
连续时间信号的卷积及信号的频域分析实验报告(1)连续时间信号的卷积及信号的频域分析实验报告一、实验目的本实验的主要目的是通过对于两个时间域信号的卷积运算,掌握信号卷积运算的基本原理及操作方法;同时,利用MATLAB软件完成信号的傅里叶变换,了解信号在频域的频谱特征。
二、实验内容1、连续时间信号的卷积运算利用MATLAB软件中conv函数进行两个信号的卷积运算,并观察结果。
2、信号在频域的频谱特征- 利用MATLAB软件中fft函数对信号进行傅里叶变换,并获取其频域表示;- 利用MATLAB软件中ifft函数对信号进行逆傅里叶变换,恢复其原始时间域信号;- 观察不同频率成分对于信号的影响,并分析其原因。
三、实验步骤1、连续时间信号的卷积运算首先在MATLAB软件中定义两个连续时间信号,如下所示:t1 = 0:0.1:10;x1 = sin(2*pi*5*t1); % 正弦波信号t2 = 0:0.1:10;x2 = exp(-(t2-5).^2); % 高斯脉冲信号然后,使用conv函数进行卷积运算,并绘制出卷积后的信号图像。
x3 = conv(x1,x2,'same'); % 卷积运算figure; % 绘制卷积后的信号图像subplot(3,1,1);plot(t1,x1);xlabel('时间/s');ylabel('幅值');title('正弦波信号');subplot(3,1,2);plot(t2,x2);xlabel('时间/s');ylabel('幅值');title('高斯脉冲信号');subplot(3,1,3);plot(t1,x3);xlabel('时间/s');ylabel('幅值');title('卷积信号');2、信号在频域的频谱特征首先,通过fft函数对于时间域信号进行傅里叶变换,获取其频域表示。
实验三常见信号的MATLAB表示及运算一、实验目的1. 熟悉常见信号的意义、特性及波形2. 学会使用MATLAB表示信号的方法并绘制信号波形3.掌握使用MATLAB进行信号基本运算的指令4.熟悉用MATLAB实现卷积积分的方法二、实验原理根据MA TLAB的数值计算功能和符号运算功能, 在MATLAB中, 信号有两种表示方法, 一种是用向量来表示, 另一种则是用符号运算的方法。
在采用适当的MATLAB语句表示出信号后, 就可以利用MATLAB中的绘图命令绘制出直观的信号波形了。
1.连续时间信号从严格意义上讲, MATLAB并不能处理连续信号。
在MATLAB中, 是用连续信号在等时间间隔点上的样值来近似表示的, 当取样时间间隔足够小时, 这些离散的样值就能较好地近似出连续信号。
在MATLAB中连续信号可用向量或符号运算功能来表示。
⑴向量表示法对于连续时间信号, 可以用两个行向量f和t来表示, 其中向量t是用形如的命令定义的时间范围向量, 其中, 为信号起始时间, 为终止时间, p为时间间隔。
向量f为连续信号在向量t所定义的时间点上的样值。
⑵符号运算表示法如果一个信号或函数可以用符号表达式来表示, 那么我们就可以用前面介绍的符号函数专用绘图命令ezplot()等函数来绘出信号的波形。
⑶常见信号的MATLAB表示单位阶跃信号单位阶跃信号的定义为:方法一: 调用Heaviside(t)函数首先定义函数Heaviside(t) 的m函数文件,该文件名应与函数名同名即Heaviside.m。
%定义函数文件,函数名为Heaviside,输入变量为x,输出变量为yfunction y= Heaviside(t)y=(t>0); %定义函数体, 即函数所执行指令%此处定义t>0时y=1,t<=0时y=0, 注意与实际的阶跃信号定义的区别。
方法二: 数值计算法在MATLAB中, 有一个专门用于表示单位阶跃信号的函数, 即stepfun( )函数, 它是用数值计算法表示的单位阶跃函数。
第三章.连续时间系统的频域分析一、任意信号在完备正交函数系中的表示法(§)信号分解的目的:● 将任意信号分解为单元信号之和,从而考查信号的特性。
●简化电路分析与运算,总响应=单元响应之和。
1.正交函数集任意信号)(t f 可表示为n 维正交函数之和:原函数()()()t g t g t g r Λ21,相互正交:⎩⎨⎧=≠=⋅⎰nm K nm dt t g t g m t t n m ,,0)()(21()t g r 称为完备正交函数集的基底。
一个信号可用完备的正交函数集表示,.正弦函数集有许多方便之处,如易实现等,我们主要讨论如何用正弦函数集表示信号。
2.能量信号和功率和信号(§一)设()t i 为流过电阻R 的电流,瞬时功率为R t i t P )()(2=一般说来,能量总是与某一物理量的平方成正比。
令R = 1Ω,则在整时间域内,实信号()t f 的能量,平均功率为: 讨论上述两个式子,只可能出现两种情况: ✍∞<<W 0(有限值) 0=P✍∞<<P 0(有限值)∞=W满足✍式的称为能量信号,满足✍式称功率信号。
3.帕斯瓦尔定理设{})(t g r 为完备的正交函数集,即信号的能量 基底信号的能量 各分量此式称为帕斯瓦尔定理 P331 式(6-81) (P93, P350) 左边是信号能量,右边是各正交函数的能量。
物理意义:一个信号所含有的能量(功率)恒等于此信号在完备正交函数集中各分量能量(功率)之和。
二、周期信号的频谱分析——傅里叶级数(1) 周期信号傅里叶级数有两种形式三角形式: ()∑∞=++=1110sin cos )(n n nt n b t n aa t f ωω=∑∞=++110)cos(n n nt n cc ϕω指数形式:t jn n e n F t f 1)()(1ωω∑∞-∞==(2) 周期信号的频谱是离散谱,三个性质收敛性()↓↑)(,1ωn F n谐波性:(离散性)谱线只出现在1ωn 处,唯一性:)(t f 的谱线唯一(3)两种频谱图的关系● 三角形式:ω~n c ,ωφ~n 单边频谱● 指数形式:ωω~)(1n F , ωφ~n 双边频谱两者幅度关系 )(1ωn F =()021≠n c n000a c F ==● 指数形式的幅度谱为偶函数 ●指数形式的相位谱为奇函数(4) 引入负频率对于双边频谱,负频率)(1ωn ,只有数学意义,而无物理意义。
实验三 连续时间LTI 系统的频域分析一、实验目的1、掌握系统频率响应特性的概念及其物理意义;2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法,理解具有不同频率响应特性的滤波器对信号的滤波作用;3、学习和掌握幅度特性、相位特性以及群延时的物理意义;4、掌握用MATLAB 语言进行系统频响特性分析的方法。
基本要求:掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法,深刻理解LTI 系统的频率响应特性的物理意义,理解滤波和滤波器的概念,掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。
二、实验原理及方法1 连续时间LTI 系统的频率响应所谓频率特性,也称为频率响应特性,简称频率响应(Frequency response ),是指系统在正弦信号激励下的稳态响应随频率变化的情况,包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。
上图中x(t)、y(t)分别为系统的时域激励信号和响应信号,h(t)是系统的单位冲激响应,它们三者之间的关系为:)(*)()(t h t x t y =,由傅里叶变换的时域卷积定理可得到:)()()(ωωωj H j X j Y =3.1或者: )()()(ωωωj X j Y j H =3.2)(ωj H 为系统的频域数学模型,它实际上就是系统的单位冲激响应h(t)的傅里叶变换。
即⎰∞∞--=dt et h j H tj ωω)()( 3.3由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换,如果h(t)是收敛的,或者说是绝对可积(Absolutly integrabel )的话,那么H(j ω)一定存在,而且H(j ω)通常是复数,因此,也可以表示成复数的不同表达形式。
在研究系统的频率响应时,更多的是把它表示成极坐标形式:)()()(ωϕωωj ej H j H = 3.4上式中,)j (ωH 称为幅度频率相应(Magnitude response ),反映信号经过系统之后,信号各频率分量的幅度发生变化的情况,)(ωϕ称为相位特性(Phase response ),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。
课程实验报告
题目:连续时间信号的卷积
及信号的频域分析
学院通信与信息工程学院
学生姓名
班级学号
指导教师
开课学院
日期 2010.11.18
实验内容:(一)连续时间信号的卷积
问题1:用计算机算卷积是把连续信号进行采样,得到一个个离散数值,然后用数值计算代替连续信号的卷积,请推导数值计算与连续信号的卷积之间的关系。
(学生回答问题)
答:连续函数x(t)和y(t)的卷积为:τττd t h x t h t x t y )()()()()(-=*=⎰∞
∞-(F2-1)
若x(t)和h(t)分别仅在时间区间),(21t t 和),(43t t 有非零值,则ττετεττετετεεεεd t t t t t h t t x t t t t t h t t t t t x t y )]()()[()]()([)()]
()()[()]()()[()(43214321------∙---=---*---=⎰∞∞-
要使y (t )为非零值,必须有:1)()(21=---t t τετε和1)()(43=-----t t t t τετε 从而,应同时满足:21t t <<τ和43t t t +<<+ττ,即4231t t t t t +<<+
由此得出结论:若x(t)和h(t)分别仅在时间区间),(21t t 和),(43t t 有非零值,则卷积)()()(t h t x t y *=有非零值的时间区间为),(4231t t t t ++。
对卷积积分式(F2-1)进行数值计算时近似为:∆∆-∆∆=
∆∑∞
-∞=)()()(n k h n x k y n 记作∆*=∆-=
∑∞-∞=n k h k x n k h n x k y )()()()()( (F2-2)
式中,y(k)、x(k)和h(k)分别为对y(t)、x(t)和h(t)以为∆时间间隔进行采样所得的离散序列。
相应的可得出结论:若x(k)和h(k)分别心在序号区间],[21k k 和
],[43k k 有非零的值,则离散卷积(卷积和))()()(t h t x t y *=有非零值的序号区间为],[4231k k k k ++。
上机题1.已知两个信号)2()1()(1---=t t t x εε和)1()()(2--=t t t x εε,试分别画出)(),(21t x t x 和卷积)()()(21t x t x t y *=的波形。
(上机原程序及所画出的波形图)
T=0.01;
t1=1,t2=2;
t3=0,t4=1;
t=0:T:t2+t4;
x1=ones(size(t)).*((t>t1)-(t>t2));
x2=ones(size(t)).*((t>t3)-(t>t4));
y=conv(x1,x2)*T;
subplot(3,1,1),plot(t,x1);
ylabel('x1(t)');
subplot(3,1,2),plot(t,x2);
ylabel('x2(t) ');
subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1));
ylabel('y(t)=x1*x2');
xlabel('――――→t/s');
上机题2.已知两个信号)()(t e t x t ε-=和)()(2/t te t h t ε-=,试用数值计算法求卷积,并分别画出)(),(t h t x 和卷积)()()(t h t x t y *=的波形。
(上机原程序及所画出的波形图)
syms tao;
t=sym('t', 'positive');
xt=exp(-t);
ht=t*exp(-t/2);
xh_tao=subs(xt,t,tao)*subs(ht,t,t-tao);
yt=int(xh_tao,tao,0,t)
执行结果为
yt =
4*exp(-t)+2*t*exp(-1/2*t)-4*exp(-1/2*t)
t2=3;t4=11;
T=0.01;
t=0:T:t2+t4;
x=exp(-t).*((t>0)-(t>t2));
h=t.*exp(-t/2).*((t>0)-(t>t4));
y=conv(x,h)*T;
yt=4*exp(-t)+2*t.*exp(-1/2*t)-4*exp(-1/2*t);
subplot(3,1,1),plot(t,x);
ylabel('x(t)');
subplot(3,1,2),plot(t,h);
ylabel('h(t) ');
subplot(3,1,3),plot(t,y(1:(t2+t4)/T+1),t,yt, '--r');
legend('by numerical ’,’Theoretical');
ylabel('y=x*h');
xlabel('---→t/s');
实验内容:(二)信号的频域分析
上机题3.求周期矩形脉冲信号的频谱图,已知s T s A 5.0,1.0,1===τ (上机原程序及所画出的波形图)
a=1;tao=0.1;t=0.5;
n0=t/tao;
n=0:2*n0;
fn_p=a*tao/t*(sin(n*pi*tao/t+eps*(n==0)))./(n*pi*tao/t+eps*(n==0));
fn_pabs=abs(fn_p);
fn_pang=angle(fn_p);
fn_mabs=fliplr(fn_pabs(2:11));
fn_mang=-fliplr(fn_pang(2:11));
fnabs=[fn_mabs fn_pabs];
fnang=[fn_mang fn_pang];
subplot(2,1,1),stem((-2*n0:2*n0),fnabs);
text(4,0.11, 'amplitude spectrum');
subplot(2,1,2),stem((-2*n0:2*n0),fnang);
text(-2,2, 'phase spectrum');
xlabel('n');grid
问题2:改变信号的周期,比较他们的计算结果。
直接在程序中改变T会报错。
应在第4行输入:xian_no=2*no+1,而且将程序中的11改为xian_no (学生通过上机画出不同周期对应的计算结果并回答问题)
T=0.2
T=2
当T增大时,第一个零交点不变,即信号频带宽度不变;由于T增大,主峰高度减小,谱线间隔变密。
可以想见,如果周期无限增加,那么,这时周期信号就转化为非周期信号。
各条谱线间将趋向于零,那么,周期信号离散频谱将过度到非
周期信号的连续频谱,频谱包络线的形状没有改变。