Как найти угол пересечения линий

Не такая тривиальная задача, скажу я вам. Всякий раз, когда возникает необходимость посчитать координату пересечения пары прямых, каждая из которых задана парой точек, снова беру блокнот и вывожу пару формул. И всякий раз – блин, ну это уже когда-то было, опять надо что-то делать с параллельными прямыми, опять появляется пакостная строго вертикальна линия, когда на (x1-x2) никак не разделить и т.д.

Поэтому – в подборку теории и практики, пригодится, сэкономим блокнот, спасем дерево.

Коэффициенты А, B, C

Все помним со школы формулу:

Latex formula

Тоже самое, но с претензией на образование (некоторые индивидуумы утверждают, что существует такая, и только такая, и никакая другая, формулировка):

Latex formula

Те же фаберже, только сбоку.

В теории надо составить и решить систему уравнений для первой и второй линии, где переменными будут X и Y точки пересечения.

Latex formula

Загвоздка в том, что мы не знаем коэффициенты для обеих линий.

В нашем случае известны координаты двух точек, по которым проходит линия. Поэтому мне, как последователю геометрического агностицизма, более привлекательная следующая формула:

Latex formula

Путем несложных операций приходим к следующей записи:

Latex formula

Глядя на вариант в исполнении высшего образования, получаем следующие формулы для нахождения коэффициентов:

Latex formula

Пока все идет отлично, нигде вероятного деления на ноль не встретилось.

Итак, мы можем легко найти два набора коэффициентов для первой и второй прямых. Переходим к системе уравнений.

Система уравнений

Как правило, подобная система уравнений решается путем выражения одной переменной через другую, подстановкой во второе уравнение, получая таким образом уравнение одной переменной. Далее переменная находится, подставляется, решается. Или определяется, что система решения не имеет.

Но нас интересует метод Крамера. Потому что с помощью этого метода можно получить сразу значения для обеих переменных, без дополнительных телодвижений.

Сразу же запишем метод под нашу систему.

Имеем следующую систему:

Latex formula

Определители будут такими:

Latex formula

Latex formula

Latex formula

Исходя из метода, решение выглядит так:

Latex formula

Latex formula

Ага! Вот и возможное деление на ноль, скажете вы. И правильно! В этой, в высшей степени непозволительной ситуации, когда знаменатель равен нулю, решения нет, прямые либо параллельны, либо совпадают (что, впрочем, частный случай параллельности).  В коде, естественно, этот момент надо учитывать.

Практика 1

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

//*******************************************************

//  Нахождение точки пересечения прямых (p1,p2) и (p3,p4)

//  Результат — факт пересечения

//*******************************************************

function CrossLines(const p1,p2,p3,p4: TxPoint;

  var res: TxPoint): Boolean;

const

  Prec = 0.0001;

var

  a1, a2: Extended;

  b1, b2: Extended;

  c1, c2: Extended;

  v: Extended;

begin

  a1 := p2.y p1.y;

  a2 := p4.y p3.y;

  b1 := p1.x p2.x;

  b2 := p3.x p4.x;

  v := a1*b2 a2*b1;

  Result := (abs(v) > Prec);

  if Result then

  begin

    c1 := p2.x*p1.y p1.x*p2.y;

    c2 := p4.x*p3.y p3.x*p4.y;

    res.X := (c1*b2 c2*b1)/v;

    res.Y := (a1*c2 a2*c1)/v;

  end;

end;

Пересечение прямых, построенных по двум точкам. Точка пересечения принадлежеит обоим отрезкам.

Рис.1. Пересечение прямых

Частные случаи

  • Прямые параллельны: ∆ab = 0
    • (A1B2 – B1A2 = 0);
  • Прямые совпадают:  ∆ab = ∆X = ∆Y = 0 
    • (A1B2 – B1A2 = 0) И (A1C2 — A2C1 = 0) И (C1B2 -B1C2 = 0);
  • Прямые перпендикулярны:
    • (A1 A2 + B1 B2 = 0).

Пересечение перпендикулярных прямых

Рис.2. Пересечение перпендикулярных прямых
Параллельные прямые
Рис.3. Параллельные прямые не пересекаются

Принадлежность точки отрезку

В общем случае, чтобы определить принадлежность точки отрезку, надо установить две вещи:

  1. Точка принадлежит прямой, проходящей через конечные точки отрезка. Для этого достаточно подставить значение X и Y в уравнение прямой и проверить получившееся равенство. В нашем случае, этот пункт уже выполнен, т.к. точка пересечения априори принадлежит обеим прямым.
  2. Проверить факт нахождения точки между концами отрезка.

Займемся пунктом 2. Данный факт можно установить двумя способами:

  • Логически, т.е. (x1 <= x <= x2) ИЛИ (x1 >= x >= x2). На случай «вертикальности» линии добавить проверку на Y:
    • (y1 <= y <= y2) ИЛИ (y1 >= y >= y2).
  • Арифметически. Сумма отрезков |x-x1| + |x-x2| должна быть равна длине отрезка |x1-x2|. Аналогично, на случай «вертикальности» , добавить проверку:
    • |y-y1| + |y-y2| = |y1-y2|

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

