4_15 四种信号的时域频域对应关系

实验五 信号与系统的复频域分析

实验五 信号与系统的复频域分析 王靖 08通信 12号 实验目的 （1）掌握利用MA TLAB 进行连续时间信号与系统的复频域分析。 （2）掌握利用MA TLAB 进行离散系统的复频域分析。 实验环境 安装MATLAB7.0以上版本的计算机 实验内容 1. 利用help 命令了解以下命令的基本用法 residue ，roots ，pzmap ，cart2pol ，residuez ，tf2zp ，zplane 2. 部分分式展开的MATLAB 实现 用部分分式展开法求X(s)的反变换。 2321 ()452s X s s s s +=+++ 步骤一：建立新的m 文件，保存并命名为program1.m 。 步骤二：输入以下命令，理解每条命令的含义。 %program1，部分分式展开法求反变换 [10 1];[1452];[,,](,) n u m d en r p k resid u e n u m d en === 步骤三：保存程序并运行，记录得到的结果。 如右图所示 步骤四：由得到的结果可以直接获得X(s)展开表示式 25 4 2 ()21(1)X s s s s =-++++： 步骤五：由此可得到X(s)反变换的原函数，记录。 X(t)=(5exp(-2*t)-4exp(-t)+2texp(-t)) 思考：将其转换成极坐标形式，应该如何使用cart2pol 命令？离散系统的部分分式展开，如何使用命 令residuez ，得到的结果如何利用？ 将笛卡尔坐标转化为极坐标用 [angle,mag]=cart2pol(real(r),imag(r)) [r,p,k] = residuez(nun,,den)

北京交通大学信号与系统时域分析

北京交通大学信号与系统时域分析

