من بررسی کردم. دستور fft را کنار بگذارید و باید از طریق همان سری فوریه عمل کنیم. من قبلا محاسبه سری فوریه با استفاده از عبارت تابع را در یک مبحث نوشته بودم، اما با توجه به پرسش شما، کلید زیر را نوشتم :

http://www.kelidestan.com/keys/keys.php?key=2548

در کلید فوق، با استفاده از مقادیر تابع در یک بازه (نه عبارت تابع)، مقادیر سری فوریه را به صورت تقریبی به دست می آوریم (به هر حال، همینکه داده ها گسسته هستند، خود به خود تقریب ایجاد می کند). بنابراین با استفاده از کدهای نوشته شده در کلید فوق، ضرایب سری فوریه را به دست آورید و سپس بر اساس آن، مولفه اصلی (هارمونیک اصلی) را رسم کنید.

در ضمن، لینک کلیدی که در مورد به دست آوردن ضرایب سری فوریه با استفاده از خود عبارت تابع بود رو هم براتون اینجا قرار می دهم :

http://www.kelidestan.com/keys/keys.php?key=103