//*****************************************************

//  Проверка факта нахождения точки res между

//  концами отрезка (p1,p2).

//  Решение с помощью условных операторов и

//  коэффициентов A=(y2-y1) B=(x1-x2).

//  Выступают в качестве параметров, чтобы не тратить

//  время на их подсчет, т.к. в вызывающей стороне

//  они уже посчитаны

//*****************************************************

function CheckCrossPoint(const p1, p2, res: TxPoint;

  const A,B: Extended): Boolean;

begin

  Result :=

    (((B<0) and (p1.X < res.X) and (p2.X > res.X)) or

     ((B>0) and (p1.X > res.X) and (p2.X < res.X)) or

     ((A<0) and (p1.y > res.Y) and (p2.Y < res.Y)) or

     ((A>0) and (p1.y < res.Y) and (p2.Y > res.Y)));

end;

//*****************************************************

//  Проверить факт нахождения точки res между

//  концами отрезка (p1,p2)

//  Арифметическое решение без коэффициентов

//*****************************************************

function CheckCrossPoint(const p1, p2, res: TxPoint): Boolean;

begin

  Result :=

    (abs(p2.xp1.x)>= abs(p2.xres.x) + abs(p1.xres.x)) and

    (abs(p2.yp1.y)>= abs(p2.yres.y) + abs(p1.yres.y));

end;

Практика показывает, что арифметический способ быстрее примерно в 3 раза. Когда-то я считал, что операции сравнения самые быстрые. Это давно уже не так.

Задача нахождения принадлежности точки P(x,y) отрезку, заданного двумя точками с координатами P1(x1, y1) и P2(x2, y2) подробно рассмотрена в отдельной статье.

Угол пересечения прямых

Угол пересечения прямых — это угол пересечения направляющих векторов. Т.е., взяв уже знакомые ранее точки p1 и p2, получим направляющий вектор V(p1,p2), и аналогично второй вектор M(p3,p4). В теории мы должны вычислить достаточно «затратную» функцию, с корнями, квадратами, дробями и арккосинусом.

Давайте не будем останавливаться на ней, она долгая, нудная и в нашем случае ненужная. Рассмотрим вектор:

Вектор из точки p1 в точку p2 с указанием расстояний по Y и X

Рис.4. Вектор V(p1,p2)

α — угол наклона вектора к оси X, который можно найти, как:

α = arctan (A1 / B1)

Где расстояния:

A1 = (y1 — y2)

B1 = (x2 — x1)

Что-то знакомое? Да это ни что иное, как коэффициенты в уравнении прямой от образованных фанатов. Может они и правы в своем испепеляющем фанатизме…

Одним словом, коэффициенты (расстояния) у нас уже есть по обеим прямым.

Пересекающиеся векторы

Рис.5. Пересекающиеся вектор V(p1,p2) и вектор M(p3,p4)

Судя по рисунку, угол между векторами, это сумма углов наклона векторов к оси X. Ммм… не совсем так, на самом деле это разность.

Пересекающиеся векторы

Рис.6. Пересекающиеся векторы в положительной Y

По рисунку явно видно, что угол между векторам это γ = (βα).

В предыдущем примере все правильно, просто знаки углов разные, т.к. находятся по разные стороны от оси X, а формула работает та же.

От теории к практике

Теперь в плане практического применения. Мне нужно точно знать, откуда, куда и в каком направлении этот угол. В теории, углом между прямыми считается наименьший из пары γ и (180-γ). Так вот, нам это не надо. Какой угол получится – такой нам и нужен.

Поэтому, под углом между векторами понимаем угол от вектора V(p1,p2) к вектору M(p3,p4). Если знак угла – отрицательный, понимаем, что он против часовой стрелки, иначе – по часовой стрелке.

Следует заметить, что, зная коэффициенты, для нахождения угла пересечения, координаты уже не нужны. Листинг таков:

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

//**********************************************************

// Посчитать угол пересечения векторов по коэфф-ам А и B

//**********************************************************

function CalcCrossAngle(const a1,b1: Extended;

  const a2,b2: Extended): Extended;

var

  c1, c2: Extended;

begin

  c1 := ArcTan2(a1,b1);

  c2 := ArcTan2(a2,b2);

  Result := c2c1;

  if Result < pi then

    Result := 2*pi + Result;

  if Result > pi then

    Result := Result 2*pi;

end;

Тут ситуация с вертикальной прямой, т.е. когда теоретически происходит деление на ноль, явно не обрабатывается. Она корректно обрабатывается функцией ArcTan2, которая вернет в этом случае и знак, и 90 градусов.  

Пересечение перпендикулярных векторов с верным подсчетом особого "вертикального" случая.

Рис.7. Пересечение перпендикулярных векторов

Практика 2

В дополнение к функции нахождения точки пересечения, напишем «продвинутую» функцию, которая находит эту точку, определяет нахождение на каждом из отрезков, и определяет угол между направляющими векторами. Или же определяет, что прямые параллельны/совпадают.

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

75

76

77

78

79

80

//**********************************************************

//  Тип пересечения прямых (p1,p2) и (p3,p4)

//**********************************************************

type

  TxCrossLineResult = (

    xclrEqual    = 32// эквивалентны

   ,xclrParallel = 16// параллельны

   ,xclrOk       = 0  // как минимум пересечение есть

   ,xclrFirst    = 1  // попадает в первый отрезок

   ,xclrSecond   = 2  // попадает во второй отрезок

   ,xclrBoth     = 3  // попадает в оба

   ,xclrPerpend  = 4  // перпендикулярны

   // можно найти по маске через AND, но для полноты картины

   ,xclrFirstP   = 5  // перпендикулярны и попадает в первый

   ,xclrSecondP  = 6  // перпендикулярны и попадает в второй

   ,xclrBothP    = 7  // перпендикулярны и попадает в оба

   );

//**********************************************************

//  Нахождение точки пересечения прямых (p1,p2) и (p3,p4)

//  Определяет параллельность, совпадение,

//  перпендикулярность, пересечение.

//  Определяет, каким отрезкам принадлежит.

//  Находит угол(рад.) от (p1,p2) к (p3,p4):

//    отрицательное значение — против часовой

//    положительное — по часовой

//**********************************************************

function CrossLines(const p1,p2,p3,p4: TxPoint;

  var res: TxPoint; var Angle: Extended): TxCrossLineResult;

const

  Prec = 0.0001;

var

  a1, a2: Extended;

  b1, b2: Extended;

  c1, c2: Extended;

  v: Extended;

begin

  Angle := 0;

  a1 := p2.y p1.y;

  a2 := p4.y p3.y;

  b1 := p1.x p2.x;

  b2 := p3.x p4.x;

  c1 := p2.x*p1.y p1.x*p2.y;

  c2 := p4.x*p3.y p3.x*p4.y;

  v := a1*b2 a2*b1;

  if abs(v) > Prec then

  begin

    Result := xclrOk;

    res.X := (c1*b2 c2*b1)/v;

    res.Y := (a1*c2 a2*c1)/v;

    if CheckCrossPoint(p1,p2,res) then

      Result := TxCrossLineResult(Integer(Result) +

        Integer(xclrFirst));

    if CheckCrossPoint(p3,p4,res) then

      Result := TxCrossLineResult(Integer(Result) +

        Integer(xclrSecond));

    if (abs(a1*a2 + b1*b2) < Prec) then

      Result := TxCrossLineResult(Integer(Result) +

        Integer(xclrPerpend));

    Angle := CalcCrossAngle(a1,b1,a2,b2);

  end else

  begin

    Result := xclrParallel;

    if ((abs(c1*b2 c2*b1) < Prec) and

       (abs(a1*c2 a2*c1) < Prec))

    then

      Result := xclrEqual;

  end;

end;

Исходники

Небольшие комментарии по интерфейсу.

Интерфейс программы

Рис.8. Интерфейс программы

Скачать (219 Кб): Исходники (Delphi XE 7-10)

Скачать (1.14 Мб): Исполняемый файл

При запуске генерируется случайным образом 4 точки, по две на прямую. Точки и отрезки можно перетаскивать мышкой. Также, слева присутствует панель, на которой можно ввести координаты точек или коэффициенты уравнения прямой. При нажатии «Enter» или когда элемент ввода теряет фокус, происходит перерасчет и перерисовка.

Внизу есть 4 кнопки переключения режимов отображения. Начиная со второй, помимо координат точки пересечений в верхнем левом углу будет отображаться текущий угол пересечения между направляющими векторами.

Если точка пересечений попадает в какой-либо из отрезков, соответствующий заголовок линии отрезка станет жирным. На рисунке это зеленая линия 2.

По умолчанию, рабочее поле системы координат имеет размерность [-10..10], которую можно изменить ползунком в нижнем правом углу.

Угол между двумя пересекающимися кривыми определяется как угол между двумя прямыми, касательными к кривым в точке их пересечения (рис. 1) по формуле

displaystyle tg: varphi =frac{k_{1}-k_{2}}{1+k_{1}k_{2}},; ; (2)

где displaystyle k_{1} и displaystyle k_{2} — угловые коэффициенты касательных к кривым в точке их пересечения displaystyle P(x_{0},y_{0}),
т. е. частные значения в точке displaystyle x_{0} производных от y по x из уравнений этих кривых:

displaystyle k_{1}=tg: alpha _{1}=left ( frac{dy_{1}}{dx} right )_{x=x_{0}};; k_{2}=tg: alpha _{2}=left ( frac{dy_{2}}{dx} right )_{x=x_{0}}.


Угол между двумя кривыми. Практикум по математическому анализу. Урок 36

Рис.1

Пример 1. Найти углы, под которыми пересекаются следующие линии:
1) прямая displaystyle x+y-4=0 и парабола displaystyle 2y=8-x^{2};
2) эллипс displaystyle x^{2}+4y^{2}=4 и парабола displaystyle 4y=4-5x^{2};
3) синусоида displaystyle y=sin x и косинусоида displaystyle y=cos x.
Решение.
1) Совместно решая уравнения параболы и прямой, находим, что они пересекаются в двух точках: A(0;4) и B(2;2), рис.2.
Угол между двумя кривыми. Практикум по математическому анализу. Урок 36