【研讨题目2】 信号与系统时域分析专题研讨 【目的】 1.研究用离散方法近似计算连续信号的卷积积分； 2.通过分析近似计算卷积积分过程中出现的问题，锻炼学生分析问题和解决问题的能力； 【知识点】 信号时域分析，卷积积分，卷积和 【研讨题目】连续信号卷积积分的数值近似计算 两个连续信号的卷积积分定义为 τττd )()()(-= ?∞ ∞ -t h x t y 为了能用数值方法进行计算，需对连续信号进行抽样。记x [k ]=x (k ), h [k ]=h (k ),为 进行数值计算所选定的抽样间隔，可以证明连续信号卷积积分可近似的表示为 (Δ)Δ([][])y k x k h k ≈?* (1) 由式(1)可知，可以利用Matlab 提供的conv 函数近似计算连续信号的卷积积分。 一、(*)理论分析 为了对近似计算的结果进行分析，用解析的方法计算下列卷积积分，推出卷积积分的解析表达式； (1) 时限信号卷积积分 x 1(t )=u (t )-u (t -1)，y 1(t )=x 1(t )*x 1(t ); 卷积结果为：y1(t)= x 1(t )*x 1(t )=r(t)-2*r(t-1)+r(t-2) (2) 分段常数信号卷积积分 x 2(t )= x 1(t )+2 x 1(t -1)+ x 1(t -2)，h 2(t )= x 1(t )- x 1(t -1)， y 2(t )=x 2(t )*h 2(t ); 卷积结果为：y2(t)= x 2(t )*h 2(t ) =y1(t)+y1(t-1)-y1(t-2)-y1(t-3) =r(t)-r(t-1)-2*r(t-2)+2*r(t-3)+r(t-4)-r(t-5) (3) 非时限信号卷积积分 x 3(t )=u (t )，h 3(t )=e -t u (t ), y 3(t )=x 3(t )*h 3(t ) 卷积结果为：y3= x 3(t )*h 3(t ) =[1-exp(-t)]*u(t) 二、(*)时限信号卷积积分的近似计算 取不同的△值，用Matlab 函数conv 近似计算卷积积分y 1(t )并画出其波形，讨论的取值对计算结果的影响。

信号时域频域及其转换

信号分析方法概述： 通用的基础理论是信号分析的两种方法：1 是将信号描述成时间的函数 2 是将信号描述成频率的函数。也有用时域和频率联合起来表示信号的方法。时域、频域两种分析方法提供了不同的角度，它们提供的信息都是一样，只是在不同的时候分析起来哪个方便就用哪个。 思考： 原则上时域中只有一个信号波（时域的频率实际上是开关器件转动速度或时钟循环次数，时域中只有周期的概念），而对应频域（纯数学概念）则有多个频率分量。 人们很容易认识到自己生活在时域与空间域之中（加起来构成了三维空间），所以比较好理解时域的波形（其参数有：符号周期、时钟频率、幅值、相位）、空间域的多径信号也比较好理解。 但数学告诉我们，自己生活在N维空间之中，频域就是其中一维。时域的信号在频域中会被对应到多个频率中，频域的每个信号有自己的频率、幅值、相位、周期（它们取值不同，可以表示不同的符号，所以频域中每个信号的频率围就构成了一个传输信道。 时域中波形变换速度越快（上升时间越短），对应频域的频率点越丰富。 所以：OFDM中，IFFT把频域转时域的原因是：IFFT的输入是多个频率抽样点（即各子信道的符号），而IFFT之后只有一个波形，其中即OFDM符号，只有一个周期。 时域 时域是真实世界，是惟一实际存在的域。因为我们的经历都是在时域中发展和验证的，已经习惯于事件按时间的先后顺序地发生。而评估数字产品的性能时，通常在时域中进行分析，因为产品的性能最终就是在时域中测量的。 时钟波形的两个重要参数是时钟周期和上升时间。 时钟周期就是时钟循环重复一次的时间间隔，通产用ns度量。时钟频率Fclock，即1秒钟时钟循环的次数，是时钟周期Tclock的倒数。 Fclock=1/Tclock 上升时间与信号从低电平跳变到高电平所经历的时间有关，通常有两种定义。一种是10-90上升时间，指信号从终值的10%跳变到90%所经历的时间。这通常是一种默认的表达方式，可以从波形的时域图上直接读出。第二种定义方式是20-80上升时间，这是指从终值的20%跳变到80%所经历的时间。 时域波形的下降时间也有一个相应的值。根据逻辑系列可知，下降时间通常要比上升时间短一些，这是由典型CMOS输出驱动器的设计造成的。在典型的输出驱动器中，p管和n 管在电源轨道Vcc和Vss间是串联的，输出连在这个两个管子的中间。在任一时间，只有一个晶体管导通，至于是哪一个管子导通取决于输出的高或低状态。 假设周期矩形脉冲信号f(t)的脉冲宽度为τ，脉冲幅度为E,重复周期为T，

北京理工大学信号与系统实验报告5-连续时间系统的复频域分析

北京理工大学信号与系统实验报告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 的有理函数 110 110 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++= +++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下： r=roots(c)，c 为多项式的系数向量，返回值r 为多项式的根向量。 求取零极点以及绘制系统函数的零极点分布图可以采用pzmap 函数，调用格式如下： pzmap(sys)绘出由系统模型sys 描述的系统的零极点分布图。 [p,z]=pzmap(sys)这种调用方式返回极点与零点，不绘出零极点分布图。 还有两个专用函数tf2zp 和zp2tf 可实现系统的传递函数模型和零极点增益模型的转换。调用格

离散信号与系统时域分析

目录 第1章设计任务及要求 (1) 1.1课程设计内容 (1) 1.2课程设计要求 (1) 第2章设计原理 (2) 2.1离散信号与系统的时域分析设计 (2) 2.1.1描写系统特性的方法介绍 (2) 2.1.2系统的时域特性 (2) 第3章设计实现 (3) 3.1实验内容与方法 (3) 3.1.1实验内容 (3) 第4章设计结果及分析 (3) 4.1程序设计结果及分析 (4) 总结 (7) 参考文献： (7) 附录： (8)

第1章 设计任务及要求 1.1课程设计内容 编制Matlab 程序，完成以下功能，产生系统输入信号；根据系统差分方程求解单位脉冲响应序列；根据输入信号求解输出响应；用实验方法检查系统是否稳定；绘制相关信号的波形。具体要求如下： (1) 给定一个低通滤波器的差分方程为 ()0.05()0.05(1)0.9(1)y n x n x n y n =+-+- 输入信号分别为182()=()()()x n R n x n u n =， ① 分别求出系统响应，并画出其波形。 ② 求出系统的单位脉冲响应，画出其波形。 (2) 给定系统的单位脉冲响应为1102()=()()() 2.5(1) 2.5(2)(3)h n R n h n n n n n δδδδ=+-+-+-，用线性卷积法求18()=()x n R n 分别对系统h1(n)和h2(n)的输出响应，并画出波形。 (3) 给定一谐振器的差分方程为() 1.8237(1)-0.9801(2)()(2)o o y n y n y n b x n b x n =--++-令b0=1/100.49,谐振器的谐振频率为0.4rad 。 1) 用实验方法检查系统是否稳定。输入信号为u(n)时,画出系统输出波形。 2) 给定输入信号为()=sin(0.014)sin(0.4)x n n n +求出系统的输出响应，并画出其波形。 1.2课程设计要求 1. 要求独立完成设计任务。 2. 课程设计说明书封面格式要求见《天津城市建设学院课程设计教学工作规范》附表1 3. 课程设计的说明书要求简洁、通顺，计算正确，图纸表达内容完整、清楚、规范。 4. 简述离散系统时域分析方法和通过实验判断系统稳定性的方法；完成以上设计实验并对结果进行分析和解释；打印程序清单和要求画出的信号波形；写出本次课程设计的收获和体会。 5. 课设说明书要求： 1) 说明题目的设计原理和思路、采用方法及设计流程。 2) 详细介绍运用的理论知识和主要的Matlab 程序。 3) 绘制结果图形并对仿真结果进行详细的分析。

