数字信号第三章第四章上机.docVIP

课后作业：的频谱，要求频率分辨率为0.5HZ，抽样点数必须为2的整数倍，试确定采样频率、采样长度进行分析。 解：fs=200; N=512; k=0:N-1; x=sin(30*pi*1/fs*k)+sin(80*pi*1/fs*k); Xk=fft(x); f=k/N*fs; plot(f,abs(Xk)); axis([0,100,0,200]); 2、两个数据保存在matlabdata.mat文件中，这两个数据是转子的振动信号数据，一个是正常无故障的振动数据，一个是有油膜涡动故障的数据。数据按照转速的64倍频采集，每周采集64点，采集8个周期，共计512点数据。正常无故障的振动数据进行fft分析谱图应该只有1倍频的转速谱线，油膜涡动故障信号的fft分析谱图有约0.5倍频的谱线，试用fft分析方法辨别出哪个信号是油膜涡动故障信号。 y=fft(SignalDatal); YKK=abs(y); subplot(YKK); axis([1,64,0,max(YKK)+5]); set(gca,xtick,[1:8:64]); set(gca,xticklabel,0|1x|2x|3x|4x|5x|6x|7x|8x); y=fft(SignalData2); YKK=abs(y); subpot(YKK); axis([1,64,max(YKK)+5]); set(gca,xtick,[1:8:64]); set(gca,xticklabel,0|1x|2x|3x|4x|5x|6x|7x|8x); 从图中可以看出，信号1是油膜涡动故障信号。 例3.10 N=24;M=8;m=3; n=0:N-1; x1=0.8.^n;x2=[(n=0)(nM)]; xn=x1.*x2; [xm,nm]=sigshift(xn,n,m); xc=xn(mod(n,8)+1); xcm=xn(mod(n-m,8)+1); subplot(2,2,1);stem(n,xn,.);axis([0,length(n),0,1]);title(x(n)) subplot(2,2,2);stem(nm,xm,.);axis([0,length(nm),0,1]);title(x(n-m)) subplot(2,2,3);stem(n,xc,.);axis([0,length(n),0,1]);title(x(n)的周期延拓) subplot(2,2,4);stem(n,xcm,.);axis([0,length(n),0,1]);title(x(n)的循环位移) 例3.11 N=16;N1=8; n=0:N-1;k=0:N1-1; x1n=exp(j*pi*n/8); X1k=fft(x1n,N); X2k=fft(x1n,N1); x2n=cos(pi*n/8); X3k=fft(x2n,N); X4k=fft(x2n,N1); subplot(2,2,1);stem(n,abs(X1k),.);axis([0,20,0,20]);ylabel(|X1(k)|) title(16点的DFT[x1(n)]) subplot(2,2,2);stem(n,abs(X3k),.);axis([0,20,0,20]);ylabel(|X2k)|) title(16点的DFT[x2(n)]) subplot(2,2,3);stem(k,abs(X2k),.);axis([0,20,0,20]);ylabel(|X1(k)|) title(16点的DFT[x1(n)]) subplot(2,2,4);stem(k,abs(X4k),.);axis([0,20,0,20]);ylabel(|X2(k)|) title(8点的DFT[x2(n)]) 例3.12 N1=8;N2=16; n=0;N1-1;k1=0:N1-1;k2=0:N2-1; w=2*pi*(0:2047)/2048; Xw=(1-exp(-j*4*w))./(1-exp(-j*w)); xn=[(n=0)(n4)]; X1k=fft(xn,N1); X2k=fft(xn,N2); subplot(3,2,1);plot(w/pi,abs(Xw));xlabel(w/|D) subplot(3,2,2);plot(w/pi,angle(Xw));axis([0,2,-pi,pi]);line([0,2],[0,0]); xlabel(w/|D) subplot(3,2,3);stem(k1,abs(X1k),.); xlabel(k(|?=2|Dk/N1));ylabel(|X1(k)|);hold on plot(N1/2*w/pi,abs