WikiDer > Дискретное преобразование Чебышева

эта статья нужны дополнительные цитаты для проверка. Пожалуйста помоги улучшить эту статью от добавление цитат в надежные источники. Материал, не полученный от источника, может быть оспорен и удален. Найдите источники: «Дискретное преобразование Чебышева» – Новости · газеты · книги · ученый · JSTOR (Февраль 2012 г.) (Узнайте, как и когда удалить этот шаблон сообщения)

В Прикладная математика, то дискретное преобразование Чебышева (ДКП), названный в честь Пафнутый Чебышев, является одной из двух основных разновидностей DCT: дискретное преобразование Чебышева на сетке «корней» Полиномы Чебышева первого вида ${ Displaystyle Т_ {п} (х)}$ и дискретное преобразование Чебышева на сетке «экстремумов» полиномов Чебышева первого рода.

Дискретное преобразование Чебышева на сетке корней

Дискретное преобразование Чебышева функции u (x) в точках ${ displaystyle {x_ {n}}}$ дан кем-то:

${ displaystyle a_ {m} = { frac {p_ {m}} {N}} sum _ {n = 0} ^ {N-1} u (x_ {n}) T_ {m} (x_ {n })}$

где:

${ displaystyle x_ {n} = - cos left ({ frac { pi} {N}} (n + { frac {1} {2}}) right)}$

${ displaystyle a_ {m} = { frac {p_ {m}} {N}} sum _ {n = 0} ^ {N-1} u (x_ {n}) cos left (m cos ^ {- 1} (x_ {n}) right)}$

где ${ displaystyle p_ {m} = 1 Leftrightarrow m = 0}$ и ${ displaystyle p_ {m} = 2}$ в противном случае.

Используя определение ${ displaystyle x_ {n}}$,

${ displaystyle a_ {m} = { frac {p_ {m}} {N}} sum _ {n = 0} ^ {N-1} u (x_ {n}) cos left ({ frac {m pi} {N}} (N + n + { frac {1} {2}}) right)}$

${ displaystyle a_ {m} = { frac {p_ {m}} {N}} sum _ {n = 0} ^ {N-1} u (x_ {n}) (- 1) ^ {m} cos left ({ frac {m pi} {N}} (n + { frac {1} {2}}) right)}$

и его обратное преобразование:

${ displaystyle u_ {n} = sum _ {m = 0} ^ {N-1} a_ {m} T_ {m} (x_ {n})}$

(Так происходит со стандартным рядом Чебышева, вычисленным на сетке корней.)

${ displaystyle u_ {n} = sum _ {m = 0} ^ {N-1} a_ {m} cos left ({ frac {m pi} {N}} (N + n + { frac {1} {2}}) right)}$

${ displaystyle , следовательно, u_ {n} = sum _ {m = 0} ^ {N-1} a_ {m} (- 1) ^ {m} cos left ({ frac {m pi} { N}} (n + { frac {1} {2}}) right)}$

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

Это можно продемонстрировать, используя следующие MATLAB код:

функцияа=fct(ж, л)% х = -cos (пи / N * ((0: N-1) '+ 1/2));ж = ж(конец:-1:1,:);А = размер(ж); N = А(1);если существуют ('A (3)', 'var') && A (3) ~ = 1    для я = 1: А (3)        а(:,:,я) = sqrt(2/N) * dct(ж(:,:,я));        а(1,:,я) = а(1,:,я) / sqrt(2);    конецеще    а = sqrt(2/N) * dct(ж(:,:,я));    а(1,:)=а(1,:) / sqrt(2);конец

Дискретное косинусное преобразование (dct) фактически вычисляется с использованием алгоритма быстрого преобразования Фурье в MATLAB.
И обратное преобразование дается кодом MATLAB:

функцияж=ifct(а, л)% x = -cos (pi / N * ((0: N-1) '+ 1/2)) k = размер(а); N=k(1);а = idct(sqrt(N/2) * [а(1,:) * sqrt(2); а(2:конец,:)]);конец

Дискретное преобразование Чебышева на сетке экстремумов

Это преобразование использует сетку:

${ displaystyle x_ {n} = - cos left ({ frac {n pi} {N}} right)}$

${ displaystyle T_ {n} (x_ {m}) = cos left ({ frac { pi mn} {N}} + n pi right) = (- 1) ^ {n} cos слева ({ frac { pi mn} {N}} right)}$

Это преобразование сложнее реализовать с помощью быстрого преобразования Фурье (БПФ). Однако он более широко используется, потому что он находится на сетке экстремумов, которая имеет тенденцию быть наиболее полезной для краевых задач. В основном потому, что к этой сетке проще применить граничные условия.

Существует дискретный (и на самом деле быстрый, потому что он выполняет dct с помощью быстрого преобразования Фурье), доступный при обмене файлами MATLAB, который был создан Грегом фон Винкелем. Поэтому здесь это опущено.

В этом случае преобразование и его обратное

${ displaystyle u (x_ {n}) = u_ {n} = sum _ {m = 0} ^ {N} a_ {m} T_ {m} (x_ {n})}$

${ displaystyle a_ {m} = { frac {p_ {m}} {N}} left [{ frac {1} {2}} (u_ {0} (- 1) ^ {m} + u_ { N}) + sum _ {n = 1} ^ {N-1} u_ {n} T_ {m} (x_ {n}) right]}$

где ${ displaystyle p_ {m} = 1 Leftrightarrow m = 0}$ и ${ displaystyle p_ {m} = 2}$ в противном случае.

Использование и реализации

Основное использование дискретного преобразования Чебышева - численное интегрирование, интерполяция и устойчивое численное дифференцирование.^[1]Реализация, обеспечивающая эти функции, приведена в C ++ библиотека Boost^[2]

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

• Полиномы Чебышева
• Дискретное косинусное преобразование
• Дискретное преобразование Фурье
• Список преобразований, связанных с Фурье

использованная литература