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;
        }
相关推荐
AI+程序员在路上19 分钟前
Qt6中模态与非模态对话框区别
开发语言·c++·qt
nbsaas-boot5 小时前
Java 正则表达式白皮书:语法详解、工程实践与常用表达式库
开发语言·python·mysql
chao_7895 小时前
二分查找篇——搜索旋转排序数组【LeetCode】两次二分查找
开发语言·数据结构·python·算法·leetcode
风无雨6 小时前
GO 启动 简单服务
开发语言·后端·golang
斯普信专业组6 小时前
Go语言包管理完全指南:从基础到最佳实践
开发语言·后端·golang
我是苏苏7 小时前
C#基础:Winform桌面开发中窗体之间的数据传递
开发语言·c#
斐波娜娜8 小时前
Maven详解
java·开发语言·maven
小码氓8 小时前
Java填充Word模板
java·开发语言·spring·word
暮鹤筠8 小时前
[C语言初阶]操作符
c语言·开发语言
Boilermaker199210 小时前
【Java EE】Mybatis-Plus
java·开发语言·java-ee