信号与系统实验三
信号与系统分析实验信号的频谱分析
实验三信号的频谱分析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中的周期方波信号的傅里叶级数的系数，并绘制出信号的幅度谱和相位谱的谱线图。
- 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描述的系统。