Рис.2

Далее находим производную от y по x из уравнения параболы: displaystyle 2y'=-2x,: y'=-x и определяем угловые коэффициенты касательных к параболе в точках A и B, как частные значения этой производной:

displaystyle y'_{A}=k_{A}=0;; y'_{B}=k_{B}=-2.

Угловой коэффициент прямой один и тот же во всех ее точках; у данной прямой он равен — 1.
Согласно формуле (2) получим

displaystyle textrm{tg}: A=1,: A=45^{circ};; textrm{tg}: B=frac{-1+2}{1+2}=frac{1}{3},; Bapprox 18,5^{circ}.

2) Решая совместно уравнения кривых, находим их общие точки: A(1,2;-0,8), B(0;1) и C(-1,2;-0,8) рис.3. Затем определяем угловые коэффициенты displaystyle k_{1} и displaystyle k_{2} касательных в любой точке эллипса и параболы как производные от y по x из их уравнений

displaystyle k_{1}=-frac{x}{4y};; k_{2}=-frac{5}{2}x.

Угол между двумя кривыми. Практикум по математическому анализу. Урок 36

Рис.3

Подставляя координаты точки A, получим displaystyle k_{1}=frac{3}{8} и displaystyle k_{2}=-3. Следовательно, в точке A:

displaystyle textrm{tg}: varphi =frac{frac{3}{8}+3}{1-frac{9}{8}}=-27;; varphi approx 92^{circ}.

Под таким же углом кривые пересекаются и в точке C вследствие их симметричности относительно оси Oy.
В точке B имеем: displaystyle k_{1}=k_{2}=0, следовательно, в точке кривые имеют общую касательную, т. е. касаются друг друга. В этой точке угол между кривыми равен нулю.
3) Абсциссы точек пересечения кривых (рис.4) определяются уравнением displaystyle sin x=cos x, решая которое, получим

displaystyle x=frac{pi }{4}+pi n; (n=0,pm 1,pm 2,...).

Дифференцированием находим угловые коэффициенты касательных к синусоиде и косинусоиде: displaystyle k_{1}=cos x;: k_{2}=-sin x.
Угол между двумя кривыми. Практикум по математическому анализу. Урок 36

Рис.4

Искомый угол между кривыми определяем по общей формуле (2)
displaystyle textrm{tg}: varphi =frac{cos x+sin x}{1-cos xsin x}=pm frac{frac{sqrt{2}}{2}+frac{sqrt{2}}{2}}{1-frac{1}{2}}=pm 2sqrt{2}.
Положительному знаку соответствует острый угол displaystyle varphi approx 70,5^{circ}, отрицательному — тупой, смежный с ним угол displaystyle varphi_{1} approx 109,5^{circ}.

Углом
между двумя кривыми

у
= f1(x)
и у
= f2(x)
в точке их пересечения М0(х0,
у0)
называется угол между касательными к
этим кривым в точке М0.
Этот угол определяется по формуле

=

.

Пример.
Найти угол между параболами

у
= 8 – х2
и у
= х2.

□ Для
нахождения координат точек пересечения
заданных кривых решим систему уравнений

В
результате получим А(2;
4) и В(−2;
4). Продифференцируем уравнения парабол:

= −2х,

= 2х.
Найдем значения

и

для точки А(2;
4):

= −4,

= 4. Следовательно,

=

=

и

=


.

Аналогично
определяется угол между кривыми в точке
В(−2;
4):

=


.

§ 21. Формула тейлора

Теорема.
Пусть функция f(x)
имеет в точке а
и некоторой ее окрестности производные
порядка п
+ 1. Пусть х
– любое значение аргумента из указанной
окрестности, х
а.
Тогда между точками а
и х
найдется точка

такая, что справедлива формула:

f(x)
= f(а)
+

(х
а)
+
(х
а)2+
…+

(х
а)п
+

+
(х
а)п+1.

Эту
формулу называют формулой
Тейлора
.

Выражение

Rn+1(x)
=

(х
а)п+1

называют
остаточным
членом

формулы Тейлора.

Запишем остаточный
член в другом виде:

так
как

(а,
х),
то найдется число

,
0 <

< 1, что

= а
+

(х
а)
и тогда

Rn+1(x)
=

(х
а)п+1,
0 <

< 1.

Эта
форма остаточного члена наиболее
употребительна в приложениях.

Если
в формуле Тейлора а
= 0, то получим формулу
Маклорена
:

f(x)
= f(0)
+

х
+

х2
+
… +

хп
+
Rn+1(x)

с
остаточным членом

Rn+1(x)
=

хп+1,
0 <

< 1.

Разложение
некоторых элементарных функций по
формуле Маклорена

1.
f(x)
= ех.

Так как

f(x)
=

=

= … = f
п+1(x)
= ех,

f(0)
=

=

= … = f
п+1(0)
= 1,

то
формула Маклорена имеет вид

ех
= 1 +

+

+

+…+

+ Rn+1(x),

где

Rn+1(x)
=

хп+1,
0 <

< 1.

