学习用Matlab计算离散信号的功率和能量

实验23 迭代法及离散卷积的计算 实验目的 学习用Matlab计算离散信号的功率和能量。 学习并掌握用迭代法求解差分方程的方法 。 掌握用Matlab进行离散卷积运算的数值方法和解析方法。加深对离散卷积的理解。 实验原理与说明 离散信号的能量与功率 与连续信号类似，离散信号也可分为能量信号和功率信号。对于非周期信号，信号能量定义为 对于周期的离散信号，由于其能量无限大，故常常用功率来作其测量参数。设有一周期为离散信号，其功率定义为： 实验原理与说明 能量有限的信号称为能量信号。功率有限的信号称为功率信号。所有周期信号都是功率信号。 离散序列的求和在MATLAB中可利用sum函数来实现，其调用形式为 y=sum(f(n:m)) 差分方程的迭代解法 实验原理与说明 令上式中，有 以此类推，通过反复迭代，就可以求出任意时刻的响应值。这种迭代方法最适合用计算机计算，下面我们用Matlab来实现这种计算。 实验原理与说明 迭代计算的一般规律 第一项 第二项与上类似 实验原理与说明 离散卷积的计算 Matlab信号处理工具箱提供了一个计算两个离散序列卷积和的函数conv()，其调用格式为 y=conv(f,h) 式中，f、h分别为待卷积的两序列的向量表示，y是卷积的结果。 如 >> f1=[2 2 2]; >> f2=[1 4 9]; >> y=conv(f1,f2) y = 2 10 28 26 18 实验原理与说明 对于有限长序列，我们建立一个通用函数，它可以计算并画出两个有限长序列卷积的结果和波形。能使三个波形的横坐标统一，间隔相同。卷积结果显示在横坐标的中间位置。这个函数取名为DSCONV()，程序自己阅读. 计算示例 例1 计算下列离散信号的能量或功率。 k=0:10;fk=3*(0.5).^k;E=sum(abs(fk).^2) k=0:3;fk=6*cos(0.5*pi.*k); E=sum(abs(fk).^2);P=E/4 k=0:3;fk=6*exp(j*0.5*pi.*k); E=sum(abs(fk).^2);P=E/4 计算示例 例2 求下述差分方程的解 其中输入信号 ，初始条件 ， 。 解 Matlab程序如下： % 计算例2的程序 exp23_1.m a=[-1.5 1];b=[0 0 2]; y0=[2 1];f0=[0 0]; n=0:30; f=ones(1,length(k)); y=recur(a,b,n,f,f0,y0); stem(n,y,'.'), xlabel('k'),ylabel('y(k)') 运行程序后，系统响应波形如图7-13所示。 例3 用Matlab求下列序列的卷积和。 (a) , 。 (b) ， 运行后显示的波形如图23-2(a)所示。 % 计算离散信号的卷积 exp23_2b.m n1=-2:2;f1=[1 1 1 1 1]; % 序列的起始点,序列值 n2=1:5;f2=n2; % 序列的起始点,序列值 M=6; % 将卷积值显示在中间,左右插入M点 dsconv(f1,n1,f2,n2,M) 在命令窗口显示的卷积结果 y = 1 3 6 10 15 14 12 9 5 运行后显示的波形如图23-2(b)所示。 实验内容 1、画出下列各信号的波形，求能量或功率。 (a) (b) (c) (d) 2、求下列差分方程的零输入响应、零状态响应和全响应。 (a) (b) 实验内容 3、求下列序列的卷积和。 (a)