matlab频率和幅值图,请教！傅里叶变换频率和幅值对不上！

拜托各位大哥帮忙看看，这个程序运行了之后幅值和应有频率对不上啊？自己找不到问题。而且相位图也不对。

clc;

close all;

clear;

A1=2.0; %频率F1信号的扰动幅度

A2=1.9; %频率F2信号的扰动幅度

F1=50; %信号1频率(Hz)

F2=50; %信号2频率(Hz)

Fs=512; %采样频率(Hz)

P1=30; %信号1相位(度)

P2=30.001; %信号2相位(度)

N=12800; %采样点数

t=[0:1/N:1]; %采样时刻

%信号

H=220*sqrt(2)*sin(2*pi*F1*t+pi/6)-A1*sqrt(2)*sin(2*pi*1.2*F1*t+P1/180);

H1=220*sqrt(2)*sin(2*pi*F2*t+pi/6)-A2*sqrt(2)*sin(2*pi*1.2*F2*t+P2/180);

%加窗

S=hanning(N)'.*H(1:N);

S1=hanning(N)'.*H1(1:N);

figure(1)

%显示原始信号

subplot(321);plot(t,H);

title('标准信号');

subplot(322);plot(t,H1);

title('被测信号');

%信号加窗

subplot(323);plot(t(1:N),S(1:N));

title('标准信号加窗');

subplot(324);plot(t(1:N),S1(1:N));

title('被测信号加窗');

Y = fft(S,N); %做FFT变换

Y1=fft(S1,N);

fuzhi = (abs(Y)); %取模

fuzhi1=(abs(Y1));

fuzhi=fuzhi*2/(N/2);   %换算成实际的幅度

F=([1:N]-1)*Fs/N; %换算成实际的频率值

subplot(325);stem(F(1:N/2),fuzhi(1:N/2));   %显示换算后的FFT模值结果

title('幅度-频率曲线图');

fuzhi1=fuzhi1*2/(N/2);   %换算成实际的幅度

F=([1:N]-1)*Fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/N

subplot(326);stem(F(1:N/2),fuzhi1(1:N/2));   %显示换算后的FFT模值结果

title('幅度-频率曲线图');

Pyy=[1:N/2];

for i=1:N/2

xiangwei(i)=angle(Y(i)); %计算相位

xiangwei(i)=xiangwei(i)*180/pi; %换算为角度

end;

figure(2)

Hline1=plot(F(1:N/2),xiangwei(1:N/2));   %显示相位图

hold on;

Pyy1=[1:N/2];

for i=1:N/2

xiangwei1(i)=angle(Y1(i)); %计算相位

xiangwei1(i)=xiangwei1(i)*180/pi; %换算为角度

end;

Hline2=plot(F(1:N/2),xiangwei1(1:N/2),'r');   %显示相位图

title('相位-频率曲线图');