正弦信号整周期采样

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

正余弦信号整周期采样的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;%设定余弦信号频率

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('余弦信号时域波形');

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('余弦信号幅频谱图');

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('余弦信号相频谱图');

grid;

angle(y(6))*180/pi %求信号初时相位。频率坐标f为[0 0.01 0.02 0.03 0.04 0.05 0.06 ...],所以谱线y中第6根谱线和信号x对应。

相关文档
最新文档