Có thể được hướng dẫn để chứng minh kết quả này từ các nguyên tắc đầu tiên và kết quả cơ bản , khai thác các thuộc tính của các hàm tạo tích lũy (chính xác như trong các bằng chứng tiêu chuẩn của Định lý giới hạn trung tâm). Nó đòi hỏi chúng ta phải hiểu tốc độ tăng trưởng của các số hài tổng quát cho Các tốc độ tăng trưởng này nổi tiếng và dễ dàng đạt được khi so sánh với các tích phân : chúng hội tụ cho và nếu không thì phân kỳ logarit cho . s = 1 , 2 , 210 . ∫ n 1 x - s d x s > 1 s = 1
H(n,s)=∑k=1nk−s
s=1,2,….∫n1x−sdxs>1s=1
Đặt và . Theo định nghĩa, hàm tạo tích lũy (cgf) của là1 ≤ k ≤ n ( X k - 1 / k ) / B nn≥21≤k≤n(Xk−1/k)/Bn
ψk,n(t)=logE(exp(Xk−1/kBnt))=−tkBn+log(1+−1+exp(t/Bn)k).
Sự mở rộng chuỗi của phía bên tay phải, thu được từ việc mở rộng quanh , có dạngz = 0log(1+z)z=0
ψk,n(t)=(k−1)2k2B2nt2+k2−3k+26k3B3nt3+⋯+kj−1−⋯±(j−1)!j!kjBjntj+⋯.
Các tử số của các phân số là đa thức tính bằng với số hạng đầu . Bởi vì mở rộng nhật ký hội tụ hoàn toàn cho , bản mở rộng này hội tụ hoàn toàn khik j - 1 | - 1 + điểm kinh nghiệm ( t / B n )kkj−1∣∣−1+exp(t/Bn)k∣∣<1
|exp(t/Bn)−1|<k.
(Trong trường hợp nó hội tụ ở mọi nơi.) Đối với giá trị cố định và tăng của , độ phân kỳ (rõ ràng) của ngụ ý miền hội tụ tuyệt đối tăng lên tùy ý. Do đó, đối với bất kỳ cố định và đủ lớn , sự mở rộng này hội tụ hoàn toàn.k n B n t nk=1knBntn
Do đó, đối với đủ lớn , chúng tôi có thể tính tổng theo kỳ hạn theo quyền hạn của để có được cgf của ,ψ k , n k t S n / B nnψk,nktSn/Bn
ψn(t)=∑k=1nψk,n(t)=12t2+⋯+1Bjn(∑k=1n(k−1−⋯±(j−1)!k−j))tjj+⋯.
Lấy các số hạng trong tổng số trên một lần yêu cầu chúng ta đánh giá các biểu thức tỷ lệ vớik
b(s,j)=1Bjn∑k=1nk−s
cho và . Sử dụng sự không triệu chứng của các số hài tổng quát được đề cập trong phần giới thiệu, nó dễ dàng theo saus = 1 , 2 , ... , jj≥3s=1,2,…,j
B2n=H(n,1)−H(n,2)∼log(n)
cái đó
b(1,j)∼(log(n))1−j/2→0
và (đối với )s>1
b(s,j)∼(log(n))−j/2→0
khi phát triển lớn. Do đó, tất cả các điều khoản trong việc mở rộng ngoài hội tụ về 0, khi đó hội tụ đến cho bất kỳ giá trị nào của . Do sự hội tụ của cgf ngụ ý sự hội tụ của hàm đặc trưng, chúng tôi kết luận từ Định lý liên tục Levy rằng tiếp cận một biến ngẫu nhiên có cgf là : đó là biến Bình thường tiêu chuẩn, QED .ψ n ( t ) t 2 ψ n ( t ) t 2 / 2 t S n / B n t 2 / 2nψn(t)t2ψn(t)t2/2tSn/Bnt2/2
Phân tích này chỉ ra mức độ hội tụ tinh tế như thế nào: trong khi trong nhiều phiên bản của Định lý giới hạn trung tâm, hệ số của là (đối với ), ở đây hệ số là chỉ : sự hội tụ chậm hơn nhiều. Theo nghĩa này, chuỗi các biến được tiêu chuẩn hóa "chỉ vừa đủ" trở thành Bình thường. O ( n 1 - j / 2 ) j ≥ 3 O ( ( ( log ( n ) ) 1 - j / 2 )tjO(n1−j/2)j≥3O(((log(n))1−j/2)
Chúng ta có thể thấy sự hội tụ chậm này trong một loạt mô phỏng. Các biểu đồ hiển thị lần lặp độc lập cho bốn giá trị của . Các đường cong màu đỏ là đồ thị của các hàm mật độ chuẩn thông thường để tham khảo trực quan. Mặc dù rõ ràng là có xu hướng dần dần về tính quy tắc, ngay cả ở (trong đó vẫn còn khá lớn) vẫn có sự phi quy tắc đáng chú ý, được chứng minh trong độ lệch (bằng trong mẫu này). (Không có gì ngạc nhiên khi độ lệch của biểu đồ này gần với , vì đó chính xác là thuật ngữ trong cgf là gì.) n n = 1000 ( log ( n ) ) - 1 / 2 ≈ 0,38 0,35 (105nn=1000(log(n))−1/2≈0.380.35 t 3(log(n))−1/2t3
Đây là R
mã cho những người muốn thử nghiệm thêm.
set.seed(17)
par(mfrow=c(1,4))
n.iter <- 1e5
for(n in c(30, 100, 300, 1000)) {
B.n <- sqrt(sum(rev((((1:n)-1) / (1:n)^2))))
x <- matrix(rbinom(n*n.iter, 1, 1/(1:n)), nrow=n, byrow=FALSE)
z <- colSums(x - 1/(1:n)) / B.n
hist(z, main=paste("n =", n), freq=FALSE, ylim=c(0, 1/2))
curve(dnorm(x), add=TRUE, col="Red", lwd=2)
}