?? psk_awgn_phase.m
字號:
% 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');
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -