clear all; clc;close all;
Rs=1e6;fs=20e6;phi=pi/2; Es=2;SNR=10; m=16;
v=randi(1,1000,m);
qam=qammod(v,m);
st=awgn(qam,SNR,'measured');
giaima=qamdemod(st,m);
[t,st1]=taoxung(qam,Es,phi,Rs,fs);
plot(t(1:1000),st1(1:1000));
figure;
st2=awgn(st1,SNR,'measured');
plot(t(1:1000),st2(1:1000));
figure;
sf1=fft(st1,length(st1));
f=0:1/(max(t)-min(t)):fs;
sf2=fftshift(sf1);
plot(f,sf2)
figure;
sf3=fft(st2,length(st2));
sf4=fftshift(sf3);
plot(f,sf4);
h=scatterplot(st,1,0);
hold on
scatterplot(qam,1,0,'ro',h);
eyediagram(st1,100);
eyediagram(st2,100);