Разложение временного ряда в ряд Фурье matlab

[vintaj~]

Написал разложение в ряд Фурье временного ряда, к сожалению, закралась ошибка. Помогите, пожалуйста, найти. Заранее, спасибо, всем откликнувшимся.

A[];%временной ряд(2556 значений)
To=2556; time=1:1:length(a);
a0=(1/To)*trapz(a);%коэф a0(+)
a1=zeros(length(a),length(a));
for x=1:1:length(a)%массив коэффицентов при cos
    for n=1:1:200
    a1(x,n)=(2/To)*trapz(a.*cos(2*w0*x*n/To));
    end
end
a2=zeros(length(a),length(a));
for x=1:1:length(a)%массив коэффицентов при sin
    for n=1:1:200
    a2(x,n)=(2/To)*trapz(a.*sin(2*w0*x*n/To));
    end
end
yras=zeros(1,length(a));
mas=zeros(1,length(a));
for x=1:1:length(a)
    for n=1:1:200
        mas(x)=mas(x)+a1(x,n)*cos(2*w0*x*n/To)+a2(x,n)*sin(2*w0*x*n/To);
    end
end

for k=1:1:length(a)%ряд Фурье
 yras(k)=a0+mas(k);
end
plot(time,a,time,yras);

а вот результат: