What is this curve?

What is this curve?



Lines (same angle space between) radiating outward from a point and intersecting a line:



Intersection Point Density Distribution



This is the density distribution of the points on the line:



Plot line Graph



I used a Python script to calculate this. The angular interval is 0.01 degrees. X = tan(deg) rounded to the nearest tenth, so many points will have the same x. Y is the number of points for its X. You can see on the plot that ~5,700 points were between -0.05 and 0.05 on the line.



What is this graph curve called? What's the function equation?





$begingroup$
What are the numbers in your graph supposed to represent?
$endgroup$
– Aretino
Sep 9 '18 at 20:42





$begingroup$
There are a few functions you could make look roughly like this, but I think its going to be of the form $amathrme^-bx^2$
$endgroup$
– aidangallagher4
Sep 9 '18 at 20:52





$begingroup$
en.wikipedia.org/wiki/Cauchy_distribution
$endgroup$
– Aretino
Sep 9 '18 at 20:54





$begingroup$
mathworld.wolfram.com/LorentzianFunction.html
$endgroup$
– Aretino
Sep 9 '18 at 20:58




3 Answers
3



The curve is a density function. The idea is the following. From your first picture, assume that the angles of the rays are spaced evenly, the angle between two rays being $alpha$. I.e. the nth ray has angle $n cdot alpha$. The nth ray's intersection point $x$ with the line then follows $tan (n cdot alpha) = x/h$ where h is the distance of the line to the origin. So from 0 to $x$, n rays cross the line.



Now you are interested in the density $p(x)$, i.e. how many rays intersect the line at $x$, per line interval $Delta x$. In the limit of small $alpha$, you have $int_0^x p(x') dx' =c n = fraccalphaarctan (x/h)$ and correspondingly, $p(x) = fracddxfraccalphaarctan (x/h) = fracc halpha (x^2+h^2)$. The constant $c$ is determined since the integral over the density function must be $1$ (in probability sense), hence



$p(x)= frachpi (x^2+h^2)$.



This curve is called a Cauchy distribution.



Obviously, $p(x)$ can be multiplied with a constant $K$ to give an expectation value distribution $E(x) = K p(x)$ over $x$, instead of a probability distribution. This explains the large value of $E(0) = 5700$ or so in your picture. The value $h$ is also called a scale parameter, it specifies the half-width at half-maximum (HWHM) of the curve and can be roughly read off to be $1$. If we are really "counting rays", then with angle spacing $alpha$, in total $pi/alpha$ many rays intersect the line and hence we must have
$$
pi/alpha = int_-infty^inftyE(x) dx = int_-infty^inftyK p(x) dx = K
$$
So the expectation value distribution of the number of rays intersecting one unit of the line at position $x$ is
$$
E(x) = fracpialpha p(x)= frachalpha (x^2+h^2)
$$
as we already had with the constant $c=1$. Reading off approximately $E(0) = 5700$ and $h=1$ gives $
E(x) = frac5700x^2+1
$ and $alpha = 1/5700$ (in rads), or in other words, $pi/alpha simeq 17900$ rays (lower half plane) intersect the line.



It can be shown assuming that the quantity from the ray is proportional to the angular interval that the denisty distribuiton function is in the form



$$y=fracdH+x^2$$



where $d$ is related to the intensity and $H$ is the distance of the source point from the line.



Here is the plot for the case $d=H=1$



enter image description here



This is also known, especially among we physicists as a Lorentz distribution.



We know that every point on the line has the same vertical distance from the source point. Let's call this $l_0$. We also know the ratio of the horizontal component of the distance to the vertical component is $tan(theta)$, where theta goes from $-piover2$ to $piover2$ in the increments given. So horizontal distance=$x=l_0*tan(theta)$.



To find the densities, first take the derivative of both sides, getting $dx=l_0*sec^2(theta)*dtheta$. From trig, we know $sec^2(theta)=1+tan^2(theta)$, and from above, we know that $tan(theta)$ is $xoverl_0$. So we can isolate $theta$ to get



$$dx over l_0* left( 1+left(xoverl_0 right)^2 right) =dtheta$$



We know theta is uniformly distributed, so we can divide both sides by $pi$. Now the probability of a given $theta$ falling between $theta$ and $theta+d theta$ is equal to the function of $x$ on the left.





$begingroup$
To clear up any possible confusion, the Lorentz distribution is an alternative name for the Cauchy distribution.
$endgroup$
– Ingolifs
Sep 10 '18 at 18:26



Thanks for contributing an answer to Mathematics Stack Exchange!



But avoid



Use MathJax to format equations. MathJax reference.



To learn more, see our tips on writing great answers.



Required, but never shown



Required, but never shown




By clicking "Post Your Answer", you acknowledge that you have read our updated terms of service, privacy policy and cookie policy, and that your continued use of the website is subject to these policies.

Popular posts from this blog

𛂒𛀶,𛀽𛀑𛂀𛃧𛂓𛀙𛃆𛃑𛃷𛂟𛁡𛀢𛀟𛁤𛂽𛁕𛁪𛂟𛂯,𛁞𛂧𛀴𛁄𛁠𛁼𛂿𛀤 𛂘,𛁺𛂾𛃭𛃭𛃵𛀺,𛂣𛃍𛂖𛃶 𛀸𛃀𛂖𛁶𛁏𛁚 𛂢𛂞 𛁰𛂆𛀔,𛁸𛀽𛁓𛃋𛂇𛃧𛀧𛃣𛂐𛃇,𛂂𛃻𛃲𛁬𛃞𛀧𛃃𛀅 𛂭𛁠𛁡𛃇𛀷𛃓𛁥,𛁙𛁘𛁞𛃸𛁸𛃣𛁜,𛂛,𛃿,𛁯𛂘𛂌𛃛𛁱𛃌𛂈𛂇 𛁊𛃲,𛀕𛃴𛀜 𛀶𛂆𛀶𛃟𛂉𛀣,𛂐𛁞𛁾 𛁷𛂑𛁳𛂯𛀬𛃅,𛃶𛁼

Edmonton

Crossroads (UK TV series)