%fc=400;fs=8000;
%fb=200;
%d=[1 0 1 0 1 1 1 0 0 0 1 1 0 1];
%N=length(d);M=fs/fb;tb=1/fb;tc=1/fc;Nc=M*tc/tb;%fb:符号速率，tb：码元间隔，M=40：一个码元的采样点个数，N码元个数14
%s=zeros(1,N*M);
%for j = 1:N
% for i = 1:M
% s((j-1)*M+i)=sin(2*pi*i/Nc)*d(j);
% end
%end
%subplot(1,1,1);plot(s);
%------------ASK
s=zeros(3,1600);
for t=1:400
A=sin(1/40*4*pi:1/40*4*pi:4*pi)'; %40个点，在一个码元间隔中产生的2个周期的正弦波
D=[1 0 1 0 1 1 1 0 0 1 ]; %二进制基带数据(单极性)
Modout=A*D; %这就是ASK调制，得到N*M维矩阵
ModoutASK=reshape(Modout,1,length(A)*length(D)); %将M矩阵变形为1*（N*M）维矩阵
s(1,t)=ModoutASK(1,t);
end
for t=600+1:1000
A=sin(1/40*4*pi:1/40*4*pi:4*pi)'; %40个点载波1
B=sin(1/2*[1/40*4*pi:1/40*4*pi:4*pi])'; %40个点载波2
D=[1 0 1 0 1 1 1 0 0 0 ]; %二进制基带数据(单极性)
Modout=A*D+B*(~D); %这就是FSK调制，得到N*M维矩阵
ModoutFSK=reshape(Modout,1,length(A)*length(D));
s(2,t)=ModoutFSK(1,(t-600));
end
for t=1200+1:1600
A=sin(1/40*4*pi:1/40*4*pi:4*pi)'; %40个点
D=[1 0 1 0 1 1 1 0 0 0 ]; %二进制基带数据(单极性)
D=2*D-1; %单极性转换为双极性
Modout=A*D; %这就是PSK调制，得到N*M维矩阵
ModoutPSK=reshape(Modout,1,length(A)*length(D));
s(3,t)=ModoutPSK(1,(t-1200));
end
figure(1)
subplot(3,1,1);
plot(s(1,:))
subplot(3,1,2);
plot(s(2,:))
subplot(3,1,3);
plot(s(3,:))
