实验4-LTI系统的频域分析
针对LTI系统频率响应，加深了对于基本概念的掌握和理解，学习并掌握了关于LTI系统频率特性的分析方法。
二，实验原理1.连续时间系统的频率响应调用函数freqs:[h,w]:freqs(b,a)计算默认频率范围内200个频率点上的频率响应的取样值，这200个频率点记录在w中。
h=freqs(b,a，w)b，a分别为表示H(jw）的有理多项式中分子和分母多项式的系数向量，w为频率取样点，返回值h就是频率响应在频率取样点上的数值向量。
[h,w]:freqs(b,a,n)计算默认频率范围内n个频率点上的频率响应的取样值，这n个频率点记录在w中。
freqs(b,a,…)这种调用格式不返回频率响应的取样值，而是以对数坐标的方式绘出系统的频率响应和相频响应。
2.离散时间系统的频率响应调用函数freqz:[H,w]:freqz(b,a，‘whole’)计算0~2πn个频率点上的频率响应的取样值，这n个频率点记录在w中。
H=freqz(b,a，n)b，a分别为有理多项式中分子和分母多项式的系数向量，返回值H就是频率响应在0到pi范围内n个频率等分点上的数值向量，w包含了这n个频率响应。
[H,w]:freqz(b,a,w) w为频率取样点，计算这些频率点上的频率响应的取样值。
freqz(b,a,…)这种调用格式不返回频率响应的取样值，而是直接绘出系统的频率响应和相频响应。
（1）已知一个RLC电路构造的二阶高通滤波器如图，其中①计算该电路系统的频率响应及高通截止频率。
答：②利用MATLAB绘制幅度响应和相位响应曲线，比较系统的频率特性与理论计算的结果是否一致。
MATLAB程序如下：b=[0.04 0 0]a=[0.04 0.4 2][H,w]=freqs(b,a)subplot(211)plot(w,abs(H))set(gca,'xtick')set(gca,'ytick',[0 0.4 0.707 1])xlabel('\omega(rad/s)')ylabel('Magnitude')title('|H(j\omega)|')grid onsubplot(212)plot(w,angle(H))set(gca,'xtick')xlabel('\omega(rad/s)')ylabel('Phase')title('\phi(\omega)')grid on（2）已知一个RC系统电路如图。
①对不同的RC值，用MATLAB画出系统的幅度响应曲线|H（w）|，观察实验结果，分析如图所示电路具有什么样的频率特性？系统的频率特性随着RC值的改变，有何变化规律？MATLAB程序如下：A=input('A=')b=[1]a=[A 1][H,w]=freqs(b,a)plot(w,abs(H))set(gca,'ytick',[0 0.4 0.707 1])xlabel('\omega(rad/s)')ylabel('Magnitude')title('|H(j\omega)|')grid on程序执行如下：A=0.01A=1A=100由程序执行结果可以看出，RC电路具有带通特性，随着RC值的减小，带通频率增加。
②系统输入信号x(t)=cos(100t)+cos(3000t)，t=0~0.2s，该信号包含了一个低频分量和一个高频分立，试确定适当的RC值，滤出信号中的高频分量，并绘出滤波前后的时域信号波形及系统的频率响应曲线。
MATLAB程序如下：A=input('A=');b=[1]a=[A 1][H,w]=freqs(b,a)plot(w,abs(H))set(gca,'ytick',[0 0.4 0.707 1])xlabel('\omega(rad/s)')ylabel('Magnitude')title('|H(j\omega)|')grid onMATLAB程序执行如下：A=0.0004t=0:0.0001:0.2x=cos(100*t)+cos(3000*t)plot(t,x)xlabel('t')title('x(t)')可得：输入：t=0:0.0001:0.2x=cos(100*t)plot(t,x)xlabel('t')title('x(t)')输出（3）已知离散系统的系统框图如图。
①写出M=8时系统的差分方程和系统函数。
x[n]+x[n-1]+x[n-2]+...+x[n-8]=y[n]H(z)=1+z^-1+z^-2+...+z^-8h[n]=1 1 1 1 1 1 1②利用MATLAB计算系统的单位抽样响应。
MATLAB程序如下：b=[1 1 1 1 1 1 1 1 1]a=[1]impz(b,a,0:20)程序执行如下：③试利用MATLAB绘出其系统零极点分布图、幅频和相频特性曲线，并分析该系统具有怎样的频率特性。
MATLAB程序如下：b=[1 1 1 1 1 1 1 1 1]a=[1][H,w]=freqz(b,a)subplot(211)plot(w/pi,abs(H))xlabel('\omega(\pi)')ylabel('Magnititude')title('|H(e^j^\omega)|')grid onsubplot(212)plot(w/pi,angle(H)/pi)xlabel('\omega(\pi)')ylabel('Phase(\pi)')title('\pheta(\omega)')grid on程序执行如下：幅频、相频：A=[1 1 1 1 1 1 1 1 1]r=roots(A)r =0.7660 + 0.6428i0.7660 - 0.6428i0.1736 + 0.9848i0.1736 - 0.9848i-0.9397 + 0.3420i-0.9397 - 0.3420i-0.5000 + 0.8660i-0.5000 - 0.8660iz=[0.7660+0.6428i 0.7660-0.6428i 0.1736+0.9848i 0.1736-0.9848i -0.9397+0.3420i -0.9397-0.3420i -0.5000+0.8660i -0.5000-0.8660i]'p=[0]'零极点分布图由程序执行可知，该系统具有高通的频率特性。
（4）已知一离散时间LTI系统的频率响应H(e^jΩ)如图，输入信号为x(n)=cos(0.3πn)+0.5cos(0.8πn)。
试分析正弦信号sin(Ω0t)通过频率响应为H(e^jΩ)的离散时间系统的响应，并根据分析结果计算系统对于x（n）的响应y（n），用MATLAB绘出系统输入与输出波形。
MATLAB程序如下：n=-10:10x=cos(0.3*pi*n)+0.5*cos(0.8*pi*n)stem(n,x)xlabel('n')title('x(n)')输入输出y=2*cos(0.3*pi*n)stem(n,y，‘filled’)xlabel('n')title('y(n)')观察实验结果，分析如图所示系统具有怎样的频率特性？从输入输出信号上怎么反映出系统的频率特性？答：可知，该系统具有带通的频率特性。
1.已知二阶高通滤波器（1）计算该电路系统的频率响应及高通截止Y(ω)=(X(ω))/(1/(1/R+1/jωL)+1/jωC)∙1/(1/R+1/jωL)H(ω)=Y(ω)/X(ω) =【(jω)】^2/(【(jω)】^2+10(jω)+50)令H(ω)=0.707得ω=7.07(2)利用matlab绘制幅度响应和相位响应曲线2．已知RC电路如图所示H(ω)=(1⁄RC)/((jω+1)⁄RC) （1）系统的频率响应随RC值的变化RC1=0.001;RC2=0.002;RC3=0.01;a1=1/RC1;a2=1/RC2;a3=1/RC3;A1=[1a1];A2=[1 a2];A3=[1 a3];B1=a1;B2=a2;B3=a3;w1=linspace(0,3000,200);w2=linspace(0,3000,400);w3=linspace(0,3000,500);H1=freqs(B1,A1,w1);H2=freqs(B2,A2,w2);H3=freqs(B3,A3,w3);subplot(311);plot(w1,abs(H1));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega( rad/s)');ylabel('频率响应');title('|H(j\omega)|,RC=0.001');subplot(312);plot(w2,abs(H2));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega( rad/s)');ylabel('频率响应');title('|H(j\omega)|,RC=0.002');subplot(313);plot(w3,abs(H3));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega( rad/s)');ylabel('频率响应');title('|H(j\omega)|,RC=0.0010');（2）滤波前后的时域信号波形和频率响应曲线t=0:0.0001:0.2;RC=0.01;a=1/RC;B=[a];A=[1 a];[H,w]=freqs(B,A);f=cos(100*t)+cos(3000*t);Hw1=a/(1i*100+a);Hw2=a/(1i*3000+a);y=abs(Hw1)*cos(100*t)+abs(Hw2)*cos(3000*t); subplot(411);plot(t,f);grid on;xlabel('time(sec)');ylabel('f(t)');title('f(t)=cos(100*t)+cos(3000*t) 滤波之前'); subplot(412);plot(t,y);grid on;xlabel('time(sec)');ylabel('y(t)');title('y(t) 滤波之后');subplot(413);plot(w,abs(H));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega(rad/s)');ylabel('Magnitude');title('|H(j\omega)|');subplot(414);plot(w,angle(H));grid on;xlabel('\omega(rad/s)');ylabel('Phase');title('\phi(\omega)');四，心得体会针在这次实验当中，我感受到了matlab在LTI系统分析中的强大作用，系统函数让分析过程变得异常简单。
MATLAB与信号实验 —— 连续LTI系统的频域分析
上机实验3 连续LTI 系统的频域分析一．实验目的(1).掌握连续时间信号傅立叶变换和傅立叶逆变换的实现方法，以及傅立叶变换的时移特性，傅立叶变换的频移特性的实现方法;(2).了解傅立叶变换的频移特性及其应用;(3).掌握函数fourier 和函数ifourier 的调用格式及作用;(4).掌握傅立叶变换的数值计算方法，以及绘制信号频谱图的方法。
二．实验原理1.系统的频率特性连续LTI 系统的频率特性又称频率响应特性，是指系统在正弦信号激励下的稳态响应随频率变化的情况，又称系统函数()ωH 。
对于一个零状态的线性系统，如图2.3-1所示。
其系统函数()ωH 定义为)()()(ωωωj X j Y j H =式中，()ωX 为系统激励信号的傅里叶变换，()ωY 为系统在零状态条件下输出响应的傅里叶变换。
系统函数()ωH 反映了系统内在的固有特性，它取决于系统自身的结构及组成系统元器件的参数，与外部激励无关，是描述系统特性的一个重要参数。
()ωH 是ω的复函数，可以表示为()ωH =()ωH ()e j ωϕ。
其中，()ωH 随ω变化的规律称为系统的幅频特性：()ωϕ随ω变化的规律称为系统的相频特性。
频率特性不仅可用函数表达式表示，还可以随频率f 变化的曲线来表示。
当频率特性曲线采用对数坐标时，又称为波特图。
)(ωj X )(ωj2.连续时间信号傅里叶变换的数值计算方法算法理论依据：()ττωτωω∑⎰∞∞---∞+∞-==e e n j t j n f dt t f j F )(lim )(（2.2-1）当)(t f 为时限信号时，或和近似的看做时限信号时，式（2-1）中的n 取值可认作是有限的，设为N ，则可得()()e n j N n kn f k F τωττ--=∑=10，0<=k<=N (2.2-2) 式（2.2-2）中k N k τπω2=。
编程中需要注意的是：要正确生成信号)(t f 的N 个样本)(τn f 的向量及向量e n j kτω-。
连续时间LTI系统的复频域分析
x(t )
6.3
k 0
dtk
k
0
dtk
对式6.3
两边做拉普拉斯变换，则有
N
akskY( s)
M
bkskX ( s)
k
0
k
0
M
bksk
Y( s)
k
0
6.4
即
H ( s)
N
X (s)
k
0
aksk
式6.4
告诉我们， 对于一个能够用线性常系数微分方程描述的连续时间
LTI系统，它的
系统函数是一个关于复变量s的有理多项式的分式， 其分子和分母的多项式系数与系统微分
%确定绘图区域
c = abs (c);
subplot (211)
%计算拉普拉斯变换
mesh (a,b,c);
surf (a,b,c);
%绘制曲面图
view (-60,20)
axis ([-0,5,-20,20,0,2]);
%调整观察视角
title ('The Laplace transform of the rectangular pulse');
则取决于系统极点位置的分布。
需要特别强调的是，MATLAB总是把由分子和分母多项式表示任何系统都当作是因果系
统。所以，利用impulse ()函数求得的单位冲激响应总
是因果信号。
5、系统函数的零极点分布与系统的滤波特性
系统具有何种滤波特性，主要取决于系统的零极点所
处的位置。没有零点的系统，通常是一个低通滤波器。
线，带返回值则将冲激响应值存于向量h之中。
2、系统函数的零极点分布图
系统函数的 零极点图（Zero-pole diagram）能够直观地表示系统的零点和极点在s平
实验4-LTI系统的频域分析
一,试验目标针对LTI体系频率响应,加深了对于根本概念的控制和懂得,进修并控制了关于LTI体系频率特征的剖析办法.二,试验道理挪用函数freqs:[h,w]:freqs(b,a)盘算默认频率规模内200个频率点上的频率响应的取样值,这200个频率点记载在w中.h=freqs(b,a,w)b,a分离为暗示H(jw）的有理多项式平分子和分母多项式的系数向量,w为频率取样点,返回值h就是频率响应在频率取样点上的数值向量.[h,w]:freqs(b,a,n)盘算默认频率规模内n个频率点上的频率响应的取样值,这n个频率点记载在w中.freqs(b,a,…)这种挪用格局不返回频率响应的取样值,而是以对数坐标的方法绘出体系的频率响应和相频响应.挪用函数freqz:[H,w]:freqz(b,a,‘whole’)盘算0~2πn个频率点上的频率响应的取样值,这n个频率点记载在w中.H=freqz(b,a,n)b,a分离为有理多项式平分子和分母多项式的系数向量,返回值H就是频率响应在0到pi规模内n个频率等分点上的数值向量,w包含了这n个频率响应.[H,w]:freqz(b,a,w)w为频率取样点,盘算这些频率点上的频率响应的取样值.freqz(b,a,…)这种挪用格局不返回频率响应的取样值,而是直接绘出体系的频率响应和相频响应.三,试验内容（1）已知一个RLC电路结构的二阶高通滤波器如图,个中①盘算该电路体系的频率响应及高通截止频率.答：②应用MATLAB绘制幅度响应和相位响应曲线,比较体系的频率特征与理论盘算的成果是否一致.MATLAB程序如下：b=[0.04 0 0]a=[0.04 0.4 2][H,w]=freqs(b,a)subplot(211)plot(w,abs(H))set(gca,'xtick')set(gca,'ytick',[0 0.4 0.707 1])xlabel('\omega(rad/s)')ylabel('Magnitude')title('|H(j\omega)|')grid onsubplot(212)plot(w,angle(H))set(gca,'xtick')xlabel('\omega(rad/s)')ylabel('Phase')title('\phi(\omega)')grid on（2）已知一个RC体系电路如图.①对不合的RC值,用MATLAB画出体系的幅度响应曲线|H（w）|,不雅察试验成果,剖析如图所示电路具有什么样的频率特征？体系的频率特征跟着RC值的转变,有何变更纪律？MATLAB程序如下：A=input('A=')b=[1]a=[A 1][H,w]=freqs(b,a)plot(w,abs(H))set(gca,'ytick',[0 0.4 0.707 1])xlabel('\omega(rad/s)')ylabel('Magnitude')title('|H(j\omega)|')grid on程序履行如下：A=1A=100由程序履行成果可以看出,RC电路具有带通特征,跟着RC值的减小,带通频率增长.②体系输入旌旗灯号x(t)=cos(100t)+cos(3000t),t=0~0.2s,该旌旗灯号包含了一个低频分量和一个高频分立,试肯定恰当的RC值,滤出旌旗灯号中的高频分量,并绘出滤波前后的时域旌旗灯号波形及体系的频率响应曲线.MATLAB程序如下：A=input('A=');b=[1]a=[A 1][H,w]=freqs(b,a)plot(w,abs(H))set(gca,'ytick',[0 0.4 0.707 1])xlabel('\omega(rad/s)')ylabel('Magnitude')title('|H(j\omega)|')grid onMATLAB程序履行如下：x=cos(100*t)+cos(3000*t)plot(t,x)xlabel('t')title('x(t)')可得：输入：x=cos(100*t)plot(t,x)xlabel('t')title('x(t)')输出（3）已知离散体系的体系框图如图.①写出M=8时体系的差分方程和体系函数.x[n]+x[n-1]+x[n-2]+...+x[n-8]=y[n]H(z)=1+z^-1+z^-2+...+z^-8h[n]=1 1 1 1 1 1 1②应用MATLAB盘算体系的单位抽样响应.MATLAB程序如下：b=[1 1 1 1 1 1 1 1 1]a=[1]impz(b,a,0:20)程序履行如下：③试应用MATLAB绘出其体系零顶点散布图.幅频和相频特征曲线,并剖析该体系具有如何的频率特征.MATLAB程序如下：b=[1 1 1 1 1 1 1 1 1]a=[1][H,w]=freqz(b,a)subplot(211)plot(w/pi,abs(H))xlabel('\omega(\pi)')ylabel('Magnititude')title('|H(e^j^\omega)|')grid onsubplot(212)plot(w/pi,angle(H)/pi)xlabel('\omega(\pi)')ylabel('Phase(\pi)')title('\pheta(\omega)')gridon程序履行如下：幅频.相频：A=[1 1 1 1 1 1 1 1 1]r=roots(A)r =z=[0.7660+0.6428i 0.7660-0.6428i 0.1736+0.9848i 0.1736-0.9848i -0.9397+0.3420i -0.9397-0.3420i-0.5000+0.8660i -0.5000-0.8660i]'p=[0]'零顶点散布图由程序履行可知,该体系具有高通的频率特征.（4）已知一离散时光LTI体系的频率响应H(e^jΩ)如图,输入旌旗灯号为x(n)=cos(0.3πn)+0.5cos(0.8πn).试剖析正弦旌旗灯号sin(Ω0t)经由过程频率响应为H(e^jΩ)的离散时光体系的响应,并依据剖析成果盘算体系对于x（n）的响应y（n）,用MATLAB 绘出体系输入与输出波形.MATLAB程序如下：n=-10:10x=cos(0.3*pi*n)+0.5*cos(0.8*pi*n)stem(n,x)xlabel('n')title('x(n)')输入输出y=2*cos(0.3*pi*n)stem(n,y,‘filled’)xlabel('n')title('y(n)')不雅察试验成果,剖析如图所示体系具有如何的频率特征？从输入输出旌旗灯号上怎么反应出体系的频率特征？答：可知,该体系具有带通的频率特征.（1）盘算该电路体系的频率响应及高通截止Y(ω)=(X(ω))/(1/(1/R+1/jωL)+1/jωC)∙1/(1/R+1/jωL)H(ω)=Y(ω)/X(ω) =【(jω)】^2/(【(jω)】^2+10(jω)+50)令H(ω)=0.707(2)应用matlab绘制幅度响应和相位响应曲线2．已知RC电路如图所示H(ω)=(1⁄RC)/((jω+1)⁄RC)（1）体系的频率响应随RC值的变更RC1=0.001;RC2=0.002;RC3=0.01;a1=1/RC1;a2=1/RC2;a3=1/RC3;A1=[1a1];A2=[1 a2];A3=[1 a3];B1=a1;B2=a2;B3=a3;w1=linspace(0,3000,200);w2=linspace(0,3000,400);w3=linspace(0,3000,500);H1=freqs(B1,A1,w1);H2=freqs(B2,A2,w2);H3=freqs(B3,A3,w3);subplot(311);plot(w1,abs(H1));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega( rad/s)');ylabel('频率响应');title('|H(j\omega)|,RC=0.001');subplot(312);plot(w2,abs(H2));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega( rad/s)');ylabel('频率响应');title('|H(j\omega)|,RC=0.002');subplot(313);plot(w3,abs(H3));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega( rad/s)');ylabel('频率响应');title('|H(j\omega)|,RC=0.0010');（2）滤波前后的时域旌旗灯号波形和频率响应曲线t=0:0.0001:0.2;RC=0.01;a=1/RC;B=[a];A=[1 a];[H,w]=freqs(B,A);f=cos(100*t)+cos(3000*t);Hw1=a/(1i*100+a);Hw2=a/(1i*3000+a);y=abs(Hw1)*cos(100*t)+abs(Hw2)*cos(3000*t); subplot(411);plot(t,f);grid on;xlabel('time(sec)');ylabel('f(t)');title('f(t)=cos(100*t)+cos(3000*t)滤波之前'); subplot(412);plot(t,y);grid on;xlabel('time(sec)');ylabel('y(t)');title('y(t) 滤波之后');subplot(413);plot(w,abs(H));grid on;set(gca,'ytick',[0 0.4 0.707 1]);xlabel('\omega(rad/s)');ylabel('Magnitude');title('|H(j\omega)|');subplot(414);plot(w,angle(H));grid on;xlabel('\omega(rad/s)');ylabel('Phase');title('\phi(\omega)');四,心得领会针在此次试验当中,我感触感染到了matlab在LTI体系剖析中的壮大感化,体系函数让剖析进程变得平常简略.。
信号频域分析LTI系统频域分析
信号的频域分析LTI系统频域分析学院：班级：姓名：学号：指导教师：目录一．实验目的 (1)二．实验原理与方法 (1)三．实验内容 (4)四．实验结果 (5)五．实验总结 (8)一、 实验目的：1、 掌握信号的MATLAB 表示及其可视化方法.2、 掌握信号基本时域运算的MATLAB 实现方法.3、利用MATLAB 分析常用信号，加深对信号时域特性的理解二、实验原理与方法1、连续周期信号的频谱分析实验原理：如果信号满足狄里赫力条件，就可以展开为傅立叶级数形式，即其中T 0 表示基波周期，w 0 为基波频率。
上式定义为周期信号复指数形式的傅里叶级数，系数c k 称为x(t)的傅里叶系数。
同时，傅里叶级数还可以用三角函数的线性组合来表示，即展开为三角函数形式的傅里叶级数。
可见，任何满足狄里赫力条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。
对于第一题中的周期矩形脉冲信号，其傅里叶级数展开式为：该题目中，当A=1，T=1，=0.5时X(t)=1/2 + 2/(n π)sin(n π/2)cos(2πnt);因此可以将信号的频谱化为各次谐波的幅度和相位来表示，根据该信号傅立叶展开式编写MATLAB 命令如下：（1） 该周期矩形信号的傅立叶级数表示为 ：（2）利用MATLAB 绘出前N 次谐波波形，观察N 变化合成信号的变化规律 >> t=-1.5:0.01:1.5;>> N=7; //绘出前7次谐波合成的信号波形 ⎪⎪⎩⎪⎪⎨⎧==⎰∑-∞-∞=T dt e t x T c e c t x t jk k k t jk k 000)(1)(0ωω>> x=zeros(size(t));>> for n= 1:1:Nx=x+2/(n*pi)*sin(n*pi/2)*cos(2*pi*n*t); //计算x的傅里叶级数前N次谐波end>> plot(t,x+0.5); //绘出图形>> title(['N=7']); //设定图形名称>> xlabel('Time(sec)'); //添加坐标轴标注时间单位为秒（3）利用MATLAB汇出周期矩形脉冲信号的频谱，观察参数T和τ变化时对频谱波形的影响。
[精品]连续时间LTI系统的频率特性及频域分析
![[精品]连续时间LTI系统的频率特性及频域分析](https://360docs-1393072327.cos.ap-beijing.myqcloud.com/2bdc77054a73f242336c1eb91a37f111f1850d9f.webp)
[精品]连续时间LTI系统的频率特性及频域分析连续时间LTI系统（Linear Time-Invariant System）是指可用于描述各种物理和工程系统运动规律的动态系统。
它们由一对连续时变系统（如模型、结构和控制）和一对线性运算符构成，其具有因变量（响应）和自变量（输入）之间的线性关联性、时间不变性、结构连续的性质，并且在响应上呈现出定义的平稳性，因而它们在描述众多系统运动规律中被广泛应用。
对于连续时间LTI系统的频域特性的研究，则涉及这些系统的相位特性、幅频特性、切趾特性等。
同时，也要探讨系统中不同频率分量的传输特性，因为有不同频率分量的信号既可以幅频分析也可以相位分析，可以衡量系统不同频率下的相应响应。
由于连续时间LTI系统在有限频率通道内传播信号时发生了部分信号丢失，因此我们引入了频域分析得到系统频响阻抗。
这样一来，它就可以用来测量系统频带上的增益，系统的模态表现，以及系统的传播属性和可控特性。
在频域分析过程中，由于信号可以被分解为离散频率分量，所以对于单个频率分量来说，有关连续时间LTI系统的分析可以比较容易地完成。
一般情况下，每一个频率分量的传播特性由一个线性系数连接，称之为频响函数，可以衡量一个系统的频率响应情况。
总的来说，对于连续时间LTI系统，研究其频率特性及频域分析具有重要的意义。
他可以提供一个系统的相位特性、幅频特性、切趾特性等详细的分析，而且由于信号可以分解为离散频率分量，因此可以很容易地实现频域分析，并衡量一个系统的频率响应情况。
此外，还可以利用频域分析来测量系统的增益，模态表现，以及系统的传播属性和可控特性，进而提高系统的性能，实现性能的优化。
连续时间LTI系统的频率特性及频域分析
实验报告实验项目名称：运用Matlab进行连续时间信号卷积运算（所属课程：信号与系统）学院：电子信息与电气工程学院专业： 10电气工程及其自动化姓名： xx学号： ************指导老师： xxx一、实验目的1、学会运用MATLAB 分析连续系统的频率特性。
2、掌握相关函数的调用。
二、实验原理1、一个连续LTI 系统的数学模型通常用常系数线性微分方程描述，即)()()()()()(01)(01)(t e b t e b t e b t r a t r a t r a m m n n +'++=+'++ (1) 对上式两边取傅里叶变换，并根据FT 的时域微分性质可得：)(])([)(])([0101ωωωωωωE b j b j b R a j a j a m m n n +++=+++101)()()()()(a j a j a b j b j b j E j R j H n n m m ++++++==ωωωωωωω H ( jω )称为系统的频率响应特性，简称系统频率响应或频率特性。
一般H ( jω )是复函数，可表示为：)()()(ωϕωωj e j H j H =其中， )(ωj H 称为系统的幅频响应特性，简称为幅频响应或幅频特性；)(ωϕ称为系统的相频响应特性，简称相频响应或相频特性。
H ( jω )描述了系统响应的傅里叶变换与激励的傅里叶变换间的关系。
H ( jω )只与系统本身的特性有关，与激励无关，因此它是表征系统特性的一个重要参数。
MATLAB 信号处理工具箱提供的freqs 函数可直接计算系统的频率响应的数值解，其语句格式为：H=freqs(b,a,w)其中，b 和a 表示H ( jω )的分子和分母多项式的系数向量；w 为系统频率响应的频率范围，其一般形式为w1:p:w2，w1 为频率起始值，w2 为频率终止值，p 为频率取值间隔。
H 返回w 所定义的频率点上系统频率响应的样值。
理学LTI系统的频域分析PPT学习教案
本节小结
第12页/共13页
例1 例2
第2页/共13页
例3：当激励是ej 0t时，用频域分析法求系统 响应：
Y( j) H( j)2( 0)
H( j0)2( 0)
0
y (t)= H(j 0) ej
t
第3页/共13页
上例的两个推广
☆ 周期信号通过系统的响应
对于fT周(t期) 信号有F：n e jn t n
y(t) FnH ( jn) e jn t n
θ (ω)
H ( j) g2C () e j td
第7页/共13页
• 理想低通filter的冲激响应h(t)
d
c h(t)= ℱ-1[g 2
( ) e -j t
]=
t 1
h(t )
c π
sin c (t td (t td )
)
c
Sa[c
(t
td
)]
t
可见，它实际上是不可实现的非因果 系统。
td
t
π
c
第8页/共13页
• 理想低通filter的阶跃响应g(t)
g(t)=h(t)*(t)
g t
1
1
2
ππ
c c
O
t0Байду номын сангаас
t
tr
上升时间：
tr
第9页/共13页
理想低通滤波器在物理上是不可实现 的，近 似理想 低通滤 波器的 实例
• 一种可实现的低通滤波器
L
v1(t) C
H j
1
R v2(t)
系统无失真传输条件
第5页/共13页
无失真传输系统的频响图： H(j)=Ke－jωtd
LTI系统的频域分析
*
LTI hHale Waihona Puke t)傅里叶变换法=
变换 y(t) ③傅氏 Y (j ω ) 反变换
②傅氏 H(jω )
×
＝
频率响应H(j)可定义为系统零状态响应的傅里叶变 换Y(j)与激励f(t)的傅里叶变换F(j)之比，即
Y ( j ) j[ y ( ) f ( )] Y ( j ) j ( ) H ( j ) e e H ( j ) F ( j ) F ( j )
|H(jω)| 1 ω - ωC 0 ωC θ (ω)
j t d , C e j t d H ( j ) g 2C ( ) e C 0,
(1)冲激响应
h(t ) F g 2c ( )e
1 jtd
c Sa[ ( t t )] c d
f (t )
FT FT
F ( )
1 2
相乘
1 Fs ( ) F ( ns ) Ts n
s
相卷积
T (t )
n
(t nT )
y(t ) K f (t td ) Y ( j ) Ke
其频谱关系为
jtd
F ( j)
(2)无失真传输条件：
系统要实现无失真传输，对系统h(t)，H(j)的要求是： (a)对h(t)的要求： |H(jω)| K h(t)=K(t–td) (b)对H(j)的要求： ω H(j)=Y(j)/F(j)=Ke-jtd 0 即 θ (ω)
例：某LTI系统的H(j)和θ()如图，
若f(t)= 2 + 4cos(5t) + 4cos(10t)，求系统的 响应。
θ(ω) π 1 -10 0 -π |H(jω)|
信号与系统实验报告实验三 连续时间LTI系统的频域分析
实验三 连续时间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),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。
LTI系统的频域分析
y(t ) h(t )* fT (t ) Fn [h(t )*e jnt ] Fn H ( jn) e jnt n n 若
则可推导出
A0 y(t ) H (0) An | H ( jn) | cos[nt n (n)] 2 n 1
h( ) e j d
y(t ) H ( j) e
j t
H ( j )反映了响应y(t)的幅度和相位。
二、一般信号f(t)作用于LTI系统的响应
e
1 2
j t
H(j ) ej t
1 2
齐次性
1 j t F ( j ) e d 2
F(j )H(j ) ej t d
FT [TS (t )] S
n
( n
S
)
如果f(t)是带限信号[即f(t)的频谱在（- m，m) 为有限值，而其余区间为0]。
设f(t)←→F(j)，取样信号fS(t)的频谱函数
1 FS ( j ) F ( j )* S ( nS ) 2 n 1 F[ j ( n S )] TS n
f (t )
FT FT
F ( )
1 2
相乘
1 Fs ( ) F ( ns ) Ts n
s
相卷积
T (t )
n
(t nT )
FT
FT
p( ) s
n Βιβλιοθήκη ( ns
)
在画取样信号fS(t)的频谱时，设定ωS ≥2ωm ,这时 其频谱不发生混叠，因此能设法(如利用低通滤波器)， 从FS(j)中取出F(j)，即从fS(t)中恢复原信号f(t)。 否则将发生混叠，而无法恢复原信号。