ask.m

该代码实现了 ook等不同通信通信方式的误码率的仿真比较。

clear
Fc=2;     %载频
Fs=40;     %系统采样频率
Fd=1;      %码速率
E=1;       %每个码能量
N=Fs/Fd;   
numSymb=25;   %进行仿真信息代码个数
M=2;           %进制数
snrlog=10;           %信噪比dB
snr=10^(snrlog/10);%信噪比数值
numPlot=15;

%产生25个二进制随机码
x=randsrc(numSymb,1,[0:M-1]);
figure(1)
stem([0:numPlot-1],x(1:numPlot));
title('二进制随机序列')
xlabel('Time');
ylabel('Amplitude');

%调制
y=zeros(1,Fs*numSymb/Fd);
y=ask_mod(x,Fc,Fd,Fs);
numModPlot=numPlot*Fs;
t=[0:numModPlot-1]./Fs;
figure(2)
plot(t,y(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('调制后的信号')
xlabel('Time');
ylabel('Amplitude');
%在已调信号中加入高斯白噪声
sgma=sqrt(E/snr)/2;
noise=sgma*randn(1,length(y));
r=y+noise;
figure(3)
plot(t,r(1:length(t)),'b-');
axis([min(t) max(t) -1.5 1.5]);
title('调制后加噪声的信号')
xlabel('Time');
ylabel('Amplitude');
%相干解调

%非相干解调

%t=0:0.01:8;
%y=sin(2*pi*t);   %载波
%x=[ones(1,100),zeros(1,100),ones(1,100),zeros(1,100),ones(1,100),zeros(1,100),ones(1,100),...
%       zeros(1,101)];
%z=x.* y;
%plot(t,z)


%function grsv=gngauss(m,sgma)  
%输入参数可以是平均值m，方差都有sgma，也可以只有方差sgma
%缺省值平均0，方差为1
%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;
%grsv(1)=m+z*cos(2*pi*u);   grsv(2)=m+z*cos(2*pi*u);