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):