সমাধান — অধ্যায় ৭.৬ · Independence, 0–1 Law ও SLLN¶
অধ্যায় ফাইল:
part-7-measure-theoretic/07-06-independence-zero-one-slln.md(§৭ অনুশীলনী)। গোটা অংশে \((\Omega,\mathcal F,\mathbb P)\) একটি probability space; random variable = measurable \(X:\Omega\to\mathbb R\) (7.3), \(\mathbb E[X]=\int_\Omega X\,d\mathbb P\) (7.4), \(X\in L^1\iff\mathbb E\lvert X\rvert<\infty\)। "i.o." = infinitely often; \(\limsup_n A_n=\bigcap_N\bigcup_{n\ge N}A_n=\{A_n\ \text{i.o.}\}\); tail σ-algebra \(\mathcal T=\bigcap_m\sigma(X_m,X_{m+1},\dots)\)।canonical তথ্য (সংখ্যাগত উত্তর reproducible, seed
np.random.default_rng(20260619)): - SLLN \(\text{Exp}(1)\) (\(\mathbb E[X]=1\)): \(n=10,10^2,10^3,10^4,10^5,10^6\to\bar X_n=0.5194,\,0.9710,\,0.9603,\,0.9841,\,0.9978,\,1.0007\) (→ \(1\) a.s.); \(\text{Bernoulli}(0.3)\): \(\bar X_{10^6}=0.2999\to0.3\)। - Cauchy (\(\mathbb E\lvert X\rvert=\infty\)): \(\bar X_n\) থামে না, ঘোরে — \(1.126,\,0.851,\,-0.173\) at \(n=10^2,10^4,10^6\) (চিহ্নও বদলায়)। - BC-I \(\mathbb P(A_n)=1/n^2\): \(\sum=\pi^2/6\approx1.6449<\infty\), ≈২ ঘটনা ঘটে (saturate), \(\mathbb P(\text{i.o.})=0\)। - BC-II \(\mathbb P(A_n)=1/n\): \(\sum=\infty\), ≈৮ ঘটনা (\(N=10^5\)), \(\sim\ln N\approx11.5\) হারে বাড়ে, \(\mathbb P(\text{i.o.})=1\)। - সীমারেখা \(\mathbb P(A_n)=1/(n\log n)\): \(\sum\sim\log\log N=\infty\) (অতি-ধীর divergent), \(\mathbb P(\text{i.o.})=1\)।
ক · ধারণাগত (conceptual)¶
সমাধান ১ (★)¶
(ক) pairwise independence। \(X_1,X_2\) independent fair coin, তাই \((X_1,X_2)\) চারটি মান \((0,0),(0,1),(1,0),(1,1)\) প্রতিটি probability \(\tfrac14\)-এ নেয়। \(X_3=X_1\oplus X_2\) (XOR): \(X_3=0\) যদি \(X_1=X_2\), নয়তো \(1\)। তিনটি জোড়া verify: - \(X_1,X_2\): সংজ্ঞা থেকেই independent। - \(X_1,X_3\): \(\mathbb P(X_1{=}a,X_3{=}c)\) — স্থির \(X_1=a\)-তে \(X_3=c\) মানে \(X_2=a\oplus c\), যার probability \(\tfrac12\); তাই \(\mathbb P(X_1{=}a,X_3{=}c)=\tfrac12\cdot\tfrac12=\tfrac14=\mathbb P(X_1{=}a)\mathbb P(X_3{=}c)\) (কারণ \(X_3\)-ও marginally fair: \(\mathbb P(X_3{=}c)=\tfrac12\))। independent। - \(X_2,X_3\): symmetry-তে একইভাবে independent।
তাই \(X_1,X_2,X_3\) জোড়ায়-জোড়ায় independent।
(খ) mutual independence ব্যর্থ। \(X_3\) প্রথম দুটোর deterministic ফাংশন (\(X_3=X_1\oplus X_2\)), তাই \(X_1,X_2\) জানলে \(X_3\) সম্পূর্ণ নির্ধারিত — কোনো বাড়তি randomness নেই। ফলে তিন-মুখী factorization ভাঙে; কংক্রিটভাবে:
অথচ গুণফল \(\mathbb P(X_1{=}0)\mathbb P(X_2{=}0)\mathbb P(X_3{=}1)=\tfrac12\cdot\tfrac12\cdot\tfrac12=\tfrac18\ne0\)। একটিমাত্র নির্বাচন ভাঙলেই পারস্পরিক independence যায়।
(গ) এক বাক্যে। pairwise (সব জোড়া independent) থেকে mutual (সব উপসেটে factorize) আসে না — ঠিক এই ফাঁক বন্ধ করতেই σ-algebra-র independence-সংজ্ঞা "যেকোনো \(A_i\in\mathcal F_i\) নির্বাচনে \(\mathbb P(\bigcap A_i)=\prod\mathbb P(A_i)\)" দাবি করে (যা \(\Omega\) বসিয়ে সব উপসেট ঢাকে), শুধু জোড়া নয়। (← 2.2)
সমাধান ২ (★)¶
(ক) tail event কী। একটি ঘটনা \(A\) tail event যদি \(A\in\mathcal T=\bigcap_m\sigma(X_m,X_{m+1},\dots)\) — equivalently, প্রতিটি \(m\)-এ \(A\) কেবল \(X_m,X_{m+1},\dots\)-এর উপর নির্ভর করে, অর্থাৎ প্রথম \(m-1\)টি \(X_i\) যা-খুশি বদলালেও \(A\)-র সত্য-মিথ্যা বদলায় না — সব \(m\)-এ। সংক্ষেপে: tail event "অসীম-দূরের লেজ"-এর তথ্য, কোনো সসীম উপসর্গে নির্ভরশীল নয়।
(খ) কোনগুলো tail। - (i) \(\{\sum_n X_n\ \text{converges}\}\) — tail। প্রথম \(m-1\)টি পদ বদলালে কেবল একটি সসীম যোগফল সরে, যা অভিসারিতা বদলায় না। - (ii) \(\{\limsup_n\bar X_n>c\}\) — tail। \(\limsup\) ও চলমান-গড়ের asymptotic আচরণ লেজ-নির্ভর; সসীম-সংখ্যক পদ বদলালে \(\bar X_n\)-এর limit-আচরণ অপরিবর্তিত। - (iii) \(\{X_1>0\}\) — tail নয়। শুধু \(X_1\)-এর উপর নির্ভর; \(X_1\) বদলালেই ঘটনা বদলায়, তাই \(\sigma(X_2,X_3,\dots)\)-এ নেই। - (iv) \(\{\sum_{n=1}^{100}X_n>5\}\) — tail নয়। কেবল \(X_1,\dots,X_{100}\)-এর (সসীম) ফাংশন; \(\sigma(X_{101},X_{102},\dots)\)-এ পড়ে না।
(গ) কেন \(0.5\) নয়। independent \((X_n)\)-এ Kolmogorov 0–1 law বলে \(\mathcal T\) নিজের থেকে independent (\(\mathcal T\perp\mathcal T\)), তাই যেকোনো \(A\in\mathcal T\)-এ \(\mathbb P(A)=\mathbb P(A\cap A)=\mathbb P(A)^2\), যার একমাত্র সমাধান \(\mathbb P(A)\in\{0,1\}\)। কাজেই \(0.5\) অসম্ভব — একটি সত্যিকারের asymptotic (tail) প্রশ্নের উত্তর নির্ধারিত (deterministic): হয় প্রায়-সব নমুনা-পথে ঘটে, নয় প্রায়-কোনোটিতেই নয়।
সমাধান ৩ (★★)¶
(ক) BC-I-এ independence লাগে না। চাবি-ধর্ম countable subadditivity (7.2)। \(B_N=\bigcup_{n\ge N}A_n\) লিখলে
আর \(\sum_n\mathbb P(A_n)<\infty\) বলে ডান পাশ (\(=\) converged ধারার লেজ) \(N\to\infty\)-এ \(0\)-তে নামে। তাই \(\mathbb P(A_n\ \text{i.o.})=\lim_N\mathbb P(B_N)=0\) — কেবল additivity ও convergence, কোনো independence নয়। স্বজ্ঞা: probability-এর যোগফল সসীম মানে গড়ে সসীম-সংখ্যক ঘটনা (\(\mathbb E[\sum_n\mathbf 1_{A_n}]=\sum_n\mathbb P(A_n)<\infty\)), তাই "অসীম-বার ঘটা" শূন্য-সম্ভাব্য।
(খ) BC-II-তে independence কোথায়। এখানে দেখাতে হয় "\(N\)-এর পর কিছুই ঘটে না" (\(C_N=\bigcap_{n\ge N}A_n^c\)) শূন্য-সম্ভাব্য। সসীম intersection-এ independence দিয়ে গুণফলে ভাঙা হয়:
যেখানে \(1-x\le e^{-x}\) inequality ও divergent লেজ \(\sum_{n\ge N}\mathbb P(A_n)=\infty\) ব্যবহৃত। independence ছাড়া এই factorization-টাই সম্ভব নয় — এটিই BC-I ও BC-II-এর গঠনগত পার্থক্য।
(গ) প্রতিউদাহরণ (independence ছাড়া BC-II মিথ্যা)। একই ঘটনা পুনরাবৃত্তি করুন: \(A_n=A\) সব \(n\)-এ, যেখানে \(0<\mathbb P(A)<1\) (যেমন \(\mathbb P(A)=0.1\))। তখন
কিন্তু \(\{A_n\ \text{i.o.}\}=\bigcap_N\bigcup_{n\ge N}A=A\) (প্রতিটি union ঠিক \(A\)), তাই \(\mathbb P(A_n\ \text{i.o.})=\mathbb P(A)=0.1\ne1\)। divergent যোগফল সত্ত্বেও i.o.-probability \(1\) নয় — কারণ ঘটনাগুলো পরস্পর-নির্ভরশীল (আসলে অভিন্ন)। এটি দেখায় BC-II-তে independence সত্যিই অপরিহার্য, অলংকার নয়।
খ · গণনামূলক (computational)¶
সমাধান ৪ (★)¶
প্রতিটি ক্ষেত্রে যোগফলের আচরণ verify করে Borel–Cantelli প্রয়োগ:
(ক) \(\mathbb P(A_n)=1/n^2\)। (i) \(\sum_{n\ge1}1/n^2=\pi^2/6\approx1.6449<\infty\) (Basel series — convergent)। (ii) যোগফল সসীম, তাই BC-I (independence লাগে না)। (iii) \(\mathbb P(A_n\ \text{i.o.})=0\) — a.s. কেবল সসীম-সংখ্যক \(A_n\) ঘটে।
(খ) \(\mathbb P(A_n)=1/n\)। (i) \(\sum_{n\ge1}1/n=\infty\) (harmonic series — divergent, বৃদ্ধি \(\sum_{n\le N}1/n\approx\ln N\))। (ii) divergent + independent, তাই BC-II। (iii) \(\mathbb P(A_n\ \text{i.o.})=1\) — a.s. অসীম-সংখ্যক \(A_n\) ঘটে।
(গ) \(\mathbb P(A_n)=1/(n\log n)\) (\(n\ge2\))। (i) integral test (সমাধান ৫ক): \(\int_2^\infty\frac{dx}{x\log x}=[\log\log x]_2^\infty=\infty\), তাই divergent (অতি-ধীর, আংশিক-যোগফল \(\sim\log\log N\))। (ii) divergent + independent, তাই BC-II। (iii) \(\mathbb P(A_n\ \text{i.o.})=1\)।
সারসংক্ষেপ-টেবিল:
| ক্ষেত্র | \(\mathbb P(A_n)\) | \(\sum\mathbb P(A_n)\) | lemma | \(\mathbb P(A_n\ \text{i.o.})\) |
|---|---|---|---|---|
| (ক) | \(1/n^2\) | \(\pi^2/6\approx1.6449\) (convergent) | BC-I | \(0\) |
| (খ) | \(1/n\) | \(\infty\) (\(\sim\ln N\)) | BC-II | \(1\) |
| (গ) | \(1/(n\log n)\) | \(\infty\) (\(\sim\log\log N\)) | BC-II | \(1\) |
লক্ষণীয়: \(1/n^2\) ও \(1/n\)-এর মাঝে সীমারেখা \(1/(n\log n)\) এখনও divergent দিকেই — তাই i.o. \(=1\), কিন্তু একদম গা ঘেঁষে।
সমাধান ৫ (★★)¶
(ক) integral test দিয়ে অপসারিতা। \(f(x)=\frac1{x\log x}\) (\(x\ge2\)) ধনাত্মক, হ্রাসমান, তাই integral test প্রযোজ্য। \(u=\log x\), \(du=\frac{dx}{x}\) substitution-এ:
integral divergent, তাই \(\sum_{n\ge2}\frac1{n\log n}\)-ও divergent। (তুলনায় \(\frac1{n(\log n)^2}\) convergent হত, কারণ \(\int\frac{du}{u^2}<\infty\) — exponent \(\log\)-এ ঠিক \(1\)-এর সীমারেখা।)
(খ) বৃদ্ধি-হার। integral test-এর তুলনা থেকে আংশিক-যোগফল \(\sum_{2\le n\le N}\frac1{n\log n}\approx\log\log N\) (অগ্রণী পদ)। আনুমানিক মান: - \(N=10^3\): \(\log\log 10^3=\log(3\ln 10)=\log(6.908)\approx1.93\)। - \(N=10^6\): \(\log\log 10^6=\log(6\ln 10)=\log(13.82)\approx2.63\)।
তুলনায় \(\sum 1/n\approx\ln N\) দেয় \(\ln 10^3\approx6.91\), \(\ln 10^6\approx13.82\)। অর্থাৎ \(N\) হাজার-গুণ বাড়লেও \(\log\log N\) মোটে \(1.93\to2.63\) (≈\(0.7\) বৃদ্ধি), যেখানে \(\ln N\) \(6.91\to13.82\) (দ্বিগুণ) — \(\log\log\) অসম্ভব ধীরে বাড়ে, তবু সীমাহীন।
(গ) উপসংহার। divergent যোগফল + independent ঘটনা ⇒ BC-II ⇒ \(\mathbb P(A_n\ \text{i.o.})=1\)। স্পেকট্রামে: \(1/n^2\) (যোগফল সসীম → i.o. \(0\), "থামে") আর \(1/n\) (যোগফল অসীম → i.o. \(1\), "চলে")-এর মধ্যে \(1/(n\log n)\) পড়ে "চলে" দিকে, যদিও সীমারেখার একদম গা ঘেঁষে — অপসারিতা যত ধীরই হোক, BC-II দ্বিধাহীনভাবে \(1\) দেয়।
সমাধান ৬ (★★)¶
fair-coin iid \(X_i\in\{0,1\}\), \(\mathbb P(X_i{=}1)=\tfrac12\), \(\mathbb E[X_i]=\tfrac12\), \(X_i\) আবদ্ধ তাই \(\mathbb E\lvert X_i\rvert=\tfrac12<\infty\) — সর্বত্র SLLN প্রযোজ্য।
(ক) \(\bar X_n=\frac1n\sum X_i\)। \(X_i\) iid, \(\mathbb E[X_i]=\tfrac12\), তাই সরাসরি SLLN:
(এটাই "আপেক্ষিক কম্পাঙ্ক → probability"-র সবল রূপ: মাথা-আসার হার \(\to\tfrac12\)।)
(খ) \(\frac1n\sum_{i=1}^n X_iX_{i+1}\)। \(Y_i:=X_iX_{i+1}\in\{0,1\}\), \(\mathbb E[Y_i]=\mathbb P(X_i{=}1,X_{i+1}{=}1)=\tfrac12\cdot\tfrac12=\tfrac14\) (\(X_i,X_{i+1}\) independent)। সীমা \(\tfrac14\) a.s.। সাবধানতা: \(Y_i\)-রা independent নয় (\(Y_i,Y_{i+1}\) overlap-করা \(X_{i+1}\) ভাগ করে), তাই মূল iid-SLLN সরাসরি লাগে না — কিন্তু \((Y_i)\) stationary ও ergodic, তাই ergodic theorem (SLLN-এর সাধারণ রূপ) দেয় \(\to\mathbb E[Y_1]=\tfrac14\)। প্রাথমিক যুক্তি: জোড়/বিজোড় সূচকে ভাঙুন — \(Y_1,Y_3,Y_5,\dots\) (অ-overlap) iid, তাই তাদের গড় \(\to\tfrac14\) a.s. (SLLN); একইভাবে \(Y_2,Y_4,\dots\); দুই উপ-গড়ের combination-ও \(\to\tfrac14\)। সীমা:
(গ) \(\frac1n\sum_{i=1}^n(2X_i-1)^2\)। \(X_i\in\{0,1\}\) বলে \(2X_i-1\in\{-1,+1\}\), তাই \((2X_i-1)^2=1\) সর্বদা — একটি ধ্রুবক random variable। গড় ঠিক \(1\) প্রতিটি \(n\)-এ, তাই সীমা:
(\(Z_i=1\) ধ্রুবক, \(\mathbb E[Z_i]=1<\infty\) — SLLN তুচ্ছভাবে মেটে।)
গ · প্রমাণভিত্তিক (proof-based)¶
সমাধান ৭ (★★)¶
Borel–Cantelli lemma I — পূর্ণ প্রমাণ।
(ক) monotone tail ও continuity-from-above। সংজ্ঞা থেকে \(\limsup_n A_n=\bigcap_{N\ge1}\big(\bigcup_{n\ge N}A_n\big)\)। লিখি \(B_N:=\bigcup_{n\ge N}A_n\)। বড় \(N\)-এ union-এ কম পদ, তাই
অর্থাৎ \(B_N\downarrow\limsup_n A_n\)। যেহেতু \(\mathbb P(B_1)\le1<\infty\) (probability measure সসীম), 7.2-এর continuity from above (continuity from above) বৈধ:
(খ) subadditivity দিয়ে বাঁধা। যেকোনো স্থির \(N\)-এ, 7.2-এর countable subadditivity সরাসরি দেয়:
(গ) লেজ \(\to0\) ও উপসংহার। ধরা শর্ত \(\sum_{n\ge1}\mathbb P(A_n)<\infty\) — converged অঋণাত্মক ধারা। যেকোনো converged ধারার লেজ-যোগফল শূন্যে নামে:
কারণ আংশিক-যোগফল পূর্ণ-যোগফলে converged। (2)-এ \(0\le\mathbb P(B_N)\le\sum_{n\ge N}\mathbb P(A_n)\to0\), তাই sandwich-এ \(\mathbb P(B_N)\to0\)। (1)-এ বসিয়ে \(\mathbb P(\limsup_n A_n)=0\)। \(\blacksquare\)
কোথাও independence ব্যবহৃত হয়নি — lemma সম্পূর্ণ সাধারণ, কেবল additivity ও convergence-এর উপর দাঁড়ানো।
সমাধান ৮ (★★★)¶
\(\{\sum_n X_n\ \text{converges}\}\) একটি tail event।
(ক) অভিসারিতা প্রথম \(N-1\) পদে নিরপেক্ষ। যেকোনো \(N\ge1\)-এর জন্য দুই ধারার পার্থক্য একটি সসীম যোগফল:
একটি সসীম সংখ্যা যোগ করলে অভিসারিতা বদলায় না (converged \(+\) সসীম \(=\) converged, এবং divergent \(+\) সসীম \(=\) divergent)। তাই
(খ) প্রতিটি \(\sigma(X_N,\dots)\)-এ, তাই \(\mathcal T\)-তে। আংশিক-যোগফল \(S_k^{(N)}:=\sum_{n=N}^k X_n\) (\(k\ge N\)) প্রতিটি \(X_N,\dots,X_k\)-এর measurable ফাংশন, তাই \(\sigma(X_N,X_{N+1},\dots)\)-measurable। অভিসারিতা মানে \(\lim_k S_k^{(N)}\)-এর অস্তিত্ব, যা measurable function-অনুক্রমের limit-অস্তিত্ব — একটি measurable ঘটনা \(\big(\{(S_k^{(N)})\ \text{Cauchy}\}=\bigcap_\varepsilon\bigcup_M\bigcap_{k,l\ge M}\{\lvert S_k^{(N)}-S_l^{(N)}\rvert<\varepsilon\}\big)\) — সম্পূর্ণরূপে \(\sigma(X_N,X_{N+1},\dots)\)-এ। \((\ast)\) বলছে ঘটনা \(\{\sum_n X_n\ \text{converges}\}\) ঠিক এটির সমান, তাই
সব \(N\)-এ সত্য বলে ছেদ:
অর্থাৎ এটি একটি tail event।
(গ) এক বাক্যে। \((X_n)\) independent বলে Kolmogorov 0–1 law প্রযোজ্য, তাই \(\mathbb P(\sum_n X_n\ \text{converges})\in\{0,1\}\) — অভিসারিতা একটি নির্ধারিত ব্যাপার (হয় a.s. converged, নয় a.s. নয়); একইভাবে \(\{\bar X_n\ \text{converges}\}\) tail, তাই সীমা থাকলে তা a.s. ধ্রুবক, আর SLLN শুধু সেই ধ্রুবকটি \(\mu\) বলে চিহ্নিত করে।
সমাধান ৯ (★★)¶
0–1 law-এর হৃৎপিণ্ড।
(ক) \(\mathbb P(A)=\mathbb P(A)^2\Rightarrow\mathbb P(A)\in\{0,1\}\)। \(A\) নিজের থেকে independent মানে \(\mathbb P(A\cap A)=\mathbb P(A)\mathbb P(A)\); কিন্তু \(A\cap A=A\), তাই \(\mathbb P(A)=\mathbb P(A)^2\)। ধরি \(p=\mathbb P(A)\):
(একটি বাস্তব সংখ্যা নিজের বর্গের সমান হয় কেবল \(0\) ও \(1\)-এ।)
(খ) কেন \(\mathcal T\perp\mathcal T\) (যুক্তি-শৃঙ্খল)। তিন ধাপ: 1. প্রতিটি \(n\)-এ \(\mathcal T\subseteq\sigma(X_{n+1},X_{n+2},\dots)\) (ছেদের একটি সদস্য)। আর \((X_k)\) independent বলে \(\sigma(X_1,\dots,X_n)\perp\sigma(X_{n+1},\dots)\), তাই \(\mathcal T\perp\sigma(X_1,\dots,X_n)\) — প্রতিটি \(n\)-এ। 2. \(\mathcal A:=\bigcup_n\sigma(X_1,\dots,X_n)\) একটি π-system (nested, তাই intersection-বদ্ধ) এবং \(\sigma(\mathcal A)=\sigma(X_1,X_2,\dots)\)। ধাপ-১ বলছে \(\mathcal T\) এই π-system \(\mathcal A\) থেকে independent; π-system criterion (প্রমাণ ১) দিয়ে \(\mathcal T\perp\sigma(X_1,X_2,\dots)\)। 3. কিন্তু \(\mathcal T\subseteq\sigma(X_1,X_2,\dots)\) (প্রতিটি \(\sigma(X_m,\dots)\) এর উপসেট, তাই ছেদও)। একটি σ-algebra তার ধারক থেকে independent থাকলে বিশেষত নিজের থেকেও: \(\mathcal T\perp\mathcal T\)।
(গ) tail RV a.s. ধ্রুবক। ধরা যাক \(Y\) একটি \(\mathcal T\)-measurable random variable (যেমন \(\limsup_n X_n\))। প্রতিটি \(t\)-এ \(\{Y\le t\}\in\mathcal T\), তাই (ক) দিয়ে \(F_Y(t)=\mathbb P(Y\le t)\in\{0,1\}\)। অর্থাৎ \(F_Y\) একটি অ-হ্রাসমান, ডান-continuous, \(\{0,1\}\)-মানের function — যা অবশ্যই একটি বিন্দু \(c=\inf\{t:F_Y(t)=1\}\)-তে \(0\) থেকে \(1\)-এ লাফ দেয়। তাই \(\mathbb P(Y=c)=F_Y(c)-F_Y(c^-)=1-0=1\) — \(Y\) a.s. ধ্রুবক (\(=c\))। \(\blacksquare\)
ঘ · কোডিং (coding)¶
সমাধান ১০ (★★)¶
SLLN convergence (\(\text{Exp}(1)\to1\)) বনাম Cauchy non-convergence — পাশাপাশি।
import numpy as np
# --- (ক) Exp(1): E|X| = 1 < ∞, তাই SLLN প্রযোজ্য — X̄_n → 1 ---
rng = np.random.default_rng(20260619)
X = rng.exponential(1.0, size=10**6) # সত্যিকারের গড় = 1
print("=== Exp(1): X̄_n → 1 (SLLN, E|X|=1<∞) ===")
for n in [10, 100, 10**3, 10**4, 10**5, 10**6]:
print(f" n={n:>7}: X̄_n = {X[:n].mean():.4f}")
# --- (খ) standard Cauchy: E|X| = ∞, তাই SLLN ভাঙে — X̄_n ঘোরে ---
rng = np.random.default_rng(20260619)
C = rng.standard_cauchy(size=10**6) # গড় অসংজ্ঞায়িত
print("\n=== Cauchy: X̄_n থামে না, ঘোরে (E|X|=∞) ===")
for n in [10**2, 10**4, 10**6]:
print(f" n={n:>7}: X̄_n = {C[:n].mean():.3f}")
# --- (গ) পার্থক্যের গোড়া: একটিমাত্র শর্ত E|X| < ∞ ---
# Exp-এ সত্য (গড় থামে), Cauchy-তে মিথ্যা (গড় ঘোরে)।
expected আউটপুট:
=== Exp(1): X̄_n → 1 (SLLN, E|X|=1<∞) ===
n= 10: X̄_n = 0.5194
n= 100: X̄_n = 0.9710
n= 1000: X̄_n = 0.9603
n= 10000: X̄_n = 0.9841
n= 100000: X̄_n = 0.9978
n=1000000: X̄_n = 1.0007
=== Cauchy: X̄_n থামে না, ঘোরে (E|X|=∞) ===
n= 100: X̄_n = 1.126
n= 10000: X̄_n = 0.851
n=1000000: X̄_n = -0.173
(ক) Exp(1)-এ ছোট \(n\)-এ গড় random (\(n=10\)-এ \(0.5194\)), কিন্তু \(n\) বাড়ার সঙ্গে দৃঢ়ভাবে \(1\)-এ গুটিয়ে আসে (\(n=10^6\)-এ \(1.0007\)) — convergence মসৃণ নয় (\(0.9710\to0.9603\) সামান্য নামে), কিন্তু সীমায় a.s. নিশ্চিত। (খ) Cauchy-তে কোনো convergence নেই: \(1.126\to0.851\to-0.173\) (\(n=10^2,10^4,10^6\)), এমনকি চিহ্ন বদলায় — heavy tail-এর বিশাল মান একটিমাত্র নমুনাতেই গড়কে টেনে সরায় (আসলে \(\bar X_n\) নিজেও হুবহু Cauchy থাকে, সরু হয় না)। (গ) পার্থক্যের গোড়ায় ঠিক একটি শর্ত: \(\mathbb E\lvert X\rvert<\infty\) — Exp-এ সত্য (\(=1\)), Cauchy-তে মিথ্যা (\(=\infty\))।
সমাধান ১১ (★★)¶
Borel–Cantelli I বনাম II — occurrence-count সিমুলেশন।
import numpy as np
N = 10**5
n = np.arange(1, N+1)
# --- (ক) BC-I: P(A_n)=1/n^2, probability-যোগফল সসীম → i.o. = 0 (saturate) ---
rng = np.random.default_rng(20260619)
occ1 = rng.random(N) < 1.0 / n**2
print("=== BC-I: P(A_n) = 1/n^2 (∑ = π²/6 ≈ 1.6449 < ∞) ===")
print(f" মোট ঘটনা সংখ্যা : {occ1.sum()} (≈ 2, saturate)")
print(f" শেষ ঘটনার সূচক : {np.where(occ1)[0].max()+1} (ছোট n-এ আটকে)")
# --- (খ) BC-II: P(A_n)=1/n, যোগফল অসীম, independent → i.o. = 1 (বাড়ন্ত) ---
rng = np.random.default_rng(20260619)
occ2 = rng.random(N) < 1.0 / n
print("\n=== BC-II: P(A_n) = 1/n (∑ = ∞, ~ln N) ===")
print(f" মোট ঘটনা সংখ্যা : {occ2.sum()} (≈ 8, saturate করে না)")
print(f" তাত্ত্বিক প্রবণতা ln(N) = {np.log(N):.2f}")
expected আউটপুট:
=== BC-I: P(A_n) = 1/n^2 (∑ = π²/6 ≈ 1.6449 < ∞) ===
মোট ঘটনা সংখ্যা : 2 (≈ 2, saturate)
শেষ ঘটনার সূচক : ... (ছোট n-এ আটকে)
=== BC-II: P(A_n) = 1/n (∑ = ∞, ~ln N) ===
মোট ঘটনা সংখ্যা : 8 (≈ 8, saturate করে না)
তাত্ত্বিক প্রবণতা ln(N) = 11.51
(ক) \(1/n^2\)-এ যোগফল সসীম (\(\pi^2/6\approx1.6449\)), তাই \(n\le10^5\)-এ মোটে ≈২টি \(A_n\) ঘটে, আর সব ছোট \(n\)-এ (কারণ \(1/n^2\) দ্রুত শূন্যে নামে) — গণনা একটি ছোট সংখ্যায় saturate করে, \(N\) বাড়ালেও বাড়ে না। (খ) \(1/n\)-এ যোগফল divergent, তাই গণনা saturate করে না — ইতিমধ্যে ≈৮টি, এবং expected গণনা \(\sum_{n\le N}1/n\approx\ln N\approx11.5\)-এর হারে একটানা বাড়ে। (গ) তুলনা-টেবিল:
| \(\mathbb P(A_n)\) | \(\sum\mathbb P(A_n)\) | independent? | count (\(N{=}10^5\)) | \(\mathbb P(\text{i.o.})\) | |
|---|---|---|---|---|---|
| BC-I | \(1/n^2\) | \(\pi^2/6\approx1.6449\) | লাগে না | ≈২ (saturate) | \(0\) |
| BC-II | \(1/n\) | \(\infty\) | লাগে | ≈৮ (বাড়ন্ত) | \(1\) |
দুই প্রতিবেশী ক্রম \(1/n^2\) ও \(1/n\) — একটির যোগফল সসীম, অপরটির অসীম; এই একমাত্র পার্থক্য (সাথে independence) i.o.-ভাগ্যকে \(0\) থেকে \(1\)-এ উল্টে দেয়, মাঝামাঝি নেই।
সমাধান ১২ (★★★)¶
একটি tail event-এর empirical probability \(\{0,1\}\)-এ জমাট (vs non-tail → \(0.5\))।
import numpy as np
rng = np.random.default_rng(20260619)
RUNS, N = 200, 10**5
tail_hits, nontail_vals = [], []
for _ in range(RUNS):
x = (rng.random(N) < 0.5).astype(float) # iid fair coin
# (ক) tail event C = {X̄_n → 1/2}: বড় n-এ |X̄_n - 1/2| < 0.01 ?
tail_hits.append(abs(x.mean() - 0.5) < 0.01)
# (খ) non-tail event {X_1 = 1}: কেবল প্রথম পদ
nontail_vals.append(x[0])
print("=== (ক) tail event {X̄_n → 1/2} (0–1 law ⇒ এখানে 1) ===")
print(f" সফল run-এর ভগ্নাংশ = {np.mean(tail_hits):.3f} (≈ 1.0)")
print("=== (খ) non-tail event {X_1 = 1} (0–1 law খাটে না) ===")
print(f" আপেক্ষিক কম্পাঙ্ক = {np.mean(nontail_vals):.3f} (≈ 0.5)")
expected আউটপুট:
=== (ক) tail event {X̄_n → 1/2} (0–1 law ⇒ এখানে 1) ===
সফল run-এর ভগ্নাংশ = 1.000 (≈ 1.0)
=== (খ) non-tail event {X_1 = 1} (0–1 law খাটে না) ===
আপেক্ষিক কম্পাঙ্ক = 0.495 (≈ 0.5)
(ক) \(C=\{\bar X_n\to\tfrac12\}\) একটি tail event (চলমান-গড়ের সীমা লেজ-নির্ভর, প্রথম পদে নয়); SLLN বলে \(\mathbb P(C)=1\), তাই প্রতিটি বড়-\(n\) run প্রায় নিশ্চিতভাবে \(\lvert\bar X_n-\tfrac12\rvert<0.01\) মেটায় — \(200\) run-এ সফল ভগ্নাংশ ≈১.০। (খ) \(\{X_1=1\}\) কেবল প্রথম পদের উপর নির্ভর — tail নয়, তাই 0–1 law খাটে না; এর probability ঠিক \(\tfrac12\) থাকে, empirical কম্পাঙ্ক ≈০.৫। (গ) পার্থক্যটা tail বনাম non-tail: tail event-এর probability independent অনুক্রমে \(\{0,1\}\)-এ বাঁধা (\(\mathbb P=\mathbb P^2\)), তাই এখানে \(1\)-এ জমাট; কিন্তু সসীম-উপসর্গ-নির্ভর non-tail ঘটনা সেই বাঁধনের বাইরে, ত