信号与系统Word格式文档下载.docx

1、修改Program1_1后得到的程序Q1_2如下：% Program1_1% This program is used to generate a sinusoidal signal and draw its plotclear, % Clear all variablesclose all, % Close all figure windowsdt = 0.2; % Specify the step of time variablet = -2:dt:2; % Specify the interval of time%x = sin（2*pi*t）; % Generate the signa

2、lx = input（Type in signal x（t） in closed form:）plot（t,x） % Open a figure window and draw the plot of x（t）axis（0,2,0,1）grid on %grid off title（Sinusoidal signal x（t）xlabel（Time（sec） 信号x（t）=e-2t的波形图Q1-3：仿照前面的示例程序的编写方法，编写一个MATLAB程序，以Q1_3为文件名存盘，使之能够在同一个图形窗口中的两个子图中分别绘制信号xn=0.5|n| 和x（t）=cos（2t）u（t）-u（t-3）

3、。要求选择的时间窗能够表现出信号的主要部分（或特征）。编写的程序Q1_3如下：f.delta:function y = delta（t）dt = 0.01;y = （u（t）-u（t-dt）/dt;f.u:% Unit step functionfunction y = u（t）y = （t=0）; % y = 1 for t 0, else y = 0% Program1_5% This program is used to implement the time-shift operation% on a continuous-time signal and to obtain its ti

4、me-shifted versions% and to draw their plots.clear,close all,n = -10:10;x = 0.5.abs（n）; % Generate the signal 0.5.abs（n）?0.5|n|t = -5:（0.01*pi）:5;x1 = cos（2*pi*t）.*u（t）-u（t-3）; % Shift x（t） to the left by 2 second to get x1（t）subplot（211）stem（n,x,.grid on,title （Original signal xnxlabel （Time n subp

5、lot （212）plot （t,x1） Left shifted version of x（t）Time t （sec）信号xn=0.5|n| 的波形图和信号x（t）=cos（2t）u（t）-u（t-3）的波形图Q1-4已知一个序列为编写MATLAB程序Q1_4，能够将xn以N = 8为周期进行周期延拓得到一个周期为N =8的周期序列yn，并分别绘制xn和yn图形。编写的程序Q1_4如下：clear, close all;n1=-3:7;x=n1.*（u（n1）-u（n1-5）;n2=-16:16;N=8;y=0;for k=-2: y=y+（n2-k*N）.*（u（n2-k*N）-u（n2

6、-5-k*N）;endstem（n1,x,Time n （sec）subplot（212）stem（n2,y,Original signal yn信号xn的波形图 信号yn的波形图Q1-5 仿照范例程序Program1_7，编写程序Q1_5，计算并绘制由如下微分方程表示的系统在输入信号为x（t） = （e-2t - e-3t）u（t）时的零状态响应和你手工计算得到的系统零状态响应曲线。手工计算得到的系统零状态响应的数学表达式是：编写的程序Q1_5如下：% Program1_7% This program is used to compute the impulse response h（t）

7、and the step response s（t） of a % continuous-time LTI systemclear, close all,num = input（微分方程右边的系数:）; % 8den = input（微分方程左边的系数: % 1 3 2t = 0:0.01:8;表达式的输入信号x（t）: %input （exp（-2*t）-exp（-3*t）.*u（t）subplot（221）;impulse（num,den,8）;单位冲激响应 x（t）subplot（222）;step（num,den,8）;单位阶跃响应 x（t）subplot（223）;plot（t,x）

8、;零输入响应 x（t）subplot（224）;y=lsim（num,den,x,t）;plot（t,y）;零状态响应 x（t）手工计算得到的系统零状态响应曲线代码：close all, dt=0.01;t=-5:y=（4*exp（-t）-8*t.*exp（-2*t）-4*exp（-3*t）.*u（t）;axis（-2,10,0,0.3）微分方程解的图像 y（t）Time index t用MATLAB绘制的手工计算的系统响应执行程序Q1_5得到的系统响应四、实验结论与体会 在这第一次实验中，由于对实验的不了解，导致做实验的时候不知道该怎么样去做，还是经过了一段时间的熟悉，才逐渐的明白要做什么，

9、但是已经浪费了大量时间，使得原来的任务没有做完，只好推迟做事之前应当把准备工作做好、实验二连续时间信号的频域分析实验名称 信号与系统的频域分析 评分 实验日期 2015 年 12 月 22 日 指导教师 张鏊峰 姓名 易煌轩 专业班级 通信1401 学号 201403020130 1、掌握连续时间周期信号的傅里叶级数的物理意义和分析方法；2、观察截短傅里叶级数而产生的“Gibbs现象”，了解其特点以及产生的原因；3、掌握连续时间傅里叶变换的分析方法及其物理意义；4、掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质；掌握并深刻理傅里叶变换的物理意义，掌握信号的傅里叶变换的计算方

