646-MATLAB入门及其信号处理应用.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * x1=[x1,zeros(1,N-length(x1))]; %填充序列x1(n)使其长度为N1+N2-1 % (序列x1(n)的长度为N1， 序列x2(n)的长度为N2) x2=[x2,zeros(1,N-length(x2))]; %填充序列x2(n)使其长度为N1+N2-1 n=[0:1:N-1]; x2=x2(mod(-n,N)+1); %生成序列x2((-n))N H=zeros(N,N); for n=1:1:N H(n,:)=cirshiftd(x2,n-1,N); %该矩阵的k行为x2((k-1-n))N end yc=x1*H; %计算循环卷积 function y=cirshiftd(x,m,N) %directly realize circular shift for sequence x %y=cirshiftd(x,m,N); %x:input sequence whose length is less than N %m:how much to shift %N:circular length %y:output shifted sequence if length(x)N error(the length of x must be less than N); end x=[x,zeros(1,N-length(x))]; n=[0:1:N-1]; y=x(mod(n-m,N)+1); (2)　研究两者之间的关系。 clear all; n=[0:1:11]; m=[0:1:5]; N1=length(n); N2=length(m); xn=0.8.^n; %生成x(n) hn=ones(1,N2); %生成h(n) yln=conv(xn,hn); %直接用函数conv计算线性卷积 ycn=circonv(xn,hn,N1); %用函数circonv计算N1点圆周卷积 ny1=(0:1:length(yln)-1); ny2=[0:1:length(ycn)-1]; subplot(2,1,1); stem(ny1,yln); %画图 ylabel(线性卷积) subplot(2,1,2);stem(ny2,ycn);ylabel(圆周卷积); axis([0 16 0 4]) 图 7 线性卷积和圆周卷积的比较 * * * * * * * * * * * 低级语言包括机器语言和汇编语言。机器语言就是计算机指令的集合，它与计算机同时诞生，是第一代的计算机语言；汇编语言是用符号来表示计算机指令，被称为第二代语言 * * * * * * * * * * * * * * * * * Z变换 . 在MATLAB语言中有专门对信号进行正反Z变换的函数ztrans( ) 和itrans( )。其调用格式分别如下： ? F=ztrans( f ) ?????对f(n)进行Z变换，其结果为F(z) ??F=ztrans(f,v)? ?对f(n)进行Z变换，其结果为F(v) ??F=ztrans(f,u,v) ???对f(u)进行Z变换，其结果为F(v) ? f=itrans ( F )? ????对F(z)进行Z反变换，其结果为f(n) ??f=itrans(F,u)? ?对F(z)进行Z反变换，其结果为f(u) ??f=itrans(F,v,u ) ???对F(v)进行Z反变换，其结果为f(u) 注意： 在调用函数ztran( )及iztran( )之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量 Z变换 例 求数列 fn=e-n的Z变换及其逆变换。命令如下： syms n z fn=exp(-n); Fz=ztrans(fn,n,z) %求fn的Z变换 f=iztrans(Fz,z,n) %求Fz的逆Z变换 例①．用MATLAB求出离散序列?的Z变换MATLAB程序如下： syms k z f=0.5^k;????????? %定义离散信号 Fz=ztrans(f)?????? %对离散信号进行Z变换