<?xml version="1.0" encoding="UTF-8"?>
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	xmlns:atom="http://www.w3.org/2005/Atom"
	xmlns:sy="http://purl.org/rss/1.0/modules/syndication/"
	xmlns:slash="http://purl.org/rss/1.0/modules/slash/"
	>

<channel>
	<title>Some stuff &#187; theta</title>
	<atom:link href="http://blog.yhuang.org/?feed=rss2&#038;tag=theta" rel="self" type="application/rss+xml" />
	<link>https://blog.yhuang.org</link>
	<description>here.</description>
	<lastBuildDate>Wed, 27 Aug 2025 08:50:58 +0000</lastBuildDate>
	<language>en</language>
	<sy:updatePeriod>hourly</sy:updatePeriod>
	<sy:updateFrequency>1</sy:updateFrequency>
	<generator>http://wordpress.org/?v=3.1.1</generator>
		<item>
		<title>road path problem</title>
		<link>https://blog.yhuang.org/?p=242</link>
		<comments>https://blog.yhuang.org/?p=242#comments</comments>
		<pubDate>Fri, 19 Feb 2010 21:45:17 +0000</pubDate>
		<dc:creator>admin</dc:creator>
				<category><![CDATA[Uncategorized]]></category>
		<category><![CDATA[circle]]></category>
		<category><![CDATA[continuous path]]></category>
		<category><![CDATA[line]]></category>
		<category><![CDATA[math]]></category>
		<category><![CDATA[math 1]]></category>
		<category><![CDATA[pi math]]></category>
		<category><![CDATA[road]]></category>
		<category><![CDATA[secant graph]]></category>
		<category><![CDATA[theta]]></category>
		<category><![CDATA[unit radius]]></category>

		<guid isPermaLink="false">http://scripts.mit.edu/~zong/wpress/?p=242</guid>
		<description><![CDATA[Suppose there is a straight road, infinitely long at both ends, located 1 unit from your starting location. Find the most efficient path to reach the road, and the worst-case total length of this path. The trivial but wrong way is to go for 1 unit in some direction, then trace the circumference of a [...]]]></description>
			<content:encoded><![CDATA[<p>Suppose there is a straight road, infinitely long at both ends, located 1 unit from your starting location. Find the most efficient path to reach the road, and the worst-case total length of this path.</p>
<p>The trivial but wrong way is to go for 1 unit in some direction, then trace the circumference of a unit-radius circle. The road will surely be found this way, but the path length is \(1+2\pi\), which can be improved upon.<br />
<span id="more-242"></span><br />
Every possible position for the road corresponds to a tangent line to the unit circle centered at the starting location. Therefore, this is a problem about finding a shortest path that touches all tangent lines to a unit circle. To make this more concrete, let us transform the space into polar form as below.</p>
<p><img src="wp-content/uploads/images/roadpath2.png" /></p>
<p>We see that in polar form, the unit circle becomes the line \(r=1\) and each tangent line touching the circle at angle \(\theta_0\) becomes the graph \(r=\sec(\theta-\theta_0)\). The trivial solution corresponds to a path that traces \(r=0\) to \(r=1\) at \(\theta=0\), followed by \(\theta=0\) to \(\theta=2\pi\) at \(r=1\). This is not efficient.</p>
<p>In fact, any continuous path starting at \(r=0\) &#8220;reaches the road&#8221; (i.e. touches every tangent line) as long as it touches \(r=\sec(\theta)\) somewhere on \(\theta\in [0,\pi/2)\) and somewhere on \(\theta\in (3\pi/2,2\pi]\), and it stays above \(r=1\) between those two points. And up to shifts in \(\theta\), these are the only possible solutions.* Note that the path doesn&#8217;t need to be closed, like the trivial path was. So we are left to find the shortest such path (e.g. the one in green above).</p>
<p>Whereas in rectangular form, the shortest path between two points is a straight line, in polar form, the shortest path between two points is along a secant graph. It is a fairly easy matter to show that the shortest path must take the following form:</p>
<p>1. Start from \((r,\theta)=(0,\phi)\) for some \(\phi\in [0,\pi/2)\);<br />
2. Drop onto the graph \(r=\sec(\theta-2\phi)\), between the points \((r,\theta)=(\sec(\phi),\phi)\) and \((r,\theta)=(1,2\phi)\);<br />
3. Take the path \(r=1\) from \(\theta=2\phi\) to some \(\theta=2\pi-2\psi\);<br />
4. Drop onto the graph \(r=\sec(\theta-(2\pi-2\psi))\), between the points \((r,\theta)=(1,2\pi-2\psi)\) and \((r,\theta)=(\sec(2\pi-\psi),2\pi-\psi)\).</p>
<p>Finally, we need to solve two minimizations, which will find the best \(\phi\) and \(\psi\). \(\psi\) can actually be solved by inspection, but formally:</p>
\(\min_{\psi\in[0,\pi/2)} \tan(\psi) + \pi - 2\psi\)<br />
solution at \(\sec^2{\psi} = 2 \Rightarrow \psi=\pi/4\)
<p>\(\min_{\phi\in[0,\pi/2)} \sec(\phi) + \tan(\phi) + \pi - 2\phi\)<br />
solution at \(\frac{\sin(\phi)}{\cos^2(\phi)} + \sec^2(\phi) = 2 \Rightarrow \sin(\phi)=1/2 \Rightarrow \phi=\pi/6\).</p>
<p>The solved path has total length \(\tan(\pi/4) + \pi &#8211; \pi/2 + \sec(\pi/6) + \tan(\pi/6) + \pi &#8211; \pi/3 = 1+\sqrt{3}+7\pi/6\), and is plotted below. This is about 12% shorter than the trivial path.</p>
<p><img src="wp-content/uploads/images/roadpath.png" /></p>
<hr />
* Technically, we still need to prove that paths that re-enter the unit circle (go below \(r=1\)) are not worthy. One can reason about why such paths can be improved upon by replacing the path segments that re-enter the unit circle, but it seems a proof needs some global properties of the path. With the reader comment (below) about reflecting the starting point, the best path actually becomes a function \(r(\theta)\), of the variable \(\theta\), which almost certainly is a necessary condition for efficient paths. Then we can impose other constraints, like, \(\forall \theta_0: r(\theta) \not< \sec(\theta-\theta_0)\), which says the path function should not be dominated by any secant function (equivalent to crossing all tangent lines to the unit circle). This may be a direction to a complete proof.</p>
]]></content:encoded>
			<wfw:commentRss>https://blog.yhuang.org/?feed=rss2&#038;p=242</wfw:commentRss>
		<slash:comments>4</slash:comments>
		</item>
		<item>
		<title>triangular pursuit</title>
		<link>https://blog.yhuang.org/?p=123</link>
		<comments>https://blog.yhuang.org/?p=123#comments</comments>
		<pubDate>Fri, 03 Oct 2008 13:01:53 +0000</pubDate>
		<dc:creator>admin</dc:creator>
				<category><![CDATA[Uncategorized]]></category>
		<category><![CDATA[boundary conditions]]></category>
		<category><![CDATA[delta]]></category>
		<category><![CDATA[edge length]]></category>
		<category><![CDATA[equilateral]]></category>
		<category><![CDATA[equilateral triangle abc]]></category>
		<category><![CDATA[incenter]]></category>
		<category><![CDATA[math math]]></category>
		<category><![CDATA[speed]]></category>
		<category><![CDATA[sqrt]]></category>
		<category><![CDATA[theta]]></category>

		<guid isPermaLink="false">http://scripts.mit.edu/~zong/wpress/?p=123</guid>
		<description><![CDATA[Here&#8217;s a problem posed to me by a friend: Consider an equilateral triangle ABC with edge length 1. At each vertex is an object that is capable of movement at exactly speed 1. Beginning at time 0, each of the three objects moves toward its initial adjacent neighbor object, as in a game of pursuit. [...]]]></description>
			<content:encoded><![CDATA[<p>Here&#8217;s a problem posed to me by a friend:</p>
<p>Consider an equilateral triangle ABC with edge length 1. At each vertex is an object that is capable of movement at exactly speed 1. Beginning at time 0, each of the three objects moves toward its initial adjacent neighbor object, as in a game of pursuit. Of course, by symmetry, the objects will meet at the incenter of ABC. The question: how far will they have traveled?<br />
<span id="more-123"></span><br />
<img src="wp-content/uploads/images/triangle.png" align="right" /> The paths followed by the objects at first seem non-trivial, but in fact turn out to be nice. We will come back to this, but for the problem, the exact paths are irrelevant. We just need to find the amount of time it takes the objects to meet. The key part of the solution is to realize that, due to symmetry, the three objects always define an equilateral triangle at any time. Furthermore, their direction of motion is always along the edges of such a triangle, at speed v=1.</p>
<p>What we have is then a shrinking and rotating equilateral triangle that eventually decays to a point. Let us characterize this process.</p>
<p>At time \(t\), let \(s(t)\) be the length traveled by one of the objects. Let \(l(t)\) be the edge length of the triangle. Let \(A(t)\) be the area of the triangle. Then we have several relationships:</p>
<ul>
<li>\(\frac{ds}{dt}=1\)</li>
<li>\(A(t)=\frac{\sqrt{3}}{4} l(t)^2\), and \(\frac{dA}{dt} = \frac{\sqrt{3}}{2} l(t) \frac{dl}{dt}\)</li>
<li>By geometry, \(A(t+\Delta t) &#8211; A(t) = -3 (l(t) &#8211; \frac{ds}{dt} \Delta t) \frac{\sqrt{3}}{4} \frac{ds}{dt} \Delta t + o({\Delta t}^2)\), thus \(\frac{dA}{dt} = -3 \frac{\sqrt{3}}{4} l(t)\)</li>
</ul>
<p>Combining: \(-\frac{3}{4}\sqrt{3} l(t) = \frac{\sqrt{3}}{2} l(t) \frac{dl}{dt}\)<br />
\(\frac{dl}{dt} = -\frac{3}{2}\). Integrating with boundary conditions, we get \(-\frac{3}{2}t_f = -1\), \(t_f = \frac{2}{3}\), and since speed is 1, \(s(t_f) = t_f = \frac{2}{3}\). This is how far the object travels.</p>
<p>Finally, since the process is self-similar at every step along the path, the path of the object must be a <a href="http://en.wikipedia.org/wiki/Logarithmic_spiral">logarithmic spiral</a>. Furthermore, since the linear speed on the path is constant, the angular speed must be exponentially increasing (in angle) toward the center of the spiral.</p>
<p><img src="wp-content/uploads/images/spiral.png" align="right" /> To solve for this spiral explicitly, we have \(\frac{dr}{dt}=\frac{dr}{d\theta} \frac{d\theta}{dt} = \frac{\sqrt{3}}{3} \frac{dl}{dt} = -\frac{\sqrt{3}}{2}\).<br />
\(\frac{dr}{d\theta} = -k \exp(-k \theta)\), so \(\frac{d\theta}{dt} = \frac{\sqrt{3}}{2k} \exp (k\theta)\), and \(\theta(t) = -\frac{1}{k} \ln(-\frac{\sqrt{3}}{2}t + C)\). Plugging into \(r(\theta) = \exp(-k\theta)\) and applying the boundary condition \(r(\theta_0) = \frac{\sqrt{3}}{3}\), we get \(r(t) = -\frac{\sqrt{3}}{2}t + \frac{\sqrt{3}}{3}\).<br />
Next, \(\int_{\theta_0}^\infty r(\theta) d\theta = \frac{\sqrt{3}}{3k} = \frac{2}{3}\), so \(k=\frac{\sqrt{3}}{2}\) and \(\theta(t) = -\frac{2\sqrt{3}}{3} \ln(-\frac{\sqrt{3}}{2}t + \frac{\sqrt{3}}{3})\).</p>
]]></content:encoded>
			<wfw:commentRss>https://blog.yhuang.org/?feed=rss2&#038;p=123</wfw:commentRss>
		<slash:comments>0</slash:comments>
		</item>
	</channel>
</rss>
