C# 实现傅里叶变化(DFT)

1、DFT函数类

csharp 复制代码
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Threading.Tasks;

namespace DFT_FFTApp.Utils
{
    public class DFT
    {
        /// <summary>
        /// DFT
        /// </summary>
        /// <param name="data"></param>
        /// <returns></returns>
        public static List<List<double>> DFTNative(List<double> data)
        {
            int N = data.Count;
            List<List<double>> xks = new List<List<double>>();
            for (int k = 0; k < N; k++)
            {
                List<double> xk = DFTPoint(k, data);
                xks.Add(xk);
            }
            GC.Collect();
            return xks;
        }
        /// <summary>
        /// 单点DFT
        /// </summary>
        /// <param name="k"></param>
        /// <param name="data"></param>
        /// <returns></returns>
        public static List<double> DFTPoint(int k, List<double> data)
        {
            int N = data.Count;
            double real = 0;
            double image = 0;
            for (int n = 0; n < N; n++)
            {
                real += data[n] * Math.Cos(2 * Math.PI * (k * n) / N);
                image -= data[n] * Math.Sin(2 * Math.PI * (k * n) / N);
            }
            List<double> res = new List<double>() { real, image };
            return res;
        }
        /// <summary>
        /// 频域求模
        /// </summary>
        /// <param name="data"></param>
        /// <returns></returns>
        public static List<double> GetMod(List<List<double>> data)
        {
            List<double> res= new List<double>();
            foreach(var  x in data)
            {
                res.Add(Math.Sqrt(x[0] * x[0] + x[1] * x[1]));
            }
            return res;
        }
        /// <summary>
        /// 频域求相位
        /// </summary>
        /// <param name="data"></param>
        /// <returns></returns>
        public static List<double> GetPhase(List<List<double>> data)
        {
            List<double> res = new List<double>();
            foreach (var x in data)
            {
                double phase = 0;
                if (x[0] == 0)
                {
                    phase = 90;
                }
                else
                {
                    phase = Math.Atan2(x[1], x[0]);
                }
                res.Add(phase);
            }
            return res;
        }

    }
}

2、应用

csharp 复制代码
private void DFTApp()
        {
            List<double> list = GetCaseData(128);
            List<List<double>> Xk = DFT.DFTNative(list);
            List<double> res = DFT.GetMod(Xk);
            foreach (double item in res)
            {
                Console.WriteLine(item);
            }
            easyChartX1.Plot(res, 0, 1);
        }
/// <summary>
        /// 获取时域数据
        /// </summary>
        /// <param name="points"></param>
        /// <returns></returns>
        private List<double> GetCaseData(int points)
        {
            List<double> data = new List<double>();
            double w = 2*Math.PI / points;
            for(int i = 0; i < points;i++)
            {
                double d=Math.Sin(w * i);
                //Console.WriteLine($"i={i},d={d.ToString("G2")}");
                data.Add(d);
            }
            return data;
        }
相关推荐
apihz14 分钟前
域名WHOIS信息查询免费API使用指南
android·开发语言·数据库·网络协议·tcp/ip
coding随想28 分钟前
掌控网页的魔法之书:JavaScript DOM的奇幻之旅
开发语言·javascript·ecmascript
爱吃烤鸡翅的酸菜鱼1 小时前
IDEA高效开发:Database Navigator插件安装与核心使用指南
java·开发语言·数据库·编辑器·intellij-idea·database
神仙别闹1 小时前
基于C#+SQL Server实现(Web)学生选课管理系统
前端·数据库·c#
心情好的小球藻1 小时前
Python应用进阶DAY9--类型注解Type Hinting
开发语言·python
惜.己2 小时前
使用python读取json数据,简单的处理成元组数组
开发语言·python·测试工具·json
Y4090012 小时前
C语言转Java语言,相同与相异之处
java·c语言·开发语言·笔记
向宇it2 小时前
【unity组件介绍】URP Decal Projector贴花投影器,将特定材质(贴花)投影到场景中的其他对象上。
游戏·3d·unity·c#·游戏引擎·材质
古月-一个C++方向的小白7 小时前
C++11之lambda表达式与包装器
开发语言·c++
沐知全栈开发7 小时前
Eclipse 生成 jar 包
开发语言