stamp = datetime.now().strftime("%y-%m-%d_%H-%M-%S") f1 = open("./simresults/lambda_theta_LT" + stamp + ".txt", 'w') #---------------------------------------------------------------design: C = pl.CapacityBSC(N, design_p) deratepercentage = 10 tolerable_error = -2 if designbyrate: #K=int((100-deratepercentage)*C/100) K = int(C) try: I = pcon.getGCHsim('ZK', N, design_p, K) RI = pcon.getGChsim('ZK', design_p, N, N) except: (I, E) = pcon.getGChZCK(design_p, N, K) (RI, E) = pcon.getGChZCK(design_p, N, N) else: try: I = pcon.getGCHsim('ZL', N, design_p, tolerable_error) except: (I, E) = pcon.getGChZCL(design_p, N, tolerable_error) #==========================================================Construction print "LLR Lambda Channel REPORT"
json.dump("RATE Vs FER REPORT derate pch", f1) f1.write("\n") json.dump("-----------------------------", f1) f1.write("\n") json.dump("N=" + str(N), f1) f1.write("\n") json.dump("p_decode=" + str(design_p), f1) f1.write("\n") json.dump("capacity*N:" + str(C), f1) f1.write("\n") json.dump("sim ran :" + str(runsim), f1) f1.write("\n") try: I_order = pcon.getGCHsim('ZK', N, design_p, N) except: (I_order, E) = pcon.getGChZCK(design_p, N, N) used_rate = [] FER = [] Ratelist = [i * float(C) / N for i in deratelist] for R in Ratelist: if R <= (C / N): K = int(R * N) I = I_order[:K] #print R block_error = pch.polarchannelsim(N, design_p, design_p, I, runsim, False) used_rate.append(R)
stamp = datetime.now().strftime("%y-%m-%d_%H-%M-%S") f1 = open("./simresults/lambda_avg" + stamp + ".txt", 'w') #==========================================================Construction print "LLR Lambda REPORT" print "---------------------------" print "N=" + str(N) print "design_p=" + str(p) PotGCh = int(ma.floor(pl.CapacityBSC(N, p))) print "Capacity for " + str(N) + "channels:" + str(PotGCh) print "tolerable error exponent:" + str(tolerable_error) print "sim ran :" + str(runsim) #-----------------------------------------ZC #(I,E)=pcon.getGChZCL(p,N,tolerable_error) #(I,E)=pcon.getGChZCK(p,N,K) I = pcon.getGCHsim("MK_ALL", N, p, K) print "Good Channels:" print I print "Number of good channels:" print len(I) R = float(len(I)) / N print "R=" + str(R) print "Frozen channels:" B = list(set(range(N)) - set(I)) print len(B) ##########################################file json.dump("LLR Lambda CHANNEL REPORT", f1) f1.write("\n") json.dump("---------------------------", f1) f1.write("\n") json.dump("N=" + str(N), f1)