信号时域与频域分析

信号时域与频域分析 实验报告 姓名：杨 班级：机械 学号： 213

实验数据中，电机转速为1200r/min，采样频率为1280Hz。Hz3为X位移振幅数据，Hz4为Y位移振幅数据，Hz5为速度振幅数据。 Matlab中信号特征对应函数编程 ma = max(Hz) %最大值 mi = min(Hz) %最小值 me = mean(Hz) %平均值 pk = ma-mi %峰-峰值 va = var(Hz); %方差 st = std(Hz); %标准差 ku = kurtosis(Hz); %峭度 rm = rms(Hz); %均方根 一、X轴位移测量分析 plot(Fs3,Hz3)时域图： ma =52.0261 mi =56.7010 me =1.8200 pk =108.7271 va =1.3870e+03 st =37.2431 ku =1.5462 rm =37.2693 频域图： fs=1280; x=Hz3; N=length(Hz3); df=fs/N; f=0:df:N*df-df; y=fft(x); y=abs(y)*2/N; figure(1); plot(f,y); xlabel('频率/Hz') ylabel('幅值') 频谱幅值取得最大值51.9847um，频率为20Hz，与电机转速对应频率一致，应为电机轴未动平衡所致；二倍频处有较大振幅，可能为轴承间隙过大所致。

二、Y轴位移测量分析 plot(Fs4,Hz4)时域图： ma =61.3987 mi =-74.6488 me =-1.1948 pk =136.0475 av =42.6109 va =2.2428e+03 st =47.3582 ku =1.5135 rm =47.3501 频域图： fs=1280; x=Hz4; N=length(Hz4); df=fs/N; f=0:df:N*df-df; y=fft(x); y=abs(y)*2/N; figure(1); plot(f,y); xlabel('频率/Hz') ylabel('幅值') 频谱幅值取得最大值66.6319um，频率为20Hz，与电机转速对应频率一致，应为电机轴未动平衡所致；二倍频处有较大振幅，可能为轴承间隙过大所致。

时域与频域的matlab程序

