正弦信号整周期采样
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
正xx信号整周期采样的fft变换
2010-01-28 10:53
fs=1;
N=100;%频率分辨率为fs/N＝
0.01Hz，下面信号的频率
0.05是
0.01的整数倍，即为整周期采样
n=0:N-1;
t=n/fs;
f0=
0.05;%设定xx信号频率
x=cos(2*pi*f0*t);%生成正弦信号%FFT是余弦类变换，最后得到的初始相位是余弦信号的初时相位，在这里为0。
如果信号
figure
(1); %为x=sin(2*pi*f0*t);则初时相位应该是－90度而非0度。
subplot
(311);
plot(t,x);%作余弦信号的时域波形
xlabel('t');
ylabel('y');
title('xx信号时域波形');
grid;
%进行FFT变换并做频谱图
y=fft(x,N);%进行fft变换
mag=abs(y)*2/N;%求幅值乘上后面的2/N得到正确幅值f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换subplot
(312);
stem(f(1:N/2),mag(1:N/2));%做频谱图
xlabel('频率(Hz)');
ylabel('幅值');
title('xx信号幅频谱图');
grid;
phase=angle(y);%求幅值乘上后面的2/N得到正确幅值f=(0:length(y)-1)'*fs/length(y);%进行对应的频率转换subplot
(313);
stem(f(1:N/2),phase(1:N/2));%做频谱图
xlabel('频率(Hz)');
ylabel('相位');
title('xx信号相频谱图');
grid;
angle(y
(6))*180/pi %求信号初时相位。
频率坐标f为[0 0.01
0.02
0.03
0.04
0.05
0.06 ...]，所以谱线y中第6根谱线和信号x对应。