Аналогично
можно разложить по формуле Маклорена
следующие функции:

2.
f(x)
=

.

= х

+



+
…+ (−1)т+1

+ R2т(x),

где
R2т(x)
= (−1)т
·
,
0 <

< 1.

3.
f(x)
=

.

= 1

+



+
…+ (−1)т

+ R2т+1(x),

где
R2т+1(x)
= (−1)т+1
·
,
0 <

< 1.

4.
f(x)
= (1 + х)т.

(1
+ х)т
=1+
х+
х2+

х
3+…+

+

х
п
+Rn+1(x),

где

Rn+1(x)=

х
п+1(1
+

)тп-1,
0 <

< 1.

Пример.
Вычислить число е.

□ Запишем
разложение ех
по формуле Маклорена:

ех

= 1+

+

+

+…+

+
хп+1,
0 <

< 1.

Если
заменить функцию ех
ее многочленом Тейлора степени п
(отбросим остаточный член), то получим
приближенное равенство

ех

1 +

+

+

+…+

,
(1)

абсолютная
погрешность которого

| Rn+1(x)
| =

| х
|п+1,
0 <

< 1.

Если
рассматривать функцию ех
для −1 ≤ х
≤ 1, то

|
Rn+1(x)
| ≤

<

.

Полагая
в (1) х
= 1, получаем приближенное значение числа

е
≈ 1+

+

+

+ …+

.

При
этом | Rn+1(x)
| <

.

Если
требуется вычислить значение е
с точностью

= 0,001, то число п
определяется из неравенства

< 0,001, или (п
+ 1)! > 3000,

которое
выполняется при п
= 6. Следовательно,

е
≈ 1+

+

+

+ …+

.

Вычисляя
с четырьмя знаками после запятой, получим

е
≈ 2,7180.

Три
знака после запятой гарантированы.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]

  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #
  • #

VMath

Инструменты сайта

Основное

Навигация

Информация

Действия

Содержание

Длина дуги, угол между линиями, площадь области на поверхности

Краткие теоретические сведения

Зная первую квадратичную форму поверхности, мы можем решить три задачи:

2. Найти угол между двумя линиями на поверхности в точке их пересечения:
Если две линии, лежащие на поверхности с первой квадратичной формой $I_1=E,du^2+2F,du,dv+G,dv^2$, пересекаются в некоторой точке $P$ поверхности и имеют в этой точке направления $(du:dv)$ и $(delta u:delta v)$, то косинус угла между ними определяется по формуле: begin mbox,varphi = displaystylefrac<sqrtcdotsqrt> \ mbox,varphi = displaystylefrac<sqrtcdotsqrt>. end Говорим, что кривая на поверхности $vec=vec(u,v)$ в точке $(u,v)$ имеет направление $(du:dv)$, если вектор $dvec=vec_udu+vec_vdv$ является касательным вектором кривой в этой точке.

3. Найти площадь области $Omega$ на поверхности: begin S = iintlimits_sqrtdu,dv, end где $D$ — прообраз $Omega$ на плоскости $(u,v)$.

Решение задач

Задача 1 (почти Феденко 684)

Найти длину дуги кривой, заданной уравнениями $v=3u$ на поверхности с первой квадратичной формой begin I_1=du^2+frac19,mbox^2u,dv^2 end между точками $M_1(u_1,v_1)$ и $M_2(u_2,v_2)$.

Решение задачи 1

Задача 2 (почти Феденко 682)

Под каким углом пересекаются линии $$ u+v=a, ,, u-v=a,$$ лежащие на поверхности: begin x=u,mboxv, ,, y=u,mbox,v, ,, z=au. end

Решение задачи 2

Первая квадратичная форма данной поверхности: begin I_1=(1+a^2),du^2+u^2,dv^2. end

Данные линии пересекаются в точке: begin left < beginu+v&=a,\ u-v&=a. end right. quad Rightarrow quad P(u=a,v=0). end

Направления данных линий: begin du+dv=0, ,, delta u-delta v=0,, Rightarrow end begin du = -dv, ,, delta u = delta v. end

Задача 3

Дана поверхность: $$z=axy.$$ Найти углы между координатными линиями.

Решение задачи 3

Координатные линии на данной поверхности задаются уравнениями: $x=x_0$, $y=y_0$. Запишем коэффициенты первой квадратичной формы: begin &E=1+(z_x)^2=1+a^2y^2,\ &F=z_xz_y=a^2xy, \ &G=1+(z_y)^2=1+a^2x^2. end

Направления координатных линий: begin &x=x_0 ,, Rightarrow dx=0,\ &y=y_0 ,, Rightarrow delta y=0. end

Задача 4 (Дополнение к Задаче 3)

Как мы вывели в примере выше, угол между координатными линиями равен

Из формулы следует, что координатная сеть поверхности ортогональна (координатные линии пересекаются под прямым углом), тогда и только тогда, когда $F$=0.

Задача 5 (Феденко 683)

Найти периметр и внутренние углы криволинейного треугольника $$ u=pm av^2/2,,, v=1,$$ расположенного на поверхности $$I_1=du^2+(u^2+a^2)dv^2.$$

