qam_ber.m

Hierarchical－QAM的程序

% This is the main program.

priority_vector1 = [9 4 2 1];
priority_vector2 = [9 4 2 1];

alpha1 = priority_vector1*(priority_vector1)';
alpha2 = priority_vector2*(priority_vector2)';

size_of_in_phase = length(priority_vector1);
size_of_q_phase = length(priority_vector2);

err_vector1=zeros(30,size_of_in_phase);
err_vector2=zeros(30,size_of_q_phase);     

for snr=-5:1:24
	snr_no=10^(snr/10);
	noise_qty = sqrt( (alpha1+alpha2) /(snr_no));
	for i=1:size_of_in_phase
		err_vector1(snr+6,i)=recurse_other(i,priority_vector1,noise_qty);
	end;
end;

for snr=-5:1:24
        snr_no=10^(snr/10);
        noise_qty = sqrt( (alpha1+alpha2) /(snr_no));
        for i=1:size_of_q_phase
                err_vector2(snr+6,i)=recurse_other(i,priority_vector2,noise_qty);
        end;
end;

figure;

%%% Plot the vectors now

for i=1:size_of_in_phase     
        k=-5:1:24;
        semilogy(k,err_vector1(2+(k+4),i));
        hold on;
end;


for i=1:size_of_q_phase
	k=-5:1:24;
	semilogy(k,err_vector2(1+(k+5),i));
	hold on;
end;
grid;
hold off;
xlabel('Carrier to noise ratio E_{s}/N_{0} [dB]');
ylabel('Average bit error probability');