信号与系统报告 实验5 连续系统的复频域分析实验
实验:连续系统的频域分析
实验4:连续系统的频域分析一、实验目的(1)掌握连续时间信号的傅里叶变换和傅里叶逆变换的实现方法。
(2)掌握傅里叶变换的数值计算方法和绘制信号频谱的方法。
二、实验原理 1.周期信号的分解根据傅里叶级数的原理,任何周期信号都可以分解为三角级数的组合——称为()f t 的傅里叶级数。
在误差确定的前提下,可以由一组三角函数的有限项叠加而得到。
例如一个方波信号可以分解为:11114111()sin sin 3sin 5sin 7357E f t t t t t ωωωωπ⎛⎫=++++ ⎪⎝⎭合成波形所包含的谐波分量越多,除间断点附近外,它越接近于原波形,在间断点附近,即使合成的波形所含谐波次数足够多,也任存在约9%的偏差,这就是吉布斯现象(Gibbs )。
2.连续时间信号傅里叶变换的数值计算 由傅里叶变换的公式:()()lim()j tj n n F j f t edt f n e ωωττωττ∞∞---∞→=-∞==∑⎰当()f t 为时限信号时,上式中的n 取值可以认为是有限项N,则有:()(),0k Nj n n F k f n e k N ωτττ-==≤≤∑,其中2k k N πωτ=3.系统的频率特性连续LTI 系统的频率特性称为频率响应特性,是指在正弦信号激励作用下稳态响应随激励信号频率的变化而变化的情况,表示为()()()Y H X ωωω=三、实验内容与方法 1.周期信号的分解【例1】用正弦信号的叠加近似合成一个频率为50Hz 的方波。
MA TLAB 程序如下: clear all; fs=10000; t=[0:1/fs:0.1]; f0=50;sum=0; subplot(211) for n=1:2:9plot(t,4/pi*1/n*sin(2*pi*n*f0*t),’k ’); hold on; endtitle(‘信号叠加前’); subplot(212) for n=1:2:9;sum=sum+4/pi*1/n*sin(2*pi*n*f0*t);endplot(t,sum,’k ’); title(‘信号叠加后’); 产生的波形如图所示:00.010.020.030.040.050.060.070.080.090.1-2-1012信号叠加前00.010.020.030.040.050.060.070.080.090.1-2-1012信号叠加后2.傅里叶变换和逆变换的实现求傅里叶变换,可以调用fourier 函数,调用格式为F=fourier(f,u,v),是关于u 的函数f 的傅里叶变换,返回函数F 是关于v 的函数。
北京理工大学信号与系统实验报告5 连续时间系统的复频域分析
实验5连续时间系统的复频域分析(综合型实验)一、实验目的1)掌握拉普拉斯变换及其反变换的定义并掌握MATLAB 实现方法。
2)学习和掌握连续时间系统函数的定义及复频域分析方法。
3)掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法 1.拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为(s)(t)e st X x dt +∞--∞=⎰(1)拉普拉斯反变换为1(t)(s)e 2j st j x X ds j σσπ+∞-∞=⎰ (2)MATLAB 中相应函数如下:(F)L laplace = 符号表达式F 拉氏变换,F 中时间变量为t ,返回变量为s 的结果表达式。
(F,t)L laplace =用t 替换结果中的变量s 。
()F ilaplace L =以s 为变量的符号表达式L 的拉氏反变换,返回时间变量为t 的结果表达式。
(,)F ilaplace L x =用x 替换结果中的变量t 。
拉氏变换还可采用部分分式法,当(s)X 为有理分式时,它可以表示为两个多项式之比:110110...(s)(s)(s)...M M M M N N N N b s b s b N X D a s a s a ----+++==+++ (3) 上式可以采用部分分式法展成以下形式1212(s)...N Nr r rX s p s p s p =+++--- (4) 再通过查找常用拉氏变换对易得反变换。
利用residue 函数可将X(s)展成(4)式形式,调用格式为:[r,p,k]residue(b,a)=其中b 、a 为分子和分母多项式系数向量,r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。
2.连续时间系统的系统函数连续时间系统的系统函数是指系统单位冲激响应的拉氏变换(s)(t)e stH h dt +∞--∞=⎰(5)连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。
实验5-连续时间系统的复频域分析报告
实验5-连续时间系统的复频域分析报告
本实验的目的是研究连续时间系统的复频域分析。
首先,构建了一个由推力继电器组
成的系统,其模型为图1所示。
再将此系统内建模,得到开环传递函数
G(s)=K/[(s+1)(s+1)(s+2)],其中1为系统参数,s为复频变量。
然后使用MATLAB编程,实现基于Laplace变换计算复频域函数和系统振型,并以一系列频率点绘制系统频率响应
曲线等曲线,从而评估系统性能。
实验结果表明,当系统参数K处于[6.5,9.2]中时,系统的复频响应表现出了各向同
性的性能(图2),表明系统具有更一致的响应特性,并且误差幅值在0.03以内保持稳定,说明系统具有良好的稳定性性能。
此外,系统振型(图3)也说明了系统的稳定性,振型
稳定时间较短,且交叉率较小,说明系统具有良好的稳定性能。
综上,连续时间系统的复频域分析中,MATLAB编程在系统参数K为[6.5,9.2]范围内时,运用Laplace变换和求和函数,成功绘制出系统的复频响应曲线,以及相应的系统振型,从而对系统的复频响应、稳定行为等做出定量性、全面性的评估,为系统运行提供了
可靠的参考。
实验5 连续时间系统的频域和复频域分析
实验5 连续时间系统的频域和复频域分析一.实验目的1.掌握和理解连续时间函数系统频率相应、系统函数的概念和物理意义。
2.学习和掌握连续时间系统频域、复频域的分析方法。
3.掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二.实验原理1.连续时间系统的频率响应系统的频率响应定义为:ττωωτd eh j H j -∞∞-⎰=)()(H (ωj )反映了LTI 连续时间系统对不同频率信号的相应特性,是系统内在固有的特性,与外部激励无关。
H (ωj )又可以表示为)()()(ωθωωj ej H j H =其中)(ωj H 称为系统的幅度响应,)(ωθ成为系统的相应响应。
对于由下述微分方程描述的LTI 连续时间系统∑∑===Mm m n Nn n n t xb t ya 0)(0)()()(其频率响应H (ωj )可以表示为下列式子所示的ωj 的有理多项式1110111...)()(...)()()()()(a j a j a j a b j b j b j b X Y j H N N N N M M M M ++++++++==----ωωωωωωωωωMATLAB 的信号处理工具箱提供了专门的函数freqs ,用来分析连续时间系统的频率响应,该函数有下列几种调用格式:[h,w]=freqs(b,a) 计算默认频率范围内200个频率点上的频率响应的取样值,这200个频率点记录在w 中。
h=freqs (b ,a ,w ) b 、a 分别为表示H (ωj )的有理多项式中分子和分母多项式的系数向量,w 为频率取样点,返回值h 就是频率响应在频率取样点上的数值向量。
[h ,w]=freqs (b ,a ,n) 计算默认频率范围内n 个频率点上的频率响应的取样值,这n 个频率点记录在w 中。
Freqs (b ,a ,……) 这种调用格式不返回频率响应的取样值,而是以对数坐标的方式绘出来系统的频率响应和相频响应。
连续时间信号与系统的频域分析报告
连续时间信号与系统的频域分析报告1. 引言连续时间信号与系统的频域分析是信号与系统理论中的重要分支,通过将信号和系统转换到频域,可以更好地理解和分析信号的频谱特性。
本报告将对连续时间信号与系统的频域分析进行详细介绍,并通过实例进行说明。
2. 连续时间信号的频域表示连续时间信号可以通过傅里叶变换将其转换到频域。
傅里叶变换将信号分解成一系列不同频率的正弦和余弦波的和。
具体来说,对于连续时间信号x(t),其傅里叶变换表示为X(ω),其中ω表示频率。
3. 连续时间系统的频域表示连续时间系统可以通过频域中的频率响应来描述。
频率响应是系统对不同频率输入信号的响应情况。
通过系统函数H(ω)可以计算系统的频率响应。
系统函数是频域中系统输出与输入之比的函数,也可以通过傅里叶变换来表示。
4. 连续时间信号的频域分析频域分析可以帮助我们更好地理解信号的频谱特性。
通过频域分析,我们可以获取信号的频率成分、频谱特性以及信号与系统之间的关系。
常用的频域分析方法包括功率谱密度估计、谱线估计等。
5. 连续时间系统的频域分析频域分析也可以用于系统的性能评估和系统设计。
通过分析系统的频响特性,我们可以了解系统在不同频率下的增益和相位变化情况,进而可以对系统进行优化和设计。
6. 实例分析以音频信号的频域分析为例,我们可以通过对音频信号进行傅里叶变换,将其转换到频域。
通过频域分析,我们可以获取音频信号的频谱图,从而了解音频信号的频率成分和频率能量分布情况。
进一步,我们可以对音频信号进行系统设计和处理,比如对音乐进行均衡、滤波等操作。
7. 结论连续时间信号与系统的频域分析是信号与系统理论中重要的内容,通过对信号和系统进行频域分析,可以更好地理解和分析信号的频谱特性。
频域分析也可以用于系统的性能评估和系统设计,对于音频信号的处理和优化具有重要意义。
总结:通过本报告,我们了解了连续时间信号与系统的频域分析的基本原理和方法。
频域分析可以帮助我们更好地理解信号的频谱特性和系统的频响特性,对系统设计和信号处理具有重要意义。
连续时间信号与系统的频域分析实验报告
《信号与系统》课程实验报告
一.实验原理 1、傅里叶变换 实验原理如下:
傅里叶变换的调用格式
F=fourier(f):返回关于w 的函数;
F=fourier(f ,v):返回关于符号对象v 的函数,而不是w 的函数。
傅里叶逆变换的调用格式
f=ifourier(F):它是符号函数F 的fourier 逆变换,返回关于x 的函数; f=ifourier(f,u):返回关于u 的函数。
2、连续时间信号的频谱图 实验原理如下:
符号算法求解如下:
ft=sym('4*cos(2*pi*6*t)*(heaviside(t+1/4)-heaviside(t-1/4))'); Fw=simplify(fourier(ft)) subplot(121)
ezplot(ft,[-0.5 0.5]),grid on subplot(122)
ezplot(abs(Fw),[-24*pi 24*pi]),grid on 波形图如下所示:
当信号不能用解析式表达时,无法用MATLAB 符号算法求傅里叶变换,则用MATLAB 的数值计算连续信号的傅里叶变换。
∑⎰
∞
-∞
=-→-∞∞
-==n n j t
j e
n f dt e
t f j F ττωτ
ωτω)(lim
)()(0
若信号是时限的,或当时间大于某个给定值时,信号已衰减的很厉害,可以近似地看成时限信号,设n 的取值为N ,有
1
1()
a jw
++
的分母和分子多项式的系数向量,
1、在调用函数fourier()及ifourier()之前,要用syms命令对所用到的变。
信号与系统标准实验报告连续系统的幅频特性精
信号与系统标准实验报告-连续系统的幅频特性(精)电 子 科 技 大 学实 验 报 告学生姓名:楼秋文 学号:2903102008 指导教师:张鹰 一、实验室名称:信号与系统实验室 二、实验项目名称:连续系统的幅频特性测量三、实验原理:正弦波信号)cos()(0t A t x ω=输入连续LTI 系统,输出)(t y 仍为正弦波信号。
信号输入连续LTI 系统在上图中)(cos()()(000ωωωj H t j H A t y ∠+=)通过测量输入)(t x 、输出)(t y 的正弦波信号幅度,计算输入、输出的正弦波信号幅度比值,可以得到系统的幅频特性在0ω处的测量值)(0ωj H 。
改变0ω可以测出不同频率处的系统幅频特性。
四、实验目的:进一步加深对系统的频率特性的了解五、实验内容:实验内容(一)、低通滤波器的幅频特性测量; 实验内容(二)、带通滤波器的幅频特性测量 ;记录不同频率正弦波通过低通、带通滤波器的响应波形,测量其幅度,拟合出频率响应的幅度特性;分析两个滤波器的截止频率。
六、实验器材(设备、元器件):数字信号处理实验箱、信号与系统实验板的低通滤波器模块U11、高通滤波)(ωj H )(x )(t y器模块U21、PC 机端信号与系统实验软件、+5V 电源 、连接线、计算机串口连接线。
七、实验步骤:(一)、低通滤波器的幅频特性测量:实验步骤:1、 信号选择:按实验箱键盘“3”选择“正弦波”,再按“+”或“-”依次选择表3.1中一个频率。
2、 连接接口区的“输入信号1”和“输出信号”,如下图所示。
点击SSP 软件界面上的按钮,观察输入正弦波。
将正弦波频率值和幅度值(Vpp/2, Vpp 为峰-峰值)记录于表。
接口区输入信号1输入信号2输出信号采样信号备用备用观察输入正弦波的连线示意图3、 按下图的模块连线示意图连接各模块。
接口区输入信号1输入信号2输出信号采样信号备用备用低通滤波器U11输入S11输出S12模块连线示意图14、点击SSP 软件界面上的按钮,观察输入正弦波通过连续系统的响应波形;适当调整X 、Y 轴的分辨率可得到实验所需波形。
连续时间信号与系统的频域分析实验报告
《信号与系统》课程实验报告一•实验原理 1傅里叶变换实验原理如下:傅里叶变换的调用格式F=fourier(f):返回关于 W 的函数;F=fourier(f , v):返回关于符号对象V 的函数,而不是W 的函数。
傅里叶逆变换的调用格式f=ifourier(F):它是符号函数F 的fourier 逆变换,返回关于X 的函数; f=ifourier(f,u):返回关于U 的函数。
2、连续时间信号的频谱图实验原理如下: 符号算法求解如下:ft=sym('4*cos(2*pi*6*t)*(heaviside(t+1∕4)-heaviside(t-1∕4))'); FW=SimPlify(fourier(ft))subplot(121)ezplot(ft,[-0.5 0.5]),grid Onsubplot(122) ezplot(abs(Fw),[-24*pi 24*pi]),grid On波形图如下所示:当信号不能用解析式表达时,无法用换,则用MATLAB 的数值计算连续信号的傅里叶变换。
实验步骤或实验方案MATLAB 符号算法求傅里叶变F(j )f(t)ejt dt 叫nf (n )e若信号是时限的,或当时间大于某个给定值时,信号已衰减的很厉 害,可以近似地看成时限信号,设 n 的取值为N ,有4 CO$(12 I )■) (he 如引日环-IMh heaviside(t IeIXW Sin(WM ⅛)yabS(W i -144 >2)3、 用MATLAB 分析LTl 系统的频率特性当系统的频率响应H (jw )是jw 的有理多项式时,有H(S )B(W) b M (jW)Mb Mi (jW)MIL b ι(jw) b oH (jW)NN 1A(W)a N (jw)a ” ι(jw) L α(jw) a °freqs 函数可直接计算系统的频率响应的数值解,其调用格式为H=freqs(b,a,w)其中,a 和b 分别是H(jw)的分母和分子多项式的系数向量,W 定义 了系统频率响应的频率范围,P 为频率取样间隔。
连续时间信号与系统的频域分析实验报告
实验四连续时间信号与系统的频域分析一、实验目的掌握连续时间信号的傅里叶变换及傅里叶逆变换的实现方法,掌握连续时间系统的频域分析方法,熟悉MATLAB 相应函数的调用格式和作用,掌握使用MATLAB 来分析连续时间信号与系统的频域特性及绘制信号频谱图的方法。
二、实验原理(一)连续时间信号与系统的频域分析原理1、连续时间信号的额频域分析 连续时间信号的傅里叶变换为:()()dt e t f j F t j ωω-∞∞-⎰=傅里叶逆变换为:()()ωωπωd e j F t f t j ⎰∞∞-=21()ωj F 称为频谱密度函数,简称频谱。
一般是复函数,可记为:()()()ωϕωωj e j F j F =()ωj F 反映信号各频率分量的幅度随频率ω的变化情况,称为信号幅度频谱。
()ωϕ反映信号各频率分量的相位随频率ω的变化情况,称为信号相位频谱。
2、连续时间系统的频域分析 在n 阶系统情况下,数学模型为:()()()()()()()()t f b dtt df b dt t f d b dt t f d b t y a dtt dy a dt t y d a dt t y d a o m m n m m n o n n n n n n ++++=++++------11111111 令初始条件为零,两端取傅里叶变换,得:()()[]()()()[]()ωωωωωωωωj F b j b j b j b j Y a j a j a j a m n m n n n nn01110111++++=++++----表示为()()()()ωωωωj F j b j Y j a kmk kkn k k∑∑===0则 ()()()()()()()()()∑∑==----=++++++++==nk kk mk kk n n n n m m mm j a j b a j a j a j a b j b j b j b j F j Y j H 0001110111ωωωωωωωωωωω3、系统传递函数 系统传递函数定义为:()()()ωωωj H j Y j H =系统传递函数反映了系统内在的固有的特性,它取决于系统自身的结构及参数,与外部 激励无关,是描述系统特性的一个重要参数。
连续系统复频域分析报告附MATLAB实现信号与系统实验报告
计算机与信息工程学院设计性实验报告专业:通信工程年级/班级:2011级第二学年第二学期一、实验目的1.掌握用matlab分析系统时间响应的方法2.掌握用matlab分析系统频率响应的方法3.掌握系统零、极点分布与系统稳定性关系二、实验原理1.系统函数H(s)系统函数:系统零状态响应的拉氏变换与激励的拉氏变换之比.H(s)=R(s)/E(s)在matlab中可采用多种方法描述系统,本文采用传递函数(系统函数)描述法.在matlab中, 传递函数描述法是通过传递函数分子和分母关于s降幂排列的多项式系数来表示的.例如,某系统传递函数如下则可用如下二个向量num和den来表示:num=[1,1];den=[1,1.3,0.8]2.用matlab分析系统时间响应1)脉冲响应y=impulse(num,den,T)T:为等间隔的时间向量,指明要计算响应的时间点.2)阶跃响应y=setp(num,den,T)T同上.3)对任意输入的响应y=lsim(num,den,U,T)U:任意输入信号. T同上.3.用matlab分析系统频率响应特性频响特性: 系统在正弦激励下稳态响应随信号频率变化的特性.|H(jω)|:幅频响应特性.ϕ(ω):相频响应特性(或相移特性).Matlab求系统频响特性函数freqs的调用格式:h=freqs(num,den,ω)ω:为等间隔的角频率向量,指明要计算响应的频率点.4.系统零、极点分布与系统稳定性关系系统函数H(s)集中表现了系统的性能,研究H(s)在S平面中极点分布的位置,可很方面地判断系统稳定性.1) 稳定系统: H(s)全部极点落于S左半平面(不包括虚轴),则可以满足系统是稳定的.2)不稳定系统: H(s)极点落于S右半平面,或在虚轴上具有二阶以上极点,则在足够长时间后,h(t)仍继续增长, 系统是不稳定的.3)临界稳定系统: H(s)极点落于S平面虚轴上,且只有一阶,则在足够长时间后,h(t)趋于一个非零数值或形成一个等幅振荡.系统函数H(s)的零、极点可用matlab的多项式求根函数roots()求得.极点:p=roots(den)零点:z=roots(num)根据p和z用plot()命令即可画出系统零、极点分布图,进而分析判断系统稳定性.三、实验内容设①p1=-2,p2=-30; ②p1=-2,p2=31.针对极点参数①②,画出系统零、极点分布图, 判断该系统稳定性.2.针对极点参数①②,绘出系统的脉冲响应曲线,并观察t→∞时, 脉冲响应变化趋势.3.针对极点参数①, 绘出系统的频响曲线.四、实验要求1.预习实验原理;2.对实验内容编写程序(M文件),上机运行;3.绘出实验内容的各相应曲线或图。
实验5--连续时间系统的复频域分析
实验5–连续时间系统的复频域分析实验背景在连续时间系统的频域分析中,复频域分析是非常重要的一个方法。
其可以帮助我们更直观地了解系统的频率响应,包括幅频响应和相频响应,对于系统的设计和优化都有非常实际的应用价值。
因此,在本次实验中,我们将通过对一个特定系统的复频域分析来学习这一方法的基本原理和操作流程。
实验目的1.了解连续时间系统的幅频响应和相频响应2.掌握利用MATLAB对系统进行复频域分析的方法3.学会根据复频域图像对系统进行分析和优化实验原理连续时间系统幅频响应和相频响应在连续时间系统的频域分析中,使用的是拉普拉斯变换。
通过对系统的输入信号和输出信号进行拉普拉斯变换,可以得到它们在复平面上的函数,进而求得系统的传递函数H(s):H(s)=Y(s)/X(s)其中,s为复变量。
系统的幅频响应和相频响应分别定义为:H(s)的模和相位:|H(jw)|=sqrt(H(s)H(s)*) (模) arg(H(jw))=tan^-1[Im{H(jw)}]/Re{H(jw)} (相位) 其中,w为实数,j为虚数单位。
利用MATLAB进行系统复频域分析MATLAB提供了众多用于连续时间系统复频域分析的工具。
其中,最基本的是bode命令。
它可以计算和绘制给定系统的幅频响应和相频响应曲线。
常用命令格式如下:[bode(H,w)]其中,H为系统的传递函数,w为频率范围除此之外,MATLAB还提供了很多其他的命令,如nyquist、margin、freqresp 等。
它们可以帮助我们更全面地分析系统的性能和特点。
实验步骤实验环境1.一台已安装MATLAB的计算机实验流程1.根据给定的系统传递函数H(s),利用MATLAB计算和绘制其幅频响应和相频响应曲线。
%定义系统传递函数H=tf([5+j*10 0.6+0.2*j],[1 2+j 3 4-j 5+j]);%绘制幅频响应和相频响应曲线figure(1)subplot(2,1,1)bode(H);subplot(2,1,2)nyquist(H);2.根据绘制的幅频响应和相频响应曲线,对系统进行分析和优化。
信号与系统实验实验报告
实验五连续系统分析一、实验目的深刻理解连续时间系统的系统函数在分析连续系统的时域特性、频域特性及稳定性中的重要作用及意义,掌握根据系统函数的零极点设计简单的滤波器的方法。
掌握利用MATLAB分析连续系统的时域响应、频响特性和零极点的基本方法。
二、实验原理MATLAB提供了许多可用于分析线性时不变连续系统的函数,主要包含有系统函数、系统时域响应、系统频域响应等分析函数。
三、实验内容1.已知描述连续系统的微分方程为,输入,初始状态,计算该系统的响应,并与理论结果比较,列出系统响应分析的步骤。
结果分析:理论值y(t)=0.8*exp(-10t)*u(t)+0.2程序运行出的结果与理论预期结果相差较大误差随时间增大而变小,初始值相差最大,而后两曲线基本吻合,表明该算法的系统响应在终值附近有很高的契合度,而在初值附近有较大的误差。
2.已知连续时间系统的系统函数为,求输入分别为,,时,系统地输出,并与理论结果比较。
实验结果:可见数值计算和理论计算曲线基本重合。
误差分析:可见误差小于0.001,计算值与理论值契合度很高。
3. 研究具有以下零极点的连续系统:(a) 1个极点s=—0.1,增益k=1。
(b) 1个极点s=0,增益k=1。
(c) 2个共轭极点,增益k=1。
(d) 2个共轭极点,增益k=1。
(e) 零点在,极点在,增益k=1。
(f) 零点在,极点在,增益k=1。
完成下列任务:(1) 利用zpk和tf命令建立系统的系统函数,画出系统的零极点图。
(2) 分析系统是否稳定。
若稳定,画出系统的幅频特性曲线。
(3) 画出系统的冲激响应波形。
(4) 详细列出根据零极点分析系统特性的过程。
实验代码:(a)结果分析:(a)~(e)均为因果稳定系统,他们的极点都在jw轴左侧。
当且仅当H(s)的全部极点都位于s平面的左半平面时,一个具有有理系统函数H(s)的因果系统才是稳定的。
(2023)304编号北京理工大学信号与系统实验报告5连续时间系统的复频域分析(一)
(2023)304编号北京理工大学信号与系统实验报告5连续时间系统的复频域分析(一)关于北京理工大学信号与系统实验报告5实验编号(2023)304实验名称信号与系统实验报告5:连续时间系统的复频域分析实验目的通过本实验,掌握连续时间系统的复频域分析的基本原理和方法,熟练掌握求解复频域中系统的幅度谱和相位谱的方法,提高对系统频域特性的认识。
实验内容实验内容主要分为以下几部分:1.实验仪器和元件的使用2.连续时间系统的频域分析方法3.MATLAB工具箱的应用实验过程1.使用示波器、函数发生器等实验仪器,搭建连续时间系统。
2.将系统的输入信号和输出信号从时域表示转换为复频域表示。
3.根据复频域表示求解系统的幅度谱和相位谱。
4.使用MATLAB工具箱验证实验结果。
5.分析实验结果,总结连续时间系统的频域特性。
实验结论通过本实验,我们了解了连续时间系统的复频域分析方法,熟悉了求解幅度谱和相位谱的步骤,并通过实验验证了所学内容的正确性。
同时,我们也认识到了连续时间系统在频域中的特性和应用范围。
实验感受本实验对我们的信号与系统学习提供了重要的实践环节,让我们更加深入地理解了信号与系统的频域分析方法,并对自己的专业兴趣产生了更深刻的认识。
同时,实验过程中我们也体验到了探索和解决问题的乐趣,收获了宝贵的经验。
实验注意事项1.实验中的电路连接应符合要求,注意仪器的使用和安全操作。
2.合理调整示波器、函数发生器等参数,以确保实验效果。
3.对于MATLAB工具箱的使用应具备一定的基础。
4.实验报告应准确记录实验过程中的操作、数据和可视化结果。
实验改进方向1.加强理论基础知识的学习,深入了解系统的频域特性。
2.进一步利用MATLAB及其他工具箱进行系统的分析和模拟,提高实验的精度和可靠性。
3.可尝试采用不同的连续时间系统进行分析和比较,从而更好地认识连续时间系统的特性。
实验意义通过本实验,我们对信号与系统的频域分析方法和连续时间系统的特性有了更深入的了解和认识。
连续时间信号与系统的频域分析实验报告(共9篇)
连续时间信号与系统的频域分析实验报告(共9篇)信号与系统实验五__连续时间信号的频域分析实验名称:连续时间信号的频域分析报告人:姓名班级学号一、实验目的1、熟悉傅里叶变换的性质;2、熟悉常见信号的傅里叶变换;3、了解傅里叶变换的MATLAB实现方法。
二、实验内容及运行结果1、编程实现下列信号的幅度频谱:(1)求出f(t)=u(2t+1)-u(2t-1)的频谱函数F(w);请与f1(t) u(2t+1)-u(2t-1)的频谱函数F1(w)进行比较,说明两者的关系。
%(1)f(t)=u(2t+1)-u(2t-1)与f(t)=u(t+1)-u(t-1) syms t w t1 w1Gt=sym('Heaviside(2*t+1)-Heaviside(2*t-1)');Gt1=sym('Heaviside(t1+1)-Heaviside(t1-1)');Fw=fourier(Gt,t,w);Fw1=fourier(Gt1,t1,w1);FFw=maple('convert',Fw,'piecewise');FFw1=maple('convert',Fw1,'piecewise');FFP=abs(FFw);FFP1=abs(FFw1);subplot(2,1,1);ezplot(FFP,[-10*pi 10*pi]);axis([-10*pi 10*pi 0 1.5]);subplot(2,1,2);ezplot(FFP1,[-10*pi 10*pi]);grid;axis([-10*pi 10*pi 0 2.2]);不同点:F1(w)的图像在扩展,幅值是F(w)的两倍。
(2)三角脉冲f2(t)=1-|t|;|t|=1;ft=sym('(1+t)*Heaviside(t+1)-2*t*Heaviside(t)+(t-1)*Heaviside( t-1)');Fw=fourier(ft);subplot(211)ezplot(abs(Fw)); g2)');ft=ifourier(Fw,w,t)ft =exp(-4*t)*heaviside(t)-exp(4*t)*heaviside(-t)(2)F(w)=((i*w)+5*i*w-8)/((i*w)+6*i*w+5)syms t wFw=sym('((i*w)+5*i*w-8)/((i*w)+6*i*w+5)');ft=ifourier(Fw,w,t)ft =dirac(t)+(2*exp(-5*t)-3*exp(-t))*heaviside(t)三、讨论与总论通过本实验,掌握了信号的傅里叶变换的性质以及方法,对傅里叶变换的性质有进一步的提高。
连续时间系统的复频域分析
信号与系统实验报告实验题目: 实验三:连续时间系统的复频域分析实验仪器: 计算机,MATLAB 软件101b s b a s a ++++++称为系统的特征多项式,征根,也称为系统的固有频率(或自然频率)。
为将个特征根,这些特征根称为()F s 极点。
根据求函数21()(1)F s s s =-的拉氏逆变换。
源代码:num = [1]; 结果为:r =-1 1 1 a=conv([1 -1],[1 -1]);den = conv([1 0], a); p =1 1 0 [r,p,k] = residue(num, den); k=03.示例3:求函数2224()(4)s F s s -=+的拉氏逆变换源代码:num = [1 0 -4];den = conv([1 0 4], [1 0 4]); [r,p,k] = residue(num, den);结果为:r =-0.0000-0.0000i 0.5000+0.0000i -0.0000+0.0000i 0.5000-0.0000ip =-0.0000+2.0000i -0.0000+2.0000i -0.0000-2.0000i -0.0000-2.0000i k=04.示例4:已知系统函数为:321()221H s s s s =+++,利用Matlab 画出该系统的零极点分布图,分析系统的稳定性,并求出该系统的单位冲激响应和幅频响应。
源代码: num=[1];den=[1 2 2 1]; sys=tf(num,den); poles=roots(den); figure(1);pzmap(sys);xlabel('Re(s)');ylabel(' Im(s)');title('zero-pole map'); t=0:0.02:10;h=impulse(num,den,t); figure(2);plot(t,h);xlabel('t(s)');ylabel('h(t)');title('Impulse Response'); [H,w]=freqs(num,den);figure(3);plot(w,abs(H));xlabel('\omega(rad/s)');ylabel('|H(j\omega)|');title('Magenitude Response'); 结果为:poles =-1.0000 -0.5000 + 0.8660i -0.5000 - 0.8660i (2) 已知象函数,试调用residue 函数完成部分分式分解,并写出逆变换。
系统的复频域分析实验报告
实验六 系统的复频域分析信号)(t x 的拉普拉斯变换⎰∞∞--=dt e t x s X st )()( (6.1) 是连续时间傅立叶变换地推广。
连续时间傅立叶变换在研究连续时间信号与系统中是很有用的。
然而,许多信号不存在傅立叶变换而存在拉普拉斯变换,这使得拉普拉斯变换成为线性时不变系统分析的一种有用方法。
对一大类信号来说,它们的拉普拉斯变换可以表示为s 的多项式之比,即)()()(s D s N s X =这里)(s N 和)(s D 分别称作分子和分母多项式。
能表示成多项式之比的变换称为有理变换,这里作为满足线性常系数微分方程的LTI 系统的系统函数中常常出现。
除了一个标量因子外,有理变换是完全由多项式)(s N 和)(s D 的根决定的,这些根分别称为零点和极点。
由于这些根在LTI 系统的研究中起着重要的作用,所以它们以零极点图的方式展现出来的是很方便的。
这一章将用拉普拉斯变换在复频域研究LTI 系统的一些性质。
§6.1 MATLAB 函数lsim (用于系统函数)目的用lsim 仿真由系统函数表征的因果LTI 系统的输出。
相关知识第二章所讨论的是如何用lsim 命令仿真一个输出满足一个线性常系数微分方程的因果LTI 连续时间系统。
因为系统函数唯一地表征了关联系统输入和输出的微分方程。
所以由系统函数表征的因果LTI 系统的输出也能够用lsim 仿真。
如果系统函数给出如下形式:)()1()1()()1()1()(N a s N a s a M b s M b s b s H N M +-+++-++= (6.2)那么,对输入地系统)(t x 的系统输出就能用lsim(b,a,x,t)仿真,其中MATLAB 向量b 和a 包含了分子分母s 多项式的系数。
例如,考虑系统函数221)(-+=s s s H ,其系数由向量b=[1 1/2]和a=[1 -2]定义。
命令lsim(b,a,x,t)将系统对由向量x 给出的输入,在向量t 给定的时刻上的系统的时间响应应存入向量y 中,向量y 和输入向量x 有相同数量的元素。
信号与系统实验五 连续线性时不变系统分析
信号与系统实验陈述课程名称:信号与系统实验实验项目名称:连续线性时不变系统分析专业班级:姓名:学号:完成时间:年月日一、实验目的1.掌握连续LTI系统的单位冲激响应、单位阶跃响应和任意激励对应响应的求解方法。
2.掌握连续LTI系统的频域分析方法。
3.掌握连续LTI系统的复频域分析方法。
4.掌握连续LTI系统的时域、频域和复频域分析方法的相互转换。
二、实验原理1.连续LTI系统的时域分析(1)连续线性时不变系统的描述设连续线性时不变系统的激励为,响应为,则描述系统的微分方程可暗示为为了在Matlab编程中调用有关函数,我们可以用向量和来暗示该系统,即这里要注意,向量和的元素排列是按微分方程的微分阶次降幂排列,缺项要用0补齐。
(2) 单位冲激响应单位冲激响应是指连续LTI系统在单位冲激信号激励下的零状态响应,因此满足线性常系数微分方程(5.1)及零初始状态,即,依照定义,它也可暗示为对于连续LTI系统,若其输入信号为,冲激响应为,则其零状态响应为可见,能够刻画和表征系统的固有特性,与何种激励无关。
一旦知道了系统的冲激响应,就可求得系统对任何输入信号所发生的零状态响应。
Matlab提供了专门用于求连续系统冲激响应的函数impulse(),该函数还能绘制其时域波形。
(3)单位阶跃响应单位阶跃响应是指连续LTI系统在单位阶跃信号激励下的零状态响应,它可以暗示为Matlab提供了专门用于求连续系统单位阶跃响应的函数step( ),该函数还能绘制其时域波形。
(4)任意激励下的零状态响应已经知道,连续LTI系统可用常系数线性微分方程(5.1)式来描述,Matlab提供的函数lsim( )能对上述微分方程描述的连续LTI系统的响应进行仿真,该函数不但能绘制指定时间范围内的系统响应波形图,而且还能求出系统响应的数值解。
其调用格式有lsim(b,a,x,t)y=lsim(b,a,x,t) :只求出系统的零状态响应的数值解,而不绘制响应曲线需要特别强调的是,Matlab总是把由分子和分母多项式暗示任何系统都当作是因果系统。
实验五 信号与系统复频域分析
信号与系统实验报告实验五信号与系统复频域分析学院专业班级姓名学号指导教师实验报告评分:_______实验五 信号与系统复频域分析一、 实验目的1.学会用MATLAB 进行部分分式展开;2.学会用MATLAB 分析LTI 系统的特性;3.学会用MATLAB 进行Laplace 正、反变换。
二、实验内容与步骤1、求信号)()(3t u te t f t -=的拉普拉斯变换。
f=sym(' t*exp(-3*t)'); F=laplace(f) 实验结果: F =1/(s+3)^22、求函数23795)(223+++++=s s s s s s F 的部分分式展开式,以及反变换。
format rat; num=[1,5,9,7]; den=[1, 3,2];[r,p]=residue(num,den) 实验结果: r =-1 2 p =-2 -1 反变换:F=sym('(s^3+5*s^2+9*s+7)/(s^2+3*s+2)'); ft=ilaplace(F)ft =dirac(1,t)+2*dirac(t)-exp(-2*t)+2*exp(-t)4、已知一个因果系统的系统函数为)2)(1()(++=s s ss H ,该系统的零点和极点分别位于?从时域和零极点分布特征两个方面说明该系统是否是稳定的系统?从频率响应特性上看,该系统具有何种滤波特性? num=[1,0]; den=[1,3,2]; sys=tf(num,den); figure(1);pzmap(sys); title('零极点分布图') t=-10:0.02:10;f=t./(t.^2 + 3* t + 2); figure(2);plot(t,f) xlabel('s') ylabel('H(s)')title('H(s)的时域波形图') [H,w]=freqs(num,den); figure(3);plot(w,abs(H)) xlabel('\omega') title('频率响应图')5、输入因果的系统函数12211)(232++++=s s s s a s H ① 此处a 取1,执行程序。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
信号与系统
实验报告
实验五连续系统的复频域分析
实验五连续系统的复频域分析
一、实验目的
1. 深刻理解拉普拉斯变换、逆变换的定义,掌握用MATLAB实现拉普拉斯变换、逆变换的方法。
2会求几种基本信号的拉氏变换。
3 掌握用MATLAB绘制连续系统零、极点的方法。
4 求解系统函数H(s)。
二
1已知连续时间信号f(t)=sin(t)u(t)、求出该信号的拉普拉斯变换,并用MATLAB
绘制拉普拉斯变换的曲面图。
syms t;
ft=sin(t)*heaviside(t);
Fs=Laplace(ft);
a=-0.5:0.08:0.5;
b=-2:0.08:2;
[a,b]=meshgrid(a,b);
c=a+i*b;
d=ones(size(a));
c=c.*c;
c=c+d;
c=1./c;
c=abs(c);
mesh(a,b,c);
surf(a,b,c)
axis([-0.5,0.5,-2,2,0,10])
colormap(hsv
)
2求[(1-e^(-at))]/t的拉氏变换。
syms t s a
f1=(1-exp(-a*t))/t;
F=laplace(f1,t,s)
F =
log(s+a)-log(s)
3求F(s)=-log(s)+ log(s+a)的拉氏逆变换syms t s a
F =log(s+a)-log(s);
f1=ilaplace(F,s,t)
f1 =
(1-exp(-a*t))/t
4已知某连续系统的系统函数为:
H(s)=(s^2+3s+2)/(8s^4+2s^3+3s^2+5)试用MATLAB求出该系统的零极点,画出零极点分布图。
b=[1 3 2];
a=[8 2 3 0 5];
zs=roots(b);
ps=roots(a);
hold on
plot(real(zs),imag(zs),'o');
plot(real(ps),imag(ps),'x');
grid
axis([-2.5,1,-1,1])
5已知H(s)=(s+1)/(s^2+s+1),绘制阶跃响应图形,冲激响应图形,频率激响应图形。
syms t s
H=(s+1)/(s^2+s+1);
f1=ilaplace(H,s,t);
f2=heaviside(t);
f3=dirac(t);
H1=laplace(f2,t,s);
H2=laplace(f3,t,s);
H4=H1*H;
H5=H2*H;
f4=ilaplace(H4,s,t);
f5=ilaplace(H5,s,t);
subplot(2,1,1);ezplot(f4);grid on
subplot(2,1,2);ezplot(f5);grid on
得到曲线图如下第一个位阶跃响应图,第二个为冲激响应图
在求频率响应,分为幅频和相频,输入以下程序
b=[1 0];
a=[1 1 1];
p=roots(a);
pxm=max(real(p));
if pxm>=0
xitongbuwending
else
freqs(b,a)
end
得到如下图形,第一个为幅频图,第二个为相频图
三、实验仪器
计算机一台,软件MATLAB7.6
四、心得体会
1要熟悉拉普拉斯变换、逆变换的定义,掌握用MATLAB实现拉普拉斯变换、逆变换的方法。
2会用MA TLAB绘制连续系统零、极点的方法。