clc;
clear?all;
close?all;

%發送端產生QPSK信號
len?=?2000;
xI?=?randint（len1）;
xQ?=?randint（len1）;
for?i?=?1:len
????if（xI（i）?==?0）
????????xI（i）?=?-1;
????end
????if（xQ（i）?==?0）
????????xQ（i）?=?-1;
????end
end
x?=?（xI?+?1i*xQ）/sqrt（2）;
scatterplot（x）;
title（‘發送端QPSK信號星座圖‘）;

hA?=?[0.04-0.050.07-0.21-0.50.720.3600.210.030.07];
hB?=?[0.4070.8150.407];
hC?=?[0.2270.460.68880.460.227];

%QPSK信號過信道
h?=?hA;?????????????????????????????????%實驗中可調整
xk?=?conv（x?h‘same‘）;

%對過信道的QPSK信號加噪
SNR?=?20;???????????????????????????????%實驗中可調整
yI?=?awgn（real（xk）?SNR?‘measured‘）;
yQ?=?awgn（imag（xk）?SNR?‘measured‘）;
y?=?yI?+?1i*yQ;
scatterplot（y）;
title（‘過信道加噪后的QPSK信號星座圖‘）;

%利用LMS算法對系統進行逆識別
len_train?=?500;
N?=?13;
xn?=?real（y（1:len_train））;
dn?=?real（x（1:len_train））;
%?diag_R?=?eig（xn*xn‘）;
%?u1?=?2