离散时间系统的时域特性分析.pdfVIP

数字信号处理实验报告 学生姓名：孙奇 学生学号 学生班级 所属专业：通信工程 实验日期：2012-11-6 实验一：离散时间系统的时域特性分析 实验目的 线性时不变离散时间系统在时域中可以通过常系数线性差 分方程来描述 冲激响应序列可以刻画其时域特性。本实验 通过使用MATLAB 函数研究离散时间系统的时域特性 以加 深对离散时间系统的差分方程、冲激响应系统的线性和时不 特性的理解。 基本原理 一个离散时间系统是将输入序列变换成输出序列的一种运 算。离散时间系统 离散时间系统最重要的最常用的是“线 性时不变系统 实验内容 程序一 clf; n=0:100; x=cos(20*pi*n/256)+cos(200*pi*n/256); subplot(3,1,1); stem(n,x); %输入信号的图形 xlabel(时间信号n); ylabel(信号幅度); title(输入信号); den1=[1]; %对应系统一的差分方程系数 num1=[0.5 0.27 0.77]; den2=[1 -0.53 0.46]; %对应系统二的差分方程系数 num2=[0.45 0.5 0.45]; y1=filter(num1,den1,x); subplot(3,1,2); stem(n,y1); %系统一输出信号的图形 y2=filter(num2,den2,x); subplot(3,1,3); stem(n,y2); %系统二输出信号的图形 3 程序二 n=40; %取冲击响应的前40 个样本 num1=[0.5 0.27 0.77]; %对应系统一的差分方程系数 den1=[1]; num2=[0.45 0.5 0.45]; %对应系统二的差分方程系数 den2=[1 -0.53 0.46]; y1=impz(num1,den1,n); %系统一的冲击响应 subplot(2,1,1); stem(y1); y2=impz(num2,den2,n); %系统二的冲击响应 subplot(2,1,2); stem(y2); 判断是否为线性 程序三（1） n=0:40; a=2; %任取两个系数 b=3; x1=cos(2*pi*0.3*n); x2=cos(2*pi*0.5*n); x=a*x1+b*x2; num=[0.45 0.5 0.45]; %对应系统二的差分方程系数 den=[1 -0.53 0.46]; y1=filter(num,den,x1); %计算出y1 （n） y2=filter(num,den,x2); %计算出y2 （n） y=filter(num,den,x); %计算出y （n） subplot(2,1,1); stem(n,y); ylabel(信号幅度); yt=a*y1+b*y2; %计算出yt （n）=a y1 （n）+b y2 （n） subplot(2,1,2); stem(n,yt); ylabel(信号幅度); 从图中可知，上下两个图完全一