第八讲MATLAB在信号处理中的应用(一)剖析.pptVIP

8.1数字信号处理模块库 信号处理模块库 信号处理输入模块库（Signal Processing Sources） ：输出一个由Delay变量指定延时时间的离散单位脉冲信号。 ：输出在连续的采样时间内从MATLAB工作空间中获取指定的信号样本。 ：输出一个由变量Constant Value指定的相同维数的常数信号。 信号处理模块库 信号处理输出模块库(Signal Processing Sinks) ：显示一个时域或者频域的矢量信号、矩阵信号、自定义的数据信号。 ：把输入的数据写入MATLAB工作空间的指定阵列中，数据要直到仿真结束或暂停的时候才能获得。 ：用来显示在仿真时产生的信号。 信号处理模块库 变换模块库（Transforms） ：当输入信号是实信号或复数信号时，输出是以2为基（序列点数N为2的整数幂，如果按时域抽选(Decimation-In-Time，DIT) （库利-图基算法）或按频域抽选（Decimation-In-Frequency，DIF）（桑德-图基算法））的快速傅里叶变换（Fast Fourier Transform,FFT）的结果。 ：输出为输入信号的快速离散傅里叶反变换。 信号处理模块库 管理模块库（Signal Management） Buffers子模块库 ：缓冲模块，用于将输入的标量信号转换成基于帧的信号，或者改变输入帧信号的大小。 本章数字信号处理仿真与建模除了以上模块要掌握外还要掌握P49离散系统模块库中的2.整数延迟模块Integer Delay(整数延迟)模块。 8.2离散时间信号 8.2.1常见的离散序列及其表示 1.单位样本序列（类似于连续时间的冲激函数） 其中箭头指出在n=0的样本。 算法思想：用函数zeros（1，N）产生N个零的行向量 利用它可以实现在一个有限区间上的 。 非零值我们可以利用逻辑关系n==0实现样本序列 例如要实现在区间 上的样本序列，可表示为 用函数文件描述样本序列： function[x,n]=imseq(n0,n1,n2) n=n1:n2; x=[n-n0==0]； [x,n]=imseq(3,0,5) stem(n,x) 2.单位阶跃序列（类似于连续时间的单位阶跃函数） 算法思想：用函数ones（1，N）产生N个1的行向量 。利用它可以产生在一个有限区间上的 。零值部 分可以利用逻辑关系n=0实现。 另外我们可以直接利用逻辑关系n=0,实现在 区间上的 ，其表达式为 利用MATLAB函数文件表示为： function [x,n]=stepseq(n0,n1,n2) n=n1:n2; x=[(n-n0)=0]; [x,n]=stepseq(3,0,5) stem(n,x) 1.信号相加、减 在MATLAB中可以用符号“+-”来实现离散信号加减。 但是 和 的长度必须相同，如果序列长度 不相等，或即使长度相等而样值位置不同，也不能 直接用运算符+。所以首先必须对 和 扩大 或延长以使它们具有相同的位置向量，然后才能进 行加、减运算。 例8-1P121例3-1 法一： n=[0:10]; x=impseq(1,0,10)-impseq(3,0,10); stem(n,x); xlabel(n); ylabel(f(n)) 方法二：自定义函数sigadd function [y,k]=sigadd(f1,k1,f2,k2) k=min(min(k1),min(k2)):max(max(k1),max(k2)); y1=zeros(1,length(k));s2=s1 y1(find((k=min(k1))(k=max(k1))==1))=f1; %find函数表示寻找非零元素的下标 y2(find((k=min(k2))(k=max(k2))==1))=f2; y=y1+y2 保存为sigadd.m文件 [x1,n1]=impseq(1,0,10) [x2,n2]=impseq(3,0,10) [y,n]=sigadd(x1,n1,-x2,n2); stem(n,y) 法三：P121例3-1仿真模型的建立 2.信号相乘 function [y,k]=sigmul