Реализация импульсного отклика с двойным суммированием

Question

Я пытаюсь реализовать следующую импульсную характеристику для фильтра.

Когда t постоянно, выходной сигнал верный, но когда он переменный, как вimage (от nT до (n + 1) T) я не могу получить правильный вывод.

Это код, который я написал -

a=3;
b=4;

H = zeros(1,14);

c=[9 10];

for m=0:2
 for n=0:1
   for t=1:14
    H(t)=H(t)+(c(n+1)*((a*((t-(n*T))/T)+b).^m));
   end
 end  
end

disp(H(4));