10、法，掌握利用MATLAB编程完成相关的傅里叶变换的计算。 学会利用MATLAB求解信号与系统的相关问题；三、实验内容Q2-1 编写程序Q2_1，绘制下面的信号的波形图：其中，0 = 0.5，要求将一个图形窗口分割成四个子图，分别绘制cos（0t）、cos（30t）、cos（50t） 和x（t） 的波形图，给图形加title，网格线和x坐标标签，并且程序能够接受从键盘输入的和式中的项数。抄写程序Q2_1如下：clear,dt = 0.001;w0 = 0.5*pi;i = input（请输入和式中的项数：x = 0;for r = 1:i; x = x+1/（r*2-1）*sin（2*r-1）*

11、pi/2）*cos（2*r-1）*w0*t）;end;subplot（2,2,4）;grid on;Signal x（t）Time t（sec）y1 = cos（w0*t）;subplot（2,2,1）;plot（t,y1）;Signal y = cos（w0t）y2 = cos（3*w0*t）;subplot（2,2,2）;plot（t,y2）;Signal y = cos（3w0t）y3 = cos（5*w0*t）;subplot（2,2,3）;plot（t,y3）;Signal y = cos（5w0t）执行程序Q2_1所得到的图形如下：在键盘输入：Q2-2 给程序Program2_1增

12、加适当的语句，并以Q2_2存盘，使之能够计算例题2-1中的周期方波信号的傅里叶级数的系数，并绘制出信号的幅度谱和相位谱的谱线图。通过增加适当的语句修改Program2_1而成的程序Q2_2抄写如下：close all;T = 2;dt = 0.00001;x1 = u（t） - u（t-1-dt）;for m = -1:1; x = x + u（t-m*T） - u（t-1-m*T-dt）;w0 = 2*pi/T;N = 10;L = 2*N+1;for k = -N: N; ak（N+1+k） = （1/T）*x1*exp（-j*k*w0*t）*dt;phi = angle（ak）;subp

13、lot（211）;k = -N:stem（k,abs（ak）,幅度谱Frequency index kylabel（|Fn|subplot（212）;stem（k,phi,相位谱n执行程序Q2_2得到的图形Q2-3 反复执行程序Program2_2，每次执行该程序时，输入不同的N值，并观察所合成的周期方波信号。通过观察，你了解的吉伯斯现象的特点是： N=10 N=100答:总是第一个波峰的值不变 且随着n次谐波的增大 从第二个开始波峰、波谷开始接近于一条直线。1、周期信号的傅里叶级数与GIBBS现象给定如下两个周期信号：Q2-4 分别手工计算x1（t） 和x2（t） 的傅里叶级数的系数。信号x

14、1（t） 在其主周期内的数学表达式为：x（t）=（t+1）*（u（t+1） - u（t） （t-1）*（u（t） - u（t-1）计算x1（t） 的傅里叶级数的系数的计算过程如下：= 1/24/通过计算得到的x1（t）的傅里叶级数的系数的数学表达式是：信号x2（t） 在其主周期内的数学表达式为： -0.2 t 0.2计算x2（t） 的傅里叶级数的系数的计算过程如下：用MATLAB帮助你计算出你手工计算的傅里叶级数的系数ak从-10到10共21个系数。从命令窗口上抄写x1（t）的21个系数如下：ak = Columns 1 through 5 -0.5000 - 0.0000i 0.5000 +

15、 0.0000i -0.5000 - 0.0000i 0.5000 + 0.0000i -0.5000 - 0.0000i Columns 6 through 10 0.5000 + 0.0000i -0.5000 - 0.0000i 0.5000 + 0.0000i -0.5000 - 0.0000i 0.5000 + 0.0000i Columns 11 through 15 -0.5000 0.5000 - 0.0000i -0.5000 + 0.0000i 0.5000 - 0.0000i -0.5000 + 0.0000i Columns 16 through 20 0.5000 -

16、 0.0000i -0.5000 + 0.0000i 0.5000 - 0.0000i -0.5000 + 0.0000i 0.5000 - 0.0000i Column 21 -0.5000 + 0.0000i从命令窗口上抄写x2（t）的21个系数如下：Columns 1 through 40.0000 + 0.0000i 0.0000 + 0.0354i 0.0000 - 0.0000i 0.0000 + 0.0455iColumns 5 through 80.0000 - 0.0000i 0.0000 + 0.0637i 0.0000 - 0.0000i 0.0000 + 0.1061i