Вершины треугольника: begin &A(u=0,, v=0),\ &B(u=-frac<2>,, v=1), \ &C(u=frac<2>,, v=1). end

Зная координаты вершин и уравнения сторон, найдем длины дуг, составляющих стороны треугольника $ABC$, и углы между линиями в точках их пересечения, то есть в вершинах треугольника: begin &s_1 = |BC| = a,\ &s_2 = |AC| = frac76 a,\ &s_3 = |BC| = frac76 a,\ &P_<triangle ABC>=s_1+s_2+s_3=frac<10><3>a. end begin &mbox,A = 1, ,, mbox,B=mbox,C=frac23. end

Пересечение прямых, угол и координаты пересечения

IP76 > Пересечение прямых, угол и координаты пересечения

Не такая тривиальная задача, скажу я вам. Всякий раз, когда возникает необходимость посчитать координату пересечения пары прямых, каждая из которых задана парой точек, снова беру блокнот и вывожу пару формул. И всякий раз – блин, ну это уже когда-то было, опять надо что-то делать с параллельными прямыми, опять появляется пакостная строго вертикальна линия, когда на (x1-x2) никак не разделить и т.д.

Поэтому – в подборку теории и практики, пригодится, сэкономим блокнот, спасем дерево.

Коэффициенты А, B, C

Все помним со школы формулу:

Тоже самое, но с претензией на образование (некоторые индивидуумы утверждают, что существует такая, и только такая, и никакая другая, формулировка):

Те же фаберже, только сбоку.

В теории надо составить и решить систему уравнений для первой и второй линии, где переменными будут X и Y точки пересечения.

Загвоздка в том, что мы не знаем коэффициенты для обеих линий.

В нашем случае известны координаты двух точек, по которым проходит линия. Поэтому мне, как последователю геометрического агностицизма, более привлекательная следующая формула:

Путем несложных операций приходим к следующей записи:

Глядя на вариант в исполнении высшего образования, получаем следующие формулы для нахождения коэффициентов:

Пока все идет отлично, нигде вероятного деления на ноль не встретилось.

Итак, мы можем легко найти два набора коэффициентов для первой и второй прямых. Переходим к системе уравнений.

Система уравнений

Как правило, подобная система уравнений решается путем выражения одной переменной через другую, подстановкой во второе уравнение, получая таким образом уравнение одной переменной. Далее переменная находится, подставляется, решается. Или определяется, что система решения не имеет.

Но нас интересует метод Крамера. Потому что с помощью этого метода можно получить сразу значения для обеих переменных, без дополнительных телодвижений.

Сразу же запишем метод под нашу систему.

Имеем следующую систему:

Исходя из метода, решение выглядит так:

Ага! Вот и возможное деление на ноль, скажете вы. И правильно! В этой, в высшей степени непозволительной ситуации, когда знаменатель равен нулю, решения нет, прямые либо параллельны, либо совпадают (что, впрочем, частный случай параллельности). В коде, естественно, этот момент надо учитывать.

Практика 1

Частные случаи

  • Прямые параллельны: ∆ab = 0
    • (A1B2 – B1A2 = 0);
  • Прямые совпадают: ∆ab = ∆X = ∆Y = 0
    • (A1B2 – B1A2 = 0) И (A1C2 — A2C1 = 0) И (C1B2 -B1C2 = 0);
  • Прямые перпендикулярны:
    • (A1 A2 + B1 B2 = 0).

Рис.2. Пересечение перпендикулярных прямых Рис.3. Параллельные прямые не пересекаются

Принадлежность точки отрезку

В общем случае, чтобы определить принадлежность точки отрезку, надо установить две вещи:

  1. Точка принадлежит прямой, проходящей через конечные точки отрезка. Для этого достаточно подставить значение X и Y в уравнение прямой и проверить получившееся равенство. В нашем случае, этот пункт уже выполнен, т.к. точка пересечения априори принадлежит обеим прямым.
  2. Проверить факт нахождения точки между концами отрезка.

Займемся пунктом 2. Данный факт можно установить двумя способами:

  • Логически, т.е. (x1 = x >= x2). На случай «вертикальности» линии добавить проверку на Y:
    • (y1 = y >= y2).
  • Арифметически. Сумма отрезков |x-x1| + |x-x2| должна быть равна длине отрезка |x1-x2|. Аналогично, на случай «вертикальности» , добавить проверку:
    • |y-y1| + |y-y2| = |y1-y2|

Практика показывает, что арифметический способ быстрее примерно в 3 раза. Когда-то я считал, что операции сравнения самые быстрые. Это давно уже не так.

Задача нахождения принадлежности точки P(x,y) отрезку, заданного двумя точками с координатами P1(x1, y1) и P2(x2, y2) подробно рассмотрена в отдельной статье.

Угол пересечения прямых

Угол пересечения прямых — это угол пересечения направляющих векторов. Т.е., взяв уже знакомые ранее точки p1 и p2, получим направляющий вектор V(p1,p2), и аналогично второй вектор M(p3,p4). В теории мы должны вычислить достаточно «затратную» функцию, с корнями, квадратами, дробями и арккосинусом.

