深入淺出通信原理二

1、連載 51：利用旋轉向量理解正交解調（一）3連載 52：利用旋轉向量理解正交解調（二）10連載 53：利用旋轉向量理解正交解調（三）12連載 54：PSK/QAM 調制僅僅是指部分嗎？13連載 55：調制解調與級數(shù)展開的關系14連載 56：利用求復系數(shù)的方法實現(xiàn)解調15連載 57：如何求復系數(shù)？16連載 58：OFDM 與級數(shù)展開18連載 59：如何求系數(shù)？20連載 60：周期方波信號的復系數(shù)22連載 61：sinc 函數(shù)23連載 62：周期方波信號的頻譜24連載 63：周期矩形波的頻譜26連載 64：周期矩形波頻譜對比（一）30連載 65：周期矩形波的頻譜對比（二）32連載 66：非周期矩形

2、信號的頻譜35連載 67：連續(xù)型頻譜36連載 68：周期矩形波的連續(xù)譜38連載 69：周期矩形波的連續(xù)譜和離散譜對比44連載 70：非周期矩形信號的連續(xù)譜46連載 71：非周期信號的連續(xù)譜（一）48連載 72：非周期信號的連續(xù)譜（二）491連載 73：非周期信號的連續(xù)譜（三）50連載 74：非周期信號的連續(xù)譜（四）52連載 75：已知頻譜求非周期信號54連載 76：變換56連載 77：對調制過程進行深入的頻譜分析57連載 78：信號與復指數(shù)信號相乘的頻譜變化59連載 79：矩形脈沖信號調制余弦載波（時域）61連載 80：矩形脈沖信號調制余弦載波（頻域）64連載 81：矩形脈沖調制余弦載波（三）

3、66連載 82：矩形脈沖信號調制余弦載波（四）68連載 83：正負脈沖的幅度譜和相位譜70連載 84：采用對數(shù)坐標的矩形脈沖幅度譜74連載 85：BPSK 調制信號的頻譜（一）75連載 86：BPSK 調制信號的頻譜（二）77連載 87：調制正弦載波前后的信號頻譜變化80連載 88：矩形脈沖調制余弦和正弦載波的頻譜對比82連載 89：QPSK 調制信號的頻譜（一）84連載 90：QPSK 調制信號的頻譜（二）86連載 91：BPSK 解調的頻域分析（一）89連載 92：BPSK 解調的頻域分析（二）93連載 93：在時域進行 BPSK 解調95連載 94：在時域進行 QPSK 解調1002連載

4、 95：QPSK 解調的頻域分析104連載 96：信號的頻譜分析方法可否統(tǒng)一？113連載 97：沖激函數(shù)114連載 98：周期信號的變換115連載 99：復指數(shù)信號的變換117連載 100：余弦信號的變換119連載 47：利用旋轉向量理解 BPSK 解調（三）34連載 48：用復數(shù)運算實現(xiàn) BPSK 調制和解調5連載 49：利用實數(shù)運算實現(xiàn) BPSK 調制和解調下面這段代碼可以畫出一個旋轉向量：>> t=-5:0.001:5;>> x=cos(2*pi *t);>> y=sin(2*pi *t);>> plot3(x,t,y);>>

5、xlabel('x');>> ylabel('t');>> zlabel('y');>> set(gca,'YDir','reverse');>> grid on;直接拷貝到中運行即可。6連載 50：利用旋轉向量理解正交調制78連載 51：利用旋轉向量理解正交解調（一）9連載 52：利用旋轉向量理解正交解調（二）1011連載 53：利用旋轉向量理解正交解調（三）12連載 54：PSK/QAM 調制僅僅是指嗎？部分13連載 55：調制解調與系級數(shù)展開的關14連載 56：

6、利用求復解調系數(shù)的方法實現(xiàn)如何求復系數(shù)呢？15連載 57：如何求復系數(shù)？1617連載 58：OFDM 與級數(shù)展開1819連載 59：如何求系數(shù)？2021連載 60：周期方波信號的復系數(shù)22連載 61：sinc 函數(shù)23連載 62：周期方波信號的頻譜2425連載 63：周期矩形波的頻譜26272829連載 64：周期矩形波頻譜對比（一）>> subplot(2,1,1)>> x=-8:1:8;>> stem(x,0.5.*sinc(0.5.*x)>> subplot(2,1,2)>> x=-16:1:16;>> stem(x

7、,0.25.*sinc(0.25.*x)>> axis(-16 16 -0.3 +0.3);>> grid on;3031連載 65：周期矩形波的頻譜對比（二）>> x=-8:0.001:8;>> subplot(3,1,1)>>plot(x,square(2*pi*x+0.5*pi,50)>> axis(-8 8 -1.5 +1.5);>> subplot(3,1,2)>> plot(x,square(pi*x+0.25*pi,25)>> axis(-8 8 -1.5 +1.5);&g

8、t;> subplot(3,1,3)32>> plot(x,square(0.5*pi*x+0.125*pi,12.5)>> axis(-8 8 -1.5 +1.5);>> subplot(3,1,1)>> x=-8:1:8;>> stem(x,0.5.*sinc(0.5.*x)>> subplot(3,1,2)>> x=-16:1:16;>> stem(x,0.25.*sinc(0.25.*x)>> axis(-16 16 -0.3 +0.3);>> subplot(

9、3,1,3)>> x=-32:1:32;>> stem(x,0.125.*sinc(0.125.*x)>> axis(-32 32 -0.2 +0.2);將縱坐標調整到一致：3334連載 66：非周期矩形信號的頻譜>> subplot(1,2,1)>> x=-8:0.001:8;>> plot(x,square(0.5/4*pi*x+0.125/4*pi,12.5/4)>> axis(-8 8 -1.5 +1.5);>> subplot(1,2,2)>> x=-128:1:128;>

10、;> stem(x,0.125/4.*sinc(0.125/4.*x)>> axis(-128 128 -0.02 +0.05);>> box on35連載 67：連續(xù)型頻譜3637連載 68：周期矩形波的連續(xù)譜383940414243連載 69：周期矩形波的連續(xù)譜和離散譜對比>> subplot(3,2,1)>> x=-8:1:8;>> stem(x,0.5.*sinc(0.5.*x)>> axis(-8 8 -0.2 +0.6);>> title('T=T_0');>> s

11、ubplot(3,2,3)>> x=-16:1:16;>> stem(x,0.25.*sinc(0.25.*x)>> axis(-16 16 -0.2 +0.6);>> title('T=2T_0');44>> subplot(3,2,5)>> x=-32:1:32;>> stem(x,0.125.*sinc(0.125.*x)>> axis(-32 32 -0.2 +0.6);>> title('T=4T_0');>> subplot(3,2

12、,2)>> d=zeros(16,2);>> x=-8:1:8;>> for i=1:16>> d(i)=x(i);>> d(16+i)= 0.5.*sinc(0.5.*x(i);>> end>> x=-8:0.001:8;>> plot(x,0.5.*sinc(0.5*x), 'r:');>> hold on;>> s=pulstran(x-0.5,d,'rectpuls',1);plot(x,s) ;>> axis(-8 8 -

13、0.2 +0.6);>> subplot(3,2,4)>> d=zeros(32,2);>> x=-16:1:16;>> for i=1:32>> d(i)=x(i);>> d(32+i)= 0.5.*sinc(0.25*x(i);>> end>> x=-16:0.001:16;>> plot(x,0.5.*sinc(0.25*x), 'r:');>> hold on;>> s=pulstran(x-0.5,d,'rectpuls'

14、,1);plot(x,s) ;>> axis(-16 16 -0.2 +0.6);>> subplot(3,2,6)>> d=zeros(64,2);45>> x=-32:1:32;>> for i=1:64>> d(i)=x(i);>> d(64+i)= 0.5.*sinc(0.125*x(i);>> end>> x=-32:0.001:32;>> plot(x,0.5.*sinc(0.125*x), 'r:');>> hold on;>&

15、gt; s=pulstran(x-0.5,d,'rectpuls',1);plot(x,s) ;>> axis(-32 32 -0.2 +0.6);連載 70：非周期矩形信號的連續(xù)譜46>> x=-8:0.001:8;>> y=sinc(x);>> plot(x,y)>> grid on47連載 71：非周期信號的連續(xù)譜（一）48連載 72：非周期信號的連續(xù)譜（二）49連載 73：非周期信號的連續(xù)譜（三）5051連載 74：非周期信號的連續(xù)譜（四）5253連載 75：已知頻譜求非周期信號5455連載 76：變換56連載

16、 77：對調制過程進行深入的頻譜分析5758連載 78：信號與復指數(shù)信號相乘的頻譜變化5960連載 79：矩形脈沖信號調制余弦載波（時域）>> t=-1.5:0.001:1.5;>> y=rectpuls(t,1);>> plot(t,y)>> axis(-1.5 1.5 -0.5 1.5);>> grid on61>> t=-1.5:0.001:1.5;>> z=cos(2*pi*5*t);>> plot(t,z)>> axis(-1.5 1.5 -1.5 1.5);>>

17、grid on62>> t=-1.5:0.001:1.5;>> y=rectpuls(t,1);>> z=cos(2*pi*5*t);>> plot(t,y.*z)>> axis(-1.5 1.5 -1.5 1.5);>> grid on63連載 80：矩形脈沖信號調制余弦載波（頻域）64>> f=-4:0.001:4;>> X=sinc(f);>> plot(f,X)>> grid on>> f=-10:0.001:10;>> X=0.5.*sinc

18、(f-5)+ 0.5.*sinc(f+5);>> plot(f,X)>> grid on65連載 81：矩形脈沖調制余弦載波（三）>> f=-1200:0.001:1200;>> X=0.5.*sinc(f-1000)+ 0.5.*sinc(f>> plot(f,X)>> axis(-1200 1200 -0.2 0.6);>> grid on);6667連載 82：矩形脈沖信號調制余弦載波（四）>> t=-1.5:0.001:1.5;>> y=-rectpuls(t,1);>&g

19、t; plot(t,y)>> axis(-1.5 1.5-1.50.5);>> grid on68>> f=-4:0.001:4;>> X=-sinc(f);>> plot(f,X)>> grid on69連載 83：正負脈沖的幅度譜和相位譜>> f=-4:0.001:4;>> X=abs(sinc(f);>> subplot(2,1,1);plot(f,X)>> grid on>> f= -4,-3,-2,-1,0,1,2,3,4;>> P= pi,

20、0,pi,0,0,-pi,0,-pi,0;>> subplot(2,1,2);stairs(f,P)>> grid on70>> f=-4:0.001:4;>> X=abs(-sinc(f);>> subplot(2,1,1);plot(f,X)>> grid on>> f= -4,-3,-2,-1,0,1,2,3,4;>> P=0,pi,0,pi,-pi,0,-pi,0,-pi;>> subplot(2,1,2);stairs(f,P)>> grid on71%正脈沖%&g

21、t;> t=-1.5:0.001:1.5;>> y=rectpuls(t,1);>> subplot(3,2,1);plot(t,y)>> axis(-1.5 1.5 -0.5 1.5);>> grid on%負脈沖%>> t=-1.5:0.001:1.5;>> y=-rectpuls(t,1);>> subplot(3,2,2);plot(t,y)>> axis(-1.5 1.5 -1.5 0.5);>> grid on72%正脈沖幅度頻譜%>> f=-4:0.001

22、:4;>> X=abs(sinc(f);>> subplot(3,2,3);plot(f,X)>> grid on%正脈沖相位頻譜%>> f= -4,-3,-2,-1,0,1,2,3,4;>> P= pi,0,pi,0,0,-pi,0,-pi,0;>> subplot(3,2,5);stairs(f,P)>> grid on%負脈沖幅度頻譜%>> f=-4:0.001:4;>> X=abs(-sinc(f);>> subplot(3,2,4);plot(f,X)>>

23、; grid on%負脈沖相位頻譜%>> f= -4,-3,-2,-1,0,1,2,3,4;>> P=0,pi,0,pi,-pi,0,-pi,0,-pi;>> subplot(3,2,6);stairs(f,P)>> grid on很容易看出：正負脈沖的幅度頻譜是相同的，相位相差 （反相）。73連載 84：采用對數(shù)坐標的矩形脈沖幅度譜%正脈沖幅度頻譜%>> f=-4:0.001:4;>> X=abs(sinc(f);>> semilogy(f,X);>> axis(-4 4 0.001 1.1);&

24、gt;> set(gca,'YTickLabel','-60' '-40' '-20' '0')>> xlabel('Frequency(Hz)');>> ylabel('Amplitude(dB)');>> grid on74連載 85：BPSK 調制信號的頻譜（一）%輸入信號>> subplot(3,1,1);>> t=0:0.001:8;>> d=0 0 ;1 1 ;2 1 ;3 0 ;4 1 ;5

25、1 ;6 0 ;7 0;%時延 1 幅度 1 ;時延 2 幅度 2>> s=pulstran(t-0.5,d,'rectpuls');plot(t,s) ;>> axis(0 8 -0.5 1.5);>> text(0.5,1.2,'0') ; text(1.5,1.2,'1') ; text(2.5,1.2,'1') ; text(3.5,1.2,'0') ;>> text(4.5,1.2,'1') ; text(5.5,1.2,'1'

26、;) ; text(6.5,1.2,'0') ; text(7.5,1.2,'0') ;75%BPSK 基帶調制信號>> subplot(3,1,2);>> t=0:0.001:8;>> d=0 +1 ;1 -1;2 -1;3 +1;4 -1;5 -1;6 +1;7 +1;>> s=pulstran(t-0.5,d,'rectpuls');plot(t,s) ;>> axis(0 8 -1.5 1.5);>> text(0.5,1.2,'+1') ; text

27、(1.5,1.2,'-1') ; text(2.5,1.2,'-1') ; text(3.5,1.2,'+1') ;>> text(4.5,1.2,'-1') ; text(5.5,1.2,'-1') ; text(6.5,1.2,'+1') ; text(7.5,1.2,'+1') ;%BPSK 頻帶調制信號>> subplot(3,1,3);>> t=0:0.001:8;>> d=0 +1 ;1 -1;2 -1;3 +1;4 -1;

28、5 -1;6 +1;7 +1;>> s=pulstran(t-0.5,d,'rectpuls').*cos(2*pi*5*t);plot(t,s) ;>> axis(0 8 -1.5 1.5);>> text(0.5,1.2,'0') ; text(1.5,1.2, 'pi') ; text(2.5,1.2,'pi') ; text(3.5,1.2,'0') ;>> text(4.5,1.2, 'pi') ; text(5.5,1.2, 'pi

29、') ; text(6.5,1.2,'0') ; text(7.5,1.2,'0')；76連載 86：BPSK 調制信號的頻譜（二）>> f=-9:0.001:-1;>> x=0.5*sinc(f+5);>> fill3(x,f,0*f,'r');>> hold on>> f=1:0.001:9;>> x=0.5*sinc(f-5);>> fill3(x,f,0*f,'r');>> grid on>> ylabel(

30、'f');77>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0);>> f=-9:0.001:-1;>> x=-0.5*sinc(f+5);>> fill3(x,f,0*f,'r');>> hold on>&

31、gt; f=1:0.001:9;>> x=-0.5*sinc(f-5);78>> fill3(x,f,0*f,'r');>> grid on>> ylabel('f');>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0

32、 0);>> f=-9:0.001:-1;>> X=abs(0.5*sinc(f+5);>> fill(f,X,'r');>> hold on>> f=1:0.001:9;79>> X=abs(0.5*sinc(f-5);>> fill(f,X,'r');>> grid on>> xlabel('f');>> ylabel('X');>> axis(-10 10 -0.1 0.6);連載化87：調制正

33、弦載波前后的信號頻譜變8081連載 88：矩形脈沖調制余弦和正弦載波的頻譜對比>> subplot(1,2,1);>> f=-9:0.001:-1;>> x=0.5*sinc(f+5);>> fill3(x,f,0*f,'r');>> hold on>> f=1:0.001:9;>> x=0.5*sinc(f-5);82>> fill3(x,f,0*f,'r');>> grid on>> ylabel('f');>>

34、 zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0);>> subplot(1,2,2);>> f=-9:0.001:-1;>> y=0.5*sinc(f+5);>> fill3(0*f,f,y,'b');>> hold on>&

35、gt; f=1:0.001:9;>> y=-0.5*sinc(f-5);>> fill3(0*f,f,y,'b');>> grid on>> ylabel('f');>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0

36、);83連載 89：QPSK 調制信號的頻譜（一）84>> a=1/sqrt(2);>> f=-9:0.001:-1;>> x=0.5*a*sinc(f+5);>> y=0.5*a*sinc(f+5);>> fill3(x,f,y,'g');>> hold on>> f=1:0.001:9;>> x=0.5*a *sinc(f-5);>> y=-0.5*a*sinc(f-5);>> fill3(x,f,y,'g');>> grid

37、on>> ylabel('f');85>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0);連載 90：QPSK 調制信號的頻譜（二）>> %數(shù)據(jù)為 10>> subplot(2,2,1);>> a=1/sqrt(2);>&g

38、t; f=-9:0.001:-1;86>> x=0.5*a*sinc(f+5);>> y=0.5*a *sinc(f+5);>> fill3(x,f,y,'g');>> hold on>> f=1:0.001:9;>> x=0.5*a *sinc(f-5);>> y=-0.5*a *sinc(f-5);>> fill3(x,f,y,'g');>> grid on>> ylabel('f');>> zlabel(

39、9;y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0);>> %數(shù)據(jù)為 01>> subplot(2,2,2);>> a=1/sqrt(2);>> f=-9:0.001:-1;>> x=-0.5*a*sinc(f+5);>> y=0.5*a *sinc(f+5);

40、>> fill3(x,f,y,'g');>> hold on>> f=1:0.001:9;>> x=-0.5*a *sinc(f-5);>> y=-0.5*a *sinc(f-5);>> fill3(x,f,y,'g');>> grid on>> ylabel('f');>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6

41、 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0);87>> %數(shù)據(jù)為 11>> subplot(2,2,3);>> a=1/sqrt(2);>> f=-9:0.001:-1;>> x=-0.5*a*sinc(f+5);>> y=-0.5*a *sinc(f+5);>> fill3(x,f,y,'g');>> hold on>> f=1:0.001:9

42、;>> x=-0.5*a *sinc(f-5);>> y=+0.5*a *sinc(f-5);>> fill3(x,f,y,'g');>> grid on>> ylabel('f');>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line

43、(0 0,-10 10,0 0);>> %數(shù)據(jù)為 00>> subplot(2,2,4);>> a=1/sqrt(2);>> f=-9:0.001:-1;>> x=0.5*a*sinc(f+5);>> y=-0.5*a *sinc(f+5);>> fill3(x,f,y,'g');>> hold on>> f=1:0.001:9;>> x=0.5*a *sinc(f-5);>> y=0.5*a *sinc(f-5);>> fill3(x

44、,f,y,'g');>> grid on>> ylabel('f');88>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -10 10 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-10 10,0 0);連載 91：BPSK 解調的頻域分析（一）89909192連載 92：BPSK 解調的頻域分析（二）9394連載

45、 93：在時域進行 BPSK 解調95%端輸入信號>> subplot(4,1,1);>> t=0:0.001:8;>> d=0 0 ;1 1 ;2 1 ;3 0 ;4 1 ;5 1 ;6 0 ;7 0;%時延 1 幅度 1 ;時延 2 幅度 2>> s=pulstran(t-0.5,d,'rectpuls');plot(t,s) ;>> axis(0 8 -0.5 1.5);>> text(0.5,1.2,'0') ; text(1.5,1.2,'1') ; text(2.

46、5,1.2,'1') ; text(3.5,1.2,'0') ;>> text(4.5,1.2,'1') ; text(5.5,1.2,'1') ; text(6.5,1.2,'0') ; text(7.5,1.2,'0') ;%端 BPSK 基帶調制信號>> subplot(4,1,2);>> t=0:0.001:8;>> d=0 +1 ;1 -1;2 -1;3 +1;4 -1;5 -1;6 +1;7 +1;>> s=pulstran(t

47、-0.5,d,'rectpuls');plot(t,s) ;>> axis(0 8 -1.5 1.5);>> text(0.5,1.2,'+1') ; text(1.5,1.2,'-1') ; text(2.5,1.2,'-1') ; text(3.5,1.2,'+1') ;>> text(4.5,1.2,'-1') ; text(5.5,1.2,'-1') ; text(6.5,1.2,'+1') ; text(7.5,1.2,&

48、#39;+1') ;96%端和接收端 BPSK 頻帶調制信號>> subplot(4,1,3);>> t=0:0.001:8;>> d=0 +1 ;1 -1;2 -1;3 +1;4 -1;5 -1;6 +1;7 +1;>> s=pulstran(t-0.5,d,'rectpuls').*cos(2*pi*5*t);plot(t,s) ;>> axis(0 8 -1.5 1.5);>> text(0.5,1.2,'0') ; text(1.5,1.2, 'pi') ;

49、text(2.5,1.2,'pi') ; text(3.5,1.2,'0') ;>> text(4.5,1.2, 'pi') ; text(5.5,1.2, 'pi') ; text(6.5,1.2,'0') ; text(7.5,1.2,'0') ;%接收端 BPSK 頻帶調制信號*cos0t>> subplot(4,1,4);>> t=0:0.001:8;>> d=0 +1 ;1 -1;2 -1;3 +1;4 -1;5 -1;6 +1;7 +1;&

50、gt;> s=pulstran(t-0.5,d,'rectpuls').*cos(2*pi*5*t) .*cos(2*pi*5*t);plot(t,s) ;>> axis(0 8 -1.5 1.5);>> text(0.5,1.2,'+1') ; text(1.5,1.2,'-1') ; text(2.5,1.2,'-1') ; text(3.5,1.2,'+1') ;>> text(4.5,1.2,'-1') ; text(5.5,1.2,'-1&

51、#39;) ; text(6.5,1.2,'+1') ; text(7.5,1.2,'+1') ;97>> subplot(3,1,1);>> t=0:0.001:8;>> d=0 +1 ;1 -1;2 -1;3 +1;4 -1;5 -1;6 +1;7 +1;>> s1=pulstran(t-0.5,d,'rectpuls').*cos(2*pi*5*t) .*cos(2*pi*5*t);>> plot(t,s1) ;>> line(0,8,0,0)>> axis

52、(0 8 -1.5 1.5);>> subplot(3,1,2);>> t=0:0.001:8;>> d=0 0.5 ;1 -0.5;2 -0.5;3 +0.5;4 -0.5;5 -0.5;6 +0.5;7 +0.5;>> s2=pulstran(t-0.5,d,'rectpuls');>> plot(t,s2) ;>> line(0,8,0,0)>> axis(0 8 -1.5 1.5);>> subplot(3,1,3);>> t=1:1:8;>> d=+

53、1,-1,-1,+1,-1,-1,+1,+1;>> stem(t,d) ;>> line(0,8,0,0)>> axis(0 8 -1.5 1.5);9899連載 94：在時域進行 QPSK 解調100%端輸入信號>> subplot(6,1,1);>> t=0:0.001:8;>> d=0 0 ;0.5 1;1 1;1.5 0;2 1 ;2.5 1;3 0;3.5 0;4 0;4.5 1 ;5 1 ;5.5 0 ;6 1 ;6.5 1 ;7 0 ;7.5 0;>> s=pulstran(t-0.25,d,&#

54、39;rectpuls',0.5);plot(t,s) ;>> axis(0 8 -0.5 1.5);>> text(0.25,1.2,'0') ; text(0.75,1.2,'1') ; text(1.25,1.2,'1') ; text(1.75,1.2,'0') ;>> text(2.25,1.2,'1') ; text(2.75,1.2,'1') ; text(3.25,1.2,'0') ; text(3.75,1.2,'

55、0') ;>> text(4.25,1.2,'0') ; text(4.75,1.2,'1') ; text(5.25,1.2,'1') ; text(5.75,1.2,'0') ;>> text(6.25,1.2,'1') ; text(6.75,1.2,'1') ; text(7.25,1.2,'0') ; text(7.75,1.2,'0') ;%端 I 路信號>> subplot(6,1,2);>> t=

56、0:0.001:8;>> a=1/sqrt(2);>> d=0 -a ;1 +a;2 -a;3 +a; 4 -a ;5 +a;6 -a;7 +a;101>> s=pulstran(t-0.5,d,'rectpuls');plot(t,s) ;>> axis(0 8 -2 2);>> text(0.5,1.5,'-0.7') ; text(1.5,1.5,'+0.7') ;text(2.5,1.5,'-0.7') ;text(3.5,1.5,'+0.7');

57、>> text(4.5,1.5,'-0.7') ; text(5.5,1.5,'+0.7') ;text(6.5,1.5,'-0.7') ;text(7.5,1.5,'+0.7');%端 Q 路信號>> subplot(6,1,3);>> t=0:0.001:8;>> d=0 +a;1 -a;2 -a;3 +a; 4 +a;5 -a;6 -a;7 +a;>> s=pulstran(t-0.5,d,'rectpuls');plot(t,s) ;>>

58、; axis(0 8 -2 2);>> text(0.5,1.5,'+0.7') ; text(1.5,1.5,'-0.7') ; text(2.5,1.5,'-0.7') ; text(3.5,1.5,'+0.7')>> text(4.5,1.5,'+0.7') ; text(5.5,1.5,'-0.7') ; text(6.5,1.5,'-0.7') ; text(7.5,1.5,'+0.7')%端和接收端 QPSK 頻帶調制信號>&

59、gt; subplot(6,1,4);>> t=0:0.001:8;>> d1=0 -a ;1 +a;2 -a;3 +a; 4 -a ;5 +a;6 -a;7 +a;>> s1=pulstran(t-0.5,d1,'rectpuls').*cos(2*pi*5*t) ;>> d2=0 +a;1 -a;2 -a;3 +a; 4 +a;5 -a;6 -a;7 +a;>> s2=pulstran(t-0.5,d2,'rectpuls').*sin(2*pi*5*t);>> plot(t,s1-s2

60、) ;>> axis(0 8 -2 2);>> text(0.3,1.5,'3pi/4') ; text(1.3,1.5, '5pi/4') ; text(2.3,1.5,'7pi/4') ; text(3.3,1.5,'pi/4') ;>> text(4.3,1.5, '3pi/4') ; text(5.3,1.5, '5pi/4') ; text(6.3,1.5,'7pi/4') ; text(7.3,1.5,'pi/4') ;

61、%接收端 QPSK 頻帶調制信號乘以 cos0t 后的波形>> subplot(6,1,5);>> t=0:0.001:8;>> d1=0 -a ;1 +a;2 -a;3 +a; 4 -a ;5 +a;6 -a;7 +a;>> s1=pulstran(t-0.5,d1,'rectpuls').*cos(2*pi*5*t) ;>> d2=0 +a;1 -a;2 -a;3 +a; 4 +a;5 -a;6 -a;7 +a;>> s2=pulstran(t-0.5,d2,'rectpuls').*s

62、in(2*pi*5*t);>> s=(s1-s2).* cos(2*pi*5*t) ;102>> plot(t,s) ;>> axis(0 8 -2 2);>> line(0,8,0,0)%接收端 QPSK 頻帶調制信號乘以-sin0t 后的波形>> subplot(6,1,6);>> t=0:0.001:8;>> d1=0 -a ;1 +a;2 -a;3 +a; 4 -a ;5 +a;6 -a;7 +a;>> s1=pulstran(t-0.5,d1,'rectpuls').*co

63、s(2*pi*5*t) ;>> d2=0 +a;1 -a;2 -a;3 +a; 4 +a;5 -a;6 -a;7 +a;>> s2=pulstran(t-0.5,d2,'rectpuls').*sin(2*pi*5*t);>> s=-(s1-s2).* sin(2*pi*5*t) ;>> plot(t,s) ;>> axis(0 8 -2 2);>> line(0,8,0,0)103連載 95：QPSK 解調的頻域分析104>> a=1/sqrt(2);>> f=-6:0.001:-

64、4;>> x=0.5*a*sinc(f+5);>> y=0.5*a *sinc(f+5);>> fill3(x,f,y,'b');>> hold on>> f=4:0.001:6;>> x=0.5*a *sinc(f-5);>> y=-0.5*a *sinc(f-5);>> fill3(x,f,y,'b');>> grid on>> ylabel('f');>> zlabel('y');>>

65、; xlabel('x');105>> axis(-0.6 0.6 -12 12 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-12 12,0 0);106% 一分為二向左平移>> subplot(3,1,1);>> a=1/sqrt(2);>> f=-11:0.001:-9;107>> x=0.25*a*sinc(f+10);>> y=0.25*a *sinc(f+10);>> f

66、ill3(x,f,y,'b');>> hold on>> f=-1:0.001:1;>> x=0.25*a *sinc(f);>> y=-0.25*a *sinc(f);>> fill3(x,f,y,'b');>> grid on>> ylabel('f');>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -12 12 -0.6 0.6);>&g

67、t; set(gca,'YDir','reverse');>> line(0 0,-12 12,0 0);% 一分為二向右平移>> subplot(3,1,2);>> a=1/sqrt(2);>> f=9:0.001:11;>> x=0.25*a*sinc(f-10);>> y=-0.25*a *sinc(f-10);>> fill3(x,f,y,'b');>> hold on>> f=-1:0.001:1;>> x=0.25*

68、a *sinc(f);>> y=0.25*a *sinc(f);>> fill3(x,f,y,'b');>> grid on>> ylabel('f');>> zlabel('y');>> xlabel('x');>> axis(-0.6 0.6 -12 12 -0.6 0.6);>> set(gca,'YDir','reverse');>> line(0 0,-12 12,0 0);108%

69、 頻譜疊加>> subplot(3,1,3);>> a=1/sqrt(2);>> f=-11:0.001:-9;>> x=0.25*a*sinc(f+10);>> y=0.25*a *sinc(f+10);>> fill3(x,f,y,'b');>> hold on>> f=9:0.001:11;>> x=0.25*a*sinc(f-10);>> y=-0.25*a *sinc(f-10);>> fill3(x,f,y,'b');>> hold on>