Question

fft求直流分量值和一次谐波值

对信号的一个周期分别求w=0和w=4.29 （f=0和f=0.68）时的直流分量值和一次谐波的值。
信号为周期T=1.46 振幅A=1的矩形波
闭关两天终于弄出了点头绪，但结果不对。求大神指出程序错误。小女子不胜感激！

>> %%%频率为0时%%%%
>> clear
>> f=0; %角频率为零时的频率
>> N=32; %采样点数
>> fs=0; %采样频率=f*N
>> t=[0:1/fs:N/fs]; %采样时刻
>> u=square(2*pi*fs*t); %振幅为1的矩形波
>> u1=fft(u,N) %FFT变换

u1 =
Columns 1 through 9
-1 -1 -1 -1 -1 -1 -1 -1 -1
Columns 10 through 18
-1 -1 -1 -1 -1 -1 -1 -1 -1
Columns 19 through 27
-1 -1 -1 -1 -1 -1 -1 -1 -1
Columns 28 through 32
-1 -1 -1 -1 -1

>> %%%频率为0.68时%%%%
>> clear
>> f=0.68; %信号频率
>> N=32; %采样点数
>> fs=21.76; %采样频率
>> t=[0:1/fs:N/fs]; %采样时刻
>> y=square(2*pi*fs*t); %信号
>> y1=fft(y,N) %fft变换

y1 =
Columns 1 through 9
32 0 0 0 0 0 0 0 0
Columns 10 through 18
0 0 0 0 0 0 0 0 0
Columns 19 through 27
0 0 0 0 0 0 0 0 0
Columns 28 through 32
0 0 0 0 0

Answer 1

才看到问题，希望已经解决。

如果还不太清楚，首先需要知道，square产生的矩形波，本身不包含直流分量，是由基频正弦波和其奇次谐波组成的。所以你第一部分没有意义。

我对你问题的理解是：有一个信号，周期1.46s，幅值为正负1的矩形波。那么这个矩形波一定没有直流分量，如要求其交流分量，对其进行FFT看个频率信号的幅值。但是你把几个名词弄混了，

1. 信号时信号，采样是采样，在你的代码中，y=square(2*pi*fs*t); ，不应该是fs而应改为f
   

2. FFT得到的结果是对应频域的信息，需要继续分析，详细物理意义推荐你看http://insomina.blog.163.com/blog/static/18315823320117165032823/

3. 如果相加直流分量，在信号表达式中加入，如y=1+square (2*pi*f*t)
   

Answer 2

都是胡说八道，佛说一切都是虚妄之相，汝等错误看待此世界，没有好好理解，傅里叶变换，离散傅里叶变换的基本概念，故而有次迷惑。还有一堆迷惑众生还给你作答，尔如同蚕，自作自缚。诱导别人犯错。阿弥陀佛

Answer 3

但结果不对