tx.m

来自「卫星信道中的调制识别的matlab仿真分析」· M 代码 · 共 28 行

function TX(y)
 z=hilbert(y);
   for i=1:5000
       ph1(i)=angle(z(i));
   end
   subplot(312);plot(ph1);
   for j=1:4900
      ph1(j) =abs(angle(z(j)*conj(z(j+100))));%求相位差
   end
   subplot(313);plot(ph1(1:4900));
   m=100;
   for i=0:48
       for j=1:m
       ph(i+1,j)=ph1(i*m+j);    
       end
   end
     r=zeros(1,3);
             for i=1:49;
           v(i)=mean(ph(i,:));%求出每个码元内的均值，即台阶的值
           if(v(i)>1.40&v(i)<=1.70)   r(1)=1;end
           if(v(i)>2.60&v(i)<=3.20)   r(2)=1;end
           if(v(i)>4.60&v(i)<=4.90)   r(3)=1;end      
             end
             r;
             canshu4=sum(r);
     % if(canshu4==1) type=menu('输入信号是','2PSK信号');
     % else type=menu('输入信号是','QPSK信号');
     % end