1.3: Кривизна
Поки що, коли ми хотіли наблизити складну криву простою кривою біля якоїсь точки, ми намалювали дотичну лінію до кривої в точці. Це досить сире. Зокрема, дотичні лінії прямі - вони не криві. Ми отримаємо набагато краще уявлення про те, як виглядає складна крива, якщо ми наблизимо її локально дуже простою «кривою», а не прямою лінією. Напевно, найпростіша «крива крива» - це коло 1 і саме це ми і будемо використовувати.
- Коло, яке найкраще наближає задану криву поблизу заданої точки, називається окружністю кривизни або оскуляційним колом 2 у точці.
- Радіус окружності кривизни називається радіусом кривизни в точці і нормально позначаєтьсяρ.
- Кривизна в точціκ=1ρ.
- Центр окружності кривизни називається центром кривизни в точці.
Ці визначення проілюстровані на малюнку нижче. Він показує (частина) оскулюючий коло вP. точці ТочкаC є центром кривизни.
Зверніть увагу, що коли кривизнаκ велика, радіус кривизниρ невеликий і у нас дуже пишна крива. З іншого боку, колиκ кривизна мала, радіус кривизниρ великий, а наша крива майже пряма. Зокрема, прямі мають кривизну рівно нулю.
Ми зараз визначимося з тим, як знайти коло кривизни, починаючи з з'ясування, яким повинен бути її радіус. Спочатку ми розглянемо криві 3, які лежать вxy -площині, а потім перейти до кривих в 3d. Розглянемо чорну криву на малюнку нижче.
Ця фігура також містить (частину) червоного кола, яке дуже добре підходить до кривої між двома радіальними лініями, які є (дуже маленьким) кутомdθ один від одного. Отже, довжинаds дуги частини чорної кривої між двома радіальними лініями повинна бути (по суті) такою ж, як довжина дуги кола між двома радіальними лініями,ρ|dθ|, деρ знаходиться радіус кола. (Ми ставимо абсолютні значення, щоб врахувати можливість, якаdθ може бути негативною.) Такимdθ чиномds=ρ|dθ|. Коли макроскопічний кут, це, звичайно, наближення. Але в межі, якdθ→0, ми повинні закінчитися
ρ=|dsdθ|
Тепер у нас є формула радіуса кривизни, але не в дуже зручній формі, тому що для її оцінки нам потрібно було б знати довжину дуги по кривій як функцію кута наθ крайньому правому малюнку нижче. Тепер ми витратимо деякий час на розробку більш зручних формул дляρ. Спочатку розглянемо три цифри нижче. Всі вони показують ту ж криву, що і на останньому малюнку. Крайня ліва цифра якраз показує
- крива інтересу, яка є чорною кривою, і
- (синя) точка інтересу на чорній кривій. Ми хочемо знайти кривизну в цій точці.
Середня цифра показує ту саму криву і точку інтересу, а також показує
- червоне коло кривизни (тобто найкраще підходить коло) для чорної кривої на синій крапці.
- Червона крапка є центром кривизни.
На крайньому правому малюнку показана така ж чорна крива, блакитна точка інтересу і червона окружність кривизни (хоча б її частина) дещо збільшена.
- Кутθ - це кут між^ıı і вектором радіуса від червоної точки (центру кривизни) до синьої точки (точки інтересу).
- ˆTдотичний вектор до чорної кривої на синій крапці.
- Кутϕ - це кут між^ıı іˆT.ˆT Вектор також дотичний до червоного кола. Оскільки тангенс і радіус вектори для кіл перпендикулярні один одному 4, у нас є щоϕ=θ+π2 і, отже,ρ=|dsdϕ| теж.
Зараз ми в змозі розробити купу формул для радіуса кривизниρ і кривизниκ=1ρ,, які є більш зручними, ніжκ=|dsdϕ|−1. Ці формули будуть використовувати
Якщо⇀r(t) параметризована крива, то
- ⇀v(t)=d⇀rdt(t)вектор швидкості при⇀r(t)
- a(t)=d2⇀rdt2(t)вектор прискорення при⇀r(t)
- ˆT(t)є одиничним дотичним вектором до кривої в⇀r(t) цих точках у напрямку збільшенняt.
- ˆN(t)це одиничний вектор нормалі до кривої в⇀r(t) цих точках до центру кривизни.
- κ(t)це викривлення при⇀r(t)
- ρ(t)радіус кривизни при⇀r(t)
- Задано 5,s(ϕ), тобто якщо ми знаємо довжину дуги вздовж кривої як функцію кута 6,ϕ=∡(^ıı,ˆT), то
ρ=|dsdϕ|κ=|dsdϕ|−1κ=|dϕds|
- Задано,⇀r(s), тобто якщо ми маємо параметризацію кривої через довжину дуги, то
dˆTds(s)=κ(s)ˆN(s)
деˆN(s) - одиничний вектор нормалі до кривої в⇀r(s) цих точках до центру кривизни. - Задано⇀r(t), тобто якщо ми маємо загальну параметризовану криву, то
dˆTdt=κdsdtˆN⇀v(t)=dsdt(t)ˆT(t)a(t)=d2sdt2ˆT+κ(dsdt)2ˆN
- Задано(x(t),y(t)), (для кривих уxy -площині)
κ=|⇀v(t)×a(t)(dsdt)3|=|dxdtd2ydt2−dydtd2xdt2|[(dxdt)2+(dydt)2]3/2
- Заданоy(x), (знову ж таки для кривих уxy -площині)
κ=|d2ydx2|[1+(dydx)2]3/2
-
(а) Враховуючиs(ϕ), тоді
ρ=|dsdϕ|κ=|dsdϕ|−1
Як ми припускаємо, що0<ρ=|dsdϕ|<∞, обернена теорема функції говорить, що ми можемо інвертувати функціюs(ϕ) (принаймні локально), щоб отриматиϕ як функціюs, і що
κ=|dϕds|
(b) Враховуючи⇀r(s), тоді, Лема 1.1.4.c,ˆT(s)=⇀r′(s) є одиницею дотичної до кривої в⇀r(s) і
dˆTds=dˆTdϕdϕds
Тепер до знакаdϕds єκ, і просто тому, щоϕ=∡(^ıı,ˆT), зˆT одиничним вектором,
ˆT=cosϕ^ıı+sinϕ^ȷȷ⟹dˆTdϕ=−sinϕ^ıı+cosϕ^ȷȷ
ТакdˆTdϕ є одиничний вектор, який перпендикулярний 7 доˆT, і, отже, до кривої в⇀r(s), і
dˆTds(s)=κ(s)ˆN(s)
зˆN(s) одиничним вектором нормалі до кривої⇀r(s). в Фактично,ˆN(s) це одиничний нормальний вектор до кривої в⇀r(s) цих точках до центру кривизни.
Щоб побачити це, подивіться на цифри нижче 8, і зверніть увагу, що підставляючи інформацію про знак з кожної цифри в (∗) дає (†). Наприклад,
розглянемо цифру зліва внизу. На цьому малюнку
- xкомпонентˆT є негативним (ˆTліворуч вказує на малюнку),
- що робитьcosϕ негативним (див. (∗∗)),
- що робитьy компонентdˆTdϕ негативним (див. (∗∗) знову),
- такdˆTdϕ і спрямований вниз,
такdˆTdϕ=−ˆN (центр кривизни - червона точка над кривою) і
- якs збільшується (тобто при русі в напрямку стрілки на кривій),ϕ зменшується (на крайній правій частині кривої вϕ≈3π2, той час як на крайній лівій частині кривоїϕ≈π), такdϕds<0 іκ=|dϕds|=−dϕds.
- Отже, по (∗),dˆTds=dˆTdϕdϕds=(−ˆN)(−κ)=κˆN.
У кожній з трьох інших фігур ми також закінчуємоdˆTds=κ(s)ˆN(s).
Зверніть увагу, що якщоκ(s)=0, тодіˆN(s) не визначено. Це має сенс: якщо крива є (локально) прямою лінією, немає «найкращого облягаючого кола».
(c) Задано,⇀r(t), тобто якщо у нас є загальна параметризована крива, ми можемо визначити одиничний дотичний вектор за допомогою Lemma 1.1.4:
⇀v(t)=d⇀rdt(t)=dsdt(t)ˆT(t)⟹ˆT(t)=⇀r′(t)|⇀r′(t)|
Тоді ми можемо визначитиκ іˆN шляхом диференціаціїˆT(t) та використання правила ланцюга:
dˆTdt=dˆTdsdsdt=κdsdtˆN⟹κ(t)=|ˆT′(t)||⇀r′(t)|
Крім того, якщо ми диференціюємо,⇀v(t)=dsdtˆT(t), ми отримаємо, що прискорення
a(t)=d2⇀rdt2t=d2sdt2ˆT+dsdtdˆTdt=d2sdt2ˆT+κ(dsdt)2ˆN
(d) Задано(x(t),y(t)), (для кривих уxy -площині), ми можемо зчитувати кривизну з
⇀v(t)×a(t)=(dsdt(t)ˆT(t))×(d2sdt2ˆT+κ(dsdt)2ˆN)=κ(dsdt)3ˆT׈N(since ˆT׈T=⇀0
ПодумайтеˆT іˆN як 3d вектори, що чиїz -компоненти трапляються нульовими. ЯкˆT іˆN є взаємно перпендикулярними одиничними векторами вxy -площині, перехресний добутокˆT׈N буде+ˆk або або−ˆk. В обох випадках,|⇀v(t)×a(t)|=κ|dsdt|3. Так
κ=|⇀v(t)×a(t)(dsdt)3|=|[dxdt^ıı+dydt^ȷȷ]×[d2xdt2^ıı+d2ydt2^ȷȷ](dsdt)3|=|[dxdtd2ydt2−dydtd2xdt2]ˆk(dsdt)3|=|dxdtd2ydt2−dydtd2xdt2|[(dxdt)2+(dydt)2]3/2
(e) Заданоy(x), знову для кривих уxy -площині, ми можемо параметризувати кривуx, використовуючи як параметр:
⇀r(t)=(X(t),Y(t))with X(t)=t and Y(t)=y(t)
Тоді
dXdt=1d2Xdt2=0dYdt=dydxd2Ydt2=d2ydt2
і
κ=|dXdtd2Ydt2−dYdtd2Xdt2|[(dXdt)2+(dYdt)2]3/2=|d2ydt2|[1+(dydx)2]3/2
Погляньте ще раз на теорему 1.3.3 і зауважте, що
- тангенціальна складова прискорення, тобтоd2sdt2, виникає чисто від зміни швидкості при
- нормальна складова прискорення, тобтоκ(dsdt)2, виникає через кривизну і пропорційна квадрату швидкостіdsdt. Подумайте про те, що відчуваєте, коли їдете. Саме тому велодроми і (автомобільні) гоночні траси часто мають банальні повороти.
- xкомпонентˆT є негативним (ˆTліворуч вказує на малюнку),
Як приклад розминки, а також перевірка того, що наші формули мають сенс, ми знайдемоκ, радіус кривизни кривизни,ρ, одиничний тангенс вектор,ˆT, одиничний нормальний векторˆN, та центр кривизни параметризованої кривої
⇀r(t)=acost^ıı+asint^ȷȷ
зa>0. константою Це, звичайно, коло радіуса,a зосереджене на початку. Як
⇀v(t)=d⇀rdt(t)=−asint^ıı+acost^ȷȷ⟹dsdt(t)=|⇀v(t)|=a
ми маємо, що вектор дотичної одиниці
⇀T(t)=⇀v(t)|⇀v(t)|=−sint^ıı+cost^ȷȷ
Зверніть увагу, як перевірка, що це дійсно вектор довжини один і перпендикулярний вектору радіуса (як очікувалося — крива є колом). Як
dˆTdt(t)=−cost^ıı−sint^ȷȷ
ми маємо, що
ˆN(t)=dˆTdt(t)|dˆTdt(t)|=−cost^ıı−sint^ȷȷκ(t)=|dˆTdt(t)|dsdt(t)=1aρ(t)=1κ(t)=a
Тепер подивіться на малюнок.
Щоб дістатися до центру кривизни, ми повинні почати з⇀r(t) і пройти відстань,ρ(t), яка врешті-решт є радіусом кривизни, у напрямку,ˆN(T), який спрямований до центру кривизни. Таким чином, центр кривизни
⇀r(t)+ρ(t)ˆN(t)=[acost^ıı+asint^ȷȷ]+a[−cost^ıı−sint^ȷȷ]=⇀0
Це має цілком гарний сенс — радіус кривизни є радіусом вихідного кола, а центр кривизни - центр вихідного кола.
Одним з альтернативних варіантів розрахунку кривизни, використовуючиx(t)=acost,y(t)=asint, є
κ(t)=|dxdt(t)d2ydt2(t)−dydt(t)d2xdt2(t)|[(dxdt(t))2+(dydt(t))2]3/2=|−asint(−asint)−acost(−acost)|[(−asint)2+(acost)2]3/2=1a
Ще один альтернативний розрахунок кривизни, використовуючиy(x)=√a2−x2 (для частини кола сy>0),
y′(x)=−x√a2−x2=−xy(x)y″(x)=−y(x)−xy′(x)y(x)2=−y(x)2+x2y(x)3=−a2y(x)3
є
κ(x)=|d2ydx2(x)|[1+(dydx(x))2]3/2=a2y(x)3[1+x2y(x)2]3/2=a2[y(x)2+x2]3/2=1a
Як більш обчислювальний приклад, ми проаналізуємо
⇀r(t)=(cost+tsint)^ıı+(sint−tcost)^ȷȷt>0⇀v(t)=tcost^ıı+tsint^ȷȷa(t)=(cost−tsint)^ıı+(sint+tcost)^ȷȷ
Ми можемо прочитати з⇀v(t), цього
dsdt(t)=|⇀v(t)|=td2sdt2(t)=1⇀T(t)=⇀v(t)|⇀v(t)|=cost^ıı+sint^ȷȷ
Далі, відa(t), ми зачитали, що
(за теоремою 1.3.3 .c)
\ почати {вирівнювати*} &=\ cos t\,\ hat {\ pmb {\ imath}} +\ sin t\,\ hat {\ pmb {\ jmath}} +t^2\ каппа (t)\ hat {\ textbf {N}} (t)\\ означає t^2\ каппа (t)\ hat {\ textbf {N} t) &= - t\ sin t\,\ hat {\ pmb {\ imath}} + т\ cos t\ cos t\,\ hat {\ pmb {\ jmath}}\ end {align*}такt2κ(t) що довжина−tsint^ıı+tcost^ȷȷ, якої єt. Таким чином
κ(t)=1tandˆN(t)=−tsint^ıı+tcost^ȷȷt2κ(t)=−sint^ıı+cost^ȷȷ
В якості альтернативного розрахунку кривизни ми маємо
κ(t)=|⇀v(t)×a(t)|(dsdt(t))3=|[tcost^ıı+tsint^ȷȷ]×[(cost−tsint)^ıı+(sint+tcost)^ȷȷ]|(dsdt(t))3=|[tcost(sint+tcost)−tsint(cost−tsint)]ˆk|(dsdt(t))3=|t2ˆk|t3=1t
Варто подумати, перш ніж розрахувати!
Вправи
Етап 1
У цій главі є багато констант, які можуть бути новими для вас. До них може знадобитися трохи звикання. Питання 1.3.1.1-1.3.1.5 дають практику роботи з цими константами і інтерпретації їх відносин між собою.
Намалюйте криву⇀r(t)=(3sint,3cost). У(0,3), мітці точкиˆT іˆN. Дайте значенняκ іρ в цій точці також.
Розглянемо коло⇀r(t)=(3sint,3cost). ЗнайтиˆT(t) іˆT(s). Потім, використовувати частини (b) і (c) теореми 1.3.3, щоб знайтиˆN(t) іˆN(s).
Функція⇀r(t)=(tcost,tsint),t≥0, визначає спіраль, орієнтовану на початок. Використовуючи тільки геометричну інтуїцію (без розрахунку), прогнозуйтеlim
Нехай\vecs{r} (t)=(e^t,3t,\sin t)\text{.} Що таке\dfrac{ds}{dt}\text{?}
У питанні 1.2.1.5 розділу 1.2 ми виявили, що спіраль
\vecs{r} (t) = e^t (\cos t, \sin t) \nonumber
параметризований в плані довжини дуги
\textbf{R}(s)=\frac{s}{\sqrt{2}}\left(\cos\Big(\ln\Big(\frac{s}{\sqrt{2}}\Big)\Big)\,,\, \sin\Big(\ln\Big(\frac{s}{\sqrt{2}}\Big)\Big)\right). \nonumber
Знайдіть\dfrac{d\hat{\textbf{T}}}{ds} і\dfrac{d\hat{\textbf{T}}}{dt} для цієї кривої.
У цій вправі ми уточнюємо сенс, в якому оскулюлюючий коло - це коло, яке найкраще наближає плоску криву в точці.
- Переводячи і обертаючи нашу систему координат, ми завжди можемо організувати, що точка є(0,0) і що криваy=f(x) зf'(0)=0 іf''(0) \gt 0\text{.} (Ми припускаємо, що кривизна в точці ненульова.)
- y=g(x)Дозволяти бути нижньою половиною кола радіусаr, який знаходиться в центрі(0,r)\text{.}
Показати, що якщоf(x) іg(x) мають однаковий другий порядок наближення Тейлора,x=0\text{,}r то радіус кривизниy=f(x) atx=0\text{.}
Етап 2
За заданою кривою\vecs{r} (t)=(e^t,t^2+t)\text{,} обчислити такі величини:
- \displaystyle \vecs{v} (t)
- \displaystyle \textbf{a}(t)
- \displaystyle \dfrac{ds}{dt}
- \displaystyle \hat{\textbf{T}}(t)
- \displaystyle \kappa(t)
Знайти\kappa(t) кривизну\vecs{r} (t)=(\cos t+\sin t , \sin t - \cos t)\text{.}
Знайти мінімальне і максимальне значення кривизни еліпсаx(t)= a \cos t\text{,}y(t)=b\sin t\text{.} Тутa \gt b \gt 0\text{.}
- Знайти кривизнуy=e^x at(0,1)\text{.}
- Знайдіть рівняння кола, яке найкраще підходитьy=e^x на(0,1)\text{.}
Розглянемо рух пальця, що застряг в протекторі шини, яка знаходиться на велосипеді, що рухається з постійною швидкістю. Цей рух задається параметризованою кривою
\vecs{r} (t) = \big(t - \sin t\,,\, 1 - \cos t\big) \nonumber
зt \gt 0\text{.}
- Намалюйте криву вxy -площині для0 \lt t \lt 4\pi\text{.}
- Знайти і спростити формулу кривизни\kappa(t)\text{.}
- Знайти радіус кривизни оскуляційної окружності до\vecs{r} (t) att = \pi\text{.}
- Знайти рівняння оскулюючого кола до\vecs{r} (t) att = \pi\text{.}
Етап 3
Знайти кривизну\kappa як функцію довжини дугиs (вимірюється від(0,0)) для кривої
x(\theta)=\int_0^\theta \cos\big(\frac{1}{2}\pi t^2\big)dt\quad \quad y(\theta)=\int_0^\theta \sin\big(\frac{1}{2}\pi t^2\big)dt \nonumber
CДозволяти крива в\mathbb{R}^2 заданій графіком функціїy=\frac{x^3}{3}\text{.}\kappa(x) Дозволяти кривизнуC в точці(x, x^3/3)\text{.} Знайти всі точки, де\kappa(x) досягають своїх максимальних значень, або ж пояснити, чому таких точок не існує. Які межі\kappa(x) якx \rightarrow \infty іx \rightarrow -\infty\text{?}
- Кола хороші для вивчення «кривизни», тому що, на відміну від парабол наприклад, швидкість, з якою окружність крива рівномірна по всьому колу.
- «Osculare» - латинське дієслово «цілуватися». Німецький математик Готфрід Вільгельм (фон) Лейбніц (1646—1716) назвав коло «circulus osculans».
- Ми також припустимо, що криві, що цікавлять, плавні, наприклад, без порізів, а не прямі, так що радіус кривизни0 \lt \rho \lt \infty\text{.}
- Ми бачили, що в прикладі 1.1.6.
-
Рівнянняs=s(\phi) називається «внутрішнім рівнянням кривої».
-
Позначення\measuredangle(\hat{\pmb{\imath}},\hat{\textbf{T}}) означає «кут між\hat{\pmb{\imath}} і\hat{\textbf{T}}».
-
Подумайте, чому так повинно бути. Зокрема, накидайте\hat{\textbf{T}}\phi і подумайте про те, про що говорить ескіз.\dfrac{d\hat{\textbf{T}} }{d\phi}\text{.}
-
На кожній з чотирьох фігур стрілка на кривій визначає напрямок збільшення довжини дуги,s а червона крапка - центр кривої у синьої точки.