2022年信号与系统MATLAB实验报告.pdf

《2022年信号与系统MATLAB实验报告.pdf》由会员分享，可在线阅读，更多相关《2022年信号与系统MATLAB实验报告.pdf（33页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、信号与系统 MATLAB实验报告信号与系统 MATLAB 实验报告院系 : 专业 : 年级 : 班号 : 姓名 : 学号 : 实验时间 : 实验地点 :实验一连续时间信号的表示及可视化实验题目 : )()(ttf;)()(ttf;atetf)(分别取00aa及); )()(tRtf;)()(tSatf;)2()(ftSintf(分别画出不同周期个数的波形)。解题分析 : 以上各类连续函数 ,先运用 t = t1: p:t2 的命令定义时间范围向量,然后调用对应的函数 ,建立 f 与 t 的关系 ,最后调用 plot()函数绘制图像 ,并用 axis()函数限制其坐标范围。实验程序 : (1)(

2、)(ttft=-1:0、01:3 %设定时间变量 t 的范围及步长f=dirac(t) %调用冲激函数 dirac() plot(t,f) %用plot函数绘制连续函数axis(-1,3,-0、5,1 、5) %用axis 函数规定横纵坐标的范围(2)()(ttft=-1:0、01:3 %设定时间变量 t 的范围及步长f=heaviside(t) %调用阶跃函数 heaviside() plot(t,f) %用plot函数绘制连续函数精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 33 页 -

3、 - - - - - - - - - 信号与系统 MATLAB实验报告title(f(t)=heaviside(t) %用title函数设置图形的名称axis(-1,3,-0、5,1 、5) %用axis 函数规定横纵坐标的范围(3)atetf)(a=1 时: t=-5:0、01:5 %设定时间变量 t 的范围及步长f=exp(t) %调用指数函数 exp() plot(t,f) %用plot函数绘制连续函数title(f=exp(t) %用title函数设置图形的名称axis(-5,5,-1,100) %用axis 函数规定横纵坐标的范围a=2 时: t=-5:0、01:5 f=exp(2*

4、t) %调用指数函数 exp() plot(t,f) title(f=exp(2*t) axis(-5,5,-1,100) a=-2 时: t=-5:0、01:5 f=exp(-2*t) plot(t,f) title(f=exp(-2*t) axis(-5,5,-1,100) (4)()(tRtft=-5:0、01:5 f=rectpuls(t,2) %用rectpuls(t,a)表示门函数 , 默认以零点为中心 , 宽度为 a plot(t,f) title(f=R(t) axis(-5 5 -0、5 1 、5) (5)()(tSatf=1 时: 精品资料 - - - 欢迎下载 - - -

5、 - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告t=-20:0、01:20 f=sin(t)、/t %调用正弦函数 sin(),并用sin(t)、/t 实现抽样函数plot(t,f) title(f(t)=Sa(t) axis(-20,-20,-0、5,1 、1) =5 时: t=-20:0、01:20 f=sin(5*t)、/(5*t) plot(t,f) title(f(t)=Sa(5*t) axis(-20,-20,-0、5,1 、1) (6)2(

6、)(ftSintf=1 时: t=-10:0、01:10 f=sin(t) %调用正弦函数 sin() plot(t,f); title(f=sin(t) axis(-10,10,-2,2) =5 时: t=-10:0、01:10 f=sin(5*t) plot(t,f); title(f=sin(5*t) axis(-10,10,-2,2) 实验结果 ; (1) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告

7、-1-0.500.511.522.53-0.500.511.5(2)-1-0.500.511.522.53-0.500.511.5f(t)=heaviside(t)(3) a=1 时: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告-5-4-3-2-10123450102030405060708090100f=exp(t)a=2 时: -5-4-3-2-10123450102030405060708090100

8、f=exp(2*t)a=-2时: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 5 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告-5-4-3-2-10123450102030405060708090100f=exp(-2*t)(4)-5-4-3-2-1012345-0.500.511.5f=R(t)(5) =1 时: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6

9、页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告-20-15-10-505101520-0.4-0.200.20.40.60.81f(t)=Sa(t)=5 时: -20-15-10-505101520-0.4-0.200.20.40.60.81f(t)=Sa(5*t)(6) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告=1 时: -10-8-6-4-20246810-

10、2-1.5-1-0.500.511.52f=sin(t)=5 时: -10-8-6-4-20246810-2-1.5-1-0.500.511.52f=sin(5*t)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告实验心得体会 : (1)在 MATLAB中, 就是用连续信号在等时间间隔点的样值来近似地表示连续信号的, 当取样时间间隔足够小时 , 这些离散的样值就能较好地近似出连续信号。在MATLAB 中 t =

11、t1: p: t2的命令定义时间范围向量 ,t1 为信号起始时间 ,t2 为终止时间 ,p 为时间间隔。(2)plot( )函数可用于连续函数的绘制。(3)用 axis() 函数限制坐标范围 , 可使图像更加匀称美观。改进想法 : 本题中函数的表示方法都不只一种。 如阶跃函数可以借助符号函数来实现可视化。其程序与结果如下: t=-5:0、05:5 f=sign(t) %调用符号函数 sign() axis(-5,5,-1、1,1 、1) ff=1/2+1/2*f %运用阶跃函数与符号函数的关系, 表示出阶跃函数 ff plot(t,ff) axis(-5,5,-0、1,1 、1) -5-4-3

12、-2-101234500.20.40.60.81f=heaviside(t)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告实验二离散时间信号的表示及可视化实验题目 : )()(nnf;)()(nnf;anenf)(分别取00aa及); )()(nRnfN(分别取不同的 N 值);)()(nSanf; )()(nSinnf(分别取不同的值); 解题分析 : 以上各类离散函数 ,可仿照连续函数的可视化 ,先运用 n

13、 =n1: p: n2 的命令定义自变量的范围及步长 ,然后调用对应的函数 ,建立 f 与 t 的关系 ,最后调用 stem()函数绘制图像 ,并用 axis()函数限制其坐标范围。实验程序 : (1)()(nnfn=-5:0、5:5 %设定时间变量 n的范围及步长f=dirac(n) stem(n,f) %调用stem()绘制离散函数title(f=dirac(t) axis(-5,5,-3,10) %用axis 函数规定横纵坐标的范围(2)()(nnfn=-5:0、5:5 f=heaviside(n) stem(n,f) title(f=Heaviside(t) axis(-5,5,-0、

14、5,1 、5) (3)anenf)(a=1 时: n=-5:0、5:5 f=exp(n) stem(n,f) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 10 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告title(f=exp(n) a=2 时: n=-5:0、5:5 f=exp(2*n) stem(n,f) title(f=exp(2*n) a=-2 时: n=-5:0、5:5 f=exp(-2*n) stem(n,f) title(f=exp(-2*n

15、) (4)()(nRnfNn=-5:0、5:5 f=rectpuls(n,2) stem(n,f) title(f=R(n) axis(-5,5,-0、5,1 、5) (5)()(nSanf=1时: n=-20:0、5:20 f=sin(n)、/(n) stem(n,f) title(f=Sa(n) axis(-20,-20,-0、5,1 、1) =5时: n=-20:0、5:20 f=sin(5*n)、/(5*n) stem(n,f) title(f=Sa(5*n) axis(-20,-20,-1,5) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名

16、师归纳 - - - - - - - - - -第 11 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告(6)()(nSinnf=1时: n=-5:0、5:5 f=sin(n) stem(n,f) title(f=sin(n) axis(-5,5,-2,2) =5时: n=-5:0、5:5 f=sin(5*n) stem(n,f) title(f=sin(5*n) axis(-5,5,-2,2) 实验结果 ; (1) -5-4-3-2-1012345-20246810f=dirac(t)(2) 精品资料 - - - 欢迎下载 - - - - - - -

17、 - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 12 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告-5-4-3-2-1012345-0.500.511.5f=Heaviside(t)(3) a=1 时: -5-4-3-2-1012345050100150f=exp(n)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 13 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告a=2 时: -5

18、-4-3-2-101234500.511.522.5x 104f=exp(2*n)a=-2时: -5-4-3-2-101234500.511.522.5x 104f=exp(-2*n)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 14 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告(4) -5-4-3-2-1012345-0.500.511.5f=R(n)(5) =1时: -20-15-10-505101520-0.4-0.200.20.40.60.81f=S

19、a(n)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 15 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告=5时: -20-15-10-505101520-0.2-0.15-0.1-0.0500.050.10.150.20.25f=Sa(5*n)(6) =1时: -5-4-3-2-1012345-2-1.5-1-0.500.511.52f=sin(n)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - -

20、 - - - - - -第 16 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告=5时: -5-4-3-2-1012345-2-1.5-1-0.500.511.52f=sin(5*n)实验心得体会 : 用 plot()函数可以绘制离散序列,但就是与连续序列有所不同,需要在括号内加上、。但就是 plot()画出来的函数图像不直观,显得很凌乱。改进想法 : (1)对于离散函数 ,如果使用 stem(n,f, 、)函数,绘图效果更好。如抽样函数的程序 : n=-20:0、5:20 f=sin(n)、/(n) stem(n,f,、) title(f=Sa(n

21、) axis(-20,-20,-0、5,1 、1) 绘图结果如下 : 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 17 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告-20-15-10-505101520-0.4-0.200.20.40.60.81f=Sa(n)对比可知此法做出的图像更加清晰美观。(2)MATLAB 可以自动地根据曲线数据的范围选择合适的坐标系,从而使得曲线尽可能清晰地显示出来,一般情况下不必选择坐标系。但就是,如果对 MATLAB自动产生的

22、坐标轴不满意,可以利用axis 命令对坐标轴进行调整。实验三系统的时域求解实验题目 : 1、设)10()()(),()9 .0()(nununxnunhn,求)(*)()(nhnxny,并画出)(nx、)(nh、)(ny波形。2、求因果线性移不变系统)2()()2(81.0)(nxnxnyny的单位抽样响应)(nh,并绘出)(jeH的幅频及相频特性曲线。解题分析 : 1、用 heaviside()与 exp() 函数 表示出 x(n) 与 h(n), 然后调用 conv() 函数实现 x(n) 与 h(n) 的卷积 y(n) 。并且分别将三个函数图像绘出。2、通过给矩阵 a,b 赋值, 建立系

23、统差分方程 , 然后调用 impz() 函数求系统的冲激精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 18 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告响应, 再用函数 freqs(b,a)进行系统频率响应的分析。实验程序 : (1) n=-10:20 %设置变量范围 , 默认步长为 1 f=heaviside(n) x=heaviside(n)-heaviside(n-10) %阶跃函数直接相减figure(1) %产生图像窗口 1 stem(n,x) %

24、绘制函数 x title(x(n) h=0 、9、n 、*f %函数h的表达式figure(2) %产生图像窗口 2 stem(n,h) %绘制函数 h title(h(n) n1=-20:40 y=conv(h,x) %调用conv()函数求 h与x的卷积figure(3) %产生图像窗口 3 stem(y) %绘制函数 y title(y(n)=x(n)*h(n) (2) a=1 0 -0、81 %描述系统的差分方程的系数b=1 0 -1 %描述系统的差分方程的系数figure(1) h=impz(n,m,-10:10) %调用impz()函数求系统的冲激响应stem(h) %绘制函数 h

25、的离散序列title(h(n) figure(2) freqs(b,a) %对连续系统频率响应 H(jw)进行分析的函数freqs() 实验结果 ; (1) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 19 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告-10-50510152000.10.20.30.40.50.60.70.80.91x(n)-10-50510152000.10.20.30.40.50.60.70.80.9h(n)精品资料 - - - 欢迎

26、下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 20 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告0102030405060700123456y(n)=x(n)*h(n)(2) 0510152025-0.200.20.40.60.811.2h(n)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 21 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告10-2

27、10-1100101-1-0.500.51Frequency (rad/s)Phase(degrees)10-210-1100101100.01100.05100.09Frequency (rad/s)Magnitude实验心得体会 : (1)计算离散序列的卷积时 ,应考虑其结果的横坐标范围的改变。(2)向量相乘时 ,注意用、 。(3)借助 MATLAB 的内部函数 conv()可以很容易地完成两个信号的卷积运算,并且其完成的就是两个多项式的乘法运算,在 MATLAB中它们的系数构成一个行向量来表示。(3)表示系统的方法就是用系统函数分子与分母多项式系数行向量来表示。改进想法 : (1)n=-

28、10:20 %设置变量范围 , 默认步长为 1 f=heaviside(n) x=heaviside(n)-heaviside(n-10) %阶跃函数直接相减figure(1) %产生图像窗口 1 axis(-10,20,0,1) stem(n,x) %绘制函数 x title(x(n) h=0 、9、n 、*f %函数h的表达式figure(2) %产生图像窗口 2 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 22 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实

29、验报告stem(n,h) %绘制函数 h axis(-10,20,0,1) title(h(n) n1=-20:40 y=conv(h,x) %调用conv 函数求 h与x的卷积figure(3) %产生图像窗口 3 stem(y) %绘制函数 y axis(0,62,0,7) title(y(n)=x(n)*h(n) 运行结果 : -10-50510152000.10.20.30.40.50.60.70.80.91x(n)精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 23 页，共 33 页 - - -

30、 - - - - - - - 信号与系统 MATLAB实验报告-10-50510152000.10.20.30.40.50.60.70.80.91h(n)010203040506001234567y(n)=x(n)*h(n)实验四信号的 DFT 分析实验题目 : 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 24 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告计算余弦序列)()8cos()(nRnnxN的 DFT。 分别对 N=10、 16、 22 时计算 D

31、FT,绘出)(kX幅频特性曲线 ,分析就是否有差别及产生差别的原因。解题分析 : 用矩阵代替门函数给变量n 赋值, 并设定不同的 N值, 然后调用 fft()函数实现函数的傅里叶变换 , 然后用 subplot()与 stem()函数绘图。实验程序 : (1)N=10时: N=10 %设定N的值为 10 n=0:N-1 %用矩阵代替门函数给 n赋值x=cos(pi/8)、*n) %调用cos()函数y=fft(x) %调用fft()函数求 x的傅里叶变换subplot(2,1,1),stem(n,y) %绘制y的离散图title(DFTcos(pi/8)*n) subplot(2,1,2),s

32、tem(n,abs(y) %绘制y的幅频特性曲线title(X(k) (2)N=16时: N=16 %设定N的值为 16 n=0:N-1 %用矩阵代替门函数给 n赋值x=cos(pi/8)、*n) %调用cos()函数y=fft(x) %调用fft()函数求 x的傅里叶变换subplot(2,1,1),stem(n,y) %绘制y的离散图title(DFTcos(pi/8)*n) subplot(2,1,2),stem(n,abs(y) %绘制y的幅频特性曲线title(X(k) (3)N=22时: N=22 %设定N的值为 22 n=0:N-1 %用矩阵代替门函数给 n赋值x=cos(pi/

33、8)、*n) %调用cos()函数精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 25 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告y=fft(x) %调用fft()函数求 x的傅里叶变换subplot(2,1,1),stem(n,y) %绘制y的离散图title(DFTcos(pi/8)*n) subplot(2,1,2),stem(n,abs(y) %绘制y的幅频特性曲线title(X(k) 实验结果; (1)N=10时: 0123456789-10123

34、DFTcos(pi/8)*n01234567890246X(k)(2)N=16时: 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 26 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告051015-50510DFTcos(pi/8)*n05101502468X(k)(3)N=22时: 0510152025-20246DFTcos(pi/8)*n05101520250510X(k)实验结果分析 : 由图可知 ,不同的 N 值所对应的 DFT 序列与幅频响应不同 ,

35、就是因为 N 代表精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 27 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告DFT 的变换区间长度 ,当 N 取不同的值时 ,函数所对应的离散傅里叶变换与幅频特性曲线也不同。实验心得体会 : MATLAB就是计算机运算 ,无法实现无限时间信号与无限大数量的计算,故而周期信号只能取有限个谐波分量近似合成,即 N 值有限 ,且 N 值越大 ,仿真结果越接近。所以手工求取的傅里叶变换系数与MATLAB 求取存在差别。实验五系统

36、时域解的快速卷积求法实验题目 : 用快速卷积法计算系统响应)(*)()(nhnxny,已知: )()4.0sin()(15nRnnx,)(9.0)(20nRnhn。要求取不同的L 点数 ,并画出)(nx、)(nh、)(ny波形,分析就是否有差别及产生差别的原因。解题分析 : 根据离散序列卷积及傅里叶变换的性质,可先求出两函数x(n)与 h(n)的 L 点傅里叶变换 ,分别得到 Xk 与 Yk,然后求 Xk 与 Yk 之积 Hk 的傅里叶反变换 , 即得到了 x(n)与 h(n)的卷积 y(n)。实验程序 : L=10 时: n1=0:14 %用矩阵代替门函数给 n1赋值x=sin(0、4、*n

37、1) %写出x的表达式n2=0:19 %给n2赋值y=0 、9、n2 %写出y的表达式Xk=fft(x,10) %调用fft()函数求 x的L(=10)点傅里叶变换Yk=fft(y,10) %求y的L点傅里叶变换Hk=Xk 、*Yk %写出Hk的表达式h=ifft(Hk) %调用ifft()函数求 Hk的傅里叶反变换subplot(3,1,1),stem(x) %绘制x的离散图title(x(n) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 28 页，共 33 页 - - - - - - - - -

38、- 信号与系统 MATLAB实验报告subplot(3,1,2),stem(y) %绘制y的离散图title(y(n) subplot(3,1,3),stem(h) %绘制h的离散图title(h(n) xlabel(L=10) %横坐标处做标注(2)L=18时: n1=0:14 x=sin(0、4、*n1) n2=0:19 y=0 、9、n2 Xk=fft(x,18) Yk=fft(y,18) Hk=Xk 、*Yk h=ifft(Hk) subplot(3,1,1),stem(x) title(x(n) subplot(3,1,2),stem(y) title(y(n) subplot(3,

39、1,3),stem(h) title(h(n) xlabel(L=18) (3)L=28时: n1=0:14 x=sin(0、4、*n1) n2=0:19 y=0 、9、n2 Xk=fft(x,28) Yk=fft(y,28) Hk=Xk 、*Yk h=ifft(Hk) subplot(3,1,1),stem(x) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 29 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告title(x(n) subplot(3,1,2

40、),stem(y) title(y(n) subplot(3,1,3),stem(h) title(h(n) xlabel(L=28) (4)L=35时: n1=0:14 x=sin(0、4、*n1) n2=0:19 y=0 、9、n2 Xk=fft(x,35) Yk=fft(y,35) Hk=Xk 、*Yk h=ifft(Hk) subplot(3,1,1),stem(x) title(x(n) subplot(3,1,2),stem(y) title(y(n) subplot(3,1,3),stem(h) title(h(n) xlabel(L=35) 实验结果 ; (1)L=10时: 精

41、品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 30 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告051015-101x(n)0246810121416182000.51y(n)12345678910024h(n)L=10(2)L=18时: 051015-101x(n)0246810121416182000.51y(n)024681012141618-505h(n)L=18(3)L=28时: 精品资料 - - - 欢迎下载 - - - - - - - - - -

42、 - 欢迎下载 名师归纳 - - - - - - - - - -第 31 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告051015-101x(n)0246810121416182000.51y(n)051015202530-505h(n)L=28(4)L=35时: 051015-101x(n)0246810121416182000.51y(n)05101520253035-505h(n)L=35实验结果分析 : 由图可知 ,当 L 取不同的值时 ,对应的 y(n)波形形状相似 ,但就是有所不同 ,产精品资料 - - - 欢迎下载 - - - - -

43、 - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 32 页，共 33 页 - - - - - - - - - - 信号与系统 MATLAB实验报告生这种差别的原因就是L 代表傅里叶变换区间长度,当 L 取不同的值时 ,所对应的函数波形也有所差别。实验心得体会 : (1)计算离散序列的卷积 ,虽然本实验的快速卷积方法瞧上去多次变换了变量的域,使过程变复杂了 ,但实际上减少了计算量 ,就是一种快速而简单的方法。(2)用 subplot 绘图函数可将图形窗口分成若干等份,便于将多个图像进行分组或者比较。改进想法 : 当 L 取不同的值时 ,matlab自动生成的图像的横纵坐标范围不同,不便于相互比较,因此可以自己规定坐标轴范围,这样可以更加直观地瞧出各波形间的差别。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 33 页，共 33 页 - - - - - - - - - -