plot_snr.zip_16qam_SNR仿真_plot_snr_snr

clear; snr=1:1:11; error_theory=(1-(1-(2*(1-1/sqrt(16))*1/2*erfc(1/sqrt(2)*sqrt(3*4*10.^(snr/10)/(16-1))))).^2)/4; N=floor(1./error_theory)*100+100; N(find(N<5000))=5000; global p; for i=1:length(N); source=randsrc(1,N(i),[1,0;p,1-p]); [source1,source2]=Qam_modulation(source); sig_insert1=insert_value(source1,8); sig_insert2=insert_value(source2,8); [source1,source2]=rise_cos(sig_insert1,sig_insert2,0.25,2); [x1,x2]=generate_noise(source1',source2',snr(i)); sig_noise1=x1'; sig_noise2=x2'; [sig_noise1,sig_noise2]=rise_cos(sig_noise1,sig_noise2,0.25,2); [x1,x2]=pick_sig(sig_noise1,sig_noise2,8); sig_noise1=x1; sig_noise2=x2; signal=demodulate_sig(sig_noise1,sig_noise2); error_bit(i)=length(find(signal-source)~=0)/N(i); end; semilogy(snr,error_bit,'-b'); hold on semilogy(snr,error_theory,'-r');