% 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');