From 1ce76ff4dc312a5ea13850064670a4e4707a73a3 Mon Sep 17 00:00:00 2001 From: mayge Date: Wed, 24 Sep 2025 04:30:10 -0400 Subject: [PATCH] =?UTF-8?q?=E4=BF=AE=E5=A4=8D=E4=BA=86fit=5Fconstant=3D=3D?= =?UTF-8?q?False=E7=9A=84bug,=E4=BD=86=E6=98=AF=E4=B8=8D=E5=BB=BA=E8=AE=AE?= =?UTF-8?q?=E4=BD=BF=E7=94=A8?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- core/MultiplePortQR.py | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/core/MultiplePortQR.py b/core/MultiplePortQR.py index 8b33353..f080085 100644 --- a/core/MultiplePortQR.py +++ b/core/MultiplePortQR.py @@ -161,7 +161,7 @@ class MultiplePortQR: M_kp = None for i in range(self.ports): for j in range(self.ports): - M0 = np.zeros((K,N*ports**2),dtype=complex) + M0 = np.zeros((K,N*self.ports**2),dtype=complex) M0[:,index*N:(index+1)*N] = Phi M0 = np.hstack([M0, -(H[:,i,j].reshape(-1,1) * Phi)]).reshape((K, -1))[keep,:] # (K, 2N), complex index+=1 @@ -225,10 +225,10 @@ class MultiplePortQR: H_kp = None for i in range(self.ports): for j in range(self.ports): - H_kp0 = weights_kp @ (H[:,i,j]).reshape(1,-1)[keep,:] - H_kp = H_kp0 if H_kp is None else np.hstack([H_kp, H_kp0]) + H_0 = H[:,i,j][keep] + H_kp = H_0 if H_kp is None else np.hstack([H_kp, H_0]) assert H_kp is not None - H_kp = H_kp.reshape(-1,1) + H_kp = weights_kp @ H_kp.reshape(-1,1) b_re = np.real(d0 * H_kp) b_im = np.imag(d0 * H_kp) @@ -318,14 +318,14 @@ if __name__ == "__main__": network = rf.Network(f"/tmp/paramer/simulation/{id}/{id}.s2p") # network = rf.data.ring_slot ports = network.nports - K = 5 + K = 10 full_freqences = network.f[start_point:] noised_sampled_points = network.y[start_point:,:,:].reshape(-1,ports,ports) sampled_points = network.y[start_point:,:,:].reshape(-1,ports,ports) - # noised_sampled_points = network.y[start_point:,0,0].reshape(-1,1,1) - # sampled_points = network.y[start_point:,0,0].reshape(-1,1,1) + # noised_sampled_points = network.y[start_point:,1,0].reshape(-1,1,1) + # sampled_points = network.y[start_point:,1,0].reshape(-1,1,1) H,freqs = auto_select_multple_ports(noised_sampled_points,full_freqences,max_points=20) poles = generate_starting_poles(2,beta_min=1e4,beta_max=freqs[-1]*1.1)