ip_09_04.m

是通信类的MATLABA仿真抗正弦噪声干扰能力

% MATLAB script for Illustrative Problem 9.4.
echo on
Lc=20;			  		% number of chips per bit
A1=3;			  		% amplitude of the first sinusoidal interference
A2=7;			  		% amplitude of the second sinusoidal interference
A3=12;			  		% amplitude of the third sinusoidal interference
A4=0;			  		% fourth case: no interference
w0=1;			  		% frequency of the sinusoidal interference in radians
SNRindB=0:2:30;
for i=1:length(SNRindB),
  % measured error rates
  smld_err_prb1(i)=ss_Pe94(SNRindB(i),Lc,A1,w0);
  smld_err_prb2(i)=ss_Pe94(SNRindB(i),Lc,A2,w0);
  smld_err_prb3(i)=ss_Pe94(SNRindB(i),Lc,A3,w0);
  echo off ;
end;
echo on ;
SNRindB4=0:2:30;
for i=1:length(SNRindB4),
  % measured error rate when there is no interference
  smld_err_prb4(i)=ss_Pe94(SNRindB4(i),Lc,A4,w0); 
echo off ;
x=SNRindB;
y1=smld_err_prb1;
y2=smld_err_prb2;
y3=smld_err_prb3;
plot(x,y1,x,y2,x,y3,x,y4)
end;    	
echo on ;		          
% Plotting commands follow.