实验六-信号与系统复频域分析
- 格式:doc
- 大小:75.00 KB
- 文档页数:10
信号与系统—信号的频域分析频域分析是指将信号从时间域转换为频域的过程,并通过对信号在频域上的性质和特征进行分析与研究。
频域分析对于理解信号的频率特性、频谱分布等方面的特性有很大的帮助,是信号处理领域中不可或缺的分析工具。
频域分析的基本方法之一是傅里叶变换。
傅里叶变换可以将连续时间域中的信号转换为离散频域中的信号,也可以将离散时间域中的信号转换为连续频域中的信号。
它通过将信号分解为不同频率的正弦波的组合来分析信号的频谱分布。
傅里叶变换的基本公式为:两个公式其中,X(f)表示信号在频域中的频谱,x(t)表示信号在时间域中的波形,f表示频率。
傅里叶变换得到的频谱图可以展示信号在不同频率上的能量分布情况,从而能够更直观地了解信号的频率成分。
频谱图通常以频率为横轴,信号在该频率上的幅度或相位为纵轴,用于描述信号在频域中的变化情况。
除了傅里叶变换,还有其他一些常用的频域分析方法,如离散傅里叶变换(DFT)、快速傅里叶变换(FFT)等。
离散傅里叶变换是对离散时间域中的信号进行频域分析的方法,快速傅里叶变换是一种高效的计算离散傅里叶变换的方法。
频域分析主要包括信号的频谱分析和系统的频率响应分析两个方面。
在信号的频谱分析中,我们可以通过观察信号在频域上的能量分布情况来判断信号的频率成分、频率范围等信息。
而在系统的频率响应分析中,我们可以通过研究系统在不同频率上的响应特性来了解系统对不同频率信号的传输、增益、衰减等情况。
频域分析在实际应用中有着广泛的应用。
例如,在音频处理领域中,频域分析可以用于声音信号的频谱分析和音效处理等方面。
在通信系统中,频域分析可以用于信号的调制解调、信道估计、信号检测等。
在图像处理中,频域分析可以用于图像的锐化、降噪、压缩等方面。
总结起来,信号的频域分析是信号与系统课程中的重要内容,它通过将信号从时间域转换为频域来研究信号的频率特性和频谱分布等问题。
傅里叶变换是频域分析中常用的方法之一,它可以将信号分解为不同频率的正弦波的组合。
第12章连续信号与系统的复频域分析 (1)f(t)=(1-e^-0.5t)u(t) (6)f(t)=sin(πt)[u(t)-u(t-1)] (1)代码如下:实验结果如下: syms t f1=(1-exp(-0.5*t))*heaviside(t); F1=laplace(f1); F1 (6)代码如下:实验结果如下: syms t f6=sin(pi*t)*(heaviside(t)-heaviside(t-1)); F6=laplace(f1); F6 12.2利用MATLAB的ilaplace函数,求下列象函数F(s)的拉普拉斯逆变换。
(1))3)(2(1)(sssssF(2)862)(2ssssF (1)实验代码如下:实验结果如下: syms s F1=(s+2)/(s*(s+2)*(s+3)); F=ilaplace(F1) (2)实验代码如下:实验结果如下: syms s F2=(s+2)/(s^2+6*s+8); F=ilaplace(F2)
12.3利用MATLAB的residue函数求12.2题中(1),(2)各小题的拉普拉斯逆变换,并与ilaplace函数的计算结果进行比较。 (1)实验代码如下:实验结果如下: a=[1 1]; b=[1 5 6 0]; [k,p,c]=residue(a,b)
(2)实验代码如下:实验结果如下: a=[1 2]; b=[1 6 8]; [k,p,c]=residue(a,b)
12.8已知连续时间系统的系统函数H(s)分别如下: (1))4)(2(1)(sssF(4)44)(22sssF (1)实验代码如下:实验结果如下: a=[1]; b=[1 6 8]; H=tf(a,b); subplot(2,1,1) pzmap(H); p=pole(H) z=zero(H) subplot(2,1,2) impulse(a,b)
信号与系统时域频域复频域的异同点下载提示:该文档是本店铺精心编制而成的,希望大家下载后,能够帮助大家解决实际问题。
文档下载后可定制修改,请根据实际需要进行调整和使用,谢谢!本店铺为大家提供各种类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by this editor. I hope that after you download it, it can help you solve practical problems. The document can be customized and modified after downloading, please adjust and use it according to actual needs, thank you! In addition, this shop provides you with various types of practical materials, such as educational essays, diary appreciation, sentence excerpts, ancient poems, classic articles, topic composition, work summary, word parsing, copy excerpts, other materials and so on, want to know different data formats and writing methods, please pay attention!信号与系统:时域、频域与复频域的异同点在信号与系统理论中,时域、频域和复频域是描述信号特性的重要概念。
实验六 无失真传输系统一、实验目的1、了解无失真传输的概念。
2、了解无失真传输的条件。
二、实验内容1、观察信号在失真系统中的波形。
2、观察信号在无失真系统中的波形。
三、实验仪器1、信号与系统实验箱一台(主板);2、系统频域与复域的分析模块一块。
3、20M 双踪示波器一台。
四、实验原理1、一般情况下,系统的响应波形和激励波形不相同,信号在传输过程中将产生失真。
线性系统引起的信号失真有两方面因素造成,一是系统对信号中各频率分量幅度产生不同程度的衰减,使响应各频率分量的相对幅度产生变化,引起幅度失真。
另一是系统对各频率分量产生的相移不与频率成正比,使响应的各频率分量在时间轴上的相对位置产生变化,引起相位失真。
线性系统的幅度失真与相位失真都不产生新的频率分量。
而对于非线性系统则由于其非线性特性对于所传输信号产生非线性失真,非线性失真可能产生新的频率分量。
所谓无失真是指响应信号与激励信号相比,只是大小与出现的时间不同,而无波形上的变化。
设激励信号为)(t e ,响应信号为)(t r ,无失真传输的条件是)()(0t t Ke t r -= (4-1) 式中K 是一常数,0t 为滞后时间。
满足此条件时,)(t r 波形是)(t e 波形经0t 时间的滞后,虽然,幅度方面有系数K 倍的变化,但波形形状不变。
2、对实现无失真传输,对系统函数)(ωj H 应提出怎样的要求?设)(t r 与)(t e 的傅立叶变换式分别为)()(ωωj E j R 与。
借助傅立叶变换的延时定理,从式4-1可以写出)()(t j ej KE j R ωωω-= (4-2)此外还有 )()()(ωωωj E j H j R = (4-3) 所以,为满足无失真传输应有)(t j Kej H ωω-= (4-4)(4-4)就是对于系统的频率响应特性提出的无失真传输条件。
欲使信号在通过线性系统时不产生任何失真,必须在信号的全部频带内,要求系统频率响应的幅度特性是一常数,相位特性是一通过原点的直线。
利用matlab进行信号和系统的复频域分析心得
利用MATLAB进行信号和系统的复频域分析是非常常见的。
下面是一些心得和步骤供参考:
1. 导入信号数据:首先,你需要将信号数据导入到MATLAB中。
这可以通过多种方式实现,例如读取文件或直接生成信号矩阵。
2. 选择合适的频域分析方法:根据你的需求和信号类型,选择合适的频域分析方法。
常见的方法包括傅里叶变换、快速傅里叶变换(FFT)、离散傅里叶变换(DFT)等。
3. 执行频域分析:使用MATLAB提供的相应函数,如fft()或fftshift()来执行频域分析。
这将为你提供信号的频谱信息。
4. 绘制频谱图:使用MATLAB的绘图函数,如plot()或stem(),将频谱数据可视化为频谱图。
你可以选择线性频谱图或对数频谱图,具体取决于信号特性和需求。
5. 分析频谱信息:根据频谱图,你可以分析信号的频率分量、幅度特性以及相位特性。
对于系统分析,你还可以计算系统的传递函数。
6. 系统设计和优化:根据频域分析结果,你可以对系统进行设计和优化。
例如,你可以确定降噪滤波器的截止频率,或者针对特定的频率范围进行信号增强。
需要注意的是,MATLAB提供了丰富的信号处理和系统分析工具箱,可以帮助你更轻松地完成复频域分析任务。
同时,请确保使用合法授权的软件和工具,遵守中国法律政策。
信号与系统实验报告——连续时间系统的复频域分析班级:05911101学号:**********姓名:***实验五连续时间系统的复频域分析——1120111487 信息工程(实验班)蒋志科一、实验目的①掌握拉普拉斯变换及其反变换的定义,并掌握MA TLAB 实现方法 ②学习和掌握连续时间系统系统函数的定义及其复频域分析方法③掌握系统零极点的定义,加深理解系统零极点分布与系统特性的关系。
二、实验原理与方法 1、拉普拉斯变换连续时间信号x(t)的拉普拉斯变换定义为:X s =x (t )e −st dt +∞−∞拉普拉斯反变换为:x t =12πj X (s )e st ds σ+j ∞σ−j ∞在MA TLAB 中可以采用符号数学工具箱中的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 。
2、连续时间系统的系统函数连续时间系统的系统函数是系统单位冲激响应的拉氏变换H s =ℎ(t )e −st dt +∞−∞此外,连续时间系统的系统函数还可以由系统输入和输出信号的拉氏变换之比得到H s =Y(s)/X(s) 单位冲激响应h(t)反映了系统的固有性质,而H(s)从复频域反映了系统的固有性质。
对于H(s)描述的连续时间系统,其系统函数s 的有理函数H s =b M s M +b M−1s M−1+⋯+b 0a n s n +a n −1s M−1+⋯+a 03、连续时间系统的零极点分析系统的零点指使式H s 的分子多项式为零的点,极点指使分母多项式为零的点,零点使系统的值为零,极点使系统函数的值无穷大。
实验五 信号与系统复频域分析一、 实验目的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.学会用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()13F 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 已知系统函数为 321221s s s +++H(s)=试画出其零极点分布图,求系统的单位冲激响应h(t)和频率响应H ω(j ),并判断系统是否稳定。
解:其MATLAB 程序如下: num=[1];den=[1,2,2,1]; sys=tf(num,den); figure(1);pzmap(sys); t=0:0.02:10;h=impulse(num,den,t); figure(2);plot(t,h) title('Impulse Response') [H,w]=freqs(num,den); figure(3);plot(w,abs(H)) xlabel('\omega')title('Magnitude Response')3.用MATLAB 进行Laplace 正、反变换MATLAB 的符号数学工具箱提供了计算Laplace 正、反变换的函数Laplace 和ilaplace,其调用格式为上述两式右端的f 和F 分别为时域表示式和s 域表示式的符号表示,可以应用函数sym 实现,其调用格式为S=sym(A)式中,A 为待分析表示式的字符串,S 为符号数字或变量。
例6-3 试分别用Laplace 和ilaplace 函数求 (1)()sin()()t f t e at u t -=的Laplace 变换;(2)22()1s F s s =+的Laplace 反变换。
解:(1)其程序为f=sym('exp(-t)*sin(a*t)');F=laplace(f)或syms a tF=laplace(exp(-t)*sin(a*t))(2)其程序为F=sym('s^2/(s^2+1)');ft=ilaplace(F)或syms sft= ilaplace(s^2/(s^2+1))4.离散系统零极点图离散系统可以用下述差分方程描述:∑∑= =-= -Mm mNiimkfbikya)()(Z变换后可得系统函数:NNMMzazaazbzbbzFzYzH----++++++==......)()()(1111用MATLAB提供的root函数可分别求零点和极点,调用格式是p=[a0,a1…an],q=[b0,b1…bm,0,0…0], 补0使二者维数一样。
画零极点图的方法有多种,可以用MATLAB函数[z,p,k]=tf2zp(b,a)和zplane(q,p),也可用plot命令自编一函数ljdt.m,画图时调用。
function ljdt(A,B)% The function to draw the pole-zero diagram for discrete system p=roots(A);%求系统极点q=roots(B);%求系统零点p=p';%将极点列向量转置为行向量q=q';%将零点列向量转置为行向量x=max(abs([p q 1]));%确定纵坐标范围x=x+0.1;y=x;%确定横坐标范围clfhold onaxis([-x x -y y]) %确定坐标轴显示范围w=0:pi/300:2*pi;t=exp(i*w);plot(t) %画单位园axis('square')plot([-x x],[0 0]) %画横坐标轴plot([0 0],[-y y]) %画纵坐标轴text(0.1,x,'jIm[z]')text(y,1/10,'Re[z]') plot(real(p),imag(p),'x') %画极点 plot(real(q),imag(q),'o')%画零点 title('pole-zero diagram for discrete system') %标注标题hold off例6-4 求系统函数零极点图131)(45+-+=z z z z H a=[3 -1 0 0 0 1]; b=[1 1]; ljdt(a,b) p=roots(a) q=roots(b) pa=abs(p)5.离散系统的频率特性离散系统的频率特性可由系统函数求出,既令ωj e z =,MATLAB函数freqz 可计算频率特性,调用格式是:[H ,W]=freqz(b,a,n),b 和a 是系统函数分子分母系数,n 是π-0范围 n 个等份点,默认值512,H 是频率响应函数值,W 是相应频率点;[H ,W]=freqz(b,a,n,’whole’), n 是π2-0范围 n 个等份点; freqz(b,a,n),直接画频率响应幅频和相频曲线; 例6-5 系统函数zz z H 5.0)(-=运行如下语句,可得10个频率点的计算结果A=[1 0];B=[1 -0.5];[H,W]=freqz(B,A,10)继续运行如下语句,可将400个频率点的计算结果用plot语句画幅频和相频曲线B=[1 -0.5];A =[1 0];[H,w]=freqz(B,A,400,'whole');Hf=abs(H);Hx=angle(H);clffigure(1)plot(w,Hf)title('离散系统幅频特性曲线')figure(2)plot(w,Hx)title('离散系统相频特性曲线')还可用freqz语句直接画图,注意区别A=[1 0];B=[1 -0.5];freqz(B,A,400)例6-6 用几何矢量法,自编程序画频率响应原理:频率响应∏∏==--=Ni ijMj jjj p e qe e H 11)()()(ωωω编程流程:定义Z 平面单位圆上k 个频率等分点;求出系统函数所有零点和极点到这些等分点的距离;求出系统函数所有零点和极点到这些等分点的矢量的相角;求出单位圆上各 频率等分点的)()(ωϕω和j e H画指定范围内的幅频与相频。
若要画零极点图,可调用ljdt.m 函数。
function dplxy(k,r,A,B)%The function to draw the frequency response of discrete system p=roots(A); %求极点 q=roots(B);%求零点figure(1) ljdt(A,B)%画零极点图 w=0:l*pi/k:r*pi;y=exp(i*w);%定义单位圆上的k 个频率等分点 N=length(p); %求极点个数 M=length(q);%求零点个数 yp=ones(N,1)*y;%定义行数为极点个数的单位圆向量yq=ones(M,1)*y;%定义行数为零点个数的单位圆向量vp=yp-p*ones(1,r*k+1); %定义极点到单位圆上各点的向量vq=yq-q*ones(1,r*k+1); %定义零点到单位圆上各点的向量Ai=abs(vp); %求出极点到单位圆上各点的向量的模Bj=abs(vq); %求出零点到单位圆上各点的向量的模Ci=angle(vp); %求出极点到单位圆上各点的向量的相角Dj=angle(vq); %求出零点到单位圆上各点的向量的相角fai=sum(Dj,1)-sum(Ci,1); %求系统相频响应H=prod(Bj,1)./prod(Ai,1); %求系统幅频响应figure(2)plot(w,H); %绘制幅频特性曲线title('离散系统幅频特性曲线')xlabel('角频率')ylabel('幅度')figure(3)plot(w,fai)title('离散系统的相频特性曲线') xlabel('角频率') ylabel('相位')已知系统函数114/11)1(4/5)(----=zz z H ,画频率响应和零极点图。
A=[1 -1/4]; B=[5/4 -5/4];dplxy(500,2,A,B) %绘制系统2π频率范围内500个频率点的幅频和相频特性曲线及零极点图三、上机实验内容1.验证实验原理中所述的相关程序; 2.求信号)()(3t u te t f t -=的拉普拉斯变换3.求函数23795)(223+++++=s s s s s s F 的反变换 4.已知连续系统的系统函数如下,试用MATLAB 绘制系统的零极点图,并根据零极点图判断系统的稳定性 5.系统函数是321551---+++zzz求频率响应。