您现在的位置：毕业论文 >> 论文 >> 正文

通信原理仿真实验系统设计与实现第13页

更新时间：2009-6-5: 来源：毕业论文

通信原理仿真实验系统设计与实现第13页
z=x.*y; %幅频键控
plot(t,z)
%对二元序列１０１１００１０画出２ＦＳＫ（移频銉控）的波形
x=0:0.01:8;
t=[ones(1,100),zeros(1,100),ones(1,100),ones(1,100),zeros(1,100),zeros(1,100),ones(1,100),zeros(1,101)];
y=sin(x.*(2*pi+2*t));
plot(x,y)

%PSK包络的绘制pskbaoluo.m
echo on
T=1; M=8; Es=T/2; fc=6/T;
N=120; delta_T=T/(N-1);   t=0:delta_T:T;
u0=sqrt(2*Es/T)*cos(2*pi*fc*t);
u1=sqrt(2*Es/T)*cos(2*pi*fc*t+2*pi/M);%求出八个波形
u2=sqrt(2*Es/T)*cos(2*pi*fc*t+4*pi/M);
u3=sqrt(2*Es/T)*cos(2*pi*fc*t+6*pi/M);
u4=sqrt(2*Es/T)*cos(2*pi*fc*t+8*pi/M);
u5=sqrt(2*Es/T)*cos(2*pi*fc*t+10*pi/M);
u6=sqrt(2*Es/T)*cos(2*pi*fc*t+12*pi/M);
u7=sqrt(2*Es/T)*cos(2*pi*fc*t+14*pi/M);
subplot(8,1,1);              plot(t,u0);
subplot(8,1,2);              plot(t,u1);
subplot(8,1,3);              plot(t,u2);
subplot(8,1,4);              plot(t,u3);
subplot(8,1,5);              plot(t,u4);
subplot(8,1,6);              plot(t,u5);
subplot(8,1,7);              plot(t,u6);
subplot(8,1,8);              plot(t,u7);
%函数文件pskmoto.m进行蒙特卡洛仿真
function[pb,ps]=pskmoto(snr_in_dB)
%pskmoto求出以ｄＢ为单位的给定信噪比的比特误码率和符号误码率
N=1000;
E=1;
snr=10^(snr_in_dB/10);
sgma=sqrt(E/snr)/2;
%信号映射
n=[0 0];
s00=[1 0];     s01=[0 1];
s11=[-1 0];     s10=[0 -1];
for i=1:N,
    temp=rand; %区间（０，１）的一个均匀随机变量
    if(temp<0.25), %信源输出为００的概率为１／４
        dsource1(i)=0;
        dsource2(i)=0;
    elseif(temp<0.5)
        dsource1(i)=0;
        dsource2(i)=1;
    elseif(temp<0.75)
        dsource1(i)=1;
        dsource2(i)=0;
    else
        dsource1(i)=1;
        dsource2(i)=1;
    end;
end;
%判决，误码率的计算
numofsymbolerror=0;         numofbiterror=0;
for i=1:N,
    %在判决器的接收端的信号，对于第i个符号为：
    n=gngauss(sgma);
    if((dsource1(i)==0)&(dsource2(i)==0)),
        r=s00+n;
    elseif((dsource1(i)==0)&(dsource2(i)==1)),
        r=s01+n;
         if((dsource1(i)==1)&(dsource2(i)==0)),
        r=s10+n;
    else
        r=s11+n;
    end;
    %以下为计算互相关量度
    c00=dot(r,s00); c01=dot(r,s01);
    c10=dot(r,s10); c11=dot(r,s11);
    %第i个符号的判决如下进行
    c_max=max([c00 c01 c10 c11]);
    if(c00==c_max),
        decis1=0; decis2=0;
    elseif(c01==c_max),
        decis1=0; decis2=1;
    elseif(c10==c_max),
        decis1=1; decis2=0;
    else
        decis1=1; decis2=1;
    end;
    %若判决不对，计错器数加１
    symbolerror=0;
    if(decis1~=dsource1(i)),
        numofbiterror=numofbiterror+1;   symbolerror=1;
    end;
    if(decis2~=dsource1(i)),
        numofbiterror=numofbiterror+1;   symbolerror=1;
    end;
    if(symbolerror==1),
     numofsymbolerror=numofsymbolerror+1; %误比特率数加1
    end;
end;
ps=numofsymbolerror/N; %总共发出Ｎ个符号
pb=numofbiterror/(2*N);
end
%PSK系统的实现pskmotocalo.m
echo on
SNRindB1=0:2:10; %定义信噪比序列，共辣个值
SNRindB2=0:0.1:10; %扫描用的信噪比序列，采用ＤＢ作单位
for i=1:length(SNRindB1),
    [pb,ps]=pskmoto(SNRindB1(i)); %计算误比特率
    smld_bit_err_prb(i)=pb;
    smld_symbol_err_prb(i)=ps;
end;
for i=1:length(SNRindB2)
    SNR=exp(SNRindB2(i)*log(10)/10); %转化信噪比为数值表示
    theo_err_prb(i)=Qfunct(sqrt(2*SNR));
end;
semilogy(SNRindB1,smld_bit_err_prb,'*'); %以对数形式作Ｙ坐标绘图
hold
semilogy(SNRindB1,smld_symbol_err_prb,'o');%作出实际的信噪比
semilogy(SNRindB2,theo_err_prb); %作出理论信噪比
%函数文件qammoto.m
function[p]=qammoto(snr_in_dB)
N=1000;
d=1;
Eav=10*d^2;
snr=10^(snr_in_dB/10);
sgma=sqrt(Eav/(8*snr));
M=16;
for i=1:N
    temp=rand;
    dsource(i)=1+floor(M*temp);
end;