Давайте не будем останавливаться на ней, она долгая, нудная и в нашем случае ненужная. Рассмотрим вектор:

Рис.4. Вектор V(p1,p2)

α — угол наклона вектора к оси X, который можно найти, как:

Что-то знакомое? Да это ни что иное, как коэффициенты в уравнении прямой от образованных фанатов. Может они и правы в своем испепеляющем фанатизме…

Одним словом, коэффициенты (расстояния) у нас уже есть по обеим прямым.

Рис.5. Пересекающиеся вектор V(p1,p2) и вектор M(p3,p4)

Судя по рисунку, угол между векторами, это сумма углов наклона векторов к оси X. Ммм… не совсем так, на самом деле это разность.

Рис.6. Пересекающиеся векторы в положительной Y

По рисунку явно видно, что угол между векторам это γ = (βα).

В предыдущем примере все правильно, просто знаки углов разные, т.к. находятся по разные стороны от оси X, а формула работает та же.

От теории к практике

Теперь в плане практического применения. Мне нужно точно знать, откуда, куда и в каком направлении этот угол. В теории, углом между прямыми считается наименьший из пары γ и (180-γ). Так вот, нам это не надо. Какой угол получится – такой нам и нужен.

Поэтому, под углом между векторами понимаем угол от вектора V(p1,p2) к вектору M(p3,p4). Если знак угла – отрицательный, понимаем, что он против часовой стрелки, иначе – по часовой стрелке.

Следует заметить, что, зная коэффициенты, для нахождения угла пересечения, координаты уже не нужны. Листинг таков:

Математический портал

Nav view search

Navigation

Search

  • Вы здесь:
  • Home
  • Математический анализ
  • Геометрические применения производных.

Геометрический смысл производной.

Литература: Сборник задач по математике. Часть 1. Под ред А. В. Ефимова, Б. П. Демидовича.

Значение производной $f'(x_0)$ функции $y=f(x)$ в точке $x_0$ равно угловому коэффициенту $k=tgvarphi$ касательной $TT’$ к графику этой функции, проведенной через точку $M_0(x_0, y_0),$ где $y_0=f(x_0)$ (геометрический смысл производной).

Прямая $NN’,$ проходящая через точку касания $M_0$ перпендикулярно к касательной, называется нормалью к графику функции $y=f(x)$ в этой точке. Уравнение нормали $$(x-x_0)+f'(x_0)(y-y_0)=0.$$ Уравнение касательной $TT’$ к графику функции $y=f(x)$ в его точке $M_0(x_0, y_0)$ имеет вид $$y-y_0=f'(x_0)(x-x_0)$$

Углом $omega$ между кривыми $y=f_1(x)$ и $y=f_2(x)$ в их общей точке $M_0(x_0, y_0)$ называется угол между касательными к этим кривым в точке $M_0.$ Его можно вычислить по формуле $$tg,omega=frac<1+f’_1(x_0)f’_2(x_0)>.$$

Примеры.

Написать уравнения касательной и нормали к графику функции $y=f(x)$ в данной точке, если:

5.235. $y=x^2-5x+4,$ $x_0=-1.$

Решение.

Уравнение касательной будем искать по формуле $y-y_0=f'(x_0)(x-x_0);$ уравнение нормали — по формуле $(x-x_0)+f'(x_0)(y-y_0)=0.$

По условию, $x_0=-1. $

$y'(x)=2x-5Rightarrow y'(x_0)=y'(-1)=2cdot (-1)-5=-2-5=-7.$

Подставляем все найденные значения в уравнение касательной:

Теперь находим уравнение нормали:

Ответ: У равнение касательной: $7x+y-3=0;$ уравнение нормали: $ x-7y+71=0.$

Решение.

Уравнение касательной будем искать по формуле $y-y_0=f'(x_0)(x-x_0);$ уравнение нормали — по формуле $(x-x_0)+f'(x_0)(y-y_0)=0.$

По условию, $x_0=4. $

Подставляем все найденные значения в уравнение касательной:

$y-2=frac<1><4>(x-4)Rightarrow 4(y-2)=x-4Rightarrow 4y-8=x-4Rightarrow x-4y+4=0.$

Теперь находим уравнение нормали:

$(x-4)+frac<1><4>(y-2)=0Rightarrow 4(x-4)+(y-2)=0Rightarrow 4x+y-18=0.$

Ответ: У равнение касательной: $x-4y+4=0;$ уравнение нормали: $4x+y-18=0.$

5.241. Написать уравнения касательной и нормали в точке $M_0(2, 2)$ к кривой $x=frac<1+t>,$ $y=frac<3><2t^2>+frac<1><2t>,,, tneq 0.$

Найдем значение $t_0,$ подставляя координаты точки $M_0$ в уравнение кривой: $2=frac<1+t>,$ $2=frac<3><2t^2>+frac<1><2t>.$

$t^2+t-2=0Rightarrow t_1=1, t_2=-2.$

Подставим полученные решения в равенство $frac<1+t>=frac<3><2t^2>+frac<1><2t>:$

$t_2=-2: frac<1-2><-8>=frac<3><8>-frac<1><4>=frac<1><8>neq 2$ — не удовлетворяет нашей системе.

