MATLAB程序设计语言在信号处理中的应用.pptVIP

实现程序：  计算圆周卷积的函数。 function yc=circonv(x1,x2,N) %直接计算圆周卷积 y=circonv(x1,x2,N) %输出参数： 圆周卷积结果y %输入参数： 需要计算圆周卷积的序列x1,x2和圆周卷积的点数N if length(x1)N error(′N must not be less than length of x1′); end if length(x2)N error(′N must not be less than length of x2′); end %以上语句判断两个序列的长度是否小于N 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］); 图 6-7 线性卷积和圆周卷积的比较 6.5.2 利用离散傅里叶变换（DFT）分析信号的频谱 例 6-11 已知序列x(n)=2 sin(0.48πn)+cos(0.52πn) 0≤n100, 试绘制x(n)及它的离散傅里叶变换|X(k)|图。  MATLAB实现程序： clear all N=100; n=0:N-1; xn=2*sin(0.48*pi*n)+cos(0.52*pi*n); XK=fft(xn,N); magXK=abs(XK); phaXK=angle(XK); subplot(1,2,1) plot(n,xn) xlabel(′n′);ylabel(′x(n)′); title(′x(n) N=100′);  subplot(1,2,2) k=0:length(magXK)-1; stem(k,magXK,′.′); xlabel(′k′);ylabel(′|X(k)|′); tit