matlab仿真程序.docxVIP

附录一 AWGN信道下的BPSK系统仿真的MATLAB程序 main.m echo on SNRindB1=0:1:10; for i=1:length(SNRindB1) % simulated error rate smld_err_prb(i)=smldPe_awgn(SNRindB1(i); echo off; end; % Plotting commands follow semilogy(SNRindB1,smld_err_prb,*); smldPe_awgn.m function p=smldPe_awgn(snr_in_dB) E=1; SNR=exp(snr_in_dB*log(10)/10);%signal to noise ratio sgma=E/sqrt(2*SNR); %sigma,standard deviation of noise N=214; % generation of the binary data source for i=1:N temp=rand; if(temp0.5) dsource(i)=0; else dsource(i)=1; end end; % detection ,and probability of error calculation numoferr=0; for i=1:N % matched filter outputs if(dsource(i)=0) r=-E+gngauss(sgma); % if the source output is 0 else r=E+gngauss(sgma); % if the source output is 1 end; % detector follows if(r0) decis=0; else decis=1; end; if(decis=dsource(i) numoferr=numoferr+1; end; end; p=numoferr/N; % probability of error estimate gngauss.m function gsrv1,gsrv2=gngauss(m,sgma) if nargin =0 m=0; sgma=1; elseif nargin=1 sgma=m; m=0; end u=rand; z=sgma*(sqrt(2*log(1/(1-u); u=rand; gsrv1=m+z*cos(2*pi*u); gsrv2=m+z*sin(2*pi*u); 附录二 RAYLEIGH信道的MATALB程序 smldPe_rayleigh.m function p=smldPe_rayleigh(snr_in_dB) E=1; SNR=exp(snr_in_dB*log(10)/10);%signal to noise ratio sgma=E/sqrt(2*SNR);%sigma,standard deviation of noise N=214; % generation of the binary data source for i=1:N temp=rand; if(temp0.5) dsource(i)=0; else dsource(i)=1; end end; % detection ,and probability of error calculation numoferr=0; for i=1:N % matched filter outputs if(dsource(i)=0) a=ray; r=-E*a+gngauss(sgma);% if the source output is 0 else a=ray; r=E*a+gngauss(sgma); % if the source output is 1 end; % detector follows if(r0) decis=0; else decis=1; end; if(decis=dsource(i) numoferr=numoferr+1; end; end; p=numoferr/N; % probability of error estimate ray.m function a=ray Xc=randn; Xs=randn; a=sqrt(Xc2+Xs2); 附录三 DOPPLER信道的MATALB程序 smldPe_dpl.m function p=smldPe_dpl(snr_in_dB) E=1; SNR=exp(snr_in_dB*log(10)/10); %signal to noise ratio sgma=E/sqrt(2*SNR); %sigma,standard deviation of noise N=214; fm=10; % generation of the binary data source for i=1:N temp=rand; if(temp0.5) dsource(i)=0; else dsource(i)=1; end end; % detection ,and probability of error calculation b=dop_fading(N,fm); numoferr=0; for i=1:N-1 % matched filter outputs if(dsource(i)=0) a=ray; r=-E*b(i)+gngauss(sgma);% if the source output is 0 else a=ray; r=E*b(i)+gngauss(sgma);% if the source output is 1 end; % detector follows if(r0) decis=0; else decis=1; end; if(decis=dsource(i) numoferr=numoferr+1; end; end; p=numoferr/N; % probability of error estimate dop_fading.m function y=dop_fading(N,fm) %N = 64; %fm = 10; %unit: 1Hz % generate Doppler power spectrum deltaf = 2*fm/(N-1); T = 1/deltaf; sf0 = 1.5/(pi*fm); for n = 1:(N-2)/2 sf(n) = 1.5/(pi*fm*sqrt(1-(n*deltaf/fm)2); end classicf = fliplr(sf),sf0,sf; % generate two normally distributed random variables gaussN_re1 = randn(1,(N-2)/2); gaussN_im1 = randn(1,(N-2)/2); gaussN_pos1 = gaussN_re1 + j*gaussN_im1; gaussN_neg1 = conj(gaussN_pos1); % gaussN1 = fliplr(gaussN_pos1),0,gaussN_pos1; gaussN1 = fliplr(gaussN_neg1),0,gaussN_pos1; % gaussN_pos2 = randn(1,(N-2)/2); % gaussN2 = fliplr(gaussN_pos2),0,gaussN_pos2; gaussN_re2 = randn(1,(N-2)/2); gaussN_im2 = randn(1,(N-2)/2); gaussN_pos2 = gaussN_re2 + j*gaussN_im2; gaussN_neg2 = conj(gaussN_pos2); gaussN2 = fliplr(gaussN_neg2),0,gaussN_pos2; % generating flat Rayleigh fading channel x = ifft(sqrt(classicf).*gaussN1); y = ifft(sqrt(classicf).*gaussN2); rayleigh_amp = sqrt(abs(x).2+abs(y).2); y=rayleigh_amp/sqrt(mean(rayleigh_amp.*rayleigh_amp); 附录四 多径信道的MATLAB程序 smldPe_duojin.m function p=smldPe_duojin(snr_in_dB) E=1; SNR=exp(snr_in_dB*log(10)/10); %signal to noise ratio sgma=E/sqrt(2*SNR); %sigma,standard deviation of noise N=214; % generation of the binary data source for i=1:N temp=rand; if(temp0.5) dsource(i)=0; else dsource(i)=1; end end; % detection ,and probability of error calculation numoferr=0; for i=1:N % matched filter outputs if(dsource(i)=0) dsource1(i,:)=-1 -1 -1 -1 -1 -1 -1 -1 -1 -1/10; else(dsource(i)=1) dsource1(i,:)=1 1 1 1 1 1 1 1 1 1/10; end end temp1=reshape(dsource1,N*10,1) for i=1:10*N-3 temp2(i)=temp1(i+3);%delay end temp3=temp2 0 0 0; temp4=reshape(temp3,10,N); dsource2=temp4; for i=1:N for j=1:10 a1=ray; a2=ray; r1=0.5*a1*dsource1(i,j)+0.5*a2*dsource2(i,j); end r=r1+gngauss(sgma); if(r0) decis=0; else decis=1; end; if(decis=dsource(i) numoferr=numoferr+1; end; end; p=numoferr/N; % probability of error estimate 附录五 加入卷积编码的AWGN信道的MATLAB程序 smldPe_conv.m function p=smldPe_conv(snr_in_dB) E=1; N=214; % generation of the binary data source for i=1:N temp=rand; if(temp