Ai đó có thể đề nghị làm thế nào tôi có thể tính toán khoảnh khắc thứ hai (hoặc toàn bộ hàm tạo khoảnh khắc) của cosin của hai vectơ ngẫu nhiên gaussian , mỗi phân phối là , độc lập với nhau không? IE, thời điểm cho biến ngẫu nhiên sauN ( 0 , Σ )
Câu hỏi gần nhất là hàm tạo mô men của sản phẩm bên trong của hai vectơ ngẫu nhiên gaussian có nguồn gốc MGF cho sản phẩm bên trong. Ngoài ra còn có câu trả lời này từ mathoverflow liên kết câu hỏi này với phân phối giá trị riêng của ma trận hiệp phương sai mẫu, nhưng tôi không thấy ngay cách sử dụng chúng để tính toán giây thứ hai.
Tôi nghi ngờ rằng khoảnh khắc thứ hai tỷ lệ với một nửa giá trị bản địa của vì tôi nhận được kết quả này thông qua thao tác đại số cho 2 chiều và cũng cho 3 chiều từ đoán và kiểm tra. Đối với giá trị riêng thêm tối đa 1, giây thứ hai là:
Sử dụng như sau để kiểm tra số
val1[a_, b_, c_] := (a + b + c)/(Sqrt[a] + Sqrt[b] + Sqrt[c])^2
val2[a_, b_, c_] := Block[{},
x := {x1, x2, x3};
y := {y1, y2, y3};
normal := MultinormalDistribution[{0, 0, 0}, ( {
{a, 0, 0},
{0, b, 0},
{0, 0, c}
} )];
vars := {x \[Distributed] normal, y \[Distributed] normal};
NExpectation[(x.y/(Norm[x] Norm[y]))^2, vars]]
val1[1.5,2.5,3.5] - val2[1.5,2.5,3.5]
Kiểm tra công thức cho 4 biến (trong giới hạn số):
val1[a_, b_, c_,
d_] := (a + b + c + d)/(Sqrt[a] + Sqrt[b] + Sqrt[c] + Sqrt[d])^2
val2[a_, b_, c_, d_] := Block[{},
x := {x1, x2, x3, x4};
y := {y1, y2, y3, y4};
normal :=
MultinormalDistribution[{0, 0, 0,
0}, {{a, 0, 0, 0}, {0, b, 0, 0}, {0, 0, c, 0}, {0, 0, 0, d}}];
vars := {x \[Distributed] normal, y \[Distributed] normal};
NExpectation[(x.y/(Norm[x] Norm[y]))^2, vars]]
val1[0.5, 1.5, 2.5, 3.5] - val2[0.5, 1.5, 2.5, 3.5]