ellipj

Эллиптические функции Якоби

свернуть все на странице

Синтаксис

[SN,CN,DN] = ellipj(U,M)
[SN,CN,DN] = ellipj(U,M,tol)

Описание

пример

[SN,CN,DN] = ellipj(U,M) возвращает эллиптические функции Якоби SN, CN и DN, оцененный для соответствующих элементов аргумента U и параметра M. Входные параметры U и M должны быть одного размера, или или U или M должны быть скаляром.

пример

[SN,CN,DN] = ellipj(U,M,tol) вычисляет эллиптические функции Якоби с точностью tol. Значением по умолчанию tol является eps. Увеличьте tol для менее точного, но более быстро вычисленного ответа.

Примеры

свернуть все

Скрипт Open Live Script

Ищите эллиптические функции Якоби для U = 0.5 и M = 0.25. 

[s,c,d] = ellipj(0.5,0.25)
s = 0.4751

c = 0.8799

d = 0.9714
Скрипт Open Live Script

Постройте эллиптические функции Якоби для -5≤U≤5 и M = 0.7.

M = 0.7;
U = -5:0.01:5;
[S,C,D] = ellipj(U,M);
plot(U,S,U,C,U,D);
legend('SN','CN','DN','Location','best')
grid on
title('Jacobi Elliptic Functions sn,cn,dn')

Скрипт Open Live Script

Сгенерируйте объемную поверхностную диаграмму Якоби эллиптическая функция sn для позволенной области значений M и -5≤U≤5.

[M,U] = meshgrid(0:0.1:1,-5:0.1:5);
S = ellipj(U,M);
surf(U,M,S)
xlabel('U')
ylabel('M')
zlabel('sn')
title('Surface Plot of Jacobi Elliptic Function sn')

Скрипт Open Live Script

Значением по умолчанию tol является eps. Найдите время выполнения со значением по умолчанию для произвольного M с помощью tic и toc. Увеличьте tol фактором 1 000 и найдите время выполнения. Сравните время выполнения.

tic
ellipj(0.253,0.937)
ans = 0.2479

toc
Elapsed time is 0.024740 seconds.

tic
ellipj(0.253,0.937,eps*1000)
ans = 0.2479

toc
Elapsed time is 0.007491 seconds.

ellipj запускается значительно быстрее, когда допуск значительно увеличен.

Входные параметры

свернуть все

Входной массив, заданный как скалярный, векторный, матричный или многомерный массив. U ограничивается действительными значениями. Если U является нескалярным, M должен быть скаляром или нескалярным одного размера как U.

Типы данных: single | double

Входной массив, заданный как скалярный, векторный, матричный или многомерный массив. M может принять значения 0 ≤ m ≤1. Если M является нескалярным, U должен быть скаляром или нескалярным одного размера как M. Сопоставьте другие значения M в эту область значений с помощью преобразований, описанных в [1], уравнения 16.10 и 16.11.

Типы данных: single | double

Точность результата, заданного как неотрицательное вещественное число. Значением по умолчанию является eps.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

Выходные аргументы

свернуть все

Эллиптическая функция Якоби sn, возвращенный как скаляр, вектор, матрица или многомерный массив.

Эллиптическая функция Якоби cn, возвращенный как скаляр, вектор, матрица или многомерный массив.

Эллиптическая функция Якоби dn, возвращенный как скаляр, вектор, матрица или многомерный массив.

Больше о

свернуть все

Эллиптические функции Якоби

Эллиптические функции Якоби заданы с точки зрения интеграла

u=∫0ϕdθ1−msin2θ.

Затем

sn(u)=sinϕ, cn(u)=потому чтоϕ, dn(u)=1−msin2ϕ.

Некоторые определения эллиптических функций используют эллиптический модуль k или модульный угол α вместо параметра m. Они связаны

k2=m=sin2a.

Эллиптические функции Якоби повинуются многим математическим тождествам. Для хорошей выборки см. [1].

Алгоритмы

ellipj вычисляет эллиптические функции Якоби с помощью метода арифметически-среднегеометрического из [1]. Это запускается с триплета чисел

a0=1, b0=1−m, c0=m.

ellipj вычисляет последовательное использование итераций

ai=12(ai−1+bi−1)bi=(ai−1bi−1)12ci=12(ai−1−bi−1).

Затем, это вычисляет амплитуды в использовании радианов

sin(2ϕn−1−ϕn)=cnansin(ϕn),

стараться разворачивает фазы правильно. Эллиптические функции Якоби затем просто

sn(u)=sinϕ0cn(u)=потому чтоϕ0dn(u)=1−m⋅sn(u)2.

Ссылки

[1] Abramowitz, M. и я. А. Стегун, руководство математических функций, Дуврские публикации, 1965, 17.6.

Смотрите также

Представлено до R2006a

Документация MATLAB
Поддержка
Для просмотра документации необходимо авторизоваться на сайте