信号与系统实验题目.doc

实验一 连续时间系统的时域和频域分析 相关MATLAB函数 设描述连续时间系统的微分方程为： 则可用向量和表示该系统，即 注意，向量和的元素一定要以微分方程时间求导的降幂次序排列，且缺项要用0补齐。如微分方程 表示该系统的向量为 （1）求解冲激响应：impulse()函数 impulse()函数有以下四种调用格式： ① impulse(b,a) 该调用格式以默认方式绘制由向量和定义的连续时间系统的冲激响应的时域波形。 ② impulse(b,a,t) 该调用格式绘制由向量和定义的连续时间系统在时间范围内的冲激响应的时域波形。 ③ impulse(b,a, t1:p:t2) 该调用格式绘制由向量和定义的连续时间系统在时间范围内，且以时间间隔均匀抽样的冲激响应的时域波形。 ④ y=impulse(b,a,t1:p:t2) 该调用格式并不绘制系统冲激响应的波形，而是求出由向量和定义的连续时间系统在时间范围内以时间间隔均匀抽样的系统冲激响应的数值解。 （2）求解阶跃响应：step()函数 step()函数也有四种调用格式： ① step(b,a) ② step(b,a,t) ③ step(b,a, t1:p:t2) ④ y=step(b,a,t1:p:t2) 上述调用格式的功能与impulse()函数完全相同。 （3）求解零状态响应：lsim()函数 lsim()函数有以下二种调用格式： ① lsim(b,a,x,t) 在该调用格式中，和是描述系统的两个行向量，表示输入信号时间范围的向量，则是输入信号在向量定义的时间点上的取样值。 ② y=lsim(b,a,x,t) 该调用格式并不绘制系统的零状态响应曲线，而是求出与向量定义的时间间隔相一致的系统零状态响应的数值解。 2. 设系统的频率响应为 求解频率响应：freqs()函数 （1）h=freqs(b,a,w) 该调用格式中，w为形如w1:p:w2的冒号运算定义的系统频率响应的频率范围，w1为起始频率，w2为终止频率，p为频率取样间隔。向量h返回在向量w所定义的频率点上系统频响的样值。 （2）[h,w]=freqs(b,a) 该调用格式将计算默认频率范围内200个频率点的系统频率响应的样值，并赋值给返回变量h，200个频率点记录在w中。 （3）[h,w]=freqs(b,a,n) 该调用格式将计算默认频率范围内n个频率点的系统频率响应的样值，并赋值给返回变量h，n个频率点记录在w中。 （4）freqs(b,a) 该调用格式并不返回系统频率响应的样值，而是以对数坐标的方式绘制系统的幅频响应和相频响应。 实验内容： 已知描述系统的微分方程和激励信号如下 ， 要求：（1）从理论上求解系统的冲激响应和零状态响应，并根据求解结果用MATLAB绘制其时域波形； （2）分别用MATLAB的impulse()函数和lsim()函数绘制系统的冲激响应和零状态响应，验证（1）中的结果。 下图是用RLC元件构成的二阶低通滤波器。设，，，试用MATLAB的freqs()函数绘出该系统的幅频特性曲线和相频特性曲线。（求模：abs()函数，求相角：angle()函数） 已知某二阶系统的零极点分别为，，（二重零点），试用MATLAB绘出该系统在0 ~ 1kHz频率范围内的幅频特性曲线（要求用频率响应的几何求解方法实现），并说明该系统的作用。 实验二 离散时间系统的时域和频域分析 相关MATLAB函数 1. 求解零状态响应：filter()函数 y=filter(b,a,x) 由向量b和a组成的系统对输入x进行滤波 2. 求解单位序列响应：impz()函数 （1）h=impz(b,a) 计算单位序列响应的序列值，取样点个数由MATLAB自动选取 （2）h=impz(b,a,n) 计算指定范围内（0: n-1）的单位序列响应的序列值 （3）impz(b,a) 绘制单位序列响应的时域波形 3. 求解频率响应：freqz()函数 （1）[h,w]=freqz(b,a,n) 可得到n点频率响应，这n个点均匀地分布在上半单位圆（即），并将这n点频率记录在w中，相应的频率响应记录在h中。n最好能取2的幂次方，如果缺省，则n=512。 （2）[h,w]=freqz(b,a,n,whole) 在之间均匀选取n个点计算频率响应。 （3）[h,w]=freqz(b,a,n,Fs) Fs为采样频率（以Hz为单位），在0~Fs/2频率范围内选取n个频率点，计算相应的频率响应。 （4）[h,w]=freqz(b,a,n,whole,Fs) 在0~Fs之间均匀选取n个点计算频率响应