psk_rec_telek.m

Hierarchical－QAM的程序

% This does the plotting for 2/4/8-PSK

input_angle_vector = [pi/2 pi/3 pi/15 pi/40 ]'

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

size_of_constellation = length(angle_vector);
err_vector=zeros(30,size_of_constellation);
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)=recurse_other(i,angle_vector,noise_qty);
	end;
end;

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