clear all;
close all;
N=100;

A=1;
fc=2;
T=10;
fs=1/T;
N_sample=1000;
N=256;
dt=T/N_sample;

t=0:dt:T-dt;
st=[ones(1,N_sample/10),zeros(1,N_sample/10),ones(1,N_sample/10),zeros(1,N_sample/10),ones(1,N_sample/10),zeros(1,N_sample/10),ones(1,N_sample/10),zeros(1,N_sample/10),ones(1,N_sample/10),zeros(1,N_sample/10)];
figure(1)
subplot(211)
plot(t,st);
axis([0 10 0 1.2]);ylabel('輸入信號');


[f,Sf]=T2F(t,st);
subplot(212)
plot(f,abs(Sf).^2);
axis([-4 +4 0 30]);ylabel('輸入信號功率譜');


%%ook
figure(2)
ht=A*cos(2*pi*fc*t);
s_ook=st.*ht;
subplot(211)
plot(t,s_ook);
axis([0 10 -1.2 1.2]);ylabel('ook信號');

[f,S_ook]=T2F(t,s_ook);
subplot(212)
plot(f,abs(S_ook).^2);
axis([-fc-4 fc+4 0 8]);ylabel('ook信號功率譜');


%%2psk
v_2psk=2*st-1;
s_2psk=v_2psk.*ht;
figure(3)
subplot(211)
plot(t,s_2psk);
axis([0 10 -1.2 1.2]);ylabel('2psk信號');

[f,S_2psk]=T2F(t,s_2psk);
subplot(212)
plot(f,abs(S_2psk).^2);
axis([-fc-4 fc+4 0 15]);ylabel('2psk信號功率譜');

%%2fsk
figure(4)
fc1=2;fc2=5;
ht1=A*cos(2*pi*fc1*t);
ht2=A*cos(2*pi*fc2*t);
v_2fsk=st-1;
s_2fsk=st.*ht1-v_2fsk.*ht2;
subplot(211)
plot(t,s_2fsk);
axis([0 10 -1.2 1.2]);ylabel('2fsk信號');

[f,S_2fsk]=T2F(t,s_2fsk);
subplot(212)
plot(f,(abs(S_2fsk).^2));
axis([-fc-4 fc+4 0 8]);ylabel('2fsk信號功率譜');