Найдем производную функции, заданной параметрически $y’_x.$

Подставляем все найденные значения в уравнение касательной:

$y-y_0=f'(x_0)(x-x_0)Rightarrow$ $y-2=frac<7><10>(x-2)Rightarrow 10(y-2)=7(x-2)Rightarrow 10y-20=7x-14Rightarrow$ $7x-10y+6=0.$

Теперь находим уравнение нормали:

$(x-x_0)+f'(x_0)(y-y_0)=0Rightarrow$ $(x-2)+frac<7><10>(y-2)=0Rightarrow 10(x-2)+7(y-2)=0Rightarrow 10x+7y-34=0.$

Ответ: У равнение касательной: $7x-10y+6=0;$ уравнение нормали: $10x+7y-34=0.$

Найти углы, под которыми пересекаются заданные кривые:

5.254. $y=x^2$ и $y=x^3.$

Решение.

Найдем координаты точки пересечения заданных кривых. Решаем систему уравнений:

$left<begin y=x^2,\ y=x^3,endright.Rightarrow$ $left<begin y=x^2,\ x^2=x^3,endright.Rightarrow$ $left<begin y=x^2,\ x_1=0\x_2=1,endright.$ Таким образом, кривые пересекаются в точках $M_1(0, 0)$ и $M_2(1, 1).$

Далее найдем значения производных заданых функций в точках пересечения.

Подставляем найденные значения, в формулу нахождения угла:

Ответ: В точке $M_1(0, 0)$ угол равен 0. (т.е. касательные совпадают), в точке $M_2(1, 1)$ угол равен $arctgfrac<1><7>.$

Написать уравнения касательной и нормали к графику функции $y=f(x)$ в данной точке, если:

Ответ: У равнение касательной: $y-5=0;$ уравнение нормали: $x+2=0.$

Ответ: У равнение касательной: $y-2x=0;$ уравнение нормали: $2y+x=0.$

Ответ: У равнение касательной: $x-y-1=0;$ уравнение нормали: $x+y-1=0.$

5.242. Написать уравнения касательных к кривой $$x=tcos t, ,,, y=tsin t,,,, tin(-infty,,, +infty),$$ в начале координат и в точке $t=pi/4.$

5.244. Написать уравнения касательной к кривой $$x^5+y^5-2xy=0 в точке $M_0(1, 1).$

Найти углы,под которыми пересекаются заданные кривые:

5.260. Найти расстояние от начала координат до нормали к линии $y=e^<2x>+x^2,$ проведенной в точке с абсциссой $x=0.$

источники:

Пересечение прямых, угол и координаты пересечения

http://mathportal.net/index.php/matematicheskij-analiz/geometricheskie-primeneniya-proizvodnykh

Задача 61796 Найти углы,под которыми пересекаются…

Условие

Найти углы,под которыми пересекаются линии, заданные уравнениями y=x^2 и x^2+2*y^2=3

математика ВУЗ
1667

Решение

Находим координаты точек пересечения графиков

Решаем систему уравнений:
{x^2=y
{x^2=3-2y^2

Приравниваем правые части
y=3-2y^2

2y^2+y-3=0

D=1-4*2*(-3)=25

y_(1)=-3/2; y_(2)=1

x^2=-3/2 — уравнение не имеет решений

x^2=1 ⇒ x_(1)=-1; x_(2)=1

Составить уравнения касательных к каждой кривой в точках x_(1)=-1; x_(2)=1

y-f(x_(o))=f`(x_(o))*(x-x_(o))

1)
y=x^2

Точка [b]х=-1[/b]

f`(x)=2x

f`(-1)=2*(-1)=-2

y-1=-2*(x-(-1))

[b]y=-2x-1 [/b]

Точка [b]х=1[/b]

f`(1)=2*1=2

y-1=2*(x-1)

[b]y=2x-1[/b]

2)

x^2+2y^2=3

2x+4y*y`=0

y`=-2x/4y

[b]x=-1[/b] ; y=1

f`(-1)=-2*(-1)/4*(1)=1/2

y-1=(1/2)*(x-(-1))

[b]y=(1/2)x+(3/2) [/b]

Точка [b]х=1[/b] ; y=1

f`(1)=-2*1/4=-1/2

y-1=(-1/2)*(x-1)

[b]y=(-1/2)x+(3/2) [/b]

Угол между касательными в точке х=-1, т.е между прямыми [b]y=[red]-2[/red]x-1 [/b] и [b]y=[red](1/2)[/red]x+(3/2) [/b] равен 90 ° ,так как произведение угловых коэффициентов равно (-1)

Угол между касательными в точке х=1, т.е между прямыми [b]y=[red]2[/red]x-1 [/b] и [b]y=[red](-1/2)[/red]x+(3/2) [/b] равен 90 °,так как произведение угловых коэффициентов равно (-1)

Написать комментарий

Понравилась статья? Поделить с друзьями:
  • Как исправить ошибку в приложении произошел сбой
  • Как найти свой биткоин кошелек по емайл
  • Как найти украденное устройство
  • Как найти окно главного меню
  • Как найти среднюю скорость без пути