深入浅出通信原理385-435

1、连载385：汉明窗>> N=32;>> y=hamming(N);>> n=0:1:N-1;>> stem(n,y, '.');>> axis(0 32 0 1);>> N=32;>> y=hamming(N);>> n=0:1:N-1;>> Y=dft1(y',N);>> stem(n,abs(Y)/max(abs(Y), '.');>> axis(0 32 0 1.1);>> N=256;>> y

2、1=hamming(N/8);>> y2=zeros(1,N*7/8);>> y=y1' y2;>> n=0:1:N-1;>> stem(n,y, '.');>> axis(0 256 0 1);>> N=256;>> y1=hamming(N/8);>> y2=zeros(1,N*7/8);>> y=y1' y2;>> n=0:1:N-1;>> Y=dft1(y,N);>> stem(n,abs(Y) /max(abs

3、(Y), '.');>> axis(0 256 0 1.1);>> N=4096;>> y1=hamming(32);>> y2=zeros(1,N-32);>> y=y1' y2;>> n=0:1:N-1;>> Y=fft(y,N);>> plot(n,abs(Y) /max(abs(Y);>> axis(0 N 0 1.1);>> N=4096;>> y1=hamming(32);>> y2=zeros(1,N-32);&g

4、t;> y=y1' y2;>> n=0:1:N-1;>> Y=fft(y,N);>> semilogy(n,abs(Y) /max(abs(Y);>> axis(0 N 0.000001 1);>> set(gca,'YTickLabel', '-120' '-100' '-80' '-60' '-40' '-20' '0');连载386：布莱克曼窗>> N=32;>>

5、y=blackman(N);>> n=0:1:N-1;>> stem(n,y, '.');>> axis(0 32 0 1);>> N=32;>> y=blackman(N);>> n=0:1:N-1;>> Y=dft1(y',N);>> stem(n,abs(Y)/max(abs(Y), '.');>> axis(0 32 0 1.1);>> N=256;>> y1=blackman(N/8);>> y2=ze

6、ros(1,N*7/8);>> y=y1' y2;>> n=0:1:N-1;>> stem(n,y, '.');>> axis(0 256 0 1);>> N=256;>> y1=blackman(N/8);>> y2=zeros(1,N*7/8);>> y=y1' y2;>> n=0:1:N-1;>> Y=dft1(y,N);>> stem(n,abs(Y) /max(abs(Y), '.');>> a

7、xis(0 256 0 1.1);>> N=4096;>> y1=blackman(32);>> y2=zeros(1,N-32);>> y=y1' y2;>> n=0:1:N-1;>> Y=fft(y,N);>> plot(n,abs(Y) /max(abs(Y);>> axis(0 N 0 1.1);>> N=4096;>> y1=blackman(32);>> y2=zeros(1,N-32);>> y=y1' y2;>&g

8、t; n=0:1:N-1;>> Y=fft(y,N);>> semilogy(n,abs(Y) /max(abs(Y);>> axis(0 N 0.000001 1);>> set(gca,'YTickLabel', '-120' '-100' '-80' '-60' '-40' '-20' '0');连载387：四种窗函数的波形比较>> N=4096;>> y1=ones(1,N);>&g

9、t; y2=hann(N);>> y3=hamming(N);>> y4=blackman(N);>> n=0:1:N-1;>> plot(n,y1, 'r');hold on;>> plot(n,y2, 'b');hold on;>> plot(n,y3, 'g');hold on;>> plot(n,y4, 'm');>> axis(0 N 0 1);连载388：四种窗函数的频谱比较>> N=4096;>>

10、y1a=ones(1,32);>> y1b=hann(32);>> y1c=hamming(32);>> y1d=blackman(32);>> y2=zeros(1,N-32);>> ya=y1a y2;>> yb=y1b' y2;>> yc=y1c' y2;>> yd=y1d' y2;>> Ya=fft(ya,N);>> Yb=fft(yb,N);>> Yc=fft(yc,N);>> Yd=fft(yd,N);>>

11、; n=0:1:N-1;>> semilogy (n,abs(Ya)/max(abs(Ya) , 'r');hold on;>> semilogy (n,abs(Yb)/max(abs(Yb) , 'b');hold on;>> semilogy (n,abs(Yc)/max(abs(Yc) , 'g');hold on;>> semilogy (n,abs(Yd)/max(abs(Yd) , 'm');hold on;>> axis(0 N 0.000001 1);&g

12、t;> set(gca,'YTickLabel', '-120' '-100' '-80' '-60' '-40' '-20' '0');>> N=4096;>> y1a=ones(1,32);>> y1b=hann(32);>> y1c=hamming(32);>> y1d=blackman(32);>> y2=zeros(1,N-32);>> ya=y1a y2;>&g

13、t; yb=y1b' y2;>> yc=y1c' y2;>> yd=y1d' y2;>> Ya=fft(ya,N);>> Yb=fft(yb,N);>> Yc=fft(yc,N);>> Yd=fft(yd,N);>> n=0:1:N-1;>> semilogy (n,abs(Ya)/max(abs(Ya) , 'r');hold on;>> semilogy (n,abs(Yb)/max(abs(Yb) , 'b');hold on

14、;>> semilogy (n,abs(Yc)/max(abs(Yc) , 'g');hold on;>> semilogy (n,abs(Yd)/max(abs(Yd) , 'm');hold on;>> axis(0 N/2-1 0.000001 1);>> set(gca,'YTickLabel', '-120' '-100' '-80' '-60' '-40' '-20' '0')

15、;连载389：窗函数的应用>> 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);>> x1a=hamming(N) ;>> x1b=blackman(N) ;>> x2(1:256-N)=0;>> x=x1 x2;>> xa=x1.*x1a' x2;>> xb=x1.*x1b' x2;

16、>> X=dft1(x,256);>> Xa=dft1(xa,256);>> Xb=dft1(xb,256);>> n=0:fs/256:255*fs/256;>> subplot(311);plot(n,abs(X);>> subplot(312);plot(n,abs(Xa);>> subplot(313);plot(n,abs(Xb);连载390：利用IDFT实现OFDM调制（一）连载391：利用IDFT实现OFDM调制（二）连载392：利用IDFT实现OFDM调制（三）连载393：利用IDFT实现OFD

17、M调制（四）连载394：利用IDFT实现OFDM调制（五）连载395：利用IDFT实现OFDM调制（六）2011-11-25 22:38 上传下载附件 (24.5 KB) >> dk= 0,+1/sqrt(2)-j*1/sqrt(2), -1/sqrt(2)-j*1/sqrt(2),0,0,0, +1/sqrt(2)+j*1/sqrt(2), -1/sqrt(2)+j*1/sqrt(2) ;>> sn=ifft(dk,8);>> n=0:1:7;>> subplot(211);stem(n,real(sn);>> subplot(21

18、2);stem(n,imag(sn);连载396：利用IDFT实现OFDM调制（七）>> t=0:0.001:1;>> s=(+1/sqrt(2)+j*1/sqrt(2)*exp(-j*4*pi*t)+( -1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t)+( +1/sqrt(2)-j*1/sqrt(2) *exp(j*2*pi*t)+( -1/sqrt(2)-j*1/sqrt(2) *exp(j*4*pi*t);>> plot3(real(s),t,imag(s), 'b');grid on;>> cl

19、f;>> subplot(211);plot (t,real(s), 'b');grid on;>> subplot(212);plot(t,imag(s), 'b');grid on;连载397：利用IDFT实现OFDM调制（八）>> t=0:0.001:1;>> s=(+1/sqrt(2)+j*1/sqrt(2)*exp(-j*4*pi*t)+( -1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t)+( +1/sqrt(2)-j*1/sqrt(2) *exp(j*2*pi*t)+( -1

20、/sqrt(2)-j*1/sqrt(2) *exp(j*4*pi*t);>> subplot(211);plot (t,real(s), 'b'); hold on;>> subplot(212);plot(t,imag(s), 'b'); hold on;>> dk= 0,+1/sqrt(2)-j*1/sqrt(2), -1/sqrt(2)-j*1/sqrt(2),0,0,0, +1/sqrt(2)+j*1/sqrt(2), -1/sqrt(2)+j*1/sqrt(2) ;>> sn=8*ifft(dk,8);&

21、gt;> t=0:0.125:0.125*7;>> subplot(211);stem(t,real(sn), 'r');>> subplot(212);stem(t,imag(sn) , 'r');连载398：利用DFT实现OFDM解调（一）连载399：利用DFT实现OFDM解调（二）连载400：利用DFT实现OFDM解调（三）连载401：利用DFT实现OFDM解调（四）>> t=0:0.001:1;>> s=(+1/sqrt(2)+j*1/sqrt(2)*exp(-j*4*pi*t)+( -1/sqrt(

22、2)+j*1/sqrt(2) *exp(-j*2*pi*t)+( +1/sqrt(2)-j*1/sqrt(2) *exp(j*2*pi*t)+( -1/sqrt(2)-j*1/sqrt(2) *exp(j*4*pi*t);>> plot3(real(s),t,imag(s), 'b');grid on;>> t=0:0.125:0.125*7;>> sn=(+1/sqrt(2)+j*1/sqrt(2)*exp(-j*4*pi*t)+( -1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t)+( +1/sqrt(2)-j*

23、1/sqrt(2) *exp(j*2*pi*t)+( -1/sqrt(2)-j*1/sqrt(2) *exp(j*4*pi*t);>> dk=fft(sn,8)/8;>> n=0:1:7;>> subplot(211);stem(n,real(dk);>> subplot(212);stem(n,imag(dk);>> dkdk =Columns 1 through 6 0.0000 - 0.0000i0.7071 - 0.7071i-0.7071 - 0.7071i0.0000 - 0.0000i0.0000 + 0.0000i0.

24、0000 + 0.0000iColumns 7 through 8 0.7071 + 0.7071i-0.7071 + 0.7071i连载402：OFDM采样频率（一）>> t=0:0.001:1;>> s=(+1/sqrt(2)-j*1/sqrt(2)*exp(-j*6*pi*t)+(+1/sqrt(2)+j*1/sqrt(2)*exp(-j*4*pi*t)+( -1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t)+( +1/sqrt(2)-j*1/sqrt(2) *exp(j*2*pi*t)+( -1/sqrt(2)-j*1/sqrt(2)

25、*exp(j*4*pi*t)+ (-1/sqrt(2)+j*1/sqrt(2)*exp(j*6*pi*t) + (-1/sqrt(2)-j*1/sqrt(2)*exp(j*8*pi*t);>> plot3(real(s),t,imag(s), 'b');grid on;>> t=0:0.125:0.125*7;>> sn=(+1/sqrt(2)-j*1/sqrt(2)*exp(-j*6*pi*t)+(+1/sqrt(2)+j*1/sqrt(2)*exp(-j*4*pi*t)+( -1/sqrt(2)+j*1/sqrt(2) *exp(-j*2

26、*pi*t)+( +1/sqrt(2)-j*1/sqrt(2) *exp(j*2*pi*t)+( -1/sqrt(2)-j*1/sqrt(2) *exp(j*4*pi*t)+ (-1/sqrt(2)+j*1/sqrt(2)*exp(j*6*pi*t) + (-1/sqrt(2)-j*1/sqrt(2)*exp(j*8*pi*t);>> dk=fft(sn,8)/8;>> n=0:1:7;>> subplot(211);stem(n,real(dk);>> subplot(212);stem(n,imag(dk);>> dkdk =Co

27、lumns 1 through 6 0.0000 + 0.0000i0.7071 - 0.7071i-0.7071 - 0.7071i-0.7071 + 0.7071i-0.7071 - 0.7071i0.7071 - 0.7071iColumns 7 through 8 0.7071 + 0.7071i-0.7071 + 0.7071i连载403：OFDM采样频率之二>> x=0:0.001:2;>> y=sin(2*pi*5*x);>> plot(x,y);>> hold on;>> x=0:0.1:2;>> y=si

28、n(2*pi*5*x);>> stem(x,y, 'r');>> x=0:0.001:2;>> y=sin(4*pi*5*x);>> plot(x,y);>> hold on;>> x=0:0.1:2;>> y=sin(4*pi*5*x);>> stem(x,y, 'r');连载404：OFDM采样频率（三）>> t=0:0.125:0.125*7;>> y=exp(j*2*pi*4*t);>> stem3(real(y),imag

29、(y),t, 'r');>> hold on;>> t=0:0.001:1;>> y=exp(j*2*pi*4*t);>> plot3(real(y),imag(y),t);>> set(gca,'YDir','reverse');>> subplot(121);>> t=0:0.125:1;>> y=exp(j*2*pi*4*t);>> stem(t,real(y), 'r');hold on;>> hold

30、on;>> t=0:0.001:1;>> y=exp(j*2*pi*4*t);>> plot(t,real(y);>> subplot(122);>> t=0:0.125:1;>> y=exp(j*2*pi*4*t);>> stem(t,imag(y), 'r');hold on;>> hold on;>> t=0:0.001:1;>> y=exp(j*2*pi*4*t);>> plot(t,imag(y);连载405：OFDM采样频率（四）>

31、;> t=0:0.125:0.125*7;>> y=exp(j*2*pi*(-4)*t);>> stem3(real(y),imag(y),t, 'r');>> hold on;>> t=0:0.001:1;>> y=exp(j*2*pi*(-4)*t);>> plot3(real(y),imag(y),t);>> set(gca,'YDir','reverse');>> subplot(121);>> t=0:0.125:1;>

32、;> y=exp(j*2*pi*(-4)*t);>> stem(t,real(y), 'r');hold on;>> hold on;>> t=0:0.001:1;>> y=exp(j*2*pi*(-4)*t);>> plot(t,real(y);>> subplot(122);>> t=0:0.125:1;>> y=exp(j*2*pi*(-4)*t);>> stem(t,imag(y), 'r');hold on;>> hold on

33、;>> t=0:0.001:1;>> y=exp(j*2*pi*(-4)*t);>> plot(t,imag(y);连载406：OFDM采样频率（五）>> t=0:0.001:1;>> s=exp(-j*8*pi*t)+exp(j*8*pi*t);>> plot3(real(s),t,imag(s), 'b');grid on;>> set(gca,'YDir','reverse');>> t=0:0.125:0.125*7;>> sn=e

34、xp(-j*8*pi*t)+exp(j*8*pi*t);>> dk=fft(sn,8)/8;>> n=0:1:7;>> subplot(211);stem(n,real(dk);>> subplot(212);stem(n,imag(dk);连载407：OFDM采样频率（六）本帖最后由 chenaijun 于 2013-5-5 17:50 编辑 2011-12-15 00:08 上传下载附件 (46.13 KB) >> t=0:0.001:1;>> s=exp(-j*6*pi*t)+exp(j*6*pi*t);>&g

35、t; plot3(real(s),t,imag(s), 'b');grid on;>> set(gca,'YDir','reverse');>> t=0:0.125:0.125*7;>> sn=exp(-j*6*pi*t)+exp(j*6*pi*t);>> dk=fft(sn,8)/8;>> n=0:1:7;>> subplot(211);stem(n,real(dk);>> subplot(212);stem(n,imag(dk);连载408：OFDM采样频率（

36、七）连载409：OFDM采样频率（八）>> N=8;k=5;>> for n=0:1:N-1>>f(n+1)=exp(j*2*pi*k*n/N);>> end;>> n=0:1:N-1;>> plot3(real(f(n+1),n, imag(f(n+1) , 'or');>> set(gca,'YDir','reverse');>> axis(-1 1 0 10 -1 1);>> t=0:0.01:N;>> f=exp(-j*

37、2*pi*(N-k)*t/N);>> hold on;>> plot3(real(f),t,imag(f), 'b'); grid on;>> N=8;k=5;>> for n=0:1:N-1>>f(n+1)=exp(j*2*pi*k*n/N);>> end;>> n=0:1:N-1;>> plot3(real(f(n+1),n, imag(f(n+1) , 'or');>> set(gca,'YDir','reverse'

38、);>> axis(-1 1 0 10 -1 1);>> t=0:0.01:N;>> f=exp(j*2*pi*k*t/N);>> hold on;>> plot3(real(f),t,imag(f), 'b'); grid on;>> N=8;k=5;>> for n=0:1:N-1>>f(n+1)=exp(j*2*pi*k*n/N);>> end;>> n=0:1:N-1;>> subplot(121);plot(n, real(f(n+1)

39、, 'or');hold on;>> subplot(122);plot(n, imag(f(n+1) , 'or');hold on;>> t=0:0.01:N;>> f=exp(-j*2*pi*(N-k)*t/N);>> subplot(121); plot(t, real(f), 'b'); grid on;>> subplot(122); plot(t, imag(f), 'b'); grid on;>> N=8;k=5;>> for n=

40、0:1:N-1>>f(n+1)=exp(j*2*pi*k*n/N);>> end;>> n=0:1:N-1;>> subplot(121);plot(n, real(f(n+1) , 'or');hold on;>> subplot(122);plot(n, imag(f(n+1) , 'or');hold on;>> t=0:0.01:N;>> f=exp(j*2*pi*k*t/N);>> subplot(121); plot(t, real(f), 'b&

41、#39;); grid on;>> subplot(122); plot(t, imag(f), 'b'); grid on;连载410：OFDM信号与周期信号>> t=0:0.001:3;>> s1=(1/sqrt(2)+j*1/sqrt(2)*exp(j*2*pi*t);>> s2= ( 1/sqrt(2)+j*1/sqrt(2) *exp(j*4*pi*t);>> s3= ( 1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t);>> s=s1+s2+s3;>> pl

42、ot3(real(s),t,imag(s), 'b');grid on;set(gca,'YDir','reverse');>> subplot(221);>> plot3(real(s1),t,imag(s1), 'b');grid on;set(gca,'YDir','reverse');axis(-1 1 0 3 -1 1);>> subplot(222);>> plot3(real(s2),t,imag(s2), 'b');gr

43、id on;set(gca,'YDir','reverse');axis(-1 1 0 3 -1 1);>> subplot(223);>> plot3(real(s3),t,imag(s3), 'b');grid on;set(gca,'YDir','reverse');axis(-1 1 0 3 -1 1);>> subplot(224);>> plot3(real(s),t,imag(s), 'b');grid on;set(gca,'Y

44、Dir','reverse');axis(-4 4 0 3 -4 4);>> t=0:0.001:1;>> s1=(1/sqrt(2)+j*1/sqrt(2)*exp(j*2*pi*t);>> s2= ( 1/sqrt(2)+j*1/sqrt(2) *exp(j*4*pi*t);>> s3= ( 1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t);>> s=s1+s2+s3;>> plot3(real(s),t,imag(s), 'b');grid on;se

45、t(gca,'YDir','reverse');axis(-4 4 0 1 -4 4);hold on;>> t=0:0.25:0.75;>> sn1=(1/sqrt(2)+j*1/sqrt(2)*exp(j*2*pi*t);>> sn2= ( 1/sqrt(2)+j*1/sqrt(2) *exp(j*4*pi*t);>> sn3= ( 1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t);>> sn=sn1+sn2+sn3;>> plot3(real(sn),t,im

46、ag(sn), 'or');>> t=0:0.001:3;>> s1=(1/sqrt(2)+j*1/sqrt(2)*exp(j*2*pi*t);>> s2= ( 1/sqrt(2)+j*1/sqrt(2) *exp(j*4*pi*t);>> s3= ( 1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t);>> s=s1+s2+s3;>> plot3(real(s),t,imag(s), 'b');grid on;set(gca,'YDir','

47、reverse');axis(-4 4 0 3 -4 4);hold on;>> t=0:0.25:2.75;>> sn1=(1/sqrt(2)+j*1/sqrt(2)*exp(j*2*pi*t);>> sn2= ( 1/sqrt(2)+j*1/sqrt(2) *exp(j*4*pi*t);>> sn3= ( 1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t);>> sn=sn1+sn2+sn3;>> plot3(real(sn),t,imag(sn), 'or');>

48、> t=0:0.25:2.75;>> sn1=(1/sqrt(2)+j*1/sqrt(2)*exp(j*2*pi*t);>> sn2= ( 1/sqrt(2)+j*1/sqrt(2) *exp(j*4*pi*t);>> sn3= ( 1/sqrt(2)+j*1/sqrt(2) *exp(-j*2*pi*t);>> sn=sn1+sn2+sn3;>> dk=fft(sn,4)/4;>> n=0:1:3;>> subplot(211);stem(n,real(dk);>> subplot(212)

49、;stem(n,imag(dk);>> dkdk =-0.0000 - 0.0000i0.7071 + 0.7071i0.7071 + 0.7071i0.7071 + 0.7071i连载411：傅立叶系数与DFT的关系（一）连载412：傅立叶系数与DFT的关系（二）2011-12-21 00:25 上传下载附件 (22.08 KB) 点评连载413：傅立叶系数与DFT的关系（三）连载414：傅立叶系数与DFT的关系（四）连载415：傅立叶系数与DFT的关系（五）连载416：傅立叶系数与DFT的关系（六）2012-1-8 22:34 上传下载附件 (48.95 KB) 连载417：O

50、FDM信号的FT与DFT>> f=-2:1:2;>> X=1 1 1 1 1;>> stem(f,X);>> axis(-10 +10 -0.3 +1.1);连载418：调制技术>> f=10000;>> t=0:0.000001:0.0003;>> y=1.2+sin(2*pi*f*t);>> plot(t,y);>> axis(0 0.0003 -2.2 2.2);>> f=1

51、00000;>> t=0:0.000001:0.0003;>> y=cos(2*pi*f*t);>> plot(t,y);>> axis(0 0.0003 -1 1);>> f1=10000;>> f2=100000;>> t=0:0.000001:0.0003;>> x= 1.2+sin(2*pi*f1*t);>> c= cos(2*pi*f2*t);>> y= x.*c;>> plot(t,y);hold on ;>> plot(t,abs(x), 'r:') ;hold on ;>> plot(t,-abs(x), 'r:');>&