clear;
clc;
close?all;

n=10000;??????????????????????????%10000個信息位??????????
b=randint（1n）;
f1=1;f2=2;
t=0:1/30:1-1/30;??????????????????%30數據點每秒??????????????????

%ASK
sa1=sin（2*pi*f1*t）;
E1=sum（sa1.^2）;
sa1=sa1/sqrt（E1）;
sa0=0*sin（2*pi*f1*t）;??????????????????%產生ASK單位波形，當信息位1時有正弦波形，0時沒有，

%FSK
sf0=sin（2*pi*f1*t）;??????????????????%產生FSK單位波形
E1=sum（sf0.^2）;
sf0=sf0/sqrt（E1）;
sf1=sin（2*pi*f2*t）;
E2=sum（sf1.^2）;
sf1=sf1/sqrt（E2）;

%PSK
sp0=sf0;??????????????????%PSK單位波形
sp1=-sf0;

%modulization
ask=[];fsk=[];psk=[];
for?i=1:n
????if（b（i）==1）
????????ask=[ask?sa1];??????%產生askfskpsk波形，根據信息位補充單位波形
????????fsk=[fsk?sf1];
????????psk=[psk?sp1];
????else
????????ask=[ask?sa0];
????????fsk=[fsk?sf0];
????????psk=[psk?sp0];
????end
end
figure（1）
subplot（411）
stairs（0:10[b（1:10）?b（10）]‘linewidth‘1.5）
axis（[0?10?-0.5?0.5]）;
grid?on;

subplot（412）
tb=0:1/30:10-1/30;
plot（tbask（1:10*30）‘b‘‘linewidth‘1.5）
grid?on;

subplot（413）
plot（tbfsk（1:10*30）‘r‘‘linewidth‘1.5）
grid?on;

subplot（414）
plot（tbpsk（1:10*30）‘k‘‘linewidth‘1.5）?????????????????????%畫出產生的波形前10個信息位
grid?on;

for?snr=0:20
????A=[];F=[];P=[];??????????????????%三個解調后數組，初始化均為0
????%awgn
????askn=awgn（asksnr）;
????fskn=awgn（fsksnr）;
????pskn=awgn（psksnr）;
????for?i=1:n
????????if?sum（sa1.*askn（1+30*（i-1）:30*i））>0.5?%ASK?detection
????????????A=[A?1];
????????else
????????????A=[A?0];