RIS + MIMO Demo

from pytwinnet.ris import RISPanel, ris_link_gain, phase_opt_greedy
from pytwinnet.mimo import mimo_rayleigh
import numpy as np

N = 64
h_tr = mimo_rayleigh(nt=N, nr=1)[:,0]
h_rr = mimo_rayleigh(nt=1, nr=N)[0,:]
theta = phase_opt_greedy(h_tr, h_rr)
panel = RISPanel(N); panel.set_phases(theta)
h_eff = ris_link_gain(h_tr, h_rr, panel.theta)
print("|h_eff| =", np.abs(h_eff))