方波的ＦＦＴ变换
我想做一个周期为2*pi的方波，然后做ＦＦＴ变换，方波图形画出来后，做ＦＦＴ变换时，图象如下，请问这样的图象对么？程序如下：
f=5;
t=0:0.001:10;
x=square(f*t);
figure(1);
plot(x);
figure(2);
Y=fft(x,1024);
pyy=Y.*conj(Y)/1024;
plot(pyy(1:100));ylabel('幅值Amm');title('x的FFT频谱');
ＦＦＴ图如下：
没问题 %%%周期为2*pi的方波
w=1;
t=0:0.001:10;
x=square(w*t);
figure(1);
plot(t,x);
%%%谱分析时要提高分辨率否则看不到1/2/pi处有个峰值。
figure(2);
psd(x,1024*16,1000,[]);
[ 本帖最后由 w89986581 于 2007-6-15 14:41 编辑 ] 对于周期信号，尽量整周期采样。
因为DFT包含了时域、频域的周期延拓过程，不是整周期采样的话，延拓出来的信号与原周期信号差别可能很大，那你得到的频谱包络形状可能和原周期信号的谱线包络相去甚远，而不只是一个常数倍的差别。 我想问的是，我怎么样才能把我所采集的的信号转成FFT的谱？
[ 本帖最后由 zhlong 于 2007-6-22 18:39 编辑 ]
回复 #5 xyxiang1 的帖子
请搜索本版有关FFT的帖子。
页:
[1]