www.pudn.com > MPSK_MFSK.rar > mfsk.m

initmpsk;  
 for m=1:11; 
     counterr=0; 
     for k=1:3000 
        a1(k)=sqrt(y(m)/2)*randn; 
        b1(k)=sqrt(y(m)/2)*randn; 
        a2(k)=sqrt(y(m)/2)*randn; 
        b2(k)=sqrt(y(m)/2)*randn; 
        a3(k)=sqrt(y(m)/2)*randn; 
        b3(k)=sqrt(y(m)/2)*randn; 
        a4(k)=sqrt(y(m)/2)*randn; 
        b4(k)=sqrt(y(m)/2)*randn; 
                for j=1:3000 
            ra=rand; 
            phi=fix(ra*2*pi); 
        rmc=sqrt(2)*cos(phi)+a1(k); 
        rms=sqrt(2)*sin(phi)+b1(k); 
                 
          if  rmc.^2+rms.^2 < a2(k).^2 + b2(k).^2 | rmc.^2+rms.^2 < a3(k).^2 + b3(k).^2 | rmc.^2+rms.^2 < a4(k).^2 + b4(k).^2 
              counterr=counterr+1; 
          end 
        end 
        end 
semilogy(x(m),counterr/9000000,'r * -') 
hold on; 
 end 
 
 for z=0:0.1:10 
semilogy(z,3/2*exp(-10.^(z/10))) 
legend ('the result of simulation') 
 end 
xlabel('Eb/N0(dB)');ylabel('BER');