বিতরণ

একটি নিয়মিত অনুশীলন হিসাবে, আমি এর বিতরণ সন্ধান করার চেষ্টা করছি $\sqrt{X^2+Y^2}$ কোথায় $X$ এবং $Y$ স্বাধীন হয় $U(0,1)$ এলোমেলো ভেরিয়েবল।

এর যৌথ ঘনত্ব $(X,Y)$ হয়

f_{X, Y} (x, y) = 1_{0 < x, y < 1}

$f_{X,Y}(x,y)=\mathbf 1_{0<x,y<1}$

পোলার স্থানাঙ্কে রূপান্তর করা $(X,Y)\to(Z,\Theta)$ যেমন যে

X = Z \cos Θ and Y = Z \sin Θ

$X=Z\cos\Theta\qquad\text{ and }\qquad Y=Z\sin\Theta$

সুতরাং, $z=\sqrt{x^2+y^2}$ এবং $0< x,y<1\implies0< z<\sqrt 2$ ।

কখন $0< z<1$ , আমাদের আছে $0< \cos\theta<1,\,0<\sin\theta<1$ যাতে $0<\theta<\frac{\pi}{2}$ ।

কখন $1< z<\sqrt 2$ , আমাদের আছে $z\cos\theta<\implies\theta>\cos^{-1}\left(\frac{1}{z}\right)$ যেমন $\cos\theta$ কমছে $\theta\in\left[0,\frac{\pi}{2}\right]$ ; এবং $z\sin\theta<1\implies\theta<\sin^{-1}\left(\frac{1}{z}\right)$ যেমন $\sin\theta$ ক্রমবর্ধমান হয় $\theta\in\left[0,\frac{\pi}{2}\right]$ ।

অনেক দূরে $1< z<\sqrt 2$ , আমাদের আছে $\cos^{-1}\left(\frac{1}{z}\right)<\theta<\sin^{-1}\left(\frac{1}{z}\right)$ ।

রূপান্তরের জ্যাকোবিয়ানের পরম মান

| J | = z

$|J|=z$

সুতরাং যৌথ ঘনত্ব $(Z,\Theta)$ দেওয়া হয়

f_{Z, Θ} (z, θ) = z 1_{{z \in (0, 1), θ \in (0, π / 2)} ⋃ {z \in (1, \sqrt{2}), θ \in (\cos^{- 1} (1 / z), \sin^{- 1} (1 / z))}}

$f_{Z,\Theta}(z,\theta)=z\mathbf 1_{\{z\in(0,1),\,\theta\in\left(0,\pi/2\right)\}\bigcup\{z\in(1,\sqrt2),\,\theta\in\left(\cos^{-1}\left(1/z\right),\sin^{-1}\left(1/z\right)\right)\}}$

একীকরণ আউট $\theta$ , আমরা এর পিডিএফ প্রাপ্ত $Z$ যেমন

f_{Z} (z) = \frac{π z}{2} 1_{0 < z < 1} + (\frac{π z}{2} - 2 z \cos^{- 1} (\frac{1}{z})) 1_{1 < z < \sqrt{2}}

$f_Z(z)=\frac{\pi z}{2}\mathbf 1_{0<z<1}+\left(\frac{\pi z}{2}-2z\cos^{-1}\left(\frac{1}{z}\right)\right)\mathbf 1_{1<z<\sqrt 2}$

আমার যুক্তি কি উপরে? যাই হোক না কেন, আমি এই পদ্ধতিটি এড়াতে এবং পরিবর্তে সিডিএফ সন্ধান করার চেষ্টা করব $Z$ সরাসরি। তবে মূল্যায়ন করার সময় আমি পছন্দসই অঞ্চলগুলি খুঁজে পাইনি $\mathrm{Pr}(Y\le \sqrt{z^2-X^2})$ জ্যামিতিক।

সম্পাদনা করুন।

এর বিতরণ ফাংশনটি সন্ধান করার চেষ্টা করেছি $Z$ যেমন

\begin{aligned} F_{Z} (z) & = Pr (Z \leq z) \\ = Pr (X^{2} + Y^{2} \leq z^{2}) \\ = \iint_{x^{2} + y^{2} \leq z^{2}} 1_{0 < x, y < 1} d x d y \end{aligned}

$\begin{align} F_Z(z)&=\Pr(Z\le z) \\&=\Pr(X^2+Y^2\le z^2) \\&=\iint_{x^2+y^2\le z^2}\mathbf1_{0<x,y<1}\,\mathrm{d}x\,\mathrm{d}y \end{align}$

