psk_awgn_phase.m

来自「Hierarchical－QAM的程序」· M 代码 · 共 52 行

% This does the plotting for 2/4/8/16-PSK
%div_const=1.1;
input_angle_vector = [pi/2 pi/4]'
phase_error=pi/100;

angle_vector=input_angle_vector(1:length(input_angle_vector)-1,1);
angle_vector(1)=pi/2-angle_vector(1);
size_of_constellation = length(angle_vector);

angle_vector1=angle_vector;
angle_vector1(1)=angle_vector(1)+phase_error;
angle_vector1;
 
angle_vector2=angle_vector;
angle_vector2(1)=angle_vector(1)-phase_error; 
angle_vector2;

err_vector=zeros(30,size_of_constellation+1);
for snr=-5:1:24
    snr;
	snr_no=10^(snr/10);
	noise_qty = sqrt(1/(snr_no));
	for i=1:size_of_constellation+1
        i;
	
err_vector(snr+6,i)=0.5*(recurse_other(i,angle_vector1,angle_vector,phase_error,snr_no)+recurse_other(i,angle_vector2,angle_vector,phase_error,snr_no));
	end;
end;


fid=fopen('ber_16psk_3_15_40.m', 'a');
fprintf(fid, '\n\n err_vector = [ ');
fprintf(fid, '%12.8f %12.8f %12.8f %12.8f; \n', err_vector');
fprintf(fid, ']; \n\n' );
fprintf(fid, 'i=-5:1:24; \n');
fprintf(fid, 'semilogy(i, err_vector); ');
fclose(fid);

k=-5:1:24;
semilogy(k,err_vector);

%for i=1:size_of_constellation+1
%	k=-5:1:24;
%	semilogy(k,err_vector(k+6,i));
%	hold on;
%end;

grid;
hold off;
xlabel('Carrier to noise ratio E_{s}/N_{0} [dB]');
ylabel('Average bit error probability');