|
![](static/image/common/ico_lz.png)
楼主 |
发表于 2014-9-17 22:02:07
|
显示全部楼层
本帖最后由 imjacob 于 2014-9-17 22:10 编辑
自己参照网上弄了个程序,基本满足我的要求:
close all; %先关闭所有图片
F1=50; %信号1频率(Hz)
Fs=256; %采样频率(Hz)
N=256; %采样点数
t=[0:1/Fs:N/Fs]; %采样时刻
%信号
S=sin(2*pi*F1*t);
%显示原始信号
plot(S);
title('原始信号');
Y = fft(S,N); %做FFT变换
Ayy = (abs(Y)); %取模
figure;
Ayy=Ayy/(N/2); %换算成实际的幅度
Ayy(1)=Ayy(1)/2;
F=([1:N]-1)*Fs/N; %换算成实际的频率值
plot(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果
title('幅度-频率曲线图');
figure;
Pyy=[1:N/2];
for i=1:N/2
Pyy(i)=phase(Y(i)); %计算相位
Pyy(i)=Pyy(i)*180/pi; %换算为角度
end;
plot(F(1:N/2),Pyy(1:N/2)); %显示相位图
title('相位-频率曲线图'); |
|