1、rec_qam16=QamMod(bit_tx,16); %方形16QAM调制star_qam16=SrarQamMod(bit_tx); %星形16QAM调制base_rec=base_shape(fd,fs,f,rec_qam16); %基带成型滤波 base_star=base_shape(fd,fs,f,star_qam16);for i=1:length(SNR) %信噪比从1dB到12dB计算误码率 SNR_=i %方形映射16QAM rf_rec_qam16=CarrierMod(fc,f,base_rec); %载波调制 rf_rec_qam16_n=awgn(rf_rec_q
2、am16,SNR(i),measured); %加噪声 rec_qam16_rx base_rec_rx=CarrierDemod(fd,fs,fc,f,rf_rec_qam16_n); %载波解调 bit_rec_rx=QamDemod(rec_qam16_rx,16); %MQAM解调 num_qam16,perr_qam16_rec(i)=biterr(bit_tx,bit_rec_rx);%误码率 qam16_data_rec(i,:)=rec_qam16_rx; %scatterplot(rec_qam16_rx); %星形映射16QAM rf_star_qam16=CarrierM
3、od(fc,f,base_star); rf_star_qam16_n=awgn(rf_star_qam16,SNR(i), star_qam16_rx base_star_rx=CarrierDemod(fd,fs,fc,f,rf_star_qam16_n); bit_star_rx=StarQamDemod(star_qam16_rx); num_qam16,perr_qam16_star(i)=biterr(bit_tx,bit_star_rx); qam16_data_star(i,:)=star_qam16_rx; %scatterplot(star_qam16_rx);end% 理
4、论误码率计算SNRtheo=0:0.1:length(SNR);length(SNRtheo) SNRdec=10.(SNRtheo(i)/10); theo_perr_qam16(i)=(3/8)*erfc(sqrt(SNRdec*2/5);% 基带波形N=200;n=1:N;t=1:N*10;%发送端波形figure;stem(n,bit_tx(n); title(发送序列 %发送序列figure(Name,发送端基带信号subplot(411);plot(t,real(base_rec(t);方形映射-Q路subplot(412);plot(t,imag(base_rec(t);方形映射
5、-I路subplot(413);plot(t,real(base_star(t);星形映射-Q路subplot(414);plot(t,imag(base_star(t);星形映射-I路%接收端波形subplot(211);stem(n,bit_rec_rx(n);方形接收序列subplot(212);stem(n,bit_star_rx(n);星形接收序列接收端基带信号plot(t,real(base_rec_rx(t);plot(t,imag(base_rec_rx(t);plot(t,real(base_star_rx(t);plot(t,imag(base_star_rx(t);%基带
6、眼图N1=20000;Tn=f/fd;eye_rex=base_rec(1:N1);eyediagram(eye_rex,Tn*4,Tn);方形基带眼图eye_star=base_star(1:eyediagram(eye_star,Tn*4,Tn);星形基带眼图% 接收端星座图%scatterplot(qam16_data_rec(12,:);方形16QAM接收端星座图for i=3:3:12 subplot(2,2,i/3); plot(real(qam16_data_rec(i,:),imag(qam16_data_rec(i,:),. xmax=5;axis(-xmax xmax -x
7、max xmax) title(Snr=,num2str(SNR(i), dB);星形16QAM接收端星座图 plot(real(qam16_data_star(i,:),imag(qam16_data_star(i,: xmax=3; grid on;% 功率谱密度f_plot;% 误码率分析%16QAM误码率曲线16QAM误码性能对比% semilogy(SNRtheo,theo_perr_qam16);% hold on;semilogy(SNR,perr_qam16_rec,*hold on;semilogy(SNR,perr_qam16_star,oxlabel(SNR in dBy
8、label(Prb of Errlegend(方形16QAM星形16QAMtitle(%-% 基带成型function base_info=base_shape(fd,fs,f,seq_16QAM)%平方根升余弦滤波器,滚降系数0.5,延迟3个采样点flt=rcosine(fd,fs,sqrt,0.5); %I路和Q路seq_Q=real(seq_16QAM);seq_I=imag(seq_16QAM);%增采样R=fs/fd;up_seq_Q=upsample(seq_Q,R);up_seq_I=upsample(seq_I,R);%升余弦调制rcos_Q=conv(up_seq_Q,flt
9、);rcos_I=conv(up_seq_I,flt);%提升rcos_Q_up=interp(rcos_Q,f/fs);rcos_I_up=interp(rcos_I,f/fs);base_info=rcos_Q_up+j*rcos_I_up;%-% 载波解调function data_rx base_rx=CarrierDemod(fd,fs,fc,f,receive)% 分两路乘正交高频载波rc_length=length(receive);t=0:rc_length-1;rc_Q=receive .* sin(2*pi*fc*t/f);rc_I=receive .* cos(2*pi*
10、fc*t/f);%减采样后根升余弦匹配滤波,注意对齐采样点down_Q=downsample(0 rc_Q,f/fs);down_I=downsample(rc_I,f/fs);low_Q_rcos=conv(down_Q,flt);low_I_rcos=conv(down_I,flt);base_rx=low_Q_rcos(1:length(low_I_rcos)+j*low_I_rcos;%两次根升余弦滤波延迟,定位初始信号位置%delay+1:end-delay-1delay=3*R*2;rc_Q_seq=(downsample(low_Q_rcos(delay+1:end-delay-
11、1),R);rc_I_seq=(downsample(low_I_rcos(delay+1: %类型转换data_rx=rc_Q_seq+j*rc_I_seq;%-% 载波调制function transmit=CarrierMod(fc,f,base)%载波调制length(base)-1;high_freq_Q=real(base) .* sin(2*pi*fc*t/f);high_freq_I=imag(base) .* cos(2*pi*fc*t/f);transmit=high_freq_Q+high_freq_I;%-% QAM解调程序,将MQAM码元还原为二进制数据% 入口参数:
12、data_QAM:QAM码元数据% M_QAM:MQAM中M的大小% 出口参数:二进制数据比特流%-function data_out=QamDemod(data_QAM,M_QAM)%每个M进制码元的bit数data_len=length(data_QAM);%码元长度%-%QAM信号放缩至与发送端相同比例data_temp=data_QAM(find(real(data_QAM0);aver=mean(real(data_temp);data_Qam_temp=data_QAM/aver*(2(k/2-1);%平移到第一象限data_Qam_temp2=(data_Qam_temp+(2(
13、k/2)-1)*(1+j)/2;%-%将实部虚部分别映射为二进制数据%实部data_real=round(real(data_Qam_temp2);%实部判决data_real(find(data_real(2(k/2)-1)=2(k/2)-1;%大于2(k/2)-1的数据判定为2(k/2)-1data_real(find(data_realdata_imag(find(data_imag2)=2;%幅度大于2的判定为2am(find(am1.5)=1;%相位判定ang=angle(qam_rev)/pi*180;%将-180度到0度变为180度到360度ang(find(ang0)=ang(find(ang7)=0;data_dec=8*(round(am)-1)+ang_n;%data_dec(find(data_dec15)=15;%大于15的判定为15%data_dec(find(data_dec%小于0的判定为0data_bit_rec=dec2bin(data_dec);data_bit=reshape(data_bit_rec,1,length(qam_rev)*4)-
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1