MATLAB实习题目实验1.docVIP

建立一个基于AWGN信道的BPSK调制的通信系统。其中AWGN信道的信噪比变化范围为-5dB~10dB，最后请画出该系统的性能曲线（即信噪比-误码率曲线，信噪比步长为0.5dB） 实验原理图： SIMULINK系统 需要设置的参数： 代码： function bpskerr_rate=trng_gauss close all; clear all; clc echo on SNRindb1=-5:0.5:10; SNRindb2=-5:0.5:10; for i=1:length(SNRindb1) smld_err_prb(i)=cm_bpske(SNRindb1(i)); %simulated error rate echo off; end echo on for i=1:length(SNRindb2) SNR=exp(SNRindb2(i)*log(10)/10); %signal to noise ratio theo_err_prb(i)=erfc(sqrt(SNR/2))/2; %theoretical symbol error rate echo off; end echo on %plotting commands follows semilogy(SNRindb1,smld_err_prb,*); xlabel(SNR(dB)----------); ylabel(BER--------------); hold on semilogy(SNRindb2,theo_err_prb); xlabel(SNR(dB)----------); ylabel(BER--------------); legend( simulated result, theoretical result); hold off function [n1,n2]=gngauss(m,sgma) % [n1 n2]=gngauss(m,sgma) % [n1 n2]=gngauss(sgma) % [n1 n2]=gngauss % GNGAUSS generates two independent gaussian variables with % mean m S.D sgma % if one of the input argement missing it takes mean as zero % if neither mean nor variance is given, it generates two %standard gaussian random variables if nargin==0 m=0; sgma=1; elseif nargin==1 sgma=m; m=0; end u=rand; % uniform random variable in (0,1) z=sgma*(sqrt(2*log(1/(1-u)))); % rayleigh distributed random variable u=rand; % another uniform distributed variable in (0,1) n1=m+z*cos(2*pi*u); n2=m+z*sin(2*pi*u); function[p]=cm_bpske(snr_in_db) % [p]=cm_bpske(snr_in_db) % cm_bpske finds the probability of error for the given value of snr_in_db N=5000; E=1; %energy per symbol snr=10^(snr_in_db/10); %signal to noise ratio sgma=E/sqrt(snr*2); %noise variance %generation of data source for i=1:N temp=rand; if(temp0.5) dsource(i)=0; else dsource(i)=1; end end %de