文档介绍：云南大学信息学院数字信号处理实验报告 1 实验二离散时间系统的时域分析实验室名称: 实验时间: 姓名: 学号: 专业: 指导教师: 成绩教师签名: 年月日一、实验目的 MATLAB 仿真一些简单的离散时间系统,并研究它们的时域特性。 MATLAB 仿真滤波器,研究它们的性质。 仿真研究离散时间系统的性质。二、实验内容 M=2 运行内容一的程序代码,生成输入 x[n]=s1[n]+s2[n] 的输出信号。 若线性时不变系统由 y[n]=(x[n]+x[n-1]) 变成 y[n]=(x[n]-x[n-1]) 对输入 x[n]=s1[n]+s2[n] 的影响是什么? 修改程序 ,用一个长度为 101 、最低频率为 0 、最高频率为 的扫频正弦信号作为输入信号(见程序 ),计算其输出信号。 运行程序 ,对由加权输入得到的 y[n] 与在相同权系数下输出 y1[n] 和 y2[n] 相加得到的 yt[n] 进行比较,这两个序列是否相等?该系统是线性系统吗? 运行程序 并比较输出序列 y[n] 和yd[n-10] 。 运行程序 ,生成式( )所给离散时间系统的冲击响应。 修改程序 产生如下因果线性时不变系统的冲激响应的前 45个样本: y[n]+[n-1]-[n-2]-[n-3]=[n]-[n-1]+[n-2]+ [n-3] 。 运行程序 ,计算输出序列 y[n] 和y2[n] 以及差值信号 d[n] 。y[n] 和 y2[n] 相等吗? ,对序列h[n] 和x[n] 求卷积,生成y[n], 并用FIR滤波器h[n] 对输入 x[n] 滤波,求得 y1[n] 。三、实验器材及软件 1台 软件云南大学信息学院数字信号处理实验报告 2 四、实验原理 ; MATLAB 表示和基本运算; 处理离散时间系统的方法。五、实验步骤 ; MATLAB 编程实现题目要求; 。六、实验记录(数据、图表、波形、程序等) 程序代码: %产生输入信号 n=0:100; s1=cos(2*pi**n); %一个低频正弦 s2=cos(2*pi**n); %一个高频正弦 x=s1+s2; %滑动平均滤波器的实现 M=input('Desired length ofthe filter ='); num =ones(1,M); y=filter(num,1,x)/M; clf; %显示输入和输出信号 subplot(2,2,1); plot(n, s1); axis([0, 100, -2, 2]); xlabel('Time index n'); ylabel('Amplitude'); title('Signal #1'); subplot(2,2,2); plot(n, s2); axis([0, 100, -2, 2]); xlabel('Time index n'); ylabel('Amplitude'); title('Signal #2'); subplot(2,2,3); plot(n, x); axis([0, 100, -2, 2]); xlabel( 'Time index n');ylabel( 'Amplitude' ); title( 'Input Signal' ); subplot(2,2,4); plot(n, y); axis([0, 100, -2, 2]); xlabel( 'Time index n');ylabel( 'Amplitude' ); title( 'Output Signal' );axis; 云南大学信息学院数字信号处理实验报告 3 显示结果: 代码(将输出信号的代码做出修改,其余代码与 相同): n=0:100; s1=cos(2*pi**n); s2=cos(2*pi**n); x=s1+s2; M=input(' Desired length ofthe filter ='); num=(-1).^[0:M-1]; y=filter(num,1,x)/M; clf; %显示输入和输出信号 subplot(2,2,1); plot(n,