在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)