Fouirer katsayılarını bulma

[seyityildirim]

fourier katsayıları ak ve bk bulurken aşağıdaki kod kullanılmış daha sonra fonksiyonun kendisi ve fourier eşdeğeri çizdirilmiş. benim anlamadığım nokta şu: ak ve bk katsayılarını bulurken bir periyot üzerinden integral alınır. ancak aşağıdaki kodlarda quadl, quad ya da herhangi bir integral komutu kullanılmamış ama çizime baktığımızda hesapları doğru. integral almadan bu durum nasıl halledilmiş?

clear all
close all
N=50;       % işaretin uzunluğu
n=0:N-1;    % 50 uzunluklu dizi
M=10;        % alınacak fourier serisi katsayı sayısı
m=0:M-1;
x=[zeros(1,28),ones(1,12),zeros(1,N-40)];
a=(1/N)*x*cos(2*pi*(n')*m/N);
b=(1/N)*x*sin(2*pi*(n')*m/N);
xr=a*cos(2*pi*(m')*n/N)+b*sin(2*pi*(m')*n/N);
subplot(2,1,1), plot(x,'LineWidth',3); grid; axis([0 50 -1 +2]);
subplot(2,1,2), plot(xr,'r','LineWidth',3); grid;

Cevap verecek olanlara şimdiden teşekkür ederim

[alicavuslu]

Burada yaptığı olay sinc ifadesini elde etmek amacı ile zamanla sönümlenen sinyalleri elde etmeye çalışmış. Kare dalganın frekans donusumu sinc'dir

http://math.stackexchange.com/questions/397586/intuitive-way-to-understand-the-square-wave-spectrum

Bu kodda sanki sinüslerden direk kare dalga oluşturulmak istenmiş gibi duruyor.

[seyityildirim]

Anladım..Allah razı olsun hocam. sayenizde hergün birşey öğreniyoruz