信号与系统实验二
信号与系统
《信号与系统》仿真作业实验一：连续信号的表示及可视化：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)波形。
信号与系统仿真作业
nGDOU-B—11—112广东海洋大学学生实验报告书（学生用表）课程名称课程号学院（系)信息学院专业班级学生姓名学号实验地点04002 实验日期实验一连时间信号的MATLAB表示和连续时间LTI系统的时域分析一、实验目的1．掌握MA TLAB产生常用连续时间信号的编程方法，并熟悉常用连续时间信号的波形和特性；2．运用MATLAB符号求解连续系统的零输入响应和零状态响应；3．运用MATLAB数值求解连续系统的零状态响应;4．运用MATLAB求解连续系统的冲激响应和阶跃响应；5．运用MATLAB卷积积分法求解系统的零状态响应。
二、实验原理1. 连续信号MATLAB实现原理从严格意义上讲,MA TLAB数值计算的方法并不能处理连续时间信号.然而，可用连续信号在等时间间隔点的取样值来近似表示连续信号,即当取样时间间隔足够小时，这些离散样值能够被MATLAB处理，并且能较好地近似表示连续信号.MATLAB提供了大量生成基本信号的函数.比如常用的指数信号、正余弦信号等都是MATLAB的内部函数。
为了表示连续时间信号，需定义某一时间或自变量的范围和取样时间间隔，然后调用该函数计算这些点的函数值，最后画出其波形图.三、实验内容1．实例分析与验证根据以上典型信号的MA TLAB函数，分析与验证下列典型信号MA TLAB程序，并实现各信号波形图的显示，连续信号的图形显示使用连续二维图函数plot（).（1）正弦信号：用MA TLAB命令产生正弦信号2sin(2/4)ππ+，并会出时间0≤t≤3的波形图。
程序如下：K=2;w=2*pi ;phi=pi/4;t=0：0.01：3；ft=K*sin （w*t+phi ）；plot(t,ft ）,grid on ；axis （[0,3,-2。
2,2.2］)title （’正弦信号’)（2） 抽样信号：用MA TLAB 中的sinc(t)函数命令产生抽样信号Sa(t),并会出时间为66t ππ-≤≤的波形图。
MATLAB信号与系统实验报告19472[五篇范文]
![MATLAB信号与系统实验报告19472[五篇范文]](https://img.taocdn.com/s3/m/a72999dee109581b6bd97f19227916888486b9f2.png)
MATLAB信号与系统实验报告19472[五篇范文]第一篇：MATLAB信号与系统实验报告19472信号与系统实验陈诉（5）MATLAB 综合实验项目二连续系统的频域阐发目的：周期信号输入连续系统的响应可用傅里叶级数阐发。
由于盘算历程啰嗦，最适适用MATLAB 盘算。
通过编程实现对输入信号、输出信号的频谱和时域响应的盘算，认识盘算机在系统阐发中的作用。
任务：线性连续系统的系统函数为11)(+=ωωjj H，输入信号为周期矩形波如图 1 所示，用MATLAB 阐发系统的输入频谱、输出频谱以及系统的时域响应。
-3-2-1 0 1 2 300.511.52Time(sec)图 1要领：1、确定周期信号 f(t)的频谱nF&。
基波频率Ω。
2、确定系统函数 )(Ω jn H。
3、盘算输出信号的频谱n nF jn H Y&&)(Ω=4、系统的时域响应∑∞-∞=Ω=nt jnn eY t y&)（MATLAB 盘算为y=Y_n*exp(j*w0*n“*t);要求（画出 3 幅图）：1、在一幅图中画输入信号f(t)和输入信号幅度频谱|F(jω)|。
用两个子图画出。
2、画出系统函数的幅度频谱|H(jω)|。
3、在一幅图中画输出信号y(t)和输出信号幅度频谱|Y(jω)|。
用两个子图画出。
解:(1)阐发盘算：输入信号的频谱为(n)输入信号最小周期为=2，脉冲宽度，基波频率Ω=2π/ =π，所以(n)系统函数为因此输出信号的频谱为系统响应为(2)步伐：t=linspace(-3,3,300);tau_T=1/4;%n0=-20;n1=20;n=n0:n1;%盘算谐波次数20F_n=tau_T*Sa(tau_T*pi*n);f=2*(rectpuls(t+1.75,0.5)+rectpuls(t-0.25,0.5)+rectpuls(t-2.25,0.5));figure(1),subplot(2,1,1),line(t,f,”linewidth“,2);%输入信号的波形 axis([-3,3,-0.1,2.1]);grid onxlabel(”Time(sec)“,”fontsize“,8),title(”输入信号“,”fontweight“,”bold“)%设定字体巨细，文本字符的粗细text(-0.4,0.8,”f(t)“)subplot(2,1,2),stem(n,abs(F_n),”.“);%输入信号的幅度频谱xlabel(”n“,”fontsize“,8),title(”输入信号的幅度频谱“,”fontweight“,”bold“)text(-4.0,0.2,”|Fn|“)H_n=1./(i*n*pi+1);figure(2),stem(n,abs(H_n),”.“);%系统函数的幅度频谱xlabel(”n“,”fontsize“,8),title(”系统函数的幅度频谱“,”fontweight“,”bold“)text(-2.5,0.5,”|Hn|“)Y_n=H_n.*F_n;y=Y_n*exp(i*pi*n”*t);figure(3),subplot(2,1,1),line(t,y,“linewidth”,2);%输出信号的波形 axis([-3,3,0,0.5]);grid onxlabel(“Time(sec)”,“fontsize”,8),title(“输出信号”,“fontweight”,“bold”)text(-0.4,0.3,“y(t)”)subplot(2,1,2),stem(n,abs(Y_n),“.”);%输出信号的幅度频谱xlabel(“n”,“fontsize”,8),title(“输出信号的幅度频谱”,“fontweight”,“bold”)text(-4.0,0.2,“|Yn|”)(3)波形：-3-2-1 0 1 2 300.511.52Time(sec)输入信号f(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输入信号的幅度频谱|Fn|-20-15-10-5 0 5 10 15 2000.10.20.30.40.50.60.70.80.91n系统函数的幅度频谱|Hn|-3-2-1 0 1 2 300.10.20.30.4Time(sec)输出信号y(t)-20-15-10-5 0 5 10 15 2000.10.20.30.4n输出信号的幅度频谱|Yn| 项目三连续系统的复频域阐发目的：周期信号输入连续系统的响应也可用拉氏变更阐发。
电路、信号与系统(2)实验指导书
描述线性时不变离散系统的差分方程为
编写求解上述方程的通用程序。
[建模]
将方程变形可得（用MATLAB语言表示）
a(1)*y(n)= b(1)*u(n)+…+ b(nb)*u(n-nb+1)- a(2)*y(n-1)-…- a(na)*y(n-na+1)
令us== [u(n),…, u(n-nb+1)]; ys=[y(n-1),…, y(n-na+1)]
x(n)={2,1,-1,3,1,4,3,7}（其中加下划线的元素为第0个采样点）在MATLAB中表示为：
n=[-3,-2,-1,0,1,2,3,4]; x=[2,1,-1,3,1,4,3,7]；
当不需要采样位置信息或这个信息是多余的时候，可以只用x向量来表示。
（一）离散信号的MATLAB表述
[问题]
实验一连续时间信号与系统分析
一、实验目的
1、了解连续时间信号的特点；
2、掌握连续时间信号的MATLAB描述；
3、掌握连续LTI系统单位冲激响应的求解方法；
4、掌握连续LTI系统的零状态响应的求解方法。
二、实验内容
严格说来，只有用符号推理的方法才能分析连续系统，用数值方法是不能表示连续信号的，因为它给出的是各个样点的数据。只有当样本点取得很密时才可看成连续信号。所谓很密，是相对于信号变化的快慢而言的。以下均假定相对于采样点密度而言，信号变化足够慢。
elseif lu<lh nh=0; nu=lh-lu;
else nu=0; nh=0;
end
dt=0.1;
lt=lmax;
u=[zeros(1, lt), uls, zeros(1, nu), zeros(1, lt)];
《信号与系统》课程实验报告
《信号与系统》课程实验报告《信号与系统》课程实验报告一图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.掌握时域内信号的四则运算基本方法；2.掌握时域内信号的平移、反转、倒相、尺度变换等基本变换；3.注意连续信号与离散信号在尺度变换运算上区别。
二、 实验原理信号的时域基本运算包括信号的相加（减）和相乘（除）。
信号的时域基本变换包括信号的平移（移位）、反转、倒相以及尺度变换。
(1) 相加（减）: ()()()t x t x t x 21±= [][][]n x n x n x 21±= (2) 相乘: ()()()t x t x t x 21∙= [][][]n x n x n x 21∙=(3) 平移（移位）: ()()0t t x t x -→ 00>t 时右移，00<t 时左移[][]N n x n x -→ 0>N 时右移，0<N 时左移(4) 反转：()()t x t x -→ [][]n x n x -→ (5) 倒相：()()t x t x -→ [][]n x n x -→ (6) 尺度变换: ()()at x t x →1>a 时尺度压缩，1<a 时尺度拉伸，0<a 时还包含反转[][]mn x n x → m 取整数1>m 时只保留m 整数倍位置处的样值，1<m 时相邻两个样值间插入1-m 个0，0<m 时还包含反转三、实验结果1、连续时间信号时域的基本运算 （1） 相加减X1(t)=t+2 , x2(t)=cos(2*pi*t) , x(t)=x1(t)+x2(t).验证：由理论得x(t)=t+2+cos(2*pi*t)，而上图x（t）满足该表达式，故得证。
（2）相乘X1(t)=t+2 , x2(t)=cos(2*pi*t) , x(t)=x1(t)*x2(t).验证：由理论得x(t)=（t+2）*cos(2*pi*t)，而上图x（t）满足该表达式，故得证。
信号与系统实验报告
信号与系统实验报告
实验名称：信号与系统实验
一、实验目的：
1.了解信号与系统的基本概念
2.掌握信号的时域和频域表示方法
3.熟悉常见信号的特性及其对系统的影响
二、实验内容：
1.利用函数发生器产生不同频率的正弦信号，并通过示波器观察其时域和频域表示。
2.通过软件工具绘制不同信号的时域和频域图像。
3.利用滤波器对正弦信号进行滤波操作，并通过示波器观察滤波前后信号的变化。
三、实验结果分析：
1.通过实验仪器观察正弦信号的时域表示，可以看出信号的振幅、频率和相位信息。
2.通过实验仪器观察正弦信号的频域表示，可以看出信号的频率成分和幅度。
3.利用软件工具绘制信号的时域和频域图像，可以更直观地分析信号的特性。
4.经过滤波器处理的信号，可以通过示波器观察到滤波前后的信号波形和频谱的差异。
四、实验总结：
通过本次实验，我对信号与系统的概念有了更深入的理解，掌
握了信号的时域和频域表示方法。
通过观察实验仪器和绘制图像，我能够分析信号的特性及其对系统的影响。
此外，通过滤波器的处理，我也了解了滤波对信号的影响。
通过实验，我对信号与系统的理论知识有了更加直观的了解和应用。
信号实验报告 2
信号与系统实验报告实验一信号与系统的时域分析一、实验目的1、熟悉和掌握常用的用于信号与系统时域仿真分析的MA TLAB函数；2、学会用MA TLAB进行信号基本运算的方法；3、掌握连续时间和离散时间信号的MA TLAB产生，掌握用周期延拓的方法将一个非周期信号进行周期信号延拓形成一个周期信号的MATLAB编程。
二、实验内容Q1-1：修改程序Program1_1，将dt改为0.2，再执行该程序，保存图形，看看所得图形的效果如何？dt = 0.01时的程序clear, % Clear all variablesclose all, % Close all figure windowsdt = 0.01; % Specify the step of time variablet = -2:dt:2; % Specify the interval of timex = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.2时的程序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 = sin(2*pi*t); % Generate the signalplot(t,x) % Open a figure window and draw the plot of x(t)title('Sinusoidal signal x(t)')xlabel('Time t (sec)')dt = 0.01时的信号波形dt = 0.2时的信号波形这两幅图形有什么区别，哪一幅图形看起来与实际信号波形更像？答：dt = 0.01的图形比dt = 0.2的图形光滑，dt = 0.01看起来与实际信号波形更像。
信号系统MATLAB实验报告
信号与系统实验报告桂林理工大学信息科学与工程学院 电子信息工程实验二 信号及其表示【实验目的】了解各种常用信号的表达方式掌握部分绘图函数【实验内容】一、绘出连续时间信号x(t)=t e 707.0 sin 32t 关于t 的曲线，t 的范围为 0~30s ，并以递增。
MATLAB 源程序为：t=0::30; %对时间变量赋值x=exp*t).*sin(2/3.*t); %计算变量所对应得函数值 plot(t,x);grid; %绘制函数曲线ylabel('x(t)');xlabel('Time(sec)')二、产生周期为的方波。
MATLAB源程序为：Fs=100000;t=0:1/Fs:1;x1=square(2*pi*50*t,20);x2=square(2*pi*50*t,80);subplot(2,1,1),plot(t,x1),axis([0,,,]); subplot(2,1,2),plot(t,x2),axis([0,,,]);三、产生sinc(x)函数波形。
MATLAB源程序为：x=linspace(-4,4);y=sinc(x);plot(x,y)四、绘制离散时间信号的棒状图。
其中x(-1)=-1,x(0)=1,x(1)=2,x(2)=1,x(3)=0,x(4)=-1,其他时间x(n)=0。
MATLAB源程序为：n=-3:5; %定位时间变量x=[0,0,-1,1,2,1,-1,0,0];stem(n,x);grid; %绘制棒状图line([-3,5],[0,0]); %画X轴线xlabel('n');ylabel('x[n]')五、单位脉冲序列δ（n-0n ）={00...1...0n n n n =≠直接实现：x=zeros(1,N);x(1,n0)=1;函数实现：利用单位脉冲序列)(0n n -δ的生成函数impseq,即 function[x,n]=impseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)==0];plot(n,x);stem(n,x);输入参数：impseq(0,0,9)——连续图形012345678900.10.20.30.40.50.60.70.80.91输入参数：impseq(0,0,9)——离散图形六、单位阶跃序列ε（n-0n ）={00...1...0n n n n ≥<直接实现：n=[ns:nf];x=[(n-n0)>=0];函数实现：利用单位阶跃序列)(0n n -ε的生成函数stepseq ，即 Function[x,n]=stepseq(n0,ns,nf)n=[ns:nf];x=[(n-n0)>=0];plot(n,x);七、实指数序列=,∀)(x n∈,Ranna直接实现：n=[ns:nf]:x=a.^n;函数实现：利用实指数序列n a(的生成函数rexpseq,即n)x=Function[x,n]=rexpseq(a,ns,nf)n=[ns:nf];x=a,^n:八、复指数序列n e n x n j ∀=+,)()(ωδ直接实现：n=[ns:nf];x=exp((sigema+jw)*n);函数实现：利用复指数序列n j e n x )()(ωδ+=的生成函数cexpseq,即 Function[x,n]=cexpseq(sigema,w,ns,nf)n=[ns:nf];x=exp((sigema+j*w)*n);0123456789-3000-2000-1000100020003000400050006000九、正（余）弦序列n wn n x ∀+=),cos()(θ直接实现：n=[ns:nf];x=cos(w*n+sita);函数实现：利用正（余）弦序列x(n)=cos(wn+θ)的生成函数cosswq,即Function[x,n]=cosseq(w,ns,nf,sita)n=[ns:nf];x=cos(w*n+sita);输入参数：cosseq,0,9,30)——连续信号0123456789-0.2-0.15-0.1-0.0500.050.10.150.2输入参数：cosseq,0,9,30)——离散信号0123456789实验三信号的运算【实验目的】了解信号处理的基本操作。
信号与系统实验报告
信号与系统实验报告目录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．连续时间信号的时域基本运算两实验之一实验分析：输出信号值就等于两输入信号相加（乘）。
由于b=2，故平移量为2时，实际是右移1，符合平移性质。
两实验之二心得体会：时域中的基本运算具有连续性，当输入信号为连续时，输出信号也为连续。
平移，伸缩变化都会导致输出结果相对应的平移伸缩。
2.离散时间信号的时域基本运算两实验之一实验分析：输出信号的值是对应输入信号在每个n值所对应的运算值，当进行拉伸变化后，n值数量不会变，但范围会拉伸所输入的拉伸系数。
两实验之二心得体会：离散时间信号可以看做对连续时间信号的采样，而得到的输出信号值，也可以看成是连续信号所得之后的采样值。
二、连续信号卷积与系统的时域分析1.连续信号卷积积分两实验之一实验分析：当两相互卷积函数为冲激函数时，所卷积得到的也是一个冲激函数，且该函数的冲激t值为函数x，函数y冲激t值之和。
两实验之二心得体会：连续卷积函数每个t值所对应的卷积和可以看成其中一个在k值取得的函数与另外一个函数相乘得到的一个分量函数，并一直移动k值直至最后，最后累和出来的最终函数便是所得到的卷积函数。
3.RC电路时域积分两实验之一实验分析：全响应结果正好等于零状态响应与零输入响应之和。
两实验之二心得体会：具体学习了零状态，零输入，全响应过程的状态及变化，与之前所学的电路知识联系在一起了。
三、离散信号卷积与系统的时域分析1.离散信号卷积求和两实验之一实验分析：输出结果的n值是输入结果的k号与另一个n-k的累和两实验之二心得体会：直观地观察到卷积和的产生，可以看成连续卷积的采样形式，从这个方面去想，更能深入地理解卷积以及采样的知识。
2.离散差分方程求解两实验之一实验分析：其零状态响应序列为0 0 4 5 7.5，零输入响应序列为2 4 5 5.5 5.75，全状态响应序列为2 4 9 10.5 13.25，即全状态=零输入+零状态。
两实验之二心得体会：求差分方程时，可以根据全状态响应是由零输入输入以及零状态相加所得，分开来求，同时也加深了自己对差分方程的求解问题的理解。
实验2 信号卷积实验
实验二 信号卷积实验一、实验目的1. 理解卷积的概念及物理意义；2. 通过实验的方法加深对卷积运算的图解方法及结果的理解。
二、实验设备1. 信号与系统实验箱 1台2. 双踪示波器1台3. 铆孔连接线 若干二、实验原理说明卷积积分的物理意义是将信号分解为冲激信号之和，借助系统的冲激响应，求解系统对任意激励信号的零状态响应。
设系统的激励信号为)t (x ，冲激响应为)t (h ，则系统的零状态响应为)(*)()(t h t x t y =()()x t h t d ττ∞-∞=-⎰。
对于任意两个信号)t (f 1和)t (f 2，两者做卷积运算定义为：()()()12f t f t f t d ττ∞-∞=-⎰=)t (f 1*)t (f 2=)t (f 2*)t (f 1。
1. 两个矩形脉冲信号的卷积过程两信号)t (x 与)t (h 都为矩形脉冲信号，如图10-1所示。
下面由图解的方法（图10-1）给出两个信号的卷积过程和结果，以便与实验结果进行比较。
0≤<∞-t210≤≤t 1≤≤t 41≤≤t ∞<≤t 2124τ（b)(a)(c)(d)(e)(f)(g)(h)(i)2卷积结果图10-1 两矩形脉冲的卷积积分的运算过程与结果2. 矩形脉冲信号与锯齿波信号的卷积信号)t(f1为矩形脉冲信号，)t(f2为锯齿波信号，如图10-2所示。
根据卷积积分的运算方法得到)t(f1和)t(f2的卷积积分结果)t(f，如图10-2(c)所示。
(a)(b)(c)图10-2 矩形脉冲信号与锯齿脉冲信号的卷积积分的结果3. 本实验进行的卷积运算的实现方法在本实验装置中采用了DSP数字信号处理芯片，因此在处理模拟信号的卷积积分运算时，是先通过A/D转换器把模拟信号转换为数字信号，利用所编写的相应程序控制DSP芯片实现数字信号的卷积运算，再把运算结果通过D/A转换为模拟信号输出。
结果与模拟信号的直接运算结果是一致的。
2010信号与系统实验2
实验二：连续和离散系统的频域分析一：实验目的1：学习傅里叶正变换和逆变换，理解频谱图形的物理含义2：了解连续和离散时间系统的单位脉冲响应3：掌握连续时间系统的频率特性二：实验原理1. 傅里叶正变换和逆变换公式 正变换：()()j t F f t e dt ωω∞--∞=⎰逆变换：1()()2j t f t F e d ωωωπ∞-∞=⎰2. 频域分析t j tj e d d e t e ωωωπωωωπ⎰⎰∞∞-∞∞-E =E =)(21)(21)(将激励信号分解为无穷多个正弦分量的和。
⎰∞∞-H E =ωωωπωd e t r tj zs )()(21)(，R(ω)为)(t r zs 傅里叶变换;πωωd )(E 各频率分量的复数振幅 激励单位冲激响应时的零状态响应→ )(t δ)(t h单位阶跃响应时的零状态响应激励→)(t u )(t g3 各函数说明：(1)impulse 冲激响应函数：[Y ,X,T]=impulse(num,den);)1()2()1()1()2()1()()()(11++++++++==--n a s a s a m b s b s b s A s B s H n n m m num 分子多项式系数； num=[b(1) b(2) … b(n+1)]; den 分母多项式系数； den=[a(1) a(2) … a(n+1)];Y ,X,T 分别表示输出响应，中间状态变量和时间变量； 如：352)(2+++=s s s s H ，等价于)(2)()(3)(5)(t e t e t r t r t r +=++ 定义den=[1 5 3];num=[1 2]; [Y ,X,T]=impulse(num,den);(2)step 阶跃响应函数：[Y,X,T]=step(num,den);num 分子多项式；den 分母多项式 Y ,X,T 分别表示输出响应，中间状态变量和时间变量；如：352)(2+++=s s s s H ，den=[1 5 3];num=[1 2];[Y ,X,T]= step (num,den);（3）impz 数字滤波器的冲激响应 [h,t] = impz(b,a,n) b 分子多项式系数；a 分母多项式系数；n 采样样本h 离散系统冲激响应；t 冲激时间，其中t=[0:n-1]', n=length(t)时间样本数（4）freqs 频域响应 [h,w] = freqs(b,a,f) b,a 定义同上，f 频率点个数 h 频域响应，w 频域变量)1()2()1()1()2()1()()()(11++++++++==--m a s a s a n b s b s b s A s B s H m m n n三．实验内容1 周期信号傅里叶级数 已知连续时间信号()()2/8cos 3/4coscos )(321ππ++++=t A t A wt A t x ，其中321,,A A A 取值如下：（X 为学号的后两位）]10,1[,5.02321∈⎪⎩⎪⎨⎧===X X A X A X A ]20,11[,55321∈⎪⎩⎪⎨⎧+==-=X X A XA X A ⎪⎩⎪⎨⎧=-=-=X A X A X A 32151020,>X 要求画出信号的时域波形和频域波形（幅度谱和相位谱）。
信号与系统的实验报告
信号与系统的实验报告信号与系统的实验报告引言：信号与系统是电子工程、通信工程等领域中的重要基础学科，它研究的是信号的传输、处理和变换过程，以及系统对信号的响应和特性。
在本次实验中，我们将通过实际操作和数据分析，深入了解信号与系统的相关概念和实际应用。
实验一：信号的采集与重构在这个实验中，我们使用了示波器和函数发生器来采集和重构信号。
首先，我们通过函数发生器产生了一个正弦信号，并将其连接到示波器上进行观测。
通过调整函数发生器的频率和幅度，我们可以观察到信号的不同特性，比如频率、振幅和相位等。
然后，我们将示波器上的信号通过数据采集卡进行采集，并使用计算机软件对采集到的数据进行处理和重构。
通过对比原始信号和重构信号，我们可以验证信号的采集和重构过程是否准确。
实验二：信号的时域分析在这个实验中，我们使用了示波器和频谱分析仪来对信号进行时域分析。
首先，我们通过函数发生器产生了一个方波信号，并将其连接到示波器上进行观测。
通过调整函数发生器的频率和占空比，我们可以观察到方波信号的周期和占空比等特性。
然后，我们使用频谱分析仪对方波信号进行频谱分析，得到信号的频谱图。
通过分析频谱图，我们可以了解信号的频率成分和能量分布情况，进而对信号的特性进行深入研究。
实验三：系统的时域响应在这个实验中，我们使用了函数发生器、示波器和滤波器来研究系统的时域响应。
首先，我们通过函数发生器产生了一个正弦信号，并将其连接到滤波器上进行输入。
然后，我们通过示波器观测滤波器的输出信号，并记录下其时域波形。
通过改变滤波器的参数，比如截止频率和增益等，我们可以观察到系统对信号的响应和滤波效果。
通过对比输入信号和输出信号的波形，我们可以分析系统的时域特性和频率响应。
实验四：系统的频域响应在这个实验中，我们使用了函数发生器、示波器和频谱分析仪来研究系统的频域响应。
首先，我们通过函数发生器产生了一个正弦信号，并将其连接到系统中进行输入。
然后，我们通过示波器观测系统的输出信号，并记录下其时域波形。
《信号与系统》实验报告
信号与系统实验报告班级:姓名:信息与通信工程学院实验一 系统的卷积响应实验性质：提高性 实验级别：必做 开课单位：信息与通信工程学院 学 时：2一、实验目的：深刻理解卷积运算，利用离散卷积实现连续卷积运算；深刻理解信号与系统的关系，学习MATLAB 语言实现信号通过系统的仿真方法。
二、实验设备： 计算机，MATLAB 软件 三、实验原理： 1、 离散卷积和： 调用函数：conv （）∑∞-∞=-==i i k f i f f f conv S )()(1)2,1(为离散卷积和，其中，f1(k), f2 (k) 为离散序列，K=…-2, -1, 0 , 1, 2, …。
但是，conv 函数只给出纵轴的序列值的大小，而不能给出卷积的X 轴序号。
为得到该值，进行以下分析：对任意输入：设)(1k f 非零区间n1~n2，长度L1=n2-n1+1；)(2k f 非零区间m1~m2，长度L2=m2-m1+1。
则：)(*)()(21k f k f k s =非零区间从n1+m1开始，长度为L=L1+L2-1，所以S （K ）的非零区间为：n1+m1~ n1+m1+L-1。
2、 连续卷积和离散卷积的关系：计算机本身不能直接处理连续信号，只能由离散信号进行近似： 设一系统（LTI ）输入为)(t P ∆，输出为)(t h ∆，如图所示。
)t)()(t h t P ∆∆→)()(lim )(lim )(0t h t h t P t =→=∆→∆∆→∆δ若输入为f(t):∆∆-∆=≈∑∞-∞=∆∆)()()()(k t P k f t f t f k得输出：∆∆-∆=∑∞-∞=∆∆)()()(k t hk f t y k当0→∆时：⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==ττδτd t f k t P k f t f t f k )()()()(lim)(lim )(0⎰∑∞∞-∞-∞=∆→∆∆→∆-=∆∆-∆==τττd t h f k t hk f t y t y k )()()()(lim)(lim )(0所以：∆∆-∆=-==∑⎰→∆)()(lim)()()(*)()(212121k t f k fd t f f t f t f t s τττ如果只求离散点上的f 值)(n f ∆])[()()()()(2121∑∑∞-∞=∞-∞=∆-∆∆=∆∆-∆∆=∆k k k n f k f k n f k fn f所以，可以用离散卷积和CONV （）求连续卷积，只需∆足够小以及在卷积和的基础上乘以∆。
信号与系统的实验报告（2）
信号与系统实验报告——连续时间系统的复频域分析班级：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 的分子多项式为零的点，极点指使分母多项式为零的点，零点使系统的值为零，极点使系统函数的值无穷大。
信号与系统期末试卷及答案
读书破万卷下笔如有神实验二利用DFT分析离散信号频谱一、实验目的应用离散傅里叶变换（DFT），分析离散信号的频谱。
深刻理解DFT分析离散信号频谱的原理，掌握改善分析过程中产生的误差的方法。
二、实验原理根据信号傅里叶变换建立的时域与频域之间的对应关系，可以得到有限长序列的离散傅里叶变换（DFT）与四种确定信号傅里叶变换之间的关系（见教材），实现由DFT分析其频谱。
三、实验内容?3的频谱；1.利用FFT分析信号x(310),nn?,1,...,n)?cos(8（1）、确定DFT计算的参数；N=32;n=0:N-1;x=cos(3*pi/8*n);X=fft(x,N);subplot(2,1,1);stem(n,abs(fftshift(X)));ylabel('Magnitude');xlabel('Frequency (rad)');title('朱艺星杨婕婕'); subplot(2,1,2);stem(n,angle(fftshift(X)));ylabel('Phase');xlabel('Frequency(rad)');读书破万卷下笔如有神进行理论值与计算值比较，讨论信号频谱分析过程中误差原因及改善2）（方法。
在频谱分析过程中由于取样频率过低或者由于信号的截取长度不当将会答：产生误差。
可以适当提高取样率，增加样点数，可能会产生混频现象,取样频率过低,来减少混叠对频谱分析所造成的误差。
对于连续周期信号，其时域取样必须kfo，即（其中K≥2*N+1N为最高谐波分量）其取样点数满足时域取样定理：2fm+fo。
≥≥2Nfo+fo;fs截取信号长度不当，会产生功率泄露，对周期序列进行频谱分析时，为避免泄露应做到：截取的长度应取一个基本周期或基本周期的整数倍，若待分析的周期信号事先不知道其确切的周期，则可截取较长时间长度的样点进行分析，以减少功率泄露误差。
信号与系统 matlab实验报告
信号与系统 matlab实验报告《信号与系统 Matlab实验报告》摘要：本实验报告通过使用 Matlab 软件进行信号与系统实验，探讨了信号与系统在数字领域的应用。
实验结果表明，Matlab 软件具有强大的信号处理和系统分析功能，能够有效地进行信号与系统的模拟和分析。
引言：信号与系统是电子工程领域中的重要基础课程，它研究了信号的产生、传输和处理，以及系统对信号的响应和影响。
在数字领域，信号与系统的理论和方法也得到了广泛的应用。
Matlab 软件作为一种强大的数学计算工具，为信号与系统的模拟和分析提供了便利和高效的途径。
实验一：信号的生成与显示在本实验中，我们首先使用 Matlab 软件生成了几种常见的信号，包括正弦信号、方波信号和三角波信号。
通过调整信号的频率、幅度和相位等参数，我们观察了信号的变化，并利用 Matlab 的绘图功能将信号图形显示出来。
实验结果表明，Matlab 软件能够方便地生成各种类型的信号，并能够直观地显示信号的波形和特性。
实验二：信号的采样与重构在本实验中，我们使用 Matlab 软件对信号进行了采样和重构。
我们首先对一个连续信号进行了离散采样，然后利用 Matlab 的插值函数对采样信号进行了重构。
实验结果表明，采样和重构过程中存在信号失真和频率混叠等问题，但通过适当的采样和重构方法，我们能够有效地还原原始信号。
实验三：系统的响应与分析在本实验中，我们使用 Matlab 软件对系统的响应进行了分析。
我们构建了几种常见的系统模型，包括线性时不变系统和滤波器系统，然后利用 Matlab 的系统分析工具对系统的频率响应、相位响应和单位脉冲响应等进行了分析。
实验结果表明，Matlab 软件能够有效地进行系统的模拟和分析，为系统设计和优化提供了有力的支持。
结论：通过本实验，我们深入了解了信号与系统在数字领域的应用，并掌握了使用 Matlab 软件进行信号与系统模拟和分析的方法。
信号与系统实验一、二
chapter1实验内容：1、画出以下连续时间信号的波形1-0)f(t)=cos(2πt)代码如下：pi=3.14159;t=0:0.01:8;fa=cos(2*pi*t);plot(t,fa);1-1)f (t)=sin(2πt)代码如下：pi=3.14159;t=0:0.01:8;fa=sin(2*t*pi); plot(t,fa);2-0)f (t)=Sa(t/π) 代码如下：pi=3.14159;t=0:0.01:8;fa=sinc(t/pi); plot(t,fa);3-0)f (t)=2[u(t 3)- u(t 5)] 代码如下：t=-1:0.01:10;ft=2*((t>=3)-(t>=5)); plot(t,ft);axis([-1,10,0,3]);4-1)f (t)=e t 代码如下：t=0:0.01:10; ft=exp(t);plot(t,ft);4-2)f (t)=e-t u(t) 代码如下：t=0:0.01:10;f1=(t>=0);f2=exp(-t);plot(t,f1.*f2);5-0）f(t)=2e j(π/4)t，画出实部、虚部、模和相角的波形代码如下：t=0:0.01:10;ft=2*exp(j*(pi/4)*t);h=real(ft); %实部g=imag(ft); %虚部r=abs(ft); %模a=angle(ft); %相角subplot(2,2,1),plot(t,h),title('实部') subplot(2,2,3),plot(t,g),title('虚部') subplot(2,2,2),plot(t,r),title('模')subplot(2,2,4),plot(t,a),title('相角')7）f (t) = u(t)代码如下：t=-1:0.01:5ft=(t>=0);plot(t,ft);axis([-1,5,0,1.5]);8）f (t) =δ(t)代码如下：t=-1:0.01:5;ft=(t>=0)-(t>=0.1); plot(t,ft);axis([-1,1,0,1.1]);9)f9为周期矩形信号，其幅度从-1 到1，占空比为75% 代码如下：pi=3.14159;t=-10:0.01/pi:10;ft=square(t,75);plot(t,ft);2、信号本身运算画出f1(t)为宽度是4，高为1，斜度为0.5 的三角脉冲，然后画出f1(-t)，f1(2t)，f1(2-2t)的波形以及f1(t)的微分和积分波形。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
实验报告二
一、实验目的
1.学会用MATLAB实现连续时间信号傅里叶变换。
2.学会用MATLAB分析LTI系统的频域特性。
二、实验内容
题目一：验证实验原理中所述的相关程序；
程序1.1：
syms t;
Fw=fourier(exp(-2*abs(t)))
运行结果：
Fw =4/(w^2 + 4)
ylabel('|H(j\omega)|');
title('H(jw)的幅频特性');
subplot(2,1,2);
plot(w,angle(H));
grid on;
xlabel('\omega(rad/s)');
ylabel('\phi(\omega)');
title('H(jw)的相频特性');
Fw=fourier(Gt)
FFP=abs(Fw);
ezplot(FFP,[-10*pi 10*pi]);
grid on;
axis([-10*pi 10*pi 0 2.2]);
运行结果：
Fw =
- (cos(w)*i - sin(w))/w + (cos(w)*i + sin(w))/w
图2
程序1.6：
图片：
图片如图2所示。
图7
w=0:0.025:5;
b=[1];
a=[1,2,2,1];
H=freqs(b,a,w);
subplot(2,1,1);
plot(w,abs(H));
grid on;
xlabel('\omega(rad/s)');
ylabel('|H(j\omega)|');
title('H(jw)的幅频特性');
ylabel('|F(j\omega)|');
grid on;
axis([-15 15 0 1.8]);
图片：
图片如图2所示。
图6
题目三：利用ifourier()函数求下列频谱函数的傅式反变换
（1） （2）
答：程序：
程序如下所示。
syms t w;
Fw1=sym('-j*2*w/(16+w^2) *pi 10*pi 0 1.1])
图片：
图片如图2所示。
图5
（3）高斯信号
答：程序：
程序如下所示。
syms t w;
f3=sym('exp(-t^2)');
F3=fourier(f3);
FFP3=abs(F3);
ezplot(FFP3,[-10*pi 10*pi]);
ft1=ifourier(Fw1,w,t)
syms t w;
Fw2=sym('(j*w)^2+5*j*w-8/((j*w)^2+6*j*w+5)');
ft2=ifourier(Fw2,w,t)
运行结果：
ft1=
(j*(pi*heaviside(-t)*exp(4*t)*i+ (pi*exp(-4*t)*dirac(t)*i)/4- (pi*exp(4*t)*dirac(t)*i)/4-pi*exp(-4*t)*heaviside(t)*i))/pi
subplot(1,2,1);
ezplot(ft,[-0.5,0.5]);
grid on;
subplot(1,2,2);
ezplot(abs(Fw),[-24*pi,24*pi]);
grid on;
运行结果：
图1
程序1.5：
syms t w;
Gt=sym('heaviside(t+1)-heaviside(t-1)');
答：程序：
程序如下所示。
syms t w;
f1=sym('heaviside(2*t)-heaviside(2*t-1)');
F1=fourier(f1);
FFP1=abs(F1);
subplot(1,2,1);
ezplot(FFP1,[-10*pi 10*pi]);
ylabel('|F1(j\omega)|');
图片：
图片如图2所示。
图4
（2）单边指数信号
答：程序：
程序如下所示。
syms t w;
f2=sym('exp(-t)*heaviside(t)');
F2=fourier(f2);
FFP2=abs(F2);
ezplot(FFP2,[-10*pi 10*pi]);
ylabel('|F(j\omega)|');
ft2 =
-(2*pi*j^2*dirac(t, 2) - (2*(pi*exp(t*(2*(1/j^2)^(1/2) - 3/j)*i)*i - pi*exp(-t*(2*(1/j^2)^(1/2) + 3/j)*i)*i + pi*sign(3*imag(1/j) - 2*imag((1/j^2)^(1/2)))*exp(t*(2*(1/j^2)^(1/2) - 3/j)*i)*i - pi*sign(3*imag(1/j) + 2*imag((1/j^2)^(1/2)))*exp(-t*(2*(1/j^2)^(1/2) + 3/j)*i)*i - pi*exp(t*(2*(1/j^2)^(1/2) - 3/j)*i)*heaviside(t)*2*i + pi*exp(-t*(2*(1/j^2)^(1/2) + 3/j)*i)*heaviside(t)*2*i))/(j^2*(1/j^2)^(1/2)) + pi*j*dirac(t, 1)*10*i)/(2*pi)
subplot(2,1,2);
plot(w,angle(H));
grid on;
xlabel('\omega(rad/s)');
ylabel('\phi(\omega)');
title('H(jw)的相频特性');
运行结果：
图3
题目二：编程实现求下列信号的幅度频谱
（1）求出 的频谱函数 ，请将它与上面门宽为2的门函数 的频谱进行比较，观察两者的特点，说明两者的关系。
grid on;
Gt=sym('heaviside(t+1)-heaviside(t-1)');
Fw=fourier(Gt);
FFP=abs(Fw);
subplot(1,2,2);
ezplot(FFP,[-10*pi 10*pi]);
ylabel('|F(j\omega)|');
grid on;
axis([-10*pi 10*pi 0 2.2]);
运行结果：
ft =
(pi*exp(-t)*heaviside(t) + pi*heaviside(-t)*exp(t))/(2*pi)
程序1.4：
ft=sym('4*cos(2*pi*6*t)*(heaviside(t+1/4)-heaviside(t-1/4))');
Fw=fourier(ft);
题目四：设 ，试用MATLAB画出该系统的幅频特性 和相频特性 。
答：程序：
程序如下所示。
w=0:0.001:20;
b=[1];
a=[0.08,0.4,1];
H=freqs(b,a,w);
subplot(2,1,1);
plot(w,abs(H));
grid on;
xlabel('\omega(rad/s)');
程序1.2：
syms t w;
ft=ifourier(1/(1+w^2),t)
运行结果：
ft =
(pi*exp(-t)*heaviside(t) + pi*heaviside(-t)*exp(t))/(2*pi)
程序1.3：
syms t w;
Fw=sym('1/(1+w^2)');
ft=ifourier(Fw,w,t)