Example #1
0
        symbol_y = sym.Addnoise(tab, symbol)  # 초기 심볼 + 노이즈
        # 초기화 부분

        if QAM == 4:
            center, temp_Y = kc.Kmeans(Cluster, symbol_y, QAM, ch)
            est_ch, mse = est.Est_Ch(center, ch)
        else:
            center, temp_Y = kc.Kmeans_16QAM(Cluster, symbol_y, QAM, ch)
            hk_center = kc.Ch_Est_Center(center, QAM)
            est_ch, mse = est.Est_Ch(hk_center, ch)
        # print(snr)
        # print(ch)
        # print(est_ch)
        # print("------")
        symbol_r = sym.DividCh(symbol_y, est_ch)

        data_r = de.Demode(QAM, symbol_r)
        SER[snr] += de.Equals(data_s, data_r)
        MSE[snr] += mse

for i in range(SNR):
    SER[i] = SER[i] / (symbol_num * Count_Total)
    MSE[i] = MSE[i] / Count_Total
    print(i)
    print(SER[i])

# 심볼점 찍는 기능
# Temp_Y
# for j in range(QAM):
#     for i in range(symbol_num):