武汉理工Matlab实验报告程序.docVIP

实验一： 1、利用zeros和ones函数来产生单位阶跃序列。 函数格式： zeros(1,n) ones(1,n) n为序列长度。 程序如下： n1=-2;n2=8;k=0; n=n1:n2; nt=length(n); nk=abs(k-n1)+1; x=[zeros(1,nk-1),ones(1,nt-nk+1)]; stem(n,x) 图1 单位阶跃序列 2、已知信号 x1(n)=3e -0.25n(4n10) x2(n)=u(n+1) (-2n6) 求x(n)=x1(n)*x2(n)，并画出对应波形。 程序如下： function[x,n]=stepseq(n1,n2,n0) n=[n1:n2]; x=[(n-n0)=0]; function [y,ny]=conv_m(x,nx,h,nh,p) nyb=nx(1)+nh(1); nyc=nx(length(x))+nh(length(h)); ny=[nyb:p:nyc]; y=conv(x,h); p=1; n1=[-4:p:10]; f1=3*exp(-0.25*n1); n2=[-2:p:6]; [f2,n2]=stepseq(-2,6,-1); [y,ny]=conv_m(f1,n1,f2,n2,p); subplot(3,1,1);stem(n1,f1) subplot(3,1,2);stem(n2,f2) subplot(3,1,3);stem(ny,y) 图2 对应波形 实验二 验证线性性质： 如果有两个有限长序列分别为x1(n)和x2(n)，长度分别为N1和N2，且y(n)=ax1(n)+bx2(n)，(a，b均为常数)，则该y(n)的N点DFT为 Y(k)=DFT[y(n)]=aX1(k)+bX2(k) (0=k=N-1) 其中：N=max(N1,N2)，X1(k)和X2(k)分别为x1(n)和x2(n)的N点DFT。 序列：x1(n)=[0,1,2,4] x2(n)=[1,0,1,0,1] 程序如下： xn1=[0,1,2,4];xn2=[1,0,1,0,1]; N1=length(xn1);N2=length(xn2); N=max(N1,N2); if N1N2 xn2=[xn2,zeros(1,N1-N2)]; elseif N2N1 xn1=[xn1,zeros(1,N2-N1)]; end yn=2*xn1+3*xn2; n=0:N-1; k=0:N-1; yk1=yn*(exp(-j*2*pi/N).^(n*k)); xk1=xn1*(exp(-j*2*pi/N).^(n*k)); xk2=xn2*(exp(-j*2*pi/N).^(n*k)); yk2=2*xk1+3*xk2; subplot(2,2,1); stem(n,yk1) subplot(2,2,2); stem(n,yk2) subplot(2,2,3); stem(n,xk1) 图1 对应波形 subplot(2,2,4); stem(n,xk2) 实验三 1.用函数conv和FFT计算同一序列的卷积，比较其计算时间。 clock函数读取瞬时时钟 etime(t1,t2)函数计算时刻t1,t2间所经历的时间。 FFT和卷积运算速度 程序如下： L=5000; N=L*2-1; n=1:L; x1=0.5*n; x2=2*n; t0=clock; yc=conv(x1,x2); conv_time=etime(clock,t0) t0=clock; yf=ifft(fft(x1,N).*fft(x2,N)); fft_time=etime(clock,t0) 运算结果 conv_time = 0.0410 fft_time = 0.0040 2.DFT与FFT运算速度 fs=100; N=128; n=0:N-1; t=n/fs; x1=0.5*n; t0=clock; yd=dft(x1,N); dft_time=etime(clock,t0) t0=clock; yf=fft(x1,N); fft_time=etime(clock,t0) 运算结果 dft_time = 0.0920 fft_time = 0 3.序列长度和fft的长度对信号频谱的影响 fs=100; Ndata=32; Nfft=32; n=0:Ndata-1; t=n/fs; x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t); y=fft(x,Nfft); mag=abs(y); f=(0:le