数字信号处理实验教材（学生）.docVIP

实验1 常见离散信号产生和实现 一、实验目的 1、加深对常用离散信号的理解； 2、熟悉使用MATLAB在时域中产生一些基本的离散时间信号。 二、实验原理 在MATLAB中可以利用zeros()函数实现。x=[1 zeros(1, n-1)] 示范程序: % Program P1_1 % Generation of a Unit Sample Sequence clf; % Generate a vector from -10 to 20 n = -10:20; % Generate the unit sample sequence u = [zeros(1,10) 1 zeros(1,20)]; % Plot the unit sample sequence stem(n,u); xlabel(Time index n);ylabel(Amplitude); title(Unit Sample Sequence); axis([-10 20 0 1.2]); 如果在时间轴上延迟了k个单位，得到即： 单位阶跃序列 在MATLAB中可以利用ones()函数实现。 实指数序列 MATLAB实现： 复指数序列 MATLAB实现： 随机序列 MATLAB提供了两种随机信号： rand(1,N) 产生[0,1]上均匀分布的随机矢量。 randn(1,N) 产生均值为0，方差为1的高斯随机序列，即白噪声序列。 基本周期波形 方波 MATLAB工具箱函数square可以产生方波; t=0:0.1*pi:6*pi; y=square(t); axis([0 7*pi -1.5 1.5]); plot(t,y); xlabel(‘时间 t’); ylabel(‘幅度y’); 正弦波 在MATLAB中 锯齿波 工具箱函数sawtooth函数可以产生锯齿波 Fs=10000; t=0:1/Fs:1.5; %抽样长度1.5s, 抽样频率为10kHz x=sawtooth(2*pi*50*t); %信号频率为50Hz plot(t, x); axis()[0 0.2 -1 1]; %画出0.2秒的波形 基本非周期波形 工具箱函数chirp能产生一种扫射频率信号，其特点是信号的瞬时频率随时间按照一定规律变化 t=0:1/1000:2 %抽样频率1kHz, 抽样时间2s. x=chirp(t, 0.1, 150) %0时刻为DC信号，1s时频率为150Hz。 specgram(x, 256, 1000, 256, 250); sinc信号 MATLAB实现： t=linspace(-5,5); x=sinc(t); plot(t,x); 序列的操作 信号加 x(n)=x1(n)+x2(n) MATLAB实现： x=x1+x2; 注意：x1和x2序列应该具有相同的长度，位置对应，才能相加。 信号乘 x(n)=x1(n)*x2(n) MATLAB实现： x=x1.*x2; %数组乘法 改变比例 y(n)=k*x(n) MATLAB实现： y=k*x; 折叠 y(n)=x(-n); MATLAB实现： y=fliplr(x); 抽样和 MATLAB实现： y=sum(x(n1:n2)); 抽样积 MATLAB实现： y=prod(x(n1:n2)); 信号能量 MATLAB实现： Ex=sum(abs(x).^2); 信号功率 MATLAB实现： Px=sum(abs(x).^2)/ N; 三、预习要求 1、预先阅读（MATLAB基础介绍）；2、、实验内容 编制程序产生上述5种信号，长度可输入确定，函数需要的参数可输入确定，并绘出其图形。绘出信号，当、时、、时的信号实部和虚部图；当时呢？此时信号周期为多少？ 实验2 离散系统的时域分析 一、实验目的 1、熟悉并掌握离散系统的差分方程表示法； 2、加深对冲激响应和卷积分析方法的理解。 二、实验原理 在时域中，离散时间系统对输入信号或者延迟信号进行运算处理，生成具有所需特性的输出信号，具体框图如下: 其输入、输出关系可用以下差分方程描述： 输入信号分解为冲激信号, 记系统单位冲激响应，则系统响应为如下的卷积计算式： 当时，h[n]是有限长度的（），称系统为FIR系统；反之，称系统为IIR系统。 三、预习要求 1、在MATLAB中，熟悉利用函数实现差分方程的仿真；2、在MATLAB中，熟悉用函数 计算卷积，用求系统冲