《信号及系统》上机实验.pptVIP

(2) 创建连续时间标号向量t=[-3:0.01:3]; (3) 画出N分别为1,3,7,19,79时，有限项近似 的波形。 ω0=2*pi/T; 第8题（吉布斯现象） 命令脚本 MATLAB命令的一个文本文件，其文件名在当前工作目录或在MATLABPATH的其它地方，以.m结尾 键入这一文件名而没有.m，包含在文件中的命令将被执行。 好处 避免重新键入所有命令 使得排除错误更容易 prob1.m 画出离散时间信号cos(πn/4)，并计算在0≤n≤16内的均值 % prob1.m %引导注释行，不予执行 n=[0:16]; x1=cos(pi*n/4); y1=mean(x1); %mean(X)求向量X所有元素的平均 stem(n,x1) title(x1=cos(pi*n/4)) xlabel(n(samples)) ylabel(x1[n]) prob1 for、if语句 FOR variable = expr statement ... statement END IF expression statements ELSEIF expression statements ELSE statements END EXAMPLE n=[-3:3]; x=2*n; 等效于 x=zeros(1,7); for i=1:7 x(i)=2*(i-4); end 函数文件 以函数名为文件名的M文件，其第1个字是function，第1行的其余部分给出这个函数输出和输入的变量名 “黑匣子” 与命令文件的区别： 命令文件的变量在文件执行完后保存在内存中 函数文件内定义的变量仅在函数文件内部起作用，当函数文件执行完后，这些内部变量将被清除 foo.m 实现一个名为foo的函数，接收一个输入x，产生y和z，分别等于2*x和(5/9)*(x-32) function [y,z]=foo(x) % [y,z]=foo(x) accepts a numerical argument x and % return two arguments y and z, where y is 2*x % and z is (5/9)*(x-32) y=2*x; z=(5/9)*(x-32); 注释行作用 解释提示作用 help foo [y,z]=foo(x) accepts a numerical argument x and return two arguments y and z, where y is 2*x and z is (5/9)*(x-32) 调用函数 foo(3) ans = 6 [y,z]=foo(-40) y = -80 z = -40 实验内容 题1－7 题8（选作） 第1题 1. 在 内画出下面每一个信号： 确定每个信号的基波周期是多少？ 第2题（离散系统性质） 2. 离散时间系统往往是用几个性质来表征，如线性、时不变性、稳定性、因果性及可逆性等。MATLAB可用来构成一些反例证明某些性质不满足。 (1) 系统 不是线性的。利用信号 来证明该系统不满足线性性质。 (2) 系统 不是因果的。利用信号 证明它。定义向量x和y分别代表在 上的输入和在 上的输出。 卷积计算 MATLAB可以用conv计算卷积和： 若 区间内非零， 上非零，那么y[n]仅在 内非零 conv不产生y[n]样本的序号。 第3题（卷积计算） 3. 有限长信号 (1) 用解析方法计算 (2) 用conv计算y。步骤： a. 定义0≤n≤5区间上的向量x b. 定义0≤n≤5区间上的向量h c. 用y=conv(x,h)计算y d. 构造y的标号向量ny e. 用stem(ny,y)画出结果 f. 验证此结果与解析导出的结果是否一致？ filter filter：计算由线性常系数差分方程 表征的因果LTI系统的输出 若x是在nx≤n≤nx+Nx-1上的输入向量，而向量a和b包含系数ak和bm，那么y=filter(