кривая на эллипсе
Модератор: Admin
кривая на эллипсе
Задана кривая на эллипсоиде в эллиптической системе координат и нужно ее построить. В Maple эта система координат определена и называется confocalellip, но воспользоваться у меня ей ни как не получается. В моем случае a=3, b=2, c=1, w=t,u=t, v=0. Буду благодарен за помощь.
-
- Сообщения: 1776
- Зарегистрирован: Вс май 01, 2005 9:02 pm
Re: кривая на эллипсе
Если чисто формально, то у Вас при Ваших данных получается мнимое значение второй координаты…
-
- Сообщения: 1776
- Зарегистрирован: Вс май 01, 2005 9:02 pm
Может быть стоит самому определить эту систему координат с помощью addcoords. Только я не понимаю как указать, что переход от декартовых зависит от четверти в который мы находимся (я имею виду плюс и минус перед корнями). Кстати подобный вопрос уже обсуждался
http://forum.exponenta.ru/viewtopic.php?t=5370&sid=bed4468188e45ee88fa95492e67539d3
правда в Mathcad.[/code]
http://forum.exponenta.ru/viewtopic.php?t=5370&sid=bed4468188e45ee88fa95492e67539d3
правда в Mathcad.[/code]
-
- Сообщения: 1776
- Зарегистрирован: Вс май 01, 2005 9:02 pm
-
- Сообщения: 1366
- Зарегистрирован: Вс дек 04, 2011 11:07 pm
confocalellip
Поскольку эта тема может представлять интерес не только для вопрошателя, то дам развернутый ответ.
Представление об этой системе координат можно получить следующим образом.
restart; with(plots): infolevel[coordplot3d] := 5:
coordplot3d(confocalellip);
plots/coordplot3d: u const values: [1/2]
plots/coordplot3d: v const values: [13/72]
plots/coordplot3d: w const values: [0]
plots/coordplot3d: u range: 129/512 .. 511/512
plots/coordplot3d: v range: 521/4608 .. 127/512
plots/coordplot3d: w range: -1/9 .. 503/4608
plots/coordplot3d: view: [0 .. 1 0 .. .5 0 .. .5]
Далее, переход от этой диковинной системы координат к декартовой задается формулами, которые можно узнать из справки ?coords,
x := sqrt((a^2-u)*(a^2-v)*(a^2-w)/((a^2-b^2)*(a^2-c^2))):
y := sqrt((b^2-u)*(b^2-v)*(b^2-w)/((-a^2+b^2)*(b^2-c^2))):
z := sqrt((c^2-u)*(c^2-v)*(c^2-w)/((-a^2+c^2)*(-b^2+c^2))):
Подставляем сюда значения параметров, указанных вопрошателем, и параметрические уравнения кривой
A := simplify(eval(`<,>`(x, y, z), [a = 3, b = 2, c = 1, u = -t, v = t, w = 4]))
Vector(3, {(1) = (1/4)*sqrt(2)*sqrt(-t^2+81), (2) = 0, (3) = (1/4)*sqrt(2*t^2-2)}).
Наконец, рисуем пространственную кривую:
spacecurve(A, t = 1 .. 9, thickness = 4, color = navy);
Посоветую вопрошателю повышать грамотность. Неграмотный текст воспринимается как ненормативная лексика со всеми последствиями для репутации вопрошателя.
Представление об этой системе координат можно получить следующим образом.
restart; with(plots): infolevel[coordplot3d] := 5:
coordplot3d(confocalellip);
plots/coordplot3d: u const values: [1/2]
plots/coordplot3d: v const values: [13/72]
plots/coordplot3d: w const values: [0]
plots/coordplot3d: u range: 129/512 .. 511/512
plots/coordplot3d: v range: 521/4608 .. 127/512
plots/coordplot3d: w range: -1/9 .. 503/4608
plots/coordplot3d: view: [0 .. 1 0 .. .5 0 .. .5]

Далее, переход от этой диковинной системы координат к декартовой задается формулами, которые можно узнать из справки ?coords,
x := sqrt((a^2-u)*(a^2-v)*(a^2-w)/((a^2-b^2)*(a^2-c^2))):
y := sqrt((b^2-u)*(b^2-v)*(b^2-w)/((-a^2+b^2)*(b^2-c^2))):
z := sqrt((c^2-u)*(c^2-v)*(c^2-w)/((-a^2+c^2)*(-b^2+c^2))):
Подставляем сюда значения параметров, указанных вопрошателем, и параметрические уравнения кривой
A := simplify(eval(`<,>`(x, y, z), [a = 3, b = 2, c = 1, u = -t, v = t, w = 4]))
Vector(3, {(1) = (1/4)*sqrt(2)*sqrt(-t^2+81), (2) = 0, (3) = (1/4)*sqrt(2*t^2-2)}).
Наконец, рисуем пространственную кривую:
spacecurve(A, t = 1 .. 9, thickness = 4, color = navy);

На самом деле эллиптическая система координат определяется следующим образом:
и отсюда имеем:
x := +/- sqrt((a^2-u)*(a^2-v)*(a^2-w)/((a^2-b^2)*(a^2-c^2))):
y := +/- sqrt((b^2-u)*(b^2-v)*(b^2-w)/((-a^2+b^2)*(b^2-c^2))):
z := +/- sqrt((c^2-u)*(c^2-v)*(c^2-w)/((-a^2+c^2)*(-b^2+c^2))):
и мой вопрос по прежнему остается открытым. Как построить кривую на эллипсоиде в эллиптической системе координат в Maple ? Этот вопрос мне кажется проще всего решить с помощью addcoords, то есть самому определить данную систему координат. Но загвоздка заключается в том, что я не знаю как задать в addcoords, что переход от декартовых зависит от четверти в которой мы находимся.

x := +/- sqrt((a^2-u)*(a^2-v)*(a^2-w)/((a^2-b^2)*(a^2-c^2))):
y := +/- sqrt((b^2-u)*(b^2-v)*(b^2-w)/((-a^2+b^2)*(b^2-c^2))):
z := +/- sqrt((c^2-u)*(c^2-v)*(c^2-w)/((-a^2+c^2)*(-b^2+c^2))):
и мой вопрос по прежнему остается открытым. Как построить кривую на эллипсоиде в эллиптической системе координат в Maple ? Этот вопрос мне кажется проще всего решить с помощью addcoords, то есть самому определить данную систему координат. Но загвоздка заключается в том, что я не знаю как задать в addcoords, что переход от декартовых зависит от четверти в которой мы находимся.
-
- Сообщения: 1366
- Зарегистрирован: Вс дек 04, 2011 11:07 pm
-
- Сообщения: 1366
- Зарегистрирован: Вс дек 04, 2011 11:07 pm