len=500;
R=randint（1?len）;
I=randint（1?len）;
Real=2*R-1;
Imag=2*I-1;
X=（Real+1j*Imag）/sqrt（2）;
scatterplot（X）;
title（‘發送端QPSK信號星座圖‘）;

h=[0.04?-0.05?0.07?-0.21?-0.5?0.72?0.36?0?0.21?0.03?0.07];
X1=conv（Xh）;
X2=zeros（1len）;
for?i=1:len
????X2（i）=X1（i）;
end
scatterplot（X2）;
title（‘經過線性系統QPSK信號星座圖‘）;

SNR=20;
Y_R=awgn（real（X2）SNR）;
Y_I=awgn（imag（X2）SNR）;
Y=Y_R+1j*Y_I;
scatterplot（Y）;
title（‘疊加噪聲后QPSK信號星座圖‘）;

N=length（h）;
H1?=?zero