深入浅出通信原理8.doc

目录连载354：利用DFT进行频谱分析2连载355：如何提高频谱密度3连载356：如何提高频谱分辨率（一）5连载357：如何提高频谱分辨率（二）6连载358：泄漏效应8连载359：为什么截断会产生泄漏（一）9连载360：为什么截断会产生泄漏（二）12连载361：为什么截断会产生泄漏（三）13连载362：为什么截断会产生泄漏（四）16连载363：为什么截断会产生泄漏（五）19连载364：为什么截断会产生泄漏（六）20连载365：为什么会产生频谱泄漏（七）21连载366：为什么会产生频谱泄漏（八）24连载367：为什么会产生频谱泄漏（九）25连载368：如何减小频谱泄漏26连载369：频谱的主瓣和旁瓣28连载370：为什么会出现主瓣和旁瓣30连载354：利用DFT进行频谱分析1.jpg (53.5 KB)2011-9-12 00:212.jpg (37.66 KB)2011-9-12 00:21 fs=400;%采样频率 T=0.04;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); subplot(211);stem(t,x); n=0:fs/N:fs*(N-1)/N; X=dft1(x,N); subplot(212);stem(n,abs(X); 连载355：如何提高频谱密度1.jpg (56.95 KB)2011-9-14 23:38 fs=400;%采样频率 T=0.04;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:256-N)=0; x=x1 x2; t=0:1/fs:16*T-1/fs; subplot(211);plot(t,x); n=0:fs/256:255*fs/256; X=dft1(x,256); subplot(212);plot(n,abs(X);连载356：如何提高频谱分辨率（一）1.jpg (55.86 KB)2011-9-16 21:53 fs=400;%采样频率 T=0.04*2;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:256-N)=0; x=x1 x2; t=0:1/fs:8*T-1/fs; subplot(211);plot(t,x); n=0:fs/256:255*fs/256; X=dft1(x,256); subplot(212);plot(n,abs(X); 连载357：如何提高频谱分辨率（二）1.jpg (42.94 KB)2011-9-19 21:42 fs=400;%采样频率 T=0.04*4;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:256-N)=0; x=x1 x2; t=0:1/fs:4*T-1/fs; subplot(211);plot(t,x); n=0:fs/256:255*fs/256; X=dft1(x,256); subplot(212);plot(n,abs(X);2.jpg (48.01 KB)2011-9-19 21:42 fs=400;%采样频率 T=0.04*8;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:256-N)=0; x=x1 x2; t=0:1/fs:2*T-1/fs; subplot(211);plot(t,x); n=0:fs/256:255*fs/256; X=dft1(x,256); subplot(212);plot(n,abs(X);3.jpg (43.8 KB)2011-9-19 21:42 for m=3:-1:1 fs=400;%采样频率 T=0.04*2m;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:256-N)=0; x=x1 x2; n=0:fs/256:255*fs/256; X=dft (x,256); subplot(3,1,m);plot(n,abs(X)/max(abs(X); end; 连载358：泄漏效应连载359：为什么截断会产生泄漏（一）1.jpg (48.53 KB)2011-9-23 00:27 fs=400;%采样频率 T=0.04*8;%截取时间长度 N=T*fs; t=-T/4:1/fs:3*T/4-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); stem(t,x1, .); axis(-T/4 +3*T/4 -3.5 +3); hold on; t=-T/4:0.001:3*T/4; x2=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); plot(t,x2, r);2.jpg (21.44 KB)2011-9-23 00:27 t=-T/4:1/fs:3*T/4-1/fs; y1=zeros(1,N/4); y2=ones(1,N/2); y=y1 y2 y1; stem(t,y, .); axis(-T/4 +3*T/4 -3.5 +3);3.jpg (23.2 KB)2011-9-23 00:27 t=-T/4:1/fs:3*T/4-1/fs; s=x1.*y; stem(t,s, .); axis(-T/4 +3*T/4 -3.5 +3); 连载360：为什么截断会产生泄漏（二）1.jpg (24.95 KB)2011-9-24 23:03 fs=400;%采样频率 T=0.04*4;%截取时间长度 N=T*fs; t=0:1/fs: T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); X=dft (x1,N); n=0:fs/N:fs*(N-1)/N; stem(n,abs(X) /max(abs(X) , .);2.jpg (31.2 KB)2011-9-24 23:03 fs=400;%采样频率 T=0.04*4;%截取时间长度 N=T*fs; t=0:1/fs: T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:256-N)=0; x=x1 x2; n=0:fs/256:255*fs/256; X=dft (x,256); stem(n,abs(X) /max(abs(X) , .);连载361：为什么截断会产生泄漏（三）1.jpg (46.11 KB)2011-9-27 23:492.jpg (36.25 KB)2011-9-27 23:49 fs=800;%采样频率 T=0.01;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); stem(t,x, r); hold on; t=0:0.00001:T; y=cos(200*pi*t); plot(t,y);3.jpg (33.6 KB)2011-9-27 23:49 n=0:fs/N：N-1)*fs/N; X=dft1(x,N); stem(n,abs(X), r); axis(0 800 0 4.5);连载362：为什么截断会产生泄漏（四）1.jpg (53.22 KB)2011-9-28 23:59 fs=800;%采样频率 T=0.02;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); subplot(211);stem(t,x); n=0:fs/N：N-1)*fs/N; X=dft1(x,N); subplot(212);stem(n,abs(X);2.jpg (46.72 KB)2011-9-28 23:59 fs=800;%采样频率 T=0.04;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); subplot(211);stem(t,x); n=0:fs/N：N-1)*fs/N; X=dft1(x,N); subplot(212);stem(n,abs(X);3.jpg (7.06 KB)2011-9-28 23:59连载363：为什么截断会产生泄漏（五）1.jpg (54.29 KB)2011-9-29 23:47 fs=800;%采样频率 T=0.0425;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); subplot(211);stem(t,x); n=0:fs/N:fs*(N-1)/N; X=dft1(x,N); subplot(212);stem(n,abs(X); 连载364：为什么截断会产生泄漏（六）1.jpg (70.51 KB)2011-10-5 21:29 fs=800;%采样频率 T=0.04;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); stem(t,x); t=T:1/fs:2*T-1/fs; hold on; stem(t,x); t=2*T:1/fs:3*T-1/fs; hold on; stem(t,x);2.jpg (58.35 KB)2011-10-5 21:29 fs=800;%采样频率 T=0.0425;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); stem(t,x); t=T:1/fs:2*T-1/fs; hold on; stem(t,x); t=2*T:1/fs:3*T-1/fs; hold on; stem(t,x); 连载365：为什么会产生频谱泄漏（七）1.jpg (59.59 KB)2011-10-6 21:51 fs=800;%采样频率 T=0.02;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t); x2(1:32-N)=0; x=x1 x2; t=0:1/fs:31/fs; stem(t,x); T1=2*T; t=T1:1/fs:2*T1-1/fs; hold on; stem(t,x); t=2*T1:1/fs:3*T1-1/fs; hold on; stem(t,x);2.jpg (43.29 KB)2011-10-6 21:51 fs=800;%采样频率 T=0.02;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t); x2(1:32-N)=0; x=x1 x2; t=0:1/fs:31/fs; subplot(211);stem(t,x); n=0:fs/32:31*fs/32; X=dft (x,32); subplot(212);stem(n,abs(X) /max(abs(X) , .); 连载366：为什么会产生频谱泄漏（八）1.jpg (61.43 KB)2011-10-7 22:302.jpg (48.98 KB)2011-10-7 22:30 fs=740;%采样频率 T=0.04;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x=cos(200*pi*t); subplot(211);stem(t,x); n=0:fs/N:fs*(N-1)/N; X=dft1(x,N); subplot(212);stem(n,abs(X); axis(0 800 0 15); 连载367：为什么会产生频谱泄漏（九）连载368：如何减小频谱泄漏1.jpg (72.58 KB)2011-10-12 22:482.jpg (24.98 KB)2011-10-12 22:48 fix(clock) fs=400;%采样频率 T=0.04*80;%截取时间长度由0.04*8改为0.04*80 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t)+sin(100*pi*t)+cos(50*pi*t); x2(1:2560-N)=0; x=x1 x2; t=0:1/fs:2*T-1/fs; subplot(211);plot(t,x); n=0:fs/2560:2559*fs/2560; X=dft1(x,2560); subplot(212);plot(n,abs(X); fix(clock)连载369：频谱的主瓣和旁瓣1.jpg (31.7 KB)2011-10-18 23:33 fs=400;%采样频率 T=0.04*2;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t); x2(1:256-N)=0; x=x1 x2; t=0:1/fs:8*T-1/fs; subplot(211);plot(t,x); n=0:fs/256:255*fs/256; X=dft1(x,256); subplot(212);plot(n,abs(X);2.jpg (31.15 KB)2011-10-18 23:33 fs=400;%采样频率 T=0.04*2;%截取时间长度 N=T*fs; t=0:1/fs:T-1/fs; x1=cos(200*pi*t); x2(1:256-N)=0; x=x1 x2; n=0:fs/256:255*fs/256; X=dft1(x,256); semilogy (n,abs(X); axis(0 400 0.001 100);