function [A,P,A1,P1,A0,P0,y]=qiuzhi(S,N,num,den)
y=filter(num,den,S);
%y=awgn(y,20);
y=y+normrnd(0,0.06,1,N+1);
ham=hamming(N+1);%加窗函數
%S=S.*ham';
%y=y.*ham';
Y = fft(S,N); %做FFT
Y0=fft(y,N);
A1 =abs(Y(1:N/2)); %取模
A1=A1/(N/2);   %換算成實際的幅度加%%
A1(1)=A1(1)/2;
P1=[1:N/2];
for i=1:N/2
 P1(i)=phase(Y(i)); %計算相位
 P1(i)=P1(i)*180/pi; %換算為角度
end;
A0 = abs(Y0(1:N/2)); %取模
A0=A1/(N/2);   %換算成實際的幅度
A0(1)=A1(1)/2;
for i=1:N/2   
 P0(i)=phase(Y0(i)); %計算相位
 P0(i)=P0(i)*180/pi; %換算為角度
end;
A=[1:N/2];
P=[1:N/2];
for i=1:N/2
 Y(i)=Y0(i)/Y(i);
end;
A=abs(Y(1:N/2)); %取模
for i=1:N/2   
 P(i)=phase(Y(i)); %計算相位
 P(i)=P(i)*180/pi; %換算為角度
end;