Используйте функцию legendre, чтобы работать с вектором и затем исследовать формат вывода.

Вычислите значения Функции Лежандра второй степени вектора.

deg = 2;
x = 0:0.1:0.2;
P = legendre(deg,x)

P = 3×3

   -0.5000   -0.4850   -0.4400
         0   -0.2985   -0.5879
    3.0000    2.9700    2.8800

Формат вывода таков что:

Уравнение для связанной Функции Лежандра второй степени

Поэтому значение

Этот результат соглашается с P(1,1) = -0.5000.

Вычислите связанные значения Функции Лежандра с несколькими нормализацией.

Вычислите ненормализованные значения Функции Лежандра первой степени. Первая строка значений соответствует, и вторая строка к.

x = 0:0.2:1;
n = 1;
P_unnorm = legendre(n,x)

P_unnorm = 2×6

         0    0.2000    0.4000    0.6000    0.8000    1.0000
   -1.0000   -0.9798   -0.9165   -0.8000   -0.6000         0

Затем, вычислите полунормализованные значения функции Шмидта. По сравнению с ненормализованными значениями форма Шмидта отличается когда масштабированием

Для первой строки эти две нормализации является тем же самым с тех пор. Для второй строки масштабирующееся постоянное умножение каждого значения-1.

P_sch = legendre(n,x,'sch')

P_sch = 2×6

         0    0.2000    0.4000    0.6000    0.8000    1.0000
    1.0000    0.9798    0.9165    0.8000    0.6000         0

C1 = (-1) * sqrt(2*factorial(0)/factorial(2))

C1 = -1

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

Этот масштабный коэффициент запрашивает все значения, таким образом, первые и вторые строки имеют различные масштабные коэффициенты.

P_norm = legendre(n,x,'norm')

P_norm = 2×6

         0    0.2449    0.4899    0.7348    0.9798    1.2247
    0.8660    0.8485    0.7937    0.6928    0.5196         0

Cm0 = sqrt((3/2))

Cm0 = 1.2247

Cm1 = (-1) * sqrt((3/2)/2)

Cm1 = -0.8660

Сферические гармоники возникают в решении уравнения Лапласа и используются, чтобы представлять функции, определяемые на поверхности сферы. Используйте legendre, чтобы вычислить и визуализировать сферическую гармонику для.

Уравнение для сферических гармоник включает термин для Функции Лежандра, а также комплексный экспоненциал:

Во-первых, создайте сетку значений, чтобы представлять все комбинации (высотный угол) и (азимутальный угол).

dx = pi/60;
alt = -pi/2:dx:pi/2;
az = 0:dx:2*pi;
[phi,theta] = meshgrid(az,alt);

Вычислите на сетке для.

l = 3;
Plm = legendre(l,cos(theta));

Поскольку legendre вычисляет ответ для всех значений, Plm содержит некоторые дополнительные значения функции. Извлеките значения для и отбросьте остальных. Используйте функцию reshape, чтобы ориентировать результаты как матрицу с тем же размером как phi и theta.

m = 2;
P32 = reshape(Plm(m+1,:,:), size(phi));

Вычислите сферические гармонические значения для.

a = (2*l+1)*factorial(l-m);
b = 4*pi*factorial(l+m);
C = sqrt(a/b);
Y32 = C .* P32 .* exp(1i*m*phi);

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

[Xm,Ym,Zm] = sph2cart(phi, theta, real(Y32));
surf(Xm,Ym,Zm)
title('$Y_3^2$ spherical harmonic','interpreter','latex')