软件无线电发射机的实现与仿真（三）

1、软件无线电发射机的实现与仿真（三）     由于频率与相位有一定的关系，为便于分析，可将式(4-56)改写为=A(n)cos n+                      (4-57)式中，表示载波的角频率。所以 =A(n)cos cos( n)-A(n)sin sin( n)= cos( n)- sin( n)   

2、0;              (4-58)式中= A(n)cos                          (4-59)= A(n)sin       

3、0;                   (4-60)这就是我们希望获得的同相和正交两个分量，根据  、，就可以对各种调制样式进行解调，三大类解调的算法如下:调幅(AM)解调:A(n)=                  

4、0;         (4-61)调相(PM)解调:=                           (4-62)          (4-63)调频(FM)解调&

5、#160;         (4-64)在利用相位差分计算瞬时频率，即 = - 时，由于计算 要进行除法和反正切运算，这对于非专用数字信号处理器来说是较复杂的，在用软件实现时也可以用下面的方法来计算瞬时频率：= =                (4-65)对于调频信号，其振幅近似恒定，设 =1，则      (4-6

6、6)式(4-66)就是利用 、 直接计算的近似公式。这种方法只有乘减运算，计算比较简便。最后得到的软件无线电数字正交解调的通用模型，如图4.2所示。 SHAPE  * MERGEFORMAT                        图4.2 数字正交解调的通用模型信号表达式:s(n)=A(n)cos(  )   

7、      ( (4-67)式中， ；为调制信号； 为载波初始相位。    对信号进行正交分解，得到同相和正交分量:同相分量: =A(n)cos                            (4-68)正交分量:= A(n)cos 

8、;                            (4-69)对同相和正交分量平方之和开方:= +m(n)                  

9、   (4-70)减去直流分量就可得到调制信号m(n)。这种方法具有着较强的抗载频适配能力，即本地载波与信号载波之间允许一定得频率偏差。当由于传输信道或其他一些原因而造成本地载波与信号的载频之间存在频差和相差时，同相分量和正交分量可表示为:           (4-71)                 (4-72)

10、式中， = - ； = - ； 、 表示差频和差相可以是常量也可以是随机变量。 为本地载波的角频率:  是本地载波的初始相位。对同相与正交分量平方之和开平方得:= +m(n)                (4-70)    所以，AM信号用正交解调算法解调时，不要求载频严格的同频同相。从以上分析过程中可知，理论上失配可以任意大，但由于失配时，同相和正交分量相当于调制在以失配频率为载频的载波上，严重失配时，

11、信号会超出数字信道而发生失真。    2. DSB解调信号表达式:S(n)=m(n)cos                        (4-74)对信号进行正交分解得:同相分量:=m(n)            &#

12、160;                (4-75)正交分量:=0                                

13、 (4-76)    解调时要求本地载频与信号载频同频同相，此时，同相分量输出就是解调信号。同频同相本地载频的提取，可以利用数字科斯塔斯环获得。数字科斯塔斯环既可以用软件实现也可以利用专门的数字信号处理硬件来实现。    3. SSB解调信号表达式:s(n)=m(n)cos sin                (4-77)对信号正交分解得:同相分量:=m(n)  &

14、#160;                          (4-78)正交分量:= (n)                     &#

15、160;    (4-79)4. FM解调信号表达式:s(n)= cos +               (4-80)式中，k为比例因子， 为常数。对信号进行正交分解得;同相分量:= cos                      (4-

16、81)正交分量:= sin                   (4-82)对正交与同相分量之比值反正切运算:=arctg                          = &

17、#160;                    (4-83)然后，求相位差分，即可求得调制信号:- =m(n)                         

18、; (4-84)    为了讨论方便，这里及以下对比例因子k及常数 忽略。FM信号用正交解调方法解调时，也具有较强的抗载频失配(指失配差频和差相是常量，非随机变量)能力，本地载波与信号的载波存在频差和相差时，同相分量和正交分量可表示为:= cos +                 (4-85)= sin +          

19、   (4-86)同样对正交与同相分量之比值反正切及差分运算，就可得到调制信号:arctg -arctg      = + + - + +          (4-87)     = +m(n)当载波失配差频和差相是常量时，解调输出只不过增加了一个直流分量 ，减去直流分量 就可得到调制信号m(n)。信号表达式:s(n)= cos( + )      

20、0;             (4-88)式中，为输入码元，且=0、1；g(n一m)是幅度为1，宽度为码元传输速率倒数的矩形脉冲门函数。    ASK的解调算法与AM解调一样:对信号进行正交分解，得同相和正交分量:同相分量:= cos( )                 

21、60;   (4-89)正交分量:= sin( )                      (4-90)对同相与正交分量平方之和开方:A(n)= =           (4-91)计算A(n)后，再对A(n)进行抽样判决，就可恢复出调制码元信号。 

22、0;  ASK的正交解调性能与AM一样，具有较强的抗载频失配能力。    2. MASK解调信号表达式:s(n)=  cos( + )                   (4-92)式中， 为输入码元，且 。解调方法与ASK一样，对信号进行正交分解，得同相和正交分量:同相分量:= cos( )      &

23、#160;             (4-93)正交分量:= sin( )                   (4-94)按照式(4-91)计算瞬时幅度A(n):A(n)=          

24、60;                    (4-95)计算出A(n)后，再进行抽样多电平幅度判决，就可恢复出调制码元信号。    MASK解调性能与ASK一样，具有较强的抗载频失配能力。    3. FSK解调信号表达式:s(n)=  cos        

25、;         (4-96)式中，为载波角频率间隔， 为输入的码元，= +1,-1 。    FSK解调类似于FM解调，对信号进行正交分解，得同相和正交分量:同相分量:= cos( n)                   (4-97)正交分量:= sin( n)   

26、60;                   (4-98)按照式(4-64)计算瞬时频率f(n):n)= arctg -arctg             =             

27、;        (4-99)在计算出瞬时频率f(n)后，对f(n)经抽样门限判决，即可恢复出传输的数据。    4. MFSK解调信号表达式:s(n)= cos( + )n            (4-100)式中， 为输入码元，且 。    MFSK解调类似于FSK解调，对信号进行正交分解，得同相和正交分量:同相分量:= cos( n) 

28、60;               (4-101)正交分量:= sin( n)                       (4-102)按照式(4-99)计算瞬时频率f(n):f(n)=    &#

29、160;                        (4-103)在计算瞬时频率f(n)后，对f (n)抽样多电平门限判决，即可恢复出数据。    5. MSK解调    信号表达式:s(n)=          (

30、4-104)式中，T为码元持续时间; 为输入码元，且 =+1,-1。=          是为保证相位连续而加入的相位常数。    MSK信号的解调同FM，对信号进行正交分解，得同相和正交分量:同相分量:= cos                    (4-105)正交分量:= sin  

31、;             (4-106)按照式(4-64)计算瞬时频率f(n):f(n)= arctg -arctg                =                

32、60;       (4-107)在计算出瞬时频率f(n)后，对f(n)抽样判决，即可恢复出码元。    6. GMSK解调    GMSK信号与MSK信号相比，仅对输入数据多加了一个预调制滤波器。因此，可按MSK信号那样解调后，再经一个滤波器 =  ( 为预调制滤波器频率响应)，即可求得码元。    7. SFSK解调信号表达式:s(n)= cos       &#

33、160; (4-108)    SFSK信号解调方法同MSK解调，对信号进行正交分解后，按照式（4-107）计算瞬时频率。在计算出瞬时频率f(n)，对f(n)抽样判决，即可恢复出码元。    8. PSK解调信号表达式:s(n)= cos +                   (4-109)式中， = , 。对信号进行正交分解后，得同相和正交分量:同相分量:= co

34、s( )                    (4-110)正交分量:= sin( )                  (4-111)按照式(4-62)求得瞬时相位:=     

35、0;                     (4-112)    在计算出瞬时相位后，对 抽样判决，即可恢复数据。在解调时需要本地载波与信号载波严格的同频同相，同频同相可由数字科斯塔斯环获得。    9. MPSK解调信号表达式:s(n)= cos +  (4-113)式中， ,。    MPSK

36、信号解调方法同PSK。在计算出瞬时相位后，对抽样进行多电平门限判决，即可恢复出码元数据。    10. QPSK解调信号表达式:s(n)= cos( )+ sin( )  (4-114)式中 ，  为双极性数据。    对信号进行正交分解，得到同相和正交分量:同相分量:=                     (4-

37、115)正交分量:=                        (4-116)由信号形式可知，I, Q分量即为恢复出的并行数据，经抽样判决，恢复出码元数据后，在并串变换，就可恢复出串行码元数据。    11.QAM解调信号表达式:s(n)= cos( )+ sin( )    (4-117)式中， , =

38、 。对信号进行正交分解，得到同相和正交分量:同相分量:=                        (4-118)正交分量:=                    

39、    (4-116)    对同相、正交分量两路信号进行抽样判决，即可恢复出并行数据，经并串变换后可得所传输的数据。                                 (5-1)式中，为第i信道的归一化中心角频率，

40、I为数据内插率。    基带信号经内插低通滤波，再与复本振 相乘，可实现将第i个信道的数字谱搬移到频带的目的。经过复本振 后，信号变为复信号，故I路合成信号需取实部后再输出。为使I个采样率为 的基带信号能够压缩在实信号所表示的频谱范围内传输，内插因子取为2I。其实现结构如图5.1所示。 SHAPE  * MERGEFORMAT  图5.1实信号输出信道化发射机的直接实现  图中每个低通滤波器的带宽均不大于，并且对应的原型理想低通滤波器的频率响应为      

41、0;               (5-2) ，第i信道的归一化中心角频率公式可重写为                             (5-3)式中，  

42、60;  当    , 这是不允许的。因此，后面的 个信道的计算公式为                             (5-4)式中，   需要指出，由式(5-3 )得到的实信号信道存在着对应的镜频 ，并且信道总数受数据内插倍数I的限制。图5.2为对应4个实信道的频谱分

43、配图。 SHAPE  * MERGEFORMAT                   图5.2 实信号的信道划分示意图         注意实信号的频谱应为正值，由式( 5-3 ) , ( 5-4 )可推出真实信道的中心频率为         

44、60;  (5-5)        (5-6)进一步由式(5-5 ) , ( 5-6)容易求出相邻信道中心频率距离为.5.2基于多相滤波器的实信号信道化发射机建模 由图5.2可得: 并定义：            代入式(5-6)可得：              把 代入式(5-8)可得:

45、    (5-9)定义： = =DFT 代入式(5-9 )可得: 则有：                            =                

46、;        (5-12)最后得：                  y(n)=                       

47、       (5-13)式中， =MOD（n/I）,MOD表示取余数。整个实现过程如图5.3所示。 SHAPE  * MERGEFORMAT             图5.3   实信号信道化软件无线电发射机数学模型对于基于多相滤波器的实信号信道化发射机模型的几点说明:(1)多相滤波器的设计步骤实现多相滤波器设计的步骤是:(1)根据原型理想低通滤波器的频率响应确定所需要的滤波器类型和阶

48、数N ;  (2)求出对应的冲击响应h(n) ; (3 )由下式确定多相滤波器：        m=0,1,2,I-1                      (5-14)若根据频率响应求得的滤波器阶数N不是I的整数倍，则需要进行反向设计，即设定滤波器的阶数N为I的整数倍后再重新计算各阶系数。利用Matla

49、b中的REMEZORD函数可以方便求出采用最佳逼近最大最小准则算法所需的原型滤波器阶数N。(2) DFT可以由快速算法FFT来完成。第六章 软件无线电发射机系统仿真 本章将构建一个基于多相滤波器的实信号信道化发射机仿真系统并用Matlab软件进行仿真，以验证其可行性。6.1基于多相滤波器的信道化发射机系统仿真 在基于多相滤波器的实信号信道化发射机仿真设计中，信道数、内插倍数和信道频率的划分是密切相关的，因此，仿真设计时进行了综合考虑，且用快速傅立叶变换对信号进行处理，不断提高系统工作效率。仿真采用Matlab软件的M文件来实现。   调制模式:AM  S

50、HAPE  * MERGEFORMAT                图6.1  多信道信道化发射机仿真结构图中，I=8为输入信号对应的信道号，Y(n)为输出信号。信道化发射机主要用在对某一带宽内的所有信道进行发射的场合，所以其信道数应很大，但考虑计算机的实际运算能力，信道数不能设置过大，而且在系统仿真中信道数量的增加只会增加计算负担，对于验证系统可行性没有多大贡献。由于信道化滤波器的最先一步运算为FFT2变换，所以信道数最好为2的整数

51、次幂，这样可以提高工作效率。基于以上考虑，信道数设置为8。(2)调制波形    语音信号虽具有形象直观的优点，但它的频谱和时域波形都比较杂乱，不能清晰地反映数字信道的问题所在，所以本节不选择语音信号。AM调制对于信道衰减敏感，本节选择一些常见波形作为调制波形，可以很容易判断发射机的性能。6.2实验结果与分析 整个仿真程序（Matlab程序）如下。a=1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0;f=1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0;I=8;n1=200;Fs=25.0;fs=Fs*2*I;kf=12.5;for k=

52、1:I  for r=1:(n1+I)     m(k,r)=a(k)*(1.0+0.5*cos(2*pi*f(k)/Fs*(r-1);  endendn0,f0,m0,w=remezord(8,12.5,1 0,0.001 0.001,fs);b=remez(287,f0,m0,w);figure(1)polt(20*lof10(abs(fft(b);grid;for r=1:36  for k=1:I    h(k,r)=b(r-1)*I+k);  endendfor r=1:(n

53、1+I)  for k=1:I     mk(k)=m(k,r);  end     mfft=fft(mk);     for k=1:Ix0(k,r)=mfft(k)*exp(j*pi/(2*I)*(k-1);  endendfor r=1:(n1+I)  for k=1:I     x00(k,(2*r-1)=x0(k,r);     x00(k,2*

54、r)=0;  endendfor k=1:I  for r=1:(n1+I)     x00k(r)=x00(k,r);  end  for r=1:36  hk(r)=h(k,r);  end  y0=conv(x00k,hk);  for r=1:n1     y0(r)=y0(r)*exp(j*pi/2*(r-1);  end  for r=1:n1     y(k,r)

55、=y0(r+36);%*(-1)(k-1)*exp(j*pi/I*(k-1);  endendfor k=1:Ifor n=1:n1*I  if mod(n-1),I)=0    y00(k,n)=y(k,(n-1)/I+1);  else y00(k,n)=0.0;  endendendfor n=I:(n1*I-I)   yout(n-I+1)=y00(1,n)+y00(2,n-1)+y00(3,n-2)+y00(4,n-3)+y00(5,n-4)+y00(6,n-5)+y00(7,n-6)+y00(8,n-7);endpoint=512;yy(1:point)=yout(101:(100+point);for n=1:point   yy(n)=(y(n)+0.001*randn)*