频谱分析

n = 1024 * 1;                                  % cai yang dian shu
fs = 1000;                                      % cai yang ping lv
t = (1 : 1 : n) / fs;                           % shi jian
signalfre = 100;                                % xin hao ping lv
jiasudu = 1000 * sin(2 * pi * signalfre * t);	% jia su du
ff = fft(jiasudu, n);                           % fft
fre = ((1 : 1 : n) - 1) * fs / n;               % fre
amp = abs(ff) / (n / 2);                        % amp
xiangwei = angle(ff) * 180 / pi;
subplot(2, 1, 1);
plot(fre(1:n/2), amp(1:n/2));
subplot(2, 1, 2);
plot(fre(1:n/2), xiangwei(1:n/2));