16QAM调制解调

clc clear all close all a=rand(1,400)>0.5; Modulated_Data=shiliuQam_modulation(a); x=Modulated_Data; I_real=real(x); %将信号分成I Q两路 Q_imag=imag(x); %%%%%%%%%%%%%%%%%I路%%%%%%%%%%%%%% I_0=find(I_real>2); I_real(I_0)=ones(size(I_0))*2; %%转换成四电平的形式 I_1=find(I_real>0&I_real<2); I_real(I_1)=ones(size(I_1))*3; I_2=find(I_real>-2&I_real<0); I_real(I_2)=ones(size(I_2))*1; I_3=find(I_real<-2); I_real(I_3)=zeros(size(I_3)); %%%%%%%%%%%%% 将四电平转换成二电平%%%%%%%% Si=[0 0;0 1;1 0;1 1]'; N_I=length(I_real); Bit_I=zeros(2,length(I_real)); for i=1:N_I switch I_real(i) case 0, Bit_I(1:2,i)=Si(1:2,1); case 1, Bit_I(1:2,i)=Si(1:2,2); case 2, Bit_I(1:2,i)=Si(1:2,3); case 3, Bit_I(1:2,i)=Si(1:2,4); end end %%%%%%%%% end %%%%%%%% %%%%%%% Q路 %%%%%%%%%%%%%%%% Q_0=find(Q_imag>2); Q_imag(Q_0)=ones(size(Q_0))*2; Q_1=find(Q_imag>0&Q_imag<2); Q_imag(Q_1)=ones(size(Q_1))*3; Q_2=find(Q_imag>-2&Q_imag<0); Q_imag(Q_2)=ones(size(Q_2))*1; Q_3=find(Q_imag<-2); Q_imag(Q_3)=zeros(size(Q_3)); %%%%%%%%%%%%% 将四电平转换成二电平%%%%%%%% Si=[0 0;0 1;1 0;1 1]'; N_Q=length(Q_imag); Bit_Q=zeros(2,length(Q_imag)); for i=1:N_Q switch Q_imag(i) case 0, Bit_Q(1:2,i)=Si(1:2,1); case 1, Bit_Q(1:2,i)=Si(1:2,2); case 2, Bit_Q(1:2,i)=Si(1:2,3); case 3, Bit_Q(1:2,i)=Si(1:2,4); end end %%%%%%%%%%%% end %%%%%%%%%%%%%%% %%%%%%%%%%%%%%%将 I Q 两路信号比特合成一路信号%%%%%%%%%%%%%%%%%%% Demodulated_Data=[Bit_I ; Bit_Q]; Demodulated_Data=reshape(Demodulated_Data,1,length(x)*4); m=sum(Demodulated_Data~=a);