کد نویسی تبدیل فوریه گسسته در متلب

نیازمندم که کد fft یا تبدیل فوریه سریع گسسته در متلب نوشته شود. یک سیگنال هارمونیکی به نرخ نمونه برداری 40 کیلوهرتز(زمان نمونه برداری 25 میکرو ثانیه) در مدت 20 میلی ثانیه دریافت شده است. یعنی برداری حاوی 800 نمونه موجود است. باید به صورت کد نویسی fft گسسته پیاده سازی شود و دامنه و فاز هر کدام از هارمونیک ها استخراج شود. هارمونیک اصلی این سیگنال 50 هرتز است و هارمونیک های 150،350 ، 750 و 950 در این شکل موج قرار دارند که تبدیل باید دامنه و فاز هر فرکانس را به صورت جدا گانه استخراج کند. کد باید به صورت ماژولار نوشته شود و تغییر کد به منظور تغییر فرکانس های انتخابی به سادگی امکان پذیر باشد. نکته مهم عدم استفاده از توابع آماده داخل متلب میباشد. صرفا با اعمال عملیات ریاضی پایه مثل جذر، حمع ، تفریق، ضرب و... تبدیل باید نوشته باشد. استفاده از محاسبات اعداد موهومی نیز مجاز نمیباشد. همچنین در صورتی که نیاز به عملیات ماتریسی(ضرب ماتریس ،دترمینان یا معکوس ماتریس) است باید با استفاده از حلقه ها و بدون استفاده از توابع و سینتکس های آماده متلب محاسبه شود. نکته مهم بعدی سرعت اجرای سریع است. زمان اجرای کد نباید زیاد طول بکشد.در نتیجه استفاده از الگوریتم های بازگشتی پیچیده و زمان بر امکان پذیر نمیباشد. حداکثر نرخ نمونه برداری 40 کیلو هرتز و برداری حاوی 800 نمونه می باشد.به منظور افزایش سرعت اجرا میتوان تعداد نمونه و نرخ نمونه برداری را کاهش داد به شرطی که دقت پاسخ نهایی قابل قبول باشد. نوع الگوریتم انتخابی برای محاسبه اهمیتی ندارد میتوانید از روش بلوستین،ریدر ، کولی-توکی یا هر الگوی جدیدی استفاده کنید. مهم سرعت و دقت محاسبه دامنه و فاز میباشد. لطفا توضیحات رو کامل بخونید و در صورتی که توانایی انجام این پروژه دارین مبلغ و زمان پیشنهادی خودتون رو ارسال کنید.