ম্যাথমেটিকা বলেছেন যে এটি হ্রাস করা উচিত

F_{Z} (z) = {\begin{cases} 0 & , if z < 0 \\ \frac{π z^{2}}{4} & , if 0 < z < 1 \\ \sqrt{z^{2} - 1} + \frac{z^{2}}{2} (\sin^{- 1} (\frac{1}{z}) - \sin^{- 1} (\frac{\sqrt{z^{2} - 1}}{z})) & , if 1 < z < \sqrt{2} \\ 1 & , if z > \sqrt{2} \end{cases}

$F_Z(z)=\begin{cases}0 &,\text{ if }z<0\\ \frac{\pi z^2}{4} &,\text{ if } 0< z<1\\ \sqrt{z^2-1}+\frac{z^2}{2}\left(\sin^{-1}\left(\frac{1}{z}\right)-\sin^{-1}\left(\frac{\sqrt{z^2-1}}{z}\right)\right) &,\text{ if }1< z<\sqrt 2\\ 1 &,\text{ if }z>\sqrt 2 \end{cases}$

যা দেখতে সঠিক মত প্রকাশের মতো। মামলার ক্ষেত্রে পার্থক্য করা যদিও এমন একটি অভিব্যক্তি উপস্থিত করে যা ইতিমধ্যে আমি যে পিডিএফটি পেয়েছি তা সহজেই সরল করে না। $F_Z$ $1< z<\sqrt 2$

পরিশেষে, আমি মনে করি সিডিএফের জন্য আমার কাছে সঠিক ছবি রয়েছে:

জন্য : $0<z<1$

এবং : $1<z<\sqrt 2$

শেড অংশগুলি অঞ্চলটির ক্ষেত্রফল নির্দেশ করবে বলে মনে করা হচ্ছে

{(x, y) : 0 < x, y < 1, x^{2} + y^{2} \leq z^{2}}

$\left\{(x,y):0<x,y< 1\,,\,x^2+y^2\le z^2\right\}$

ছবিটি সঙ্গে সঙ্গে ফলন দেয়

\begin{aligned} F_{Z} (z) & = Pr (- \sqrt{z^{2} - X^{2}} \leq Y \leq \sqrt{z^{2} - X^{2}}) \\ = {\begin{cases} \frac{π z^{2}}{4} & , if 0 < z < 1 \\ \sqrt{z^{2} - 1} + \int_{\sqrt{z^{2} - 1}}^{1} \sqrt{z^{2} - x^{2}} d x & , if 1 < z < \sqrt{2} \end{cases} \end{aligned}

$\begin{align} F_Z(z)&=\Pr\left(-\sqrt{z^2-X^2}\le Y\le\sqrt{z^2-X^2}\right) \\&=\begin{cases}\frac{\pi z^2}{4} &,\text{ if } 0<z<1\\\\ \sqrt{z^2-1}+\int_{\sqrt{z^2-1}}^1 \sqrt{z^2-x^2}\,\mathrm{d}x &,\text{ if }1< z<\sqrt 2 \end{cases} \end{align}$

যেমনটি আমি আগে খুঁজে পেয়েছি।

— StubbornAtom
সূত্র

সিডিএফ সরাসরি খুঁজে পেতে, সূচক ফাংশন ব্যবহার করুন। জন্যবাকিগুলি খাঁটি বীজগণিত কারসাজি। (সম্পাদনা করুন: আমি @

z \geq 0,

$z\ge 0,$

Pr (\sqrt{X^{2} + Y^{2}} \leq z) = \int_{0}^{1} \int_{0}^{1} I (x^{2} + y^{2} \leq z^{2}) d x d y .

$\Pr(\sqrt{X^2+Y^2}\le z)=\int_0^1\int_0^1\mathcal{I}(x^2+y^2\le z^2)\mathrm{d}x\mathrm{d}y.$

— শি'ান কেবলমাত্র

সম্পাদনাটি পুনরায়: আমি বেশ কয়েকটি পৃথক এক্সপ্রেশনও পেয়েছি এবং (ব্যবহার করে FullSimplify) তারা গণিতের বিভিন্ন সূত্রে সরল করে । তবে এগুলি সমতুল্য। এটি তাদের পার্থক্যটি চক্রান্ত করে সহজেই প্রদর্শিত হয়। স্পষ্টতই গাণিতিক জানে না যে যখন ।

\tan^{- 1} (\sqrt{z^{2} - 1}) = \sec^{- 1} (z)

$\tan ^{-1}\left(\sqrt{z^2-1}\right)=\sec ^{-1}(z)$

1 < z < \sqrt{2}

$1\lt z \lt \sqrt{2}$

— হুবুহু

পৃষ্ঠের প্রান্ত, , আপনার শেষ ছবিতে কেন্দ্র (0,0) সহ একটি (আধা-) বৃত্ত হওয়া উচিত। সুতরাং (আপনার বর্তমানে টানা) উত্তল পরিবর্তে অবতল।

\sqrt{r^{2} - x^{2}}

$\sqrt{r^2-x^2}$

— সেক্সটাস এম্পেরিকাস

উত্তর:

যে পিডিএফ সঠিক তা সাধারণ সিমুলেশন দ্বারা পরীক্ষা করা যায়

samps=sqrt(runif(1e5)^2+runif(1e5)^2)
hist(samps,prob=TRUE,nclass=143,col="wheat")
df=function(x){pi*x/2-2*x*(x>1)*acos(1/(x+(1-x)*(x<1)))}
curve(df,add=TRUE,col="sienna",lwd=3)

ভেরিয়েবলের মেরু পরিবর্তন ছাড়াই সিডিএফ অনুসন্ধান করা যায়

\begin{aligned} P r (\sqrt{X^{2} + Y^{2}} \leq z) & = P r (X^{2} + Y^{2} \leq z^{2}) \\ = P r (Y^{2} \leq z^{2} - X^{2}) \\ = P r (Y \leq \sqrt{z^{2} - X^{2}}, X \leq z) \\ = E^{X} [\sqrt{z^{2} - X^{2}} I_{[0, min (1, z)]} (X)] \\ = \int_{0}^{min (1, z)} \sqrt{z^{2} - x^{2}} d x \\ = z^{2} \int_{0}^{min (1, z^{- 1})} \sqrt{1 - y^{2}} d y [x = y z, d x = z d y] \\ = z^{2} \int_{0}^{min (π / 2, \cos^{- 1} z^{- 1})} \sin^{2} θ d θ [y = \cos (θ), d y = \sin (θ) d θ] \\ = \frac{z^{2}}{2} [min (π / 2, \cos^{- 1} z^{- 1}) - \sin {min (π / 2, \cos^{- 1} z^{- 1})} \cos {min (π / 2, \cos^{- 1} z^{- 1}}] \\ = \frac{z^{2}}{2} {\begin{cases} π / 2 & if z < 1 \\ \cos^{- 1} z^{- 1} - \sin {\cos^{- 1} z^{- 1})} z^{- 1} & if z \geq 1 \end{cases} \\ = \frac{z^{2}}{2} {\begin{cases} π / 2 & if z < 1 \\ \cos^{- 1} z^{- 1} - \sqrt{1 - z^{- 2}} z^{- 1} & if z \geq 1 \end{cases} \end{aligned}

$\begin{align*} \mathrm{Pr}(\sqrt{X^2+Y^2}\le z) &= \mathrm{Pr}(X^2+Y^2\le z^2)\\ &= \mathrm{Pr}(Y^2\le z^2-X^2)\\ &=\mathrm{Pr}(Y\le \sqrt{z^2-X^2}\,,X\le z)\\ &=\mathbb{E}^X[\sqrt{z^2-X^2}\mathbb{I}_{[0,\min(1,z)]}(X)]\\ &=\int_0^{\min(1,z)} \sqrt{z^2-x^2}\,\text{d}x\\ &=z^2\int_0^{\min(1,z^{-1})} \sqrt{1-y^2}\,\text{d}y\qquad [x=yz\,,\ \text{d}x=z\text{d}y]\\ &=z^2\int_0^{\min(\pi/2,\cos^{-1} z^{-1})} \sin^2{\theta} \,\text{d}\theta\qquad [y=\cos(\theta)\,,\ \text{d}y=\sin(\theta)\text{d}\theta]\\ &=\frac{z^2}{2}\left[ \min(\pi/2,\cos^{-1} z^{-1}) - \sin\{\min(\pi/2,\cos^{-1} z^{-1})\}\cos\{\min(\pi/2,\cos^{-1} z^{-1}\}\right]\\ &=\frac{z^2}{2}\begin{cases} \pi/2 &\text{ if }z<1\\ \cos^{-1} z^{-1}-\sin\{\cos^{-1} z^{-1})\}z^{-1}&\text{ if }z\ge 1\\ \end{cases}\\ &=\frac{z^2}{2}\begin{cases} \pi/2 &\text{ if }z<1\\ \cos^{-1} z^{-1}-\sqrt{1-z^{-2}}z^{-1}&\text{ if }z\ge 1\\ \end{cases} \end{align*}$ একই জটিলতায় শেষ হয়! (পথে আমার সম্ভাব্য ভুল!)

— সিয়ান
সূত্র

কেস যেখানে এটি কিছুটা अस्पष्ट হয়ে যায়। আমার ধারণা আমি জন্য সঠিক পিডিএফ এক্সপ্রেশনকে পৃথক করে শেষ করব না ।

1 \leq z < \sqrt{2}

$1\le z<\sqrt 2$

z \geq 1

$z\ge 1$

— জেদীআটম

$f_z(z)$ :

সুতরাং, , আমাদের কাছে $1\le z<\sqrt 2$ $\cos^{-1}\left(\frac{1}{z}\right)\le\theta\le\sin^{-1}\left(\frac{1}{z}\right)$

যখন আপনি প্রতিসাম্য ব্যবহার এবং জন্য চিহ্নগুলির মান নির্ণয়ের তুমি তোমার এক্সপ্রেশন প্রক্রিয়া সহজ করতে । এইভাবে, অর্ধেক জায়গার জন্য এবং তারপরে ফলাফল দ্বিগুণ করুন। $\theta_{min} < \theta < \frac{\pi}{4}$

তারপরে আপনি পাবেন:

P (Z \leq r) = 2 \int_{0}^{r} z (\int_{θ_{m i n}}^{\frac{π}{4}} d θ) d z = \int_{0}^{r} z (\frac{π}{2} - 2 θ_{m i n}) d z

$P(Z \leq r) = 2 \int_0^r z \left(\int_{\theta_{min}}^{\frac{\pi}{4}}d\theta\right) dz = \int_0^r z \left(\frac{\pi}{2}-2\theta_{min}\right) dz$

এবং আপনার হয় $f_z(z)$

f_{z} (z) = z (\frac{π}{2} - 2 θ_{m i n}) = {\begin{cases} z (\frac{π}{2}) & if 0 \leq z \leq 1 \\ z (\frac{π}{2} - 2 \cos^{- 1} (\frac{1}{z})) & if 1 < z \leq \sqrt{2} \end{cases}

$f_z(z) = z \left(\frac{\pi}{2}-2\theta_{min}\right) = \begin{cases} z\left(\frac{\pi}{2}\right) & \text{ if } 0 \leq z \leq 1 \\ z \left(\frac{\pi}{2} - 2 \cos^{-1}\left(\frac{1}{z}\right)\right) & \text{ if } 1 < z \leq \sqrt{2} \end{cases}$

$F_z(z)$ :

আপনি অনির্দিষ্ট অবিচ্ছেদ্য ব্যবহার করতে পারেন:

\int z \cos^{- 1} (\frac{1}{z}) = \frac{1}{2} z (z \cos^{- 1} (\frac{1}{z}) - \sqrt{1 - \frac{1}{z^{2}}}) + C

$\int z \cos^{-1}\left(\frac{1}{z}\right) = \frac{1}{2} z \left( z \cos^{-1}\left(\frac{1}{z}\right) - \sqrt{1-\frac{1}{z^2}} \right) + C$

দ্রষ্টব্য $\frac{d}{du} \cos^{-1}(u) = - (1-u^2)^{-0.5}$

এই বিশালাকার জন্য Xi'ans অভিব্যক্তি হিসাবে অনুরূপ কিছু সহজবোধ্য যথা $Pr(Z \leq z)$

যদি তবে: $1 \leq z \leq \sqrt{2}$

F_{z} (z) = z^{2} (\frac{π}{4} - \cos^{- 1} (\frac{1}{z}) + z^{- 1} \sqrt{1 - \frac{1}{z^{2}}})

$F_z(z) = {z^2} \left(\frac{\pi}{4}-\cos^{-1}\left(\frac{1}{z}\right) + z^{-1}\sqrt{1-\frac{1}{z^2}} \right)$

আপনার অভিব্যক্তির সাথে সম্পর্কটি দেখা যায় যখন আমরা কে দুটি এক্সপ্রেশনে বিভক্ত করি এবং তারপরে বিভিন্ন এক্সপ্রেশনগুলিতে রূপান্তর করি । $cos^{-1}$ $cos^{-1}$ $sin^{-1}$

জন্য আছে $z>1$

\cos^{- 1} (\frac{1}{z}) = \sin^{- 1} (\sqrt{1 - \frac{1}{z^{2}}}) = \sin^{- 1} (\frac{\sqrt{z^{2} - 1}}{z})

$\cos^{-1}\left(\frac{1}{z}\right) = \sin^{-1}\left(\sqrt{1-\frac{1}{z^2}}\right) = \sin^{-1}\left(\frac{\sqrt{z^2-1}}{z}\right)$

এবং

\cos^{- 1} (\frac{1}{z}) = \frac{π}{2} - \sin^{- 1} (\frac{1}{z})

$\cos^{-1}\left(\frac{1}{z}\right) = \frac{\pi}{2} -\sin^{-1}\left(\frac{1}{z}\right)$

সুতরাং

\begin{matrix} \cos^{- 1} (\frac{1}{z}) & = 0.5 \cos^{- 1} (\frac{1}{z}) + 0.5 \cos^{- 1} (\frac{1}{z}) \\ = \frac{π}{4} - 0.5 \sin^{- 1} (\frac{1}{z}) + 0.5 \sin^{- 1} (\frac{\sqrt{z^{2} - 1}}{z}) \end{matrix}

$\begin{array}\\ \cos^{-1}\left(\frac{1}{z}\right) & = 0.5 \cos^{-1}\left(\frac{1}{z}\right) + 0.5 \cos^{-1}\left(\frac{1}{z}\right) \\ & = \frac{\pi}{4} - 0.5 \sin^{-1}\left(\frac{1}{z}\right) + 0.5 \sin^{-1}\left(\frac{\sqrt{z^2-1}}{z}\right) \end{array}$

আপনি যখন জন্য উল্লিখিত এটি প্লাগ করেন তখন আপনার অভিব্যক্তির ফলাফল হয় $F_z(z)$ $1<z<\sqrt{2}$

— সেক্সটাস এম্পেরিকাস
সূত্র

জন্য , মাত্র ব্যাসার্ধ কোয়ার্টার বৃত্তের এলাকা যা । । এটি, $0 \leq z \leq 1$ $P\left(\sqrt{X^2+Y^2} \leq z\right)$ $z$ $\frac 14 \pi z^2$

For 0 \leq z \leq 1, area of quarter-circle = \frac{π z^{2}}{4} = P (\sqrt{X^{2} + Y^{2}} \leq z) .

$\text{For }0 \leq z \leq 1, ~\text{area of quarter-circle} = \frac{\pi z^2}{4} = P\left(\sqrt{X^2+Y^2} \leq z\right).$

জন্য , অঞ্চল যার উপর আমরা এটি সংহত করতে প্রয়োজন দুই ডান ত্রিভুজ ভাগে ভাগ করা যায় একটিতে উল্লম্ব এবং রয়েছে অন্যটির এবং একসাথে ব্যাসার্ধের বৃত্তের একটি ক্ষেত্র এবং এর সাথে কোণ । এই অঞ্চলের ক্ষেত্রফল (এবং সেইজন্য মান ) সহজেই পাওয়া যায় We $1 < z \leq \sqrt{2}$ $P\left(\sqrt{X^2+Y^2} \leq z\right)$ $\big($ $(0,0), (0,1)$ $(\sqrt{z^2-1}, 1)$ $(0,0), (1,0)$ $(1, \sqrt{z^2-1})$ $\big)$ $z$ $\frac{\pi}{2}-2\arccos\left(\frac{1}{z}\right)$ $\left( P(\sqrt{X^2+Y^2} \leq z\right)$ $1 < z \leq \sqrt{2}$ , যা মার্তিজান ওয়েটারিংয়ের উত্তরের ফলাফল।

\begin{aligned} area of region & = area of two triangles plus area of sector \\ = \sqrt{z^{2} - 1} + \frac{1}{2} z^{2} (\frac{π}{2} - 2 \arccos (\frac{1}{z})) \\ = \frac{π z^{2}}{4} + \sqrt{z^{2} - 1} - z^{2} \arccos \frac{1}{z} \\ = (P (\sqrt{X^{2} + Y^{2}} \leq z) \end{aligned}

$\begin{align}\text{area of region} &= \text{area of two triangles plus area of sector}\\ &=\sqrt{z^2-1} + \frac 12 z^2\left( \frac{\pi}{2}-2\arccos \left(\frac{1}{z}\right)\right)\\ &= \frac{\pi z^2}{4} + \sqrt{z^2-1} - z^2\arccos \frac{1}{z}\\ &= \left( P(\sqrt{X^2+Y^2} \leq z\right)\end{align}$

— দিলীপ সরওয়াতে
সূত্র

বিতরণ

fz(z)fz(z)f_z(z) :

Fz(z)Fz(z)F_z(z) :

$f_z(z)$ :

$F_z(z)$ :