ft.m

在matlab环境下利用追赶法解线性方程组的算法实现

function ft(N)
%x(n)为输入信号;N为所取点的个数
y=zeros(1,N);
for n=1:N
    x(n)=cos(0.48*pi*n)+cos(0.52*pi*n);
end
for k=1:N
    for m=1:N
        y(k)=y(k)+x(m)*exp(-j*2*pi*(k-1)*(m-1)/N);
        z(k)=norm(y(k));
    end
end
plot(z)