深入浅出通信原理二

1、连载 51：利用旋转向量理解正交解调（一）3连载 52：利用旋转向量理解正交解调（二）10连载 53：利用旋转向量理解正交解调（三）12连载 54：PSK/QAM 调制仅仅是指部分吗？13连载 55：调制解调与级数展开的关系14连载 56：利用求复系数的方法实现解调15连载 57：如何求复系数？16连载 58：OFDM 与级数展开18连载 59：如何求系数？20连载 60：周期方波信号的复系数22连载 61：sinc 函数23连载 62：周期方波信号的频谱24连载 63：周期矩形波的频谱26连载 64：周期矩形波频谱对比（一）30连载 65：周期矩形波的频谱对比（二）32连载 66：非周期矩形

2、信号的频谱35连载 67：连续型频谱36连载 68：周期矩形波的连续谱38连载 69：周期矩形波的连续谱和离散谱对比44连载 70：非周期矩形信号的连续谱46连载 71：非周期信号的连续谱（一）48连载 72：非周期信号的连续谱（二）491连载 73：非周期信号的连续谱（三）50连载 74：非周期信号的连续谱（四）52连载 75：已知频谱求非周期信号54连载 76：变换56连载 77：对调制过程进行深入的频谱分析57连载 78：信号与复指数信号相乘的频谱变化59连载 79：矩形脉冲信号调制余弦载波（时域）61连载 80：矩形脉冲信号调制余弦载波（频域）64连载 81：矩形脉冲调制余弦载波（三）

3、66连载 82：矩形脉冲信号调制余弦载波（四）68连载 83：正负脉冲的幅度谱和相位谱70连载 84：采用对数坐标的矩形脉冲幅度谱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：信号的频谱分析方法可否统一？113连载 97：冲激函数114连载 98：周期信号的变换115连载 99：复指数信号的变换117连载 100：余弦信号的变换119连载 47：利用旋转向量理解 BPSK 解调（三）34连载 48：用复数运算实现 BPSK 调制和解调5连载 49：利用实数运算实现 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：调制解调与系级数展开的关14连载 56：

6、利用求复解调系数的方法实现如何求复系数呢？15连载 57：如何求复系数？1617连载 58：OFDM 与级数展开1819连载 59：如何求系数？2021连载 60：周期方波信号的复系数22连载 61：sinc 函数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：连续型频谱3637连载 68：周期矩形波的连续谱383940414243连载 69：周期矩形波的连续谱和离散谱对比>> 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：非周期矩形信号的连续谱46>> x=-8:0.001:8;>> y=sinc(x);>> plot(x,y)>> grid on47连载 71：非周期信号的连续谱（一）48连载 72：非周期信号的连续谱（二）49连载 73：非周期信号的连续谱（三）5051连载 74：非周期信号的连续谱（四）5253连载 75：已知频谱求非周期信号5455连载 76：变换56连载

16、 77：对调制过程进行深入的频谱分析5758连载 78：信号与复指数信号相乘的频谱变化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：采用对数坐标的矩形脉冲幅度谱%正脉冲幅度频谱%>> 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 调制信号的频谱（二）>> %数据为 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);>> %数据为 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>> %数据为 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);>> %数据为 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>