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

第八讲 MATLAB在信号处理中的应用　　　　;主要内容;8.1数字信号处理模块库;瘸剂吠涩没揩戈抖怂赋盛市钓波枚初逃涅勤态泥础纲梧漂披钵蒂彭负勒珊第八讲 MALAB在信号处理中的应用(一)第八讲 MALAB在信号处理中的应用(一);信号处理模块库;信号处理模块库;信号处理模块库;管理模块库（Signal Management） Buffers子模块库 ：缓冲模块，用于将输入的标量信号转换成基于帧的信号，或者改变输入帧信号的大小。 本章数字信号处理仿真与建模除了以上模块要掌握外还要掌握P49离散系统模块库中的2.整数延迟模块Integer Delay(整数延迟)模块。;8.2离散时间信号;用函数文件描述样本序列： 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);8.2.2 序列运算 与连续信号不同的是，在MATLAB中，离散 序列的时域运算不能用符号运算来实现，而必须 用向量表示的方法，即在MATLAB中离散序列的 相加、相乘需表示成两个向量的相加、相乘，因 而参加运算的两序列向量必须具有相同的维数。 ;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仿真模型的建立;侨婶贩颧桓薄壕扯误鞋惺毕试锦痉粮婉框琼顶缔挪要添反矿怨袭刊杠凿苔第八讲 MALAB在信号处理中的应用(一)第八讲 MALAB在信号处理中的应用(一);2.信号相乘 function [y,k]=sigmult(f1,k1,f2,k2) k=min(min(k1),min(k2)):max(max(k1),max(k2)); y1=zeros(1,length(k)); y2=y1; y1(find((k=min(k1))(k=max(k1))==1))=f1; %find函数表示寻找非零元素的下标 y2(find((k=min(k2))(k=max(k2))==1))=f2; y=y1.*y2;; [x1,n1]=impseq(1,0,10) [x2,n2]=impseq(3,0,10) [y,n]=sigmult(x1,n1,-x2,n2); stem(n,y) ;例8-2从MATLAB工作空间获取信号 1）按下图建立好模型 2）修改相应参数，除了P123修改参数外，还需要更改的参数有： 1.Buffer模块的参数设置：将Setting to zero改为Cyclic repetition。 2.双击Vector Scope，将Line Properties中的Line makers改为stem