推 psallen:太感謝了!我再研究看看 謝謝 02/21 23:09
http://www.clarku.edu/~djoyce/trig/identities.html
※ 引述《psallen (psallen)》之銘言:
: short FFT(short int dir,long m,double *x,double *y)
: {
: /* Compute the FFT */
: c1 = -1.0;
: c2 = 0.0;
c1 c2分別是cos(-180度)和sin(-180度)
: l2 = 1;
: for (l=0;l<m;l++) {
: l1 = l2;
: l2 <<= 1;
: u1 = 1.0;
: u2 = 0.0;
u1 u2是cos0與sin0
: }
: z = u1 * c1 - u2 * c2;
: u2 = u1 * c2 + u2 * c1;
: u1 = z;
和角公式 cos(A+B) 與 sin(A+B)
: }
: c2 = sqrt((1.0 - c1) / 2.0);
: if (dir == 1)
: c2 = -c2;
: c1 = sqrt((1.0 + c1) / 2.0);
半角公式 sin(A/2) 與 cos(A/2)
: }
: }
可能會有搞錯的地方,請你再檢查看看。不好意思了
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 61.230.126.99
※ 編輯: DJWS 來自: 61.230.126.99 (02/21 20:27)
※ 編輯: DJWS 來自: 61.230.126.99 (02/21 20:45)
其實他還是有用sin和cos
只是換成了實際數字
這是三角函數公式表 可以參考一下