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);