这个图是怎么画出来的？

lzh19840809

我现在手里有51200个样本数据，然后我要按下图那个样子画出来它FFT变换后的频谱图和相频图。频谱图我已经画好了，但是那个相频图总感觉不太对，不知道它是怎么画出来的。大家这方面的经验丰富，所以想请教大家这个相频图怎么画？
对了，它那个图不是用matlab画得，而我的是用matlab画得。

lzh19840809

有些不明白的地方，就是它的那个纵坐标好象是Phase，也就是相位，但是刻度怎么会是600000 800000 1000000呢？我的相谱图跟它的一点都不一样。

这是我的源代码
clc
clear
close all
SampleData=load('3.txt')
x=SampleData(1:51200)
x=x-mean(x)
subplot(311)
plot(x)
title('51200个样本点，fs=100')
grid on
%---------------------------------------
%FFT变换
%---------------------------------------
fs=100
N=51200
n=0:N-1;
y=fft(x,N);%进行fft变换
mag=abs(y);%求幅值

f=n*fs/length(y);%进行对应的频率转换
subplot(312);
plot(f(1:N/2),mag(1:N/2));%做频谱图
title('幅值-频率曲线图')
grid on

Pyy=[1:N/2];
for i=1:N/2
Pyy(i)=phase(y(i)); %计算相位
Pyy(i)=Pyy(i)*180/pi; %换算为角度
end;
subplot(313);
plot(f(1:N/2),Pyy(1:N/2));    %显示相位图
title('相位-频率曲线图');
grid on

wfloveu

第一个图是用origin画出来的