一．典型连续信号和离散信号的时域波形。 1．单边指数信号)()(t u Ae t y t α=； 2．单位冲激信号)()(0t t t y +=δ； 3．单位阶跃信号)()(0t t u t y +=； 4．矩形脉冲信号)]()([)(21t t u t t u A t y +-+?=； 5．正弦信号)()sin()(t u t A t y ω?=； 6．单位序列)()(0n n n y +=δ； 7．单位阶跃序列)()(0n n u n y +=； 8．单位矩形序列)()()(21n n u n n u n y +-+=； 9．指数序列)()(n u a A n y n ?=； 10．正弦序列)()sin()(n u n A n y ω?=。

单边指数信号 function zhishu(A,a,t1,t2,dt) t1=0 t2=10 A=1 A=-0.4 dt=0.01 t=t1:dt:t2; y=A*exp(a*t); plot(t,y) axis([t1,t2,0,1.2]) xlabel('t') ylabel('y(t)') title(' 单边指数信号') 单位冲激信号 function chongji(t1,t2,t0) dt=0.01; t1=10; t2=-5; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(-t0-t1)/dt+1)=1/dt; stairs(t,x); axis([t1,t2,0,1.2/dt]) xlabel('t') ylabel('y(t)') title('单位冲激信号')

实验六-信号与系统复频域分析

实验六信号与系统复频域分析 一、实验目的 1.学会用MATLAB进行部分分式展开； 2.学会用MATLAB分析LTI系统的特性； 3.学会用MATLAB进行Laplace正、反变换。 4.学会用MATLAB画离散系统零极点图； 5.学会用MATLAB分析离散系统的频率特性； 二、实验原理及内容 1．用MATLAB进行部分分式展开 用MATLAB函数residue可以得到复杂有理分式F(s)的部分分式展开式，其调用格式为 其中，num,den分别为F(s)的分子和分母多项式的系数向量，r为部分分式的系数，p为极点，k为F(s)中整式部分的系数，若F(s)为有理真分式，则k为零。 例6-1 用部分分式展开法求F(s)的反变换 解:其MATLAB程序为 format rat; num=[1,2]; den=[1,4,3,0]; [r,p]=residue(num,den) 程序中format rat是将结果数据以分数形式显示

F(s)可展开为 210.536()13 F s s s s --=++++ 所以，F(s)的反变换为 3211()()326t t f t e e u t --??=--???? 2．用MATLAB 分析LTI 系统的特性 系统函数H （s ）通常是一个有理分式，其分子和分母均为多项式。计算H （s ）的零极点可以应用MATLAB 中的roots 函数，求出分子和分母多项式的根，然后用plot 命令画图。 在MATLAB 中还有一种更简便的方法画系统函数H （s ）的零极点分布图，即用pzmap 函数画图。其调用格式为 pzmap(sys) sys 表示LTI 系统的模型，要借助tf 函数获得，其调用格式为 sys=tf(b,a) 式中，b 和a 分别为系统函数H （s ）的分子和分母多项式的系数向量。 如果已知系统函数H （s ），求系统的单位冲激响应h(t)和频 率响应H ω（j ）可以用以前介绍过的impulse 和freqs 函数。 例6-2 已知系统函数为 321221 s s s +++H(s)= 试画出其零极点分布图，求系统的单位冲激响应h(t)和频率响应H ω（j ），并判断系统是否稳定。 解：其MATLAB 程序如下： num=[1];

信号与系统实验报告实验一 信号与系统的时域分析

实验一信号与系统的时域分析 一、实验目的 1、熟悉与掌握常用的用于信号与系统时域仿真分析的MA TLAB函数; 2、掌握连续时间与离散时间信号的MA TLAB产生,掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MA TLAB编程; 3、牢固掌握系统的单位冲激响应的概念,掌握LTI系统的卷积表达式及其物理意义,掌握卷积的计算方法、卷积的基本性质; 4、掌握利用MA TLAB计算卷积的编程方法,并利用所编写的MA TLAB程序验证卷积的常用基本性质; 掌握MA TLAB描述LTI系统的常用方法及有关函数,并学会利用MATLAB求解LTI系统响应,绘制相应曲线。 基本要求:掌握用MA TLAB描述连续时间信号与离散时间信号的方法,能够编写MATLAB程序,实现各种信号的时域变换与运算,并且以图形的方式再现各种信号的波形。掌握线性时不变连续系统的时域数学模型用MA TLAB描述的方法,掌握卷积运算、线性常系数微分方程的求解编程。 二、实验原理 信号(Signal)一般都就是随某一个或某几个独立变量的变化而变化的,例如,温度、压力、 声音,还有股票市场的日收盘指数等,这些信号都就是随时间的变化而变化的,还有一些信号,例如在研究地球结构时,地下某处的密度就就是随着海拔高度的变化而变化的。一幅图片中的每一个象素点的位置取决于两个坐标轴,即横轴与纵轴,因此,图像信号具有两个或两个以上的独立变量。 在《信号与系统》课程中,我们只关注这种只有一个独立变量(Independent variable)的信号,并且把这个独立变量统称为时间变量(Time variable),不管这个独立变量就是否就是时间变量。 在自然界中,大多数信号的时间变量都就是连续变化的,因此这种信号被称为连续时间信号(Continuous-Time Signals)或模拟信号(Analog Signals),例如前面提到的温度、压力与声音 信号就就是连续时间信号的例子。但就是,还有一些信号的独立时间变量就是离散变化的,这种信号称为离散时间信号。前面提到的股票市场的日收盘指数,由于相邻两个交易日的日收盘指数相隔24小时,这意味着日收盘指数的时间变量就是不连续的,因此日收盘指数就是离散时间信号。 而系统则用于对信号进行运算或处理,或者从信号中提取有用的信息,或者滤出信号中某些无用的成分,如滤波,从而产生人们所希望的新的信号。系统通常就是由若干部件或单元组成的一个整体(Entity)。系统可分为很多不同的类型,例如,根据系统所处理的信号的不同,系统可分为连续时间系统(Continuous-time system)与离散时间系统(Discrete-time system),根

肌电信号的时域和频域分析

肌电信号的时域和频域分析 摘要：肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位在时间和空间上的叠加，反映了神经，肌肉的功能状态，在基础医学研究、临床诊断和康复工程中有广泛的应用。 其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌 电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测 方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电 极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单，病人易接受，有着广泛的应用前景。 本次设计基于matlab用小波变换对肌电信号进行消噪处理，分别选用20N 的肌电信号数据和50N的肌电数据进行对比，最后在GUI界面上完成相应的功能处理。 关键字：肌电信号 Matlab 小波去噪 GUI 第一章绪论 肌电信号是产生肌肉力的电信号根源，它是肌肉中很多运动单元动作电位在时间和空间上的叠加，反映了神经，肌肉的功能状态，在基础医学研究、临床诊断和康复工程中有广泛的应用。 其种类重要有两种：一，临床肌电图检查多采用针电极插入肌肉检测肌电图，其优点是干扰小，定位性好，易识别，但由于它是一种有创伤的检测方法，其应用收到了一定的限制。二，表面肌电则是从人体皮肤表面通过电极记录下来的神经肌肉活动时发放的生物电信号，属于无创伤性，操作简单，病人易接受，有着广泛的应用前景。 肌电信号本身是一种较微弱的电信号。检测和记录表面肌电信号，需要考虑的主要问题是尽量消除噪声和干扰的影响, 提高信号的保真度[1]。

第二章肌电信号的时域分析 2.1 肌电信号时域图的显示及比较 肌电信号采用两个不同的数据进行比较，通过比较时域图及其特性来进行分析[2]。其图像如下所示： 如上图所示：肌电数据分别是同一个体在20N的力和50N的力所反映的图像。可以看出在不同作用力时，其图像的差别很大。 2.2 时域参数 2.2.1 均值 对于一个随机变量来说，均值是一个很重要的数值特征。粗略的说，就是来描述一个群体的平均水平。其严格的数学定义非常的简单，就是一个随机变量关于概率测度的积分。这样的积分在测度轮或者实分析里是没有什么直观的解释的。而在概率论里却成为了一个群体的主要指标。在此处，均值表示肌电信号的平均水平。 2.2.2 标准差 标准差（Standard Deviation），也称均方差（mean square error），是各数据偏离平均数的距离的平均数，它是离均差平方和平均后的方根，用σ表示。

频域和时域的关系

信号的频域 在电子学、控制系统及统计学中，频域是指在对函数或信号进行分析时，分析其和频率有关部份，而不是和时间有关的部份，和时域一词相对。函数或信号可以透过一对数学的运算子在时域及频域之间转换。例如傅里叶变换可以将一个时域信号转换成在不同频率下对应的振幅及相位，其频谱就是时域信号在频域下的表现，而反傅里叶变换可以将频谱再转换回时域的信号。 以信号为例，信号在时域下的图形可以显示信号如何随着时间变化，而信号在频域下的图形（一般称为频谱）可以显示信号分布在哪些频率及其比例。频域的表示法除了有各个频率下的大小外，也会有各个频率的相位，利用大小及相位的资讯可以将各频率的弦波给予不同的大小及相位，相加以后可以还原成原始的信号。在频域的分析中，常会用频谱分析仪来将实际的信号转换为频域下的频谱。 频域，尤其在射频和通信系统中运用较多，在高速数字应用中也会遇到频域。频域最重要的性质是：它不是真实的，而是一个数学构造。时域是惟一客观存在的域，而频域是一个遵循特定规则的数学范畴。 正弦波是频域中唯一存在的波形，这是频域中最重要的规则，即正弦波是对频域的描述，因为时域中的任何波形都可用正弦波合成。这是正弦波的一个非常重要的性质。然而，它并不是正弦波的独有特性，还有许多其他的波形也有这样的性质。正弦波有四个性质使它可以有效地描述其他任一波形： （1）时域中的任何波形都可以由正弦波的组合完全且惟一地描述。 （2）任何两个频率不同的正弦波都是正交的。如果将两个正弦波相乘并在整个时间轴上求积分，则积分值为零。这说明可以将不同的频率分量相互分离开。（3）正弦波有精确的数学定义。 （4）正弦波及其微分值处处存在，没有上下边界。 使用正弦波作为频域中的函数形式有它特别的地方。若使用正弦波，则与互连线的电气效应相关的一些问题将变得更容易理解和解决。如果变换到频域并使用正弦波描述，有时会比仅仅在时域中能更快地得到答案。 而在实际中，首先建立包含电阻，电感和电容的电路，并输入任意波形。一般情况下，就会得到一个类似正弦波的波形。而且，用几个正弦波的组合就能很容易地描述这些波形。 许多物理元件的特性会随着输入讯号的频率而改变，例如电容在低频时阻抗变大，高频时阻抗变小，而电感恰好相反，高频时阻抗变大，低频时阻抗变小。一个线性非时变系统的特性也会随频率而变化，因此也有其频域下的特性，频率响应的图形即为其代表。频率响应可以视为是一个系统在输入信号振幅相同、频率不同时，其输出信号振幅的变化，可以看出系统在哪些频率的输出较大。有些系统的定义就是以频域为主，例如低通滤波器只允许低于一定频率的讯号通过。 不论是进行拉普拉斯转换、Z转换或是傅立叶变换，其产生的频谱都是一个频率的复变函数，表示一个信号（或是系统的响应）的振幅及其相位。不过在许多的应用中相位的资讯并不重要，若不考虑相位的资讯，都可以将频谱的资讯只以不同频率下的振幅（或是功率密度）来表示。 功率谱密度是一种常应用在许多非周期性也不满足平方可积性（square-integrable）讯号的频域表示法。只要一个讯号是符合广义平稳随机过程的输出，就可以计算其对应的功率谱密度。 时域分析与频域分析是对模拟信号的两个观察面。时域分析是以时间轴为坐标表示动态信号的关系；频域分析是把信号变为以频率轴为坐标表示出来。一般来说，时域的表示较为形象与直观，频域分析则更为简练，剖析问题更为深刻和方便。目前，信号分析的趋势是从时域

北京理工大学信号与系统实验实验5连续时间系统地复频域分析报告报告材料

实验5 连续时间系统的复频域分析 一、实验目的 1.掌握拉普拉斯变换及其反变换的定义，并掌握MATLAB 实现方法。 2.学习和掌握连续时间系统系统函数的定义及复频域分析方法。 3.掌握系统零极点的定义，加深理解系统零极点分布与系统特性的关系。 二、实验原理与方法 1.拉普拉斯变换 连续时间信号)(t x 的拉普拉斯变换定义为 )1.....(..........)()(dt e t x s X st ? +∞ ∞ --= 拉普拉斯反变换定义为 )2....(..........)(21)(ds e s X j t x j j st ?∞ +∞ -=σσπ 在MATLAB 中，可以采用符号数学工具箱的laplace 函数和ilaplace 函数进行拉氏变换和反拉氏变换。 L=laplace(F)符号表达式F 的拉氏变换，F 中时间变量为t ，返回变量为s 的结果表达式。 L=laplace(F,t)用t 替换结果中的变量s 。 F=ilaplace(L)以s 为变量的符号表达式L 的拉氏反变换，返回时间变量为t 的结果表达式。 F=ilaplace(L,x)用x 替换结果中的变量t 。 除了上述ilaplace 函数，还可以采用部分分式法，求解拉普拉斯逆变换，具体原理如下： 当 X (s )为有理分式时，它可以表示为两个多项式之比： )3.(..........)()()(0 110 11a s a s a b s b s b s D s N s X N N N N M M M M +?+++?++==---- 式（3）可以用部分分式法展成一下形式 )4.....(.............)(2211N N p s r p s r p s r s X -++-+-= 通过查常用拉普拉斯变换对，可以由式（1-2）求得拉普拉斯逆变换。 利用 MATLAB 的residue 函数可以将 X (s )展成式（1-2）所示的部分分式展开式，该 函数的调用格式为：[r,p,k] = residue(b,a) 其中b 、a 为分子和分母多项式系数向量，r 、p 、k 分别为上述展开式中的部分分式系数、极点和直项多项式系数。 2.连续时间系统的系统函数

用Matlab进行信号与系统的时、频域分析

课程实验报告 题目：用Matlab进行 信号与系统的时、频域分析 学院 学生姓名 班级学号 指导教师 开课学院 日期 用Matlab进行信号与系统的时、频域分析 一、实验目的 进一步了解并掌握Matlab软件的程序编写及运行； 掌握一些信号与系统的时、频域分析实例； 了解不同的实例分析方法，如：数值计算法、符号计算法； 通过使用不同的分析方法编写相应的Matlab程序； 通过上机，加深对信号与系统中的基本概念、基本理论和基本分析方法的理解。 二、实验任务 了解数值计算法编写程序，解决实例； 在Matlab上输入三道例题的程序代码，观察波形图； 通过上机实验，完成思考题； 完成实验报告。 三、主要仪器设备

硬件：微型计算机 软件：Matlab 四、 实验内容 （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 *=的波形。 程序代码： 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 numberical','Theoretical'); ylabel('y=x*h'); xlabel('----t/s'); （3）求周期矩形脉冲信号的频谱图，已知s T s A 5.0,1.0,1===τ

连续系统的复频域分析及MATLAB 实现 - 信号与系统实验报告

计算机与信息工程学院设计性实验报告 一、实验目的 1.掌握用matlab 分析系统时间响应的方法 2.掌握用matlab 分析系统频率响应的方法 3.掌握系统零、极点分布与系统稳定性关系 二、实验原理 1.系统函数H(s) 系统函数:系统零状态响应的拉氏变换与激励的拉氏变换之比. H(s)=R(s)/E(s) 在matlab 中可采用多种方法描述系统,本文采用传递函数(系统函数)描述法. 在matlab 中, 传递函数描述法是通过传递函数分子和分母关于s 降幂排列的多项式系数来表示的.例如,某系统传递函数如下 )1(8 .03.11 )(2+++=s s s s H 则可用如下二个向量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 分析系统频率响应特性 频响特性: 系统在正弦激励下稳态响应随信号频率变化的特性. ()()() ()j s j H j H s H j e φωω ωω=== |H(j ω)|:幅频响应特性. ?(ω):相频响应特性(或相移特性).

Matlab 求系统频响特性函数freqs 的调用格式: h=freqs(num,den,ω) ω:为等间隔的角频率向量,指明要计算响应的频率点. 4.系统零、极点分布与系统稳定性关系 系统函数H(s)集中表现了系统的性能,研究H(s)在S 平面中极点分布的位置,可很方面地判断系统稳定性. 1) 稳定系统: H(s)全部极点落于S 左半平面(不包括虚轴),则可以满足 0)]([lim =∞ →t h t 系统是稳定的. 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()命令即可画出系统零、极点分布图,进而分析判断系统稳定性. 三、实验内容 设()(1)(2) s H s s p s p = -- 设①p1=-2,p2=-30; ②p1=-2,p2=3 1. 针对极点参数①②, 画出系统零、极点分布图, 判断该系统稳定性. 2. 针对极点参数①②,绘出系统的脉冲响应曲线,并观察t →∞时, 脉冲响应变化趋势. 3. 针对极点参数①, 绘出系统的频响曲线. 四、实验要求 1．预习实验原理； 2．对实验内容编写程序(M 文件),上机运行； 3．绘出实验内容的各相应曲线或图。 五、实验设备 1.装MATLAB 软件的计算机 1台

实验一 信号的时域与频域分析

实验一信号的时域和频域分析 一、实验目的 1、了解SystemView图符库的分类； 2、掌握SystemView各个功能库常用图符的功能及其使用方法； 3、掌握信号的时域与频域的分析方法； 4、掌握SystemView分析窗口的使用； 5、能利用分析窗口对波形进行时域与频域的分析。 二、实验内容 1、按照实例使用图符构建简单的通信系统，并了解每个图符的功能； 2、建立简单的调制系统，并使用分析窗口对输出信号进行时域与频域的分析， 得出分析结果。 三、SystemView常用图符库 SystemView的图符库功能十分丰富，一共分为以下几个大类： 1.基本库 SystemView的基本库包括信源库、算子库、函数库、信号接收器库等，它为该系统仿真提供了最基本的工具。 （信源库）：SystemView为我们提供了16种信号源，可以用它来产生任意信号 （算子库）功能强大的算子库多达31种算子，可以满足您所有运算的要求 （函数库）32种函数尽显函数库的强大库容！ （信号接收器库）12种信号接收方式任你挑选，要做任何分析都难不倒它 2.扩展功能库 扩展功能库提供可选择的能够增加核心库功能的用于特殊应用的库。它允许通信、DSP、射频/模拟和逻辑应用。 （通信库）：包含有大量的通信系统模块的通信库，是快速设计和仿真现代通信系统的有力工具。这些模块从纠错编码、调制解调、到各种信道模型一应俱全。 （DSP库）：DSP库能够在你将要运行DSP芯片上仿真DSP系统。该库支持大多DSP芯片的算法模式。例如乘法器、加法器、除法器和反相器的图标代表真正的DSP 算法操作符。还包括高级处理工具：混合的Radix FFT、FIR和IIR滤波器以及块传输等。 （逻辑运算库）：逻辑运算自然离不开逻辑库了，它包括象与非门这样的通用器件的图标、74系列器件功能图标及用户自己的图标等。

理工大学信号与系统实验报告连续时间系统的复频域分析

理工大学信号与系统实验报告连续时间系统的 复频域分析 Document serial number【UU89WT-UU98YT-UU8CB-UUUT-UUT108】

实验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 为有理分式时，它可以表示为两个多项式之比： 110 1 10 ...(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 N r r r X 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 st H h dt +∞ --∞ = ? （5） 连续时间系统的系统函数还可以由系统输入与输出信号的拉氏变换之比得到。 (s)(s)/X(s)H Y = （6） 单位冲激响应(t)h 反映了系统的固有性质，而(s)H 从复频域反映了系统的固有性质。由（6）描述的连续时间系统，其系统函数为s 的有理函数 110 1 10 ...(s)...M M M M N N N N b s b s b H a s a s a ----+++=+++ （7） 3.连续时间系统的零极点分析 系统的零点指使式（7）的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统的值为无穷大。通常将系统函数的零极点绘在s 平面上，零点用O 表示，极点用?表示，这样得到的图形为零极点分布图。可以通过利用MATLAB 中的求多项式根的roots 函数来实现对（7）分子分母根的求解，调用格式如下：

常用信号的频谱分析及时域采样定理

常用信号的频谱分析及时域采样定理

开课学期 2016-2017 学年第 2 学期 实验课程信号与系统仿真实验 实验项目常用信号的频谱分析及时域采样定理 班级学号学生姓名 实验时间实验台号A11 操作成绩报告成绩 一、实验目的 1.掌握常用信号的频域分析方法； 2.掌握时域采样定理； 3.掌握时域采样信号恢复为原来连续信号的方法及过程。 二、实验性质 验证性 三、预习内容 1.时域采样定理的内容及信号时域采样过程； 2.连续信号经时域采样后，信号的频谱发生的变化； 3.时域采样信号恢复为原来连续信号的方法及过程。 四、实验内容(编写程序，绘制实验结果) 1.实现周期信号的频谱 f(t)=sin( 2*80t) 程序： fa='sin(2.*pi.*80.*t)';%原信号 fs0=10000; %采样频率 tp=0.1;%时间范围 t=[-tp:1/fs0:tp];%信号持续时间范围 k1=0:999;k2=-999:-1; m1=length(k1);m2=length(k2); f=[fs0*k2/m2,fs0*k1/m1];%信号频率范围 w=[-2*pi*k2/m2,2*pi*k1/m1]; fx1=eval(fa);%把文本fa赋值给信号fx1 FX1=fx1*exp(-j*[1:length(fx1)]'*w);%进行傅立叶变换 figure subplot(2,1,1),plot(t,fx1,'r'); title('原信号');xlabel('时间t(s)');%原信号的时域波形图 axis([min(t),max(t),min(fx1),max(fx1)]); subplot(212),plot(f,abs(FX1),'r'); title('原信号频谱');xlabel ('频率f(Hz)');%频域波形图 axis([-100,100,0,max(abs(FX1))+5]);

信号与系统 信号的时域分析 及Matlab实现

信号与系统 课程实验报告 实验四 实验名称 信号的时域分析 及Matlab 实现 系 别 教师姓名 实验地点 5309 实验日期 2011-06-20 学生姓名 学号 一、实验内容 1．预习实验原理。 2．对实验内容编写程序(M 文件)，上机运行。 3．记录并整理实验数据。 二、实验目的 1．掌握用Matlab 分析系统时间响应的方法。 2．掌握用Matlab 分析连续、离散系统的冲激响应的方法。 3．理解系统零、极点分布与系统稳定性关系。 三、涉及实验的相关情况介绍（包含使用软件或实验设备等情况） 计算机一台(安装MATLAB6.5版本或以上版本) 四、实验试做记录（含程序、数据记录及分析） 1．设) 2)(1()(p s p s s s H --= 设①p1=-2，p2=-30； ②p1=-2，p2=3 （1）针对极点参数①②，画出系统零、极点分布图，判断该系统稳定性。 程序： clear num=[1]; den=[1,32,60]; [r,poles,k]=residue(num,den); p=roots(den);

z=roots(num); subplot(2,2,1);plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on T=0:0.1:10; y1=impulse(num,den,T); subplot(2,2,2);plot(T,y1);grid;title('脉冲响应') 结果： 位于s左半平面，稳定 （2）针对极点参数①②，绘出系统的脉冲响应曲线，并观察t→∞时，脉冲响应变化趋势。 程序： clear num=[1]; den=[1,-1,-6]; [r,poles,k]=residue(num,den); p=roots(den); z=roots(num); subplot(2,2,1);plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on T=0:0.1:10; y1=impulse(num,den,T); subplot(2,2,2);plot(T,y1);grid;title('脉冲响应') 结果：