Поиск Полиномов Legendre для числовых и символьных входов

Найдите полином Legendre степени 3 при 5.6.

legendreP(3,5.6)

ans =
  430.6400

Найдите полином Legendre степени 2 при x.

syms x
legendreP(2,x)

ans =
(3*x^2)/2 - 1/2

Если вы не задаете числовое значение для степени n, legendreP функция не может найти явную форму полинома и возвращает вызов функции.

syms n
legendreP(n,x)

ans =
legendreP(n, x)

Поиск Полинома Legendre с Вектором и Матрицей входами

Найдите полиномы Legendre с степенями 1 и 2 путем установки n = [1 2].

syms x
legendreP([1 2],x)

ans =
[ x, (3*x^2)/2 - 1/2]

legendreP действует поэлементно на n чтобы вернуть вектор с двумя элементами.

Если несколько входов заданы в виде вектора, матрицы или многомерного массива, входы должны быть совпадающими по размеру. Найдите полиномы Legendre, где входные параметры n и x являются матрицами.

n = [2 3; 1 2];
xM = [x^2 11/7; -3.2 -x];
legendreP(n,xM)

ans =
[ (3*x^4)/2 - 1/2,        2519/343]
[           -16/5, (3*x^2)/2 - 1/2]

legendreP действует поэлементно на n и x чтобы вернуть матрицу того же размера, что и n и x.

Дифференцируйте и найдите пределы Legendre Полиномов

Использовать limit чтобы найти предел полинома Лежандра степени 3 как x имеет тенденцию к - ∞.

syms x
expr = legendreP(4,x);
limit(expr,x,-Inf)

ans =
Inf

Использовать diff найти третью производную полинома Лежандра степени 5.

syms n
expr = legendreP(5,x);
diff(expr,x,3)

ans =
(945*x^2)/2 - 105/2

Найдите расширение полинома легенды серии Тейлора

Использовать taylor чтобы найти расширение ряда Тейлора полинома Лежандра степени 2 при x = 0.

syms x
expr = legendreP(2,x);
taylor(expr,x)

ans =
(3*x^2)/2 - 1/2

Построение легенды Полиномов

Постройте полиномы Легенды порядков 1 через 4.

syms x y
fplot(legendreP(1:4, x))
axis([-1.5 1.5 -1 1])
grid on

ylabel('P_n(x)')
title('Legendre polynomials of degrees 1 through 4')
legend('1','2','3','4','Location','best')

Найти корни полинома легенды

Использовать vpasolve чтобы найти корни полинома Лежандра степени 7.

syms x
roots = vpasolve(legendreP(7,x) == 0)

roots =
 -0.94910791234275852452618968404785
 -0.74153118559939443986386477328079
 -0.40584515137739716690660641207696
                                   0
  0.40584515137739716690660641207696
  0.74153118559939443986386477328079
  0.94910791234275852452618968404785