Выполнение преобразования Лапласа над непрерывной функцией в максимумах - PullRequest
1 голос
/ 23 февраля 2020

Я вычисляю преобразование Лапласа следующей функции, но получаю странный вывод с limit и realpart et c на (%t5) и (%t6). Я ожидаю что-то вроде a_n = 0 и b_n = -1/(n*%pi).

Почему это происходит? Я неправильно определяю функцию на (%i1)? Или это ограничение максимумов?

function plot

(%i1) f(t) := mod(t, 1);
(%o1)                          f(t) := mod(t, 1)
(%i2) plot2d(f(t), [t, -2, 2]);
(%o2)                 [/tmp/maxout1266174.gnuplot_pipes]
(%i3) load(fourie);
(%o3)         /usr/share/maxima/5.43.0/share/calculus/fourie.mac
(%i4) fourier(f(t), t, 1);
                                         1
(%t4)                               a  = -
                                     0   2

                    /
                    [
(%t5) a  = limit    I cos(%pi n t) realpart(floor(t)) dt
       n   t -> - 1 ]
                    /
                                           /
                                           [
                                  - limit  I cos(%pi n t) realpart(floor(t)) dt
                                    t -> 1 ]
                                           /

                     /
                     [
(%t6) b  = (- limit  I sin(%pi n t) realpart(floor(t)) dt)
       n      t -> 1 ]
                     /
             /
             [                                      2 sin(%pi n)   2 cos(%pi n)
  + limit    I sin(%pi n t) realpart(floor(t)) dt + ------------ - ------------
    t -> - 1 ]                                           2  2         %pi n
             /                                        %pi  n

(%o6)                           [%t4, %t5, %t6]
(%i7) 
...