17、Columns 9 through 120.0000 - 0.0000i 0.0000 + 0.3183i 0.5000 0.0000 + 0.3183iColumns 13 through 160.0000 + 0.0000i 0.0000 + 0.1061i 0.0000 + 0.0000i 0.0000 - 0.0637iColumns 17 through 200.0000 + 0.0000i 0.0000 - 0.0455i 0.0000 + 0.0000i 0.0000 - 0.0354iColumns 210.0000 -0.0000i 在第一次实验的经验下，第二次实验就不像第一

18、次那样，毫无头绪了，一开始实验就直奔重点，经过实验前的预习，在实验课上做就显得没有那么难了，可见事先做准备是没有一点问题的。实验三 连续时间LTI系统的频域分析1、掌握系统频率响应特性的概念及其物理意义；2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用；3、学习和掌握幅度特性、相位特性以及群延时的物理意义；4、掌握用MATLAB语言进行系统频响特性分析的方法。掌握LTI连续和离散时间系统的频域数学模型和频域数学模型的MATLAB描述方法，深刻理LTI系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB计算和绘制LTI

19、系统频率响应特性曲线中的编程。 学会利用MAYLAB求解信号与系统的相关问题；Q3-1 修改程序Program3_1，并以Q3_1存盘，使之能够能够接受键盘方式输入的微分方程系数向量。并利用该程序计算并绘制由微分方程Eq.3.1、Eq.3.2和Eq.3.3描述的系统的幅度响应特性、相位响应特性、频率响应的实部和频率响应的虚部曲线图。抄写程序Q3_1如下：clc,clear all, b = input（微分方程左边的系数向量a = input（微分方程右边的系数向量H,w = freqs（b,a）;Hm = abs（H）;phai = angle（H）;Hr = real（H）;Hi = im

20、ag（H）;plot（w,Hm）;幅度响应特性Frequency in rad/secplot（w,phai）;grid on, 相位响应特性plot（w,Hr）;频率响应的实部plot（w,Hi）;频率响应的虚部执行程序Q3_1，绘制的系统1的频率响应特性曲线如下：从系统1的幅度频率响应曲线看，系统1是低通、高通、全通、带通还是带阻滤波器？带通滤波器执行程序Q3_1，绘制的系统2的频率响应特性曲线如下：从系统2的幅度频率响应曲线看，系统2低通、高通、全通、带通还是带阻滤波器？全通滤波器执行程序Q3_1，绘制的系统3的频率响应特性曲线如下：从系统3的幅度频率响应曲线看，系统3是低通、高通、全通

21、、带通还是带阻滤波器？低通滤波器这三个系统的幅度频率响应、相位频率相应、频率响应的实部以及频率响应的虚部分别具有何种对称关系？请根据傅里叶变换的性质说明为什么会具有这些对称关系？Q3-2 编写程序Q3_2，使之能够能够接受键盘方式输入的输入信号x（t）的数学表达式，系统微分方程的系数向量，计算输入信号的幅度频谱，系统的幅度频率响应，系统输出信号y（t）的幅度频谱，系统的单位冲激响应h（t），并按照下面的图Q3-2的布局，绘制出各个信号的时域和频域图形。图Q3-2你编写的程序Q3_2抄写如下：close allT =0.01;dw =0.1;t=0:T:40;w=-4*pi:dw:4*pi;b=

22、input（b=a=Type in the expression of the input signal x（t）X=x*exp（-j*t*w）*T;X1=abs（X）;y=lsim（b,a,x,t）;Y=y*exp（-j*tY1=abs（Y）;h=impulse（b,a,40）;H,w= freqs（b,a）;Hm= abs（H）;subplot（324）plot（w,Hm）,axis（0 4*pi -0 1.2）;Magnitude response of the systemsubplot（321）plot（t,x）,axis（0 40 -3 3）;Input sihnal x（t）subplot（323）impulse（b,a,40）,axis（0 40 -0.2 1）;Impulse response h（t） of the systemsubplot（325）lsim （b,a,x,t）,axis（0 40 -1 1）;grid on,title（Output signal y（t）subplot（322）plot（w,X1）,axis（-4*pi 4*pi 0 20）;Magnitude response of input signal x（t）subplot（326）plot（w,Y1）,Magnitude respon