From a22b0f45c90e75452e9efa50aebfb6d8c1ef79c8 Mon Sep 17 00:00:00 2001
From: Ksyer <>
Date: Tue, 16 Jul 2024 15:55:23 +0800
Subject: [PATCH] Add Analyse_of_Psi

---
 1 - experiment/FAR/Analyse_of_Psi.m | 42 +++++++++++++++++++++++++++++
 1 file changed, 42 insertions(+)
 create mode 100755 1 - experiment/FAR/Analyse_of_Psi.m

diff --git a/1 - experiment/FAR/Analyse_of_Psi.m b/1 - experiment/FAR/Analyse_of_Psi.m
new file mode 100755
index 0000000..b821b74
--- /dev/null
+++ b/1 - experiment/FAR/Analyse_of_Psi.m	
@@ -0,0 +1,42 @@
+%% 参数设置
+N = 16;
+M = 4;
+
+A = get_Psi(N, M, 0);
+
+%% 检查每一行 / 列的二范数
+col_norms = zeros(N, 1);
+for i = 1: N
+    col_norms(i) = norm(A(i, :), 2);
+end
+
+row_norms = zeros(N * M, 1);
+for i = 1: N * M
+    row_norms(i) = norm(A(:, i), 2);
+end
+
+figure(1);
+
+subplot(211);
+plot(1: N, col_norms);
+xlabel("Col");
+ylabel("L2 norm of col");
+title("L2 norm of cols");
+
+subplot(212);
+plot(1: N * M, row_norms);
+xlabel("Row");
+ylabel("L2 norm of row");
+title("L2 norm of rows");
+
+
+%%
+AH = conj(A).';
+figure(2);
+result = abs(A * AH);
+subplot(211);
+plot(diag(result));
+subplot(212);
+heatmap(result);
+a = sum(result(:)) - sum(diag(result));
+fprintf("%f", a);