信号与系统实验三
信号与系统分析实验信号的频谱分析
实验三信号的频谱分析1方波信号的分解与合成实验1实验目的1. 了解方波的傅立叶级数展开和频谱特性。
2. 掌握方波信号在时域上进行分解与合成的方法。
3. 掌握方波谐波分量的幅值和相位对信号合成的影响。
2 实验设备PC机一台，TD－SAS系列教学实验系统一套。
3 实验原理及内容1. 信号的傅立叶级数展开与频谱分析信号的时域特性和频域特性是对信号的两种不同的描述方式。
对于一个时域的周期信号f(t)，只要满足狄利克莱条件，就可以将其展开成傅立叶级数：如果将式中同频率项合并，可以写成如下形式：从式中可以看出，信号f(t)是由直流分量和许多余弦（或正弦）分量组成。
其中第一项A0/2是常数项，它是周期信号中所包含的直流分量；式中第二项A1cos(Ωt+φ1)称为基波，它的角频率与原周期信号相同，A1是基波振幅，φ1是基波初相角；式中第三项A2cos(Ωt+φ2)称为二次谐波，它的频率是基波的二倍，A2是基波振幅，φ2是基波初相角。
依此类推，还有三次、四次等高次谐波分量。
2. 方波信号的频谱将方波信号展开成傅立叶级数为：n=1，3，5…此公式说明，方波信号中只含有一、三、五等奇次谐波分量，并且其各奇次谐波分量的幅值逐渐减小，初相角为零。
图3-1-1为一个周期方波信号的组成情况，由图可见，当它包含的分量越多时，波形越接近于原来的方波信号，还可以看出频率较低的谐波分量振幅较大，它们组成方波的主体，而频率较高的谐波分量振幅较小，它们主要影响波形的细节。
（a）基波（b）基波＋三次谐波（c）基波＋三次谐波＋五次谐波（d）基波＋三次谐波＋五次谐波＋七次谐波（e）基波＋三次谐波＋五次谐波＋七次谐波＋九次谐波图3-1-1方波的合成3. 方波信号的分解方波信号的分解的基本工作原理是采用多个带通滤波器，把它们的中心频率分别调到被测信号的各个频率分量上，当被测信号同时加到多路滤波器上，中心频率与信号所包含的某次谐波分量频率一致的滤波器便有输出。
信号与系统
《信号与系统》仿真作业实验一：连续信号的表示及可视化：f(t)=δ(t); f(t)=ε(t); f(t)=e at(分别取a>0与a<0);f(t)=R(t); f(t)=Sa(wt); f(t)=sin(2πft);(分别画出不同周期个数的波形)解：（1）f(t)=δ(t)的matlab表示：程序清单如下：》t=-5:0.01:5;k=(0-(-5))/0.01+1;y=zeros(size(t));y(k)=1/(0.01-(-0.01));plot(t,y);title('冲击函数f(t)=δ(t)')画出冲击函数的图形如下：冲击函数f(t)=δ(t)t（2） f(t)=ε(t )的matlab 表示及图形： 程序清单如下： 》t=-5:0.01:5; y=heaviside(t) plot(t,y)画出阶跃函数的图形如下：(3) f(t)=e at 的matlab 表示及图形： 程序清单如下： 》t=-10:0.01:10;y1=exp(0.1*t); y2=exp(-0.1*t); plot(t,y1,'r',t,y2,'b') 画出指数函数的图形如下：tf (t )=ε(t )(4) f(t)=R(t)的matlab 表示及图形： 程序清单如下： 》t=-5:0.01:5;y=heaviside(t+2)-heaviside(t-2); plot(t,y,'b') 画出窗函数的图形如下：(5) f(t)=Sa(wt) 的matlab 表示及图形： 程序清单如下：》ezplot('sin(t)./t',[-20,20]) grid ontf (t )=e atty =R 9t )画出抽样函数的图形如下：sin(t)/tt（6）f(t)=sin(2πft)的matlab表示及图形：程序清单如下：》ezplot('sin(2*pi*50*t)',[-.02,.02])grid on画出正弦函数的图形如下：实验二：离散信号的表示及可视化：f(t)=δ(n ); f(t)=ε(n ); f(t)=e an (分别取a>0与a<0); f(t)=R N (n ); f(t)=Sa(nw); f(t)=sin(nw );(分别取不同的w 值) 解：（1） 冲击序列f(n)=δ(n )的matlab 实现： 程序清单如下： 》n0=0; ns=-10; nf=10; n=[ns:nf];y=[zeros(1,n0-ns),1,zeros(1,nf-n0)];-0.02-0.015-0.01-0.00500.0050.010.0150.02-1-0.50.51tsin(2 50 t)stem(n,y);title('冲击序列f(n)=δ(n)')画出冲击序列的图形如下：冲击序列f(n)=δ(n)n（2）阶跃序列f(n)=ε(n)的matlab实现：程序清单如下：》n0=0;ns=-10;nf=10;n=[ns:nf];y=[zeros(1,n0-ns),ones(1,nf-n0+1)];stem(n,y);title('阶跃序列f(n)=ε(n)')阶跃序列的图形如下：（3） 指数序列f(t)=e an (分别取a>0与a<0)的matlab 实现： 程序清单如下： 》n=-10:10; y1=exp(0.1*n); y2=exp(-0.1*n); plot(n,y1,'ro',n,y2,'bo') 指数序列的图形如下：（4） 门序列f(n)=R N (n )的matlab 实现：程序清单如下： 》n1=-3;n2=3;ns=-15;nf=15;阶跃序列f(n)=ε(n)nnf (t )=e a nn=[ns:nf];y=[zeros(1,n1-ns),ones(1,n2-n1+1),zeros(1,nf-n2)]; stem(n,y);title('窗序列f(n)=R N (n )') 窗序列的图形如下：（5） 抽样序列f(t)=Sa(nw)的matlab 实现： 》n=-20:0.5:20; y=sin(n)./n; plot(n,y,'o'); title('f(t)=Sa(nw)')窗函数f(n)=R N (n)n抽样序列的图形如下：（6） 正弦序列f(t)=sin(nw )(分别取不同的w 值)的matlab 实现： 》n=-0.1:0.002:0.1 w=100 y=sin(w*n) plot(n,y,'o') grid on正弦序列的图形如下：f (t)=Sa(nw)nny =s i n (w *n )实验三：系统的时域求解1、设h(n)=(0.9)n u(n),x(n)=u(n)-u(n-10),求：y(n)=x(n)*h(n),并画出x(n),h(n),y(n)波形。
《信号与系统》课程实验报告
《信号与系统》课程实验报告《信号与系统》课程实验报告一图1-1 向量表示法仿真图形2．符号运算表示法若一个连续时间信号可用一个符号表达式来表示，则可用ezplot命令来画出该信号的时域波形。
上例可用下面的命令来实现（在命令窗口中输入，每行结束按回车键）。
t=-10:0.5:10;f=sym('sin((pi/4)*t)');ezplot(f,[-16,16]);仿真图形如下：图1-2 符号运算表示法仿真图形三、实验内容利用MATLAB实现信号的时域表示。
三、实验步骤该仿真提供了7种典型连续时间信号。
用鼠标点击图0-3目录界面中的“仿真一”按钮，进入图1-3。
图1-3 “信号的时域表示”仿真界面图1-3所示的是“信号的时域表示”仿真界面。
界面的主体分为两部分：1) 两个轴组成的坐标平面（横轴是时间，纵轴是信号值）；2) 界面右侧的控制框。
控制框里主要有波形选择按钮和“返回目录”按钮，点击各波形选择按钮可选择波形，点击“返回目录”按钮可直接回到目录界面。
图1-4 峰值为8V，频率为0.5Hz，相位为180°的正弦信号图1-4所示的是正弦波的参数设置及显示界面。
在这个界面内提供了三个滑动条，改变滑块的位置，滑块上方实时显示滑块位置代表的数值，对应正弦波的三个参数：幅度、频率、相位；坐标平面内实时地显示随参数变化后的波形。
在七种信号中，除抽样函数信号外，对其它六种波形均提供了参数设置。
矩形波信号、指数函数信号、斜坡信号、阶跃信号、锯齿波信号和抽样函数信号的波形分别如图1-5～图1-10所示。
图1-5 峰值为8V，频率为1Hz，占空比为50%的矩形波信号图1-6 衰减指数为2的指数函数信号图1-7 斜率=1的斜坡信号图1-8 幅度为5V，滞后时间为5秒的阶跃信号图1-9 峰值为8V，频率为0.5Hz的锯齿波信号图1-10 抽样函数信号仿真途中，通过对滑动块的控制修改信号的幅度、频率、相位，观察波形的变化。
《信号与系统》实验三
三：
源程序：
（1）：τ/T=1/4时的周期矩形脉冲的幅度谱和相位谱：
n=-20:20;
F=zeros(size(n));
forii=-20:20
F(ii+21)= sin(ii*pi/4)/(ii*pi+eps);
end
F(21)=1/4;
实验
内容
1.求图1所示周期信号（ ， ）的傅里叶级数，用Matlab做出其前3、9、21、45项谐波的合成波形与原信号作比较，并做出其单边幅度谱和相位谱。
图1 周期为2的三角脉冲信号
2. 求图2所示的单个三角脉冲（ ）的傅里叶变换，并做出其幅度谱和相位谱。
图2 单个三角脉冲
3. 求不同占空比下周期矩形脉冲的幅度谱和相位谱，例如 、 。
y=1/4;
forn=1:m
y=y+4/(n*n*pi*pi)*(1-cos(n*pi/2)).*cos(n*pi.*t);
end
源代码：
t=-6:0.01:6;
d=-6:2:6;
fxx=pulstran(t,d,'tripuls');
f1=fourierseries(3,t);
f2=fourierseries(9,t);
n=1:10;
a=zeros(size(n));
fori=1:10
a(i)=angle(4/(i*i*pi*pi)*(1-cos(i*pi/2)))
end
n=0:pi:9*pi
stem(n,a,'fill','linewidth',2);
axis([0,9*pi,-0.2,0.2])
信号与系统实验_矩形信号的分解
学号： 姓名：实验三、矩形信号的分解一、实验目的1、分析典型的矩形脉冲信号，了解矩形脉冲信号谐波分量的构成；2、观察矩形脉冲信号分解出各谐波分量的情况。
二、预备知识1.学习“周期信号的傅里叶级数分析”一节；2．复习matlab 软件的使用方法。
3．信号的滤波知识三、实验原理1、信号的频谱与测量信号的时域特性和频域特性是对信号的两种不同的描述方式。
对于一个时域的周期信号)t (f ，只要满足狄利克莱(Dirichlet)条件，就可以将其展开成三角形式或指数形式的傅里叶级数。
例如，对于一个周期为T 的时域周期信号)t (f ，可以用三角形式的傅里叶级数求出它的各次分量，在区间)T t ,t (11+内表示为)sin cos ()(10t n b t n a a t f n n n Ω+Ω+=∑∞=即将信号分解成直流分量及许多余弦分量和正弦分量，研究其频谱分布情况。
AA（c）图3-1 信号的时域特性和频域特性信号的时域特性与频域特性之间有着密切的内在联系，这种联系可以用图3-1来形象地表示。
其中图3-1(a)是信号在幅度--时间--频率三维座标系统中的图形；图3-1(b)是信号在幅度--时间座标系统中的图形即波形图；把周期信号分解得到的各次谐波分量按频率的高低排列，就可以得到频谱图。
反映各频率分量幅度的频谱称为振幅频谱。
图3-1(c)是信号在幅度--频率座标系统中的图形即振幅频谱图。
反映各分量相位的频谱称为相位频谱。
在本实验中只研究信号振幅频谱。
周期信号的振幅频谱有三个性质：离散性、谐波性、收敛性。
测量时利用了这些性质。
从振幅频谱图上，可以直观地看出各频率分量所占的比重。
测量方法有同时分析法和顺序分析法。
2、 矩形脉冲信号的频谱一个幅度为E ，脉冲宽度为τ，重复周期为T 的矩形脉冲信号，如图10-3所示。
图3-2 周期性矩形脉冲信号其傅里叶级数为：t n Tn Sa T E T E t f n i ωπτττcos )(2)(1∑=+= 该信号第n 次谐波的振幅为：Tn T n T E T n Sa T E a n /)/sin(2)(2τπτπττπτ== 由上式可见第n 次谐波的振幅与E 、T 、τ有关。
信号与系统实验指导
二、实验仪器
1、双踪示波器 1台
2、信号源及频率计模块S2 1块
3、数字信号处理模块S4 1块
三、实验原理
卷积积分的物理意义是将信号分解为冲激信号之和，借助系统的冲激响应，求解系统对任意激励信号的零状态响应。设系统的激励信号为 ，冲激响应为 ，则系统的零状态响应为：
2、对实现无失真传输，对系统函数 应提出怎样的要求？
设 与 的傅立叶变换式分别为 。借助傅立叶变换的延时定理，从式4-1可以写出
（4-2）
此外还有 （4-3）
所以，为满足无失真传输应有
（4-4）
式（4-4）就是对于系统的频率响应特性提出的无失真传输条件。欲使信号在通过线性系统时不产生任何失真，必须在信号的全部频带内，要求系统频率响应的幅度特性是一常数，相位特性是一通过原点的直线。
线性系统的幅度失真与相位失真都不产生新的频率分量。而对于非线性系统则由于其非线性特性对于所传输信号产生非线性失真，非线性失真可能产生新的频率分量。
所谓无失真是指响应信号与激励信号相比，只是大小与出现的时间不同，而无波形上的变化。设激励信号为 ，响应信号为 ，无失真传输的条件是
（4-1）
式中 是一常数， 为滞后时间。满足此条件时， 波形是 波形经 时间的滞后，虽然，幅度方面有系数 倍的变化，但波形形状不变。
图6-1 连续信号抽样过程
将连续信号用周期性矩形脉冲抽样而得到抽样信号，可通过抽样器来实现，实验原理电路如图6-2所示。
2、连续周期信号经周期矩形脉冲抽样后，抽样信号的频谱
它包含了原信号频谱以及重复周期为fs（f s = s/2л）、幅度按 Sa（m sτ/2）规律变化的原信号频谱，即抽样信号的频谱是原信号频谱的周期性延拓。因此，抽样信号占有的频带比原信号频带宽得多。
实验三 连续信号与系统的频域分析
学号
0174280
同组人：无
实验项目
实验三连续信号与系统的频域分析
☑必修□选修
□演示性实验☑验证性实验□操作性实验□综合性实验
实验地点
H113
实验仪器台号
F0
指导教师
蒋娜
实验日期及节次
week14->2-12
一、实验目的及要求：
1、目的
1．掌握非周期信号的傅里叶变换：fourier函数和ifourier函数；
四、实验结果与数据处理：
1．利用fourier函数求下列信号的傅里叶变换F(jω)，并用ezplot函数绘出其幅度谱和相位谱。
（1）
syms t v w phase im re;%定义变量t,v,w,phase,im re
f=sym('Heaviside(t)-Heaviside(t-2)');%
Fw=fourier(f);
plot([07.0711],[0.7070.707],':');
axis([04001.1]);
grid;
xlabel('角频率(\omega)');
ylabel('幅度');
title('H(j\omega)的幅频特性');
subplot(212);
plot(w,h2*180/pi);
axis([0400200]);
（2）
syms t v w phase im re;%定义变量t,v,w,phase,im re
f=exp(-1*t)*sym('Heaviside(t)');%
Fw=fourier(f);
subplot(311);
信号与系统实验报告-实验3--周期信号的频谱分析
信号与系统实验报告-实验3--周期信号的频谱分析信号与系统实验报告实验三周期信号的频谱分析实验三周期信号的频谱分析实验目的：1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法；2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；3、掌握各种典型的连续时间非周期信号的频谱特征。
实验内容：（1）Q3-1 编写程序Q3_1，绘制下面的信号的波形图：其中，0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制cos(0t)、cos(30t)、cos(50t) 和x(t) 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。
程序如下:clear,%Clear all variablesclose all,%Close all figure windowsdt = 0.00001; %Specify the step of time variable t = -2:dt:4; %Specify the interval of timew0=0.5*pi; x1=cos(w0.*t); x2=cos(3*w0.*t);x3=cos(5*w0.*t);N=input('Type in the number of the harmonic components N=');x=0;for q=1:N;x=x+(sin(q*(pi/2)).*cos(q*w0*t))/q;endsubplot(221)plot(t,x1)%Plot x1axis([-2 4 -2 2]);grid on,title('signal cos(w0.*t)')subplot(222)plot(t,x2)%Plot x2axis([-2 4 -2 2]); grid on,title('signal cos(3*w0.*t))')subplot(223)plot(t,x3)%Plot x3axis([-2 4 -2 2])grid on,title('signal cos(5*w0.*t))')subplot(224)plot(t,x)%Plot xtaxis([-2 4 -2 2])grid on,title('signal xt')（2）给程序3_1增加适当的语句，并以Q3_2存盘，使之能够计算例题1中的周期方波信号的傅里叶级数的系数，并绘制出信号的幅度谱和相位谱的谱线图。
实验三一阶网络频响特性测量,信号与系统,南京理工大学紫金学院实验报告
信号与系统实验报告实验名称：一阶网络频响特性测量姓名：学号：班级：通信时间：2013.6南京理工大学紫金学院电光系一、 实验目的1、 掌握一阶网络的构成方法；2、 掌握一阶网络的系统响应特性；3、 了解一阶网络频响特性图的测量方法；二、实验基本原理系统响应特性是指系统在正弦信号激励下，稳态响应随信号频率变化而变化的特性，称为系统的频率响应特性（frequency response ）简称频响特性。
一阶系统是构成复杂系统的基本单元。
学习一阶系统的特点有助于对一般系统特性的了解。
一阶系统的系统函数为H(s)，表达式可以写成：γ+⋅=s k s H 1)( k 为一常数 (3-1) 激励信号x(t)为：(3-2)按照系统频响特性的定义可求得该一阶系统的稳态响应为：(3-3)其中⎣⎦00)()(|)(00ϕj j s ej H j H s H Ω=Ω=Ω=，⎣⎦)(00Ω=j H H 。
可见，当改变系统输入信号的频率时，稳态响应的幅度和相位也随之而改变。
因果系统是稳定的要求：0>γ，不失一般性可设τγ1==k 。
该系统的频响特性为：11)(+Ω=Ωτj j H (3-4)从其频响函数中可以看出系统响应呈低通方式，其3dB 带宽点τ1。
系统的频响特性图如下图：0()sin()m x t E t =Ω000()sin()ss m y t E H t ϕ=Ω+θ图1 一阶网络频响特性图一阶低通系统的单位冲击响应与单位阶跃响应如下图：图2 一阶网络单位冲击响应与单位阶跃响应图三、实验内容及结果一阶系统的幅度谱一阶系统相位谱3、用矢量作图法作出该一阶系统的幅度谱和相位谱。
一阶系统的幅度谱一阶系统的相位谱4、作出一阶网络的单位阶跃响应波形，标注在阶跃响应最大值的(1-e-1)倍处的时间t的值，与理论值R1C1是否相符。
四、实验分析1、实验所得一阶网络的频响特性图和用矢量作图法所得的频响特性图有何异同？说明原因。
信号与系统matlab实验习题3 绘制典型信号及其频谱图
绘制典型信号及其频谱图答案在下面四个常用信号及其傅里叶变换式如表1所示。
(1)绘制单边指数信号及其频谱图的MATLAB程序如下：close all;E=1;a=1;t=0:0.01:4;w=-30:0.01:30;f=E*exp(-a*t);F=1./(a+j*w);plot(t,f);xlabel('t');ylabel('f(t)');figure;plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|');figure;max_logF=max(abs(F));plot(w,20*log10(abs(F)/max_logF));xlabel('\omega');ylabel('|F(\omega)| indB');figure;plot(w,angle(F));xlabel('\omega');ylabel('\phi(\omega)');请更改参数，调试此程序，绘制单边指数信号的波形图和频谱图。
观察参数a 对信号波形及其频谱的影响。
注：题目中阴影部分是幅频特性的对数表示形式，单位是(dB)，请查阅相关资料，了解这种表示方法的意义及其典型数值对应的线性增益大小。
(2)绘制矩形脉冲信号、升余弦脉冲信号和三角脉冲信号的波形图和频谱图，观察并对比各信号的频带宽度和旁瓣的大小。
(3)更改参数，调试程序，绘制单边指数信号的波形图和频谱图。
观察参数a对信号波形及其频谱的影响。
答案附上程序代码：close all;E=1;a=1;t=0:0.01:4;w=-30:0.01:30;f=E*exp(-a*t);F=1./(a+j*w);plot(t,f);xlabel('t');ylabel('f(t)');figure;plot(w,abs(F));xlabel('\omega');ylabel('|F(\omega)|';E=1，a=1,波形图 频谱图更改参数E=2，a=1;更改参数a ，对信号波形及其频谱的影响。
信号与系统实验三 信号的采样与恢复
实验三 信号的采样与恢复一、实验目的1、了解电信号的采样方法与过程以及信号恢复的方法。
2、验证抽样定理。
二、实验仪器1、信号与系统实验箱一台（主板）。
2、系统时域与频域分析模块一块。
3、20M 双踪示波器一台。
三、实验内容、过程及结果1)实验内容：观察低中高三种频率下不混叠时（即f ≥2B ）原信号与抽样信号以及抽样恢复信号的波形然后进行对比。
2)实验步骤：1、把系统时域与频域分析模块插在主板上，用导线接通此模块“电源接入”和主板上的电源（看清标识，防止接错，带保护电路），并打开此模块的电源开关。
2、将函数信号发生器产生一正弦波（幅度（峰值）为2V 左右，为便于观察，抽样信号频率一般选择50HZ ～400HZ 的范围，抽样脉冲的频段由开关SK1000进行选择，有“高”“中”“低”档，频率则是通过电位器“频率调节”来调节的，抽样脉冲的脉宽则是由电位器“脉宽调节”进行调节的（一般取30％）），将其送入抽样器，即用导线将函数信号发生器的输出端与本实验模块的输入端相连，用示波器测试“抽样信号”的波形，观察经抽样后的正弦波。
3、改变抽样脉冲的频率为B f s 2 ，用导线将“抽样信号”和“低通输入”相连，用示波器测试测试钩“抽样恢复”，观察复原后的信号，比较其失真程度。
3)实验结果：①低频下：原信号与抽样信号 原信号与抽样恢复信号②中频下：原信号与抽样信号原信号与抽样恢复信号③高频下：原信号与抽样信号原信号与抽样恢复信号四、实验结果分析1）由原信号、抽样信号以及复原信号的波形，能得出什么结论？答：抽样信号是从原信号中获得的离散周期性的信号，其包含了部分乃至绝大部分的原信号内容，通过对这些抽样信号内容进行还原，就可以得到近似原信号波形的结果，但是不能得到跟原信号完全一致的波形，因为失真无法完全避免，只能调试到最佳结果。
2）比较三种不同抽样频率下的fs(t)的波形，能得出什么结论？答：当fs<2B时，抽样信号的频谱会发生混叠，从发生混叠后的频谱中无法用低通滤波器获得原信号频谱的全部内容，即使fs=2B，复原后的信号失真还是难免的。
信号与系统实验教程只有答案
信号与系统实验教程（只有答案））（实验报告目录实验一信号与系统的时域分析 (2)三、实验内容及步骤 (2)实验二连续时间信号的频域分析 (14)三、实验内容及步骤 (14)实验三连续时间LTI系统的频域分析 (35)三、实验内容及步骤 (35)实验四通信系统仿真 (42)三、实验内容及步骤 (42)实验五连续时间LTI系统的复频域分析 (51)三、实验内容及步骤 (51)实验一信号与系统的时域分析三、实验内容及步骤实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。
实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。
并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。
实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。
Q1-1：修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别，哪一幅图形看起来与实际信号波形更像？答：Q1-2：修改程序Program1_1，并以Q1_2为文件名存盘，产生实指数信号x(t)=e-0.5t。
要求在图形中加上网格线，并使用函数axis()控制图形的时间范围在0~2秒之间。
然后执行该程序，保存所的图形。
修改Program1_1后得到的程序Q1_2如下：信号x(t)=e-0.5t的波形图clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = exp(-0.5*t); % Generate the signalplot(t,x)grid on;axis ([0 2 0 1 ])title('Sinusoidal signal x(t)')xlabel('Time t (sec)')Q1-3：修改程序Program1_1，并以Q1_3为文件名存盘，使之能够仿真从键盘上任意输入的一个连续时间信号，并利用该程序仿真信号x(t)=e-2t。
信号与系统实验
实验一 抽样定理与信号恢复一、实验目的1. 观察离散信号频谱，了解其频谱特点；2. 验证抽样定理并恢复原信号。
二、实验原理1. 离散信号不仅可从离散信号源获得，而且也可从连续信号抽样获得。
抽样信号 Fs （t ）=F （t ）·S （t ）。
其中F （t ）为连续信号（例如三角波），S （t ）是周期为Ts 的矩形窄脉冲。
Ts 又称抽样间隔，Fs=1Ts 称抽样频率，Fs （t ）为抽样信号波形。
F （t ）、S （t ）、Fs （t ）波形如图1-1。
t-4T S -T S 0T S 4T S8T S 12T S tt02/1τ1τ2/31τ2/1τ1τ2/31τ2/1τ-(a)(b)(c)图1-1 连续信号抽样过程将连续信号用周期性矩形脉冲抽样而得到抽样信号，可通过抽样器来实现，实验原理电路如图1-2所示。
2. 连续周期信号经周期矩形脉冲抽样后，抽样信号的频谱()∑∞∞--∙=m s s m m SaTsA j )(22s F ωωπδτωτω 它包含了原信号频谱以及重复周期为fs （f s =πω2s 、幅度按ST A τSa （2τωs m ）规律变化的原信号频谱，即抽样信号的频谱是原信号频谱的周期性延拓。
因此，抽样信号占有的频带比原信号频带宽得多。
以三角波被矩形脉冲抽样为例。
三角波的频谱 F （j ω）=∑∞-∞=-K k k sa E )2()2(12τπωδππ抽样信号的频谱Fs （j ω）=式中 取三角波的有效带宽为31ω18f f s =作图，其抽样信号频谱如图1-3所示。
图1-2 信号抽样实验原理图)(2(212s m k s m k k Sa m Sa TS EA ωωωδπτωτπ--∙∙∑∞-∞=-∞=111112ττπω==f 或(b) 抽样信号频谙图1-3 抽样信号频谱图如果离散信号是由周期连续信号抽样而得，则其频谱的测量与周期连续信号方法相同，但应注意频谱的周期性延拓。
信号与系统实验报告
信号与系统实验报告目录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．熟悉一阶系统的无源和有源电路；2．研究一阶系统时间常数T的变化对系统性能的影响；3．研究一阶系统的零点对系统响应的影响。
二、实验原理1．无零点的一阶系统无零点一阶系统的有源和无源电路图如图2-1的(a)和(b)所示。
它们的传递函数均为：10.2s1G(s)＝+(a) 有源(b) 无源图2-1 无零点一阶系统有源、无源电路图2．有零点的一阶系统(|Z|<|P|)图2-2的(a)和(b)分别为有零点一阶系统的有源和无源电路图，它们的传递函数为：10.2s1)0.2(sG(s)++=，⎪⎪⎪⎪⎭⎫⎝⎛++=S611S161G(s）(a) 有源(b) 无源图2-2 有零点(|Z|<|P|)一阶系统有源、无源电路图3．有零点的一阶系统(|Z|>|P|)图2-3的(a)和(b)分别为有零点一阶系统的有源和无源电路图，它们的传递函数为：1s10.1sG(s)＝++(a) 有源(b) 无源图2-3 有零点(|Z|>|P|)一阶系统有源、无源电路图三、实验步骤1．打开THKSS-A/B/C/D/E型信号与系统实验箱，将实验模块SS02插入实验箱的固定孔中，利用该模块上的单元组成图2-1(a)(或(b))所示的一阶系统模拟电路。
2．实验线路检查无误后，打开实验箱右侧总电源开关。
3．将“阶跃信号发生器”的输出拨到“正输出”，按下“阶跃按键”按钮，调节电位器RP1，使之输出电压幅值为1V，并将“阶跃信号发生器”的“输出”端与电路的输入端“Ui”相连，电路的输出端“Uo”接到双踪示波器的输入端，然后用示波器观测系统的阶跃响应，并由曲线实测一阶系统的时间常数T。
4．再依次利用实验模块上相关的单元分别组成图2-2(a)(或(b))、2-3(a)(或(b))所示的一阶系统模拟电路，重复实验步骤3，观察并记录实验曲线。
信号与系统实验报告实验三 连续时间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),反映信号经过系统后,信号各频率分量在相位上发生变换的情况。
信号与系统Matlab实验作业
实验一典型连续时间信号和离散时间信号一、实验目的掌握利用Matlab画图函数和符号函数显示典型连续时间信号波形、典型时间离散信号、连续时间信号在时域中的自变量变换。
二、实验内容1、典型连续信号的波形表示（单边指数信号、复指数信号、抽样信号、单位阶跃信号、单位冲击信号）1）画出教材P28习题1-1(3) ()[(63)(63)]t=----的波形图。
f t e u t u t2）画出复指数信号()()j t f t e σω+=当0.4, 8σω==(0<t<10)时的实部和虚部的波形图。
t=0:0.01:10;f1='exp(0.4*t)*cos(8*t)';f2='exp(0.4*t)*sin(8*t)';figure(1)ezplot(f1,t);grid on;figure(2)ezplot(f2,t);grid on;3）画出教材P16图1-18，即抽样信号Sa(t)的波形(-20<t<20)。
t=-10:0.01:10;f='sin(t)/t';ezplot(f,t);grid on;4）用符号函数sign画出单位阶跃信号u(t-3)的波形（0<t<10）。
t=0:0.01:10;f='(sign(t-3)+1)/2';ezplot(f,t);grid on;5）单位冲击信号可看作是宽度为∆，幅度为1/∆的矩形脉冲，即t=t 1处的冲击信号为11111 ()()0 t t t x t t t otherδ∆⎧<<+∆⎪=-=∆⎨⎪⎩画出0.2∆=, t 1=1的单位冲击信号。
t=0:0.01:2;f='5*(u(t-1)-u(t-1.2))';ezplot(f,t);grid on;axis([0 2 -1 6]);2、典型离散信号的表示（单位样值序列、单位阶跃序列、实指数序列、正弦序列、复指数序列）编写函数产生下列序列：1）单位脉冲序列，起点n0，终点n f，在n s处有一单位脉冲。
《信号与系统》实验报告
《信号与系统》实验报告目录一、实验概述 (2)1. 实验目的 (2)2. 实验原理 (3)3. 实验设备与工具 (4)二、实验内容与步骤 (5)1. 实验一 (6)1.1 实验目的 (7)1.2 实验原理 (7)1.3 实验内容与步骤 (8)1.4 实验结果与分析 (9)2. 实验二 (10)2.1 实验目的 (12)2.2 实验原理 (12)2.3 实验内容与步骤 (13)2.4 实验结果与分析 (14)3. 实验三 (15)3.1 实验目的 (16)3.2 实验原理 (16)3.3 实验内容与步骤 (17)3.4 实验结果与分析 (19)4. 实验四 (20)4.1 实验目的 (20)4.2 实验原理 (21)4.3 实验内容与步骤 (22)4.4 实验结果与分析 (22)三、实验总结与体会 (24)1. 实验成果总结 (25)2. 实验中的问题与解决方法 (26)3. 对信号与系统课程的理解与认识 (27)4. 对未来学习与研究的展望 (28)一、实验概述本实验主要围绕信号与系统的相关知识展开，旨在帮助学生更好地理解信号与系统的基本概念、性质和应用。
通过本实验，学生将能够掌握信号与系统的基本操作，如傅里叶变换、拉普拉斯变换等，并能够运用这些方法分析和处理实际问题。
本实验还将培养学生的动手能力和团队协作能力，使学生能够在实际工程中灵活运用所学知识。
本实验共分为五个子实验，分别是：信号的基本属性测量、信号的频谱分析、信号的时域分析、信号的频域分析以及信号的采样与重构。
每个子实验都有明确的目标和要求，学生需要根据实验要求完成相应的实验内容，并撰写实验报告。
在实验过程中，学生将通过理论学习和实际操作相结合的方式，逐步深入了解信号与系统的知识体系，提高自己的综合素质。
1. 实验目的本次实验旨在通过实践操作，使学生深入理解信号与系统的基本原理和概念。
通过具体的实验操作和数据分析，掌握信号与系统分析的基本方法，提高解决实际问题的能力。
信号与系统 matlab实验报告
信号与系统 matlab实验报告《信号与系统 Matlab实验报告》摘要：本实验报告通过使用 Matlab 软件进行信号与系统实验，探讨了信号与系统在数字领域的应用。
实验结果表明，Matlab 软件具有强大的信号处理和系统分析功能，能够有效地进行信号与系统的模拟和分析。
引言：信号与系统是电子工程领域中的重要基础课程，它研究了信号的产生、传输和处理，以及系统对信号的响应和影响。
在数字领域，信号与系统的理论和方法也得到了广泛的应用。
Matlab 软件作为一种强大的数学计算工具，为信号与系统的模拟和分析提供了便利和高效的途径。
实验一：信号的生成与显示在本实验中，我们首先使用 Matlab 软件生成了几种常见的信号，包括正弦信号、方波信号和三角波信号。
通过调整信号的频率、幅度和相位等参数，我们观察了信号的变化，并利用 Matlab 的绘图功能将信号图形显示出来。
实验结果表明，Matlab 软件能够方便地生成各种类型的信号，并能够直观地显示信号的波形和特性。
实验二：信号的采样与重构在本实验中，我们使用 Matlab 软件对信号进行了采样和重构。
我们首先对一个连续信号进行了离散采样，然后利用 Matlab 的插值函数对采样信号进行了重构。
实验结果表明，采样和重构过程中存在信号失真和频率混叠等问题，但通过适当的采样和重构方法，我们能够有效地还原原始信号。
实验三：系统的响应与分析在本实验中，我们使用 Matlab 软件对系统的响应进行了分析。
我们构建了几种常见的系统模型，包括线性时不变系统和滤波器系统，然后利用 Matlab 的系统分析工具对系统的频率响应、相位响应和单位脉冲响应等进行了分析。
实验结果表明，Matlab 软件能够有效地进行系统的模拟和分析，为系统设计和优化提供了有力的支持。
结论：通过本实验，我们深入了解了信号与系统在数字领域的应用，并掌握了使用 Matlab 软件进行信号与系统模拟和分析的方法。
信号与系统实验报告实验三连续时间LTI系统的频域分析
实验三 连续时间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 ），反映信号经过系统后，信号各频率分量在相位上发生变换的情况。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
信号与系统实验报告实验三连续时间LTI系统的频域分析学院：信息工程学院班级：2012级电子信息工程三班姓名：学号：2012550711指导老师：苏永新一、实验目的1、掌握系统频率响应特性的概念及其物理意义；2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用；3、学习和掌握幅度特性、相位特性以及群延时的物理意义；4、掌握用MATLAB 语言进行系统频响特性分析的方法。
基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。
二、实验内容及步骤实验前，必须首先阅读本实验原理，了解所给的MATLAB 相关函数，读懂所给出的全部范例程序。
实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。
并结合范例程序所完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。
实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。
给定三个连续时间LTI 系统，它们的微分方程分别为系统1： dt t dx t y dt t dy dt t y d )()(25)(1)(22=++ Eq.3.1系统2： )()()()(t x dtt dx t y dt t dy -=+ Eq.3.2 系统3：)(262)(262)(401)(306)(148)(48)(10)(2233445566t x t y dt t dy dtt y d dt t y d dt t y d dt t y d dt t y d =++++++ Eq.3.3Q3-1 修改程序Program3_1，并以Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。
并利用该程序计算并绘制由微分方程Eq.3.1、Eq.3.2和Eq.3.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。
抄写程序Q3_1如下： % Program Q3_1% This Program is used to compute and draw the plots of the frequency response % of a continuous-time system clc,clear,close all ;b = input('Type in the right coefficient vector of differential equation b£º'); % The coefficient vector of the right side of the differential equationa = input('Type in the left coefficient vector of differential equation a£º'); % The coefficient vector of the left side of the differential equation [H,w] = freqs(b,a); % Compute the frequency response H Hm = abs(H); % Compute the magnitude response Hm phai = angle(H); % Compute the phase response phaiHr = real(H); % Compute the real part of the frequency responseHi = imag(H); % Compute the imaginary part of the frequency responsesubplot(221)plot(w,Hm), grid on , title('Magnitude response'), xlabel('Frequency in rad/sec') subplot(223)plot(w,phai), grid on , title('Phase response'), xlabel('Frequency in rad/sec') subplot(222)plot(w,Hr), grid on , title('Real part of frequency response'), xlabel('Frequency in rad/sec') subplot(224)plot(w,Hi), grid on , title('Imaginary part of frequency response'), xlabel('Frequency in rad/sec') %Eq.3.1 b=[1 0];a=[1 1 25]; %Eq.3.2 b=[1 -1];a=[1 1];%Eq.3.3 b=[262];a=[1 10 48 148 306 401 262];执行程序Q3_1，绘制的系统1的频率响应特性曲线如下：51000.10.20.30.4Magnitude responseFrequency in rad/sec 0510-4-3-2-10Phase responseFrequency in rad/sec510-0.100.10.20.3Real part of frequency responseFrequency in rad/sec510-0.3-0.2-0.10Imaginary part of frequency response Frequency in rad/sec从系统1的幅度频率响应曲线看，系统1是低通、高通、全通、带通还是带阻滤波器？ 答：是带通的滤波器。
执行程序Q3_1，绘制的系统2的频率响应特性曲线如下：5101111Magnitude responseFrequency in rad/sec01234Phase responseFrequency in rad/sec510-1-0.500.51Real part of frequency responseFrequency in rad/sec0.51Imaginary part of frequency response Frequency in rad/sec从系统2的幅度频率响应曲线看，系统2低通、高通、全通、带通还是带阻滤波器？ 答：是低通的滤波器。
执行程序Q3_1，绘制的系统3的频率响应特性曲线如下：Magnitude responseFrequency in rad/sec510Phase responseFrequency in rad/sec-1-0.500.51Real part of frequency responseFrequency in rad/sec510-1-0.500.51Imaginary part of frequency response Frequency in rad/sec从系统3的幅度频率响应曲线看，系统3是低通、高通、全通、带通还是带阻滤波器？ 答：是带阻滤波器。
这三个系统的幅度频率响应、相位频率相应、频率响应的实部以及频率响应的虚部分别具有何种对称关系？请根据傅里叶变换的性质说明为什么会具有这些对称关系？答：H （j ） 是 ω 的偶函数 θ（ ω ）是 ω 的奇函数，系统 1 的幅度频率响应和频率响应的实部是关于 w 对称的偶函数，相位频率响应和频率响应的虚部是关于原点对称的奇函数；系统 2 的幅度频率响应是一条关于 w 一值对称的直线，频率响应的实部，频率响应的虚部以及相位频率响应都没有一定的对称关系；系统 3 的幅度频率响应和频率响应的实部没有一定的对称关系，而相位频率响应和频率响应的虚部在低频段是关于 w 一值的奇函数，随着频率的逐渐变大，没有这种对称关系。
根据傅里叶变换的性质，如果 f（t）为实的奇函数，则对应的傅里叶变换 F（w）为虚的奇函数，如果 f（t）为虚的奇函数，则对应的傅里叶变换 F（w）为虚的奇函数，如果 ft为实的偶函数，则对应的傅里叶变换 F（w）为实的偶函数，如果 f（t）为实的奇函数，则对应的傅里叶变换 F（w）为虚的奇函数Q3-2编写程序Q3_2，使之能够能够接受键盘方式输入的输入信号x(t)的数学表达式，系统微分方程的系数向量，计算输入信号的幅度频谱，系统的幅度频率响应，系统输出信号y(t)的幅度频谱，系统的单位冲激响应h(t)，并按照下面的图Q3-2的布局，绘制出各个信号的时域和频域图形。
图Q3-2你编写的程序Q3_2抄写如下：% Program Q3_2% This Program is used to compute and draw the plots of the frequency response % of a continuous-time systemclc,clear,close allT = 0.01; dw = 0.1;t = 0:T:40;w=-4*pi:dw:4*pi;b = input('Type in the right coefficient vector of differential equation b£º'); % The coefficient vector of the right side of the differential equationa = input('Type in the left coefficient vector of differential equation a£º'); % The coefficient vector of the left side of the differential equationx = input('Type in the expression of the input signal x(t)£º');X=x*exp(-j*t'*w)*T; %CTFTX1=abs(X);y=lsim(b,a,x,t);Y=y'*exp(-j*t'*w)*T;Y1=abs(Y);h=impulse(b,a,40);[H,w] = freqs(b,a); % Compute the frequency response HHm = abs(H); % Compute the magnitude response Hmsubplot(324)plot(w,Hm), axis([0 4*pi -0 1.2]);grid on, title('Magnitude response of the system')subplot(321)plot(t,x),axis([0 40 -3 3]); grid on,title('Input sihnal x(t)')subplot(323)impulse(b,a,40), axis([0 40 -0.2 1]);grid on, title('Impluse response h(t) of the system')subplot(325)lsim(b,a,x,t),axis([0 40 -1 1]); grid on, title('Output signal y(t)')w=-4*pi:dw:4*pi;subplot(322)plot(w,X1),axis([-4*pi 4*pi 0 20]); grid on, title('Magnitude response of input signal x(t)')subplot(326)plot(w,Y1),axis([-4*pi 4*pi 0 20]) ;grid on, title('Magnitude response of output signal y(t)'); xlabel('Frequency in rad/sec')%Eq.3.3 b=[262];a=[1 10 48 148 306 401 262];% x(t)=sin(t)+sin(8*t);% , xlabel('Time t ')执行程序Q3_2，输入信号x(t) = sin(t) + sin(8t)，输入由Eq.3.3描述的系统。