FFT得到的是什么

whutcl8110

我想对发动机的示功图进行简谐分析，
Pg=[3.93        3.75        3.34        2.82        2.34        1.89        1.49        1.2        0.97        0.78        0.64        0.55        0.47        0.4        0.35        0.31        0.29        0.26        0.24        0.23        0.225        0.22        0.215        0.21        0.205        0.2        0.195        0.19        0.185        0.18        0.17        0.165        0.16        0.155        0.15        0.145        0.14        0.135        0.13        0.125        0.12        0.115        0.11        0.105        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.1        0.115        0.13        0.145        0.16        0.175        0.19        0.205        0.22        0.235        0.25        0.26        0.27        0.28        0.29        0.3        0.31        0.32        0.33        0.34        0.35        0.355        0.36        0.365        0.37        0.375        0.38        0.385        0.39        0.395        0.4        0.405        0.41        0.415        0.42        0.425        0.43        0.435        0.44        0.445        0.45        0.45        0.59        0.75        0.93        1.15        1.41        1.76        2.22        3.41        3.87];
Pg是气缸压力在曲轴720度内对应的值.
Pt=Pg.*(sin(afa)+0.25/2*sin(2*afa));  afa是角度
Pt是压力作用到曲柄上的值

P=fft(Pt);
P是什么？

我想得到：
每一个谐次的幅值C1，C2……  (Pt=P0+C1sin(wt+a)+C2sin(2wt+b)+……)

请指教！

whutcl8110

请各位高手多多指教！
感谢感谢！

gh688

气缸压力在曲轴720度内对应的值 Pg     size(Pg)=160
压力作用到曲柄上的值           Pt     Pt=Pg.*(sin(afa)+0.25/2*sin(2*afa));  
这里你的正弦函数sin(afa)和sin(2*afa))的角度afa就相当于我们普通FFT中的时间，由于你有Pt=Pg.*(sin(afa)+0.25/2*sin(2*afa));  这个点乘的限制所以afa=0:4*pi/(160-1):4*pi这样size(afa)=160。
这样压力作用到曲柄上的值就求出来了（在720度范围内，采样间隔4*pi/(160-1)）。如果你要批P=FFT(Pt)，那你得到的P就是角度域afa变换到频域的一系列数组，这里你的角度就相当于时间了。

我想得到：
每一个谐次的幅值C1，C2……  (Pt=P0+C1sin(wt+a)+C2sin(2wt+b)+……)

这个地方符号也没有解释一下没太看懂，请解释一下