यह पहले सिद्धांतों और बुनियादी परिणामों से इस परिणाम को प्रदर्शित करने के लिए शिक्षाप्रद हो सकता है , क्यूम्युलेंट जनरेटिंग फ़ंक्शंस के गुणों का शोषण करता है (बिल्कुल केंद्रीय सीमा प्रमेय के मानक प्रमाणों में)। हमें सामान्यीकृत हार्मोनिक संख्या लिए की वृद्धि की दर को समझना होगा ये विकास दर अच्छी तरह से ज्ञात हैं और इंटीग्रल्स तुलना में आसानी से प्राप्त की जाती हैं : वे लिए अभिसरण करती हैं और अन्यथा लिए लघुगणकीय रूप से विचलन करती हैं ।
H(n,s)=∑k=1nk−s
s=1,2,….∫n1x−sdxs>1s=1
चलो और । परिभाषा के अनुसार, का उत्पादक कार्य (cgf) हैn≥21≤k≤n(Xk−1/k)/Bn
ψk,n(t)=logE(exp(Xk−1/kBnt))=−tkBn+log(1+−1+exp(t/Bn)k).
दाहिने हाथ की ओर की श्रृंखला का विस्तार, आसपास के विस्तार से प्राप्त होता है , जो रूप लेता हैlog(1+z)z=0
ψk,n(t)=(k−1)2k2B2nt2+k2−3k+26k3B3nt3+⋯+kj−1−⋯±(j−1)!j!kjBjntj+⋯.
अंशों के अंशों में बहुपद होते हैं साथ अग्रणी शब्द । क्योंकि लॉग का विस्तार पूरी तरह से लिए परिवर्तित होता है , यह विस्तार पूर्ण रूप से परिवर्तित होता है जबkkj−1∣∣−1+exp(t/Bn)k∣∣<1
|exp(t/Bn)−1|<k.
(केस यह हर जगह होता है।) के निश्चित और बढ़ते मूल्यों के लिए , का स्पष्ट (स्पष्ट) निरपेक्ष रूप से बड़े अभिसरण के डोमेन को बड़े पैमाने पर बढ़ता है। इस प्रकार, किसी भी निश्चित और पर्याप्त रूप से बड़े , यह विस्तार बिल्कुल परिवर्तित होता है।k=1knBntn
पर्याप्त रूप से बड़े , तो, इसलिए हम , का प्राप्त करने के लिए शक्तियों में व्यक्तिगत रूप से को टर्म से अधिक कर ,nψk,nktSn/Bn
ψn(t)=∑k=1nψk,n(t)=12t2+⋯+1Bjn(∑k=1n(k−1−⋯±(j−1)!k−j))tjj+⋯.
एक समय में पर शर्तों को लेने के लिए हमें आनुपातिक अभिव्यक्ति का मूल्यांकन करना होगाk
b(s,j)=1Bjn∑k=1nk−s
के लिए और । परिचय में वर्णित सामान्यीकृत हार्मोनिक संख्याओं के एसिम्पोटिक्स का उपयोग करना, यह आसानी से निम्नानुसार हैj≥3s=1,2,…,j
B2n=H(n,1)−H(n,2)∼log(n)
उस
b(1,j)∼(log(n))1−j/2→0
और ( )s>1
b(s,j)∼(log(n))−j/2→0
बड़े रूप में बढ़ता है। नतीजतन, से परे के विस्तार में सभी शब्द शून्य में परिवर्तित हो जाते हैं, के किसी भी मान के लिए परिवर्तित हो जाता । चूँकि cgf के अभिसरण से अभिप्राय फ़ंक्शन के अभिसरण से है, हम लेवी निरंतरता प्रमेय से यह निष्कर्ष निकालते हैं कि एक यादृच्छिक चर के पास आता है जिसका cgf 2/2 है : जो कि मानक सामान्य चर, QED है ।nψn(t)t2ψn(t)t2/2tSn/Bnt2/2
यह विश्लेषण केवल यह बताता है कि अभिसरण कितना नाज़ुक है: जबकि सेंट्रल लिमिट के कई संस्करणों में गुणांक का गुणांक ( ) है, यहाँ गुणांक है केवल : अभिसरण बहुत धीमा है। इस अर्थ में मानकीकृत चर "बस मुश्किल से" का क्रम सामान्य हो जाता है।tjO(n1−j/2)j≥3O(((log(n))1−j/2)
हम इस धीमे अभिसरण को सिमुलेशन की एक श्रृंखला में देख सकते हैं। हिस्टोग्राम्स चार मूल्यों के लिए स्वतंत्र पुनरावृत्तियों को प्रदर्शित करते हैं । लाल वक्र दृश्य संदर्भ के लिए मानक सामान्य घनत्व कार्यों के ग्राफ हैं। यद्यपि स्पष्ट रूप से सामान्यता की ओर एक क्रमिक प्रवृत्ति है, यहां तक कि (जहां अभी भी बड़े आकार का है) वहाँ प्रशंसनीय गैर-सामान्यता बनी हुई है, जैसा कि तिरछेपन में स्पष्ट है। ( इस नमूने में बराबर )। (यह कोई आश्चर्य की बात नहीं है कि इस हिस्टोग्राम का तिरछापन करीब है , क्योंकि यह ठीक है कि cgf में टर्म क्या है।)105nn=1000(log(n))−1/2≈0.380.35(log(n))−1/2t3
यहाँ R
उन लोगों के लिए कोड है जो आगे प्रयोग करना चाहते हैं।
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)
}