প্রথম নীতি এবং মৌলিক ফলাফলগুলি থেকে এই ফলাফলটি প্রদর্শন করা নির্দেশক হতে পারে , ক্রমবর্ধমান উত্পাদনের কার্যাদি (কেন্দ্রীয় সীমাবদ্ধ তত্ত্বের আদর্শ প্রমাণ হিসাবে) ব্যবহার করে properties এটা তোলে আয় বৃদ্ধির হার বুঝতে আমাদের প্রয়োজন সাধারণ সুরেলা সংখ্যার জন্য এই বৃদ্ধির হারগুলি সুপরিচিত এবং সহজেই অখণ্ড তুলনা করে প্রাপ্ত হয় : তারা জন্য রূপান্তর করে এবং অন্যথায় জন্য ।
H(n,s)=∑k=1nk−s
s=1,2,….∫n1x−sdxs>1s=1
আসুন এবং । সংজ্ঞা অনুসারে, (এক্স_ এর কোমল্যান্ট জেনারেটিং ফাংশন (সিজিএফ) হয়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+⋯.
ভগ্নাংশের numerators এবং মধ্যে polynomials হয় নেতৃস্থানীয় শব্দটি সঙ্গে । কারণ লগ সম্প্রসারণ সম্পূর্ণরূপে , যখন এই সম্প্রসারণটি একেবারে রূপান্তরিত হয়kkj−1∣∣−1+exp(t/Bn)k∣∣<1
|exp(t/Bn)−1|<k.
( ক্ষেত্রে এটি সর্বত্র রূপান্তরিত হয়।) স্থির এবং ক্রমবর্ধমান মানের জন্য , এর (স্পষ্ট) বিচ্যুতিটি বোঝায় নিরঙ্কুশ রূপান্তরটির ডোমেনটি নির্বিচারে বড় হয়। সুতরাং, কোনও স্থির এবং পর্যাপ্ত পরিমাণে বড় এই সম্প্রসারণটি একেবারে রূপান্তরিত হয়।k=1knBntn
ভালোই বড় জন্য , তারপর, সুতরাং আমরা পৃথক যোগফল পারে উপর শব্দটি শব্দটি দ্বারা ক্ষমতা মধ্যে এর cgf প্রাপ্ত ,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
হিসাবে বড় হয়। ফলে সম্প্রসারণ সমস্ত পদ পরলোক শুন্যতে মিলিত, কোথা এগোয় করার কোন মানের জন্য । যেহেতু সিজিএফের রূপান্তরটি বৈশিষ্ট্যযুক্ত ফাংশনটির রূপান্তরকে বোঝায়, আমরা লেভি কন্টিনিউটি উপপাদ্য থেকে উপসংহারে যে একটি র্যান্ডম ভেরিয়েবলের কাছে পৌঁছায় যার সিজিএফ 2/2 : এটি স্ট্যান্ডার্ড নরমাল ভেরিয়েবল, কিউইডি ।nψn(t)t2ψn(t)t2/2tSn/Bnt2/2
এই বিশ্লেষণের খুঁজে পান ঠিক কিভাবে উপাদেয় অভিসৃতি হল: যেহেতু সহগ কেন্দ্রীয় সীমা উপপাদ্য অনেক সংস্করণ মধ্যে হয় (জন্য ), এখানে সহগ হয় কেবল : রূপান্তরটি অনেক ধীর গতি সম্পন্ন। এই অর্থে মানকযুক্ত ভেরিয়েবলগুলির ক্রম "সবেমাত্র সবে" স্বাভাবিক হয়ে যায়।tjO(n1−j/2)j≥3O(((log(n))1−j/2)
আমরা এই ধীর অভিযাত্রাটি বিভিন্ন সিমুলেশনগুলিতে দেখতে পাই। হিস্টোগ্রামগুলি চারটি মানের জন্য স্বতন্ত্র পুনরাবৃত্তি প্রদর্শন করে । লাল রেখাচিত্রগুলি ভিজ্যুয়াল রেফারেন্সের জন্য স্ট্যান্ডার্ড সাধারণ ঘনত্ব ফাংশনের গ্রাফ। যদিও স্পষ্টতই স্বাভাবিকতার দিকে ধীরে ধীরে প্রবণতা রয়েছে, এমনকি (যেখানে এখনও বড় আকারের) সেখানে প্রশংসনীয় অ-স্বাভাবিকতা রয়ে গেছে, যেমন স্কিউনেসে প্রমাণিত ( এই নমুনায় সমান )। (এই হিস্টগ্রামের সঙ্কোচনের to এর খুব কাছাকাছি হওয়া কোনও আশ্চর্যের বিষয় নয় , কারণ এটি সিজিএফ-তে term টার্মটি হ'ল তাই)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)
}