求相干法的DQPSK调制解调的Matlab仿真程序

kylebit

毕设中涉及DQPSK调制解调的Matlab仿真程序，必须要用相干法做，其中有提取载波和码元定时之类，不知哪位达人编过此类程序？请赐教

happy

DQPSk调制过程 解调过程 系统误码率 matlab程序
提 供 者:         海涛

1. N_Trials=1000;
   
2. N_number=100;
   
3. N_snr=10;
   
4. Q=16;
   
5. 
   
6. BER_MATRIX=[];
   
7. SER_MATRIX=[];
   
8. for trials=1:N_Trials
   
9.     trials
   
10. noise=randn(1,N_number/2)+j.*randn(1,N_number/2);
    
11. s10=round(rand(1,N_number));
    
12. s=(s10*2-1)/sqrt(2);
    
13. sreal=s(1:2:N_number);
    
14. simage=s(2:2:N_number);
    
15. sc=sreal+j.*simage;
    
16. sgma=1;
    
17. Error_ber=[];
    
18. Error_ser=[];
    
19.     for snr_db=0:1:N_snr
    
20.         snr=10.^(snr_db./10);
    
21. 
    
22.         N0=2*sgma.^2;
    
23.         Eb=snr*N0;
    
24.         Es=2*Eb;
    
25.         yy=sqrt(Es).*sc+noise;
    
26.         y_real=sign(real(yy))./sqrt(2);
    
27.         y_image=sign(imag(yy))./sqrt(2);
    
28.         
    
29.         d_bit=[];
    
30.         for k=1:length(y_real)
    
31.             d_bit=[d_bit,[y_real(k),y_image(k)]];
    
32.         end
    
33.         d_symbol=y_real+j.*y_image;
    
34.         
    
35.           dif_bit=s-d_bit;  
    
36.           dif_symbol=sc-d_symbol;
    
37.          
    
38.           ber_snr=0;
    
39.           for k=1:N_number
    
40.           if dif_bit(k)~=0
    
41.               ber_snr=ber_snr+1;
    
42.           end
    
43.       end
    
44.       Error_ber=[Error_ber,ber_snr];
    
45.       
    
46.       ser_snr=0;
    
47.       for k=1:N_number/2
    
48.           if dif_symbol(k)~=0
    
49.               ser_snr=ser_snr+1;
    
50.           end
    
51.       end
    
52.       Error_ser=[Error_ser,ser_snr];
    
53.       
    
54.         %s_e=sign(y_real);
    
55.         %s_e10=(s_e+1)./2;
    
56.         %Error_snr=sum(abs(s10-s_e10));
    
57.         %Error_v=[Error_v,Error_snr];
    
58.     end
    
59.    BER_MATRIX=[BER_MATRIX;Error_ber] ;
    
60.    SER_MATRIX=[SER_MATRIX;Error_ser];
    
61. end
    
62. 
    
63.    BER_sum=mean(BER_MATRIX);
    
64.    BER=BER_sum./N_number;
    
65.    SER_sum=mean(SER_MATRIX);
    
66.    SER=SER_sum./(N_number./2);
    
67.    
    
68. BER_T=[];
    
69. SER_T=[];
    
70.    for snr_db=0:1:N_snr
    
71.       snr=10.^(snr_db./10);
    
72.        BER_THEROY=Qfunct(sqrt(2.*snr));
    
73.        SER_THEROY=1-(1-0.5.*erfc(sqrt(snr))).^2;
    
74.        BER_T=[BER_T,BER_THEROY];
    
75.        SER_T=[SER_T,SER_THEROY];
    
76.    end
    
77.    i=0:1:10;
    
78.    semilogy(i,BER,'-r',i,BER_T ,'*g',i,SER,'.c',i,SER_T,':k');

复制代码

1. 
   
2. function[y]=Qfunct(x);
   
3. y=(1/2)*erfc(x/sqrt(2));

复制代码

liyanping2007

请问一下谁有pi/4-dqpsk的调制与解调的MATLAB仿真程序，谢拉！

dwardancer

pi/4 qpsk 不就是dqpsk吗？
谁有这个仿真，要求有频谱图，谢谢啦:@Q

200810100659

专业的happy

Happy99

2F提供的程序, 专业不懂, 但感觉好像fortran似的:@)
直觉应该可以优化, 看谁愿意练习下