离散哈特莱变换(DHT)用C语言实现

2024-03-02 22:58

本文主要是介绍离散哈特莱变换(DHT)用C语言实现,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

离散哈特莱变换(DHT)

摘 要

离散哈特莱变换(DHT)是一种与傅里叶变换相关的转换。类似于离散傅里叶变换。与傅里叶变换在信号处理及其他相关领域有相似的应用。
本设计介绍了DHT的定义以及使用C语言实现其算法。

关键字:傅里叶变换 哈特莱变换

二、设计平台

Linux平台、GCC编译器、VIM、windows、Visual Studio 2017

三、设计原理

DHT的定义
设x(n),n=0,1,…,N−1x(n),n=0,1,…,N−1,为一实序列,其DHT定义为:

avatar

式中cas(a)=cos(a)+sin(a)cas(a)=cos(a)+sin(a)逆变换(IDHT)为:

avatar

DHT的正交证明:

avatar

DHT和DFT关系
用X(k)表示实序列x(n)的,用XH(k)表示x(n)的DHT,分别用XHe(k), XHo(k)表示XH(k)的偶对称分量与奇对称分量,即:

avatar

其中:

avatar

DHT的优点

DHT为实值,避免了复数运算DHT正反变换形式基本一致DHT与DFT的转换容易实现

DHT的性质

DHT的性质与DFT的性质类似,但由于DHT是实序列间的变换,有些性质有具体的表达形式。这里只给出结论。
设x(n)、y(n)的DHT分别为Xh(k)、Yh(k)。用符号x(n)↔Xh(k)表示Xh(k)=DHT(x(n))。

1.线性性
avatar

2.逆序列x(N-n)的DHT
avatar


avatar

当k=0时,可得Xh(N)=Xh(0)。

3.循环位移的性质
avatar

四、实现代码

【C语言】

#include <stdio.h>
#include <stdlib.h>
#include <math.h>const int N = 1024;
const float PI = 3.1416;inline void swap(float &a, float &b)
{float t;t = a;a = b;b = t;
}void bitrp(float xreal[], float ximag[], int n)
{// 位反转置换 Bit-reversal Permutationint i, j, a, b, p;for (i = 1, p = 0; i < n; i *= 2){p++;}for (i = 0; i < n; i++){a = i;b = 0;for (j = 0; j < p; j++){b = (b << 1) + (a & 1);    // b = b * 2 + a % 2;a >>= 1;        // a = a / 2;}if (b > i){swap(xreal[i], xreal[b]);swap(ximag[i], ximag[b]);}}
}void FFT(float xreal[], float ximag[], int n)
{// 快速傅立叶变换,将复数 x 变换后仍保存在 x 中,xreal, ximag 分		//	别是 x 的实部和虚部float wreal[N / 2], wimag[N / 2], treal, timag, ureal, uimag, arg;int m, k, j, t, index1, index2;bitrp(xreal, ximag, n);// 计算 1 的前 n / 2 个 n 次方根的共轭复数 W'j = wreal [j] + i * 		wimag [j] , j = 0, 1, ... , n / 2 - 1arg = -2 * PI / n;treal = cos(arg);timag = sin(arg);wreal[0] = 1.0;wimag[0] = 0.0;for (j = 1; j < n / 2; j++){wreal[j] = wreal[j - 1] * treal - wimag[j - 1] * timag;wimag[j] = wreal[j - 1] * timag + wimag[j - 1] * treal;}for (m = 2; m <= n; m *= 2){for (k = 0; k < n; k += m){for (j = 0; j < m / 2; j++){index1 = k + j;index2 = index1 + m / 2;t = n * j / m;    // 旋转因子 w 的实部在 wreal [] 中					//的下标为 ttreal = wreal[t] * xreal[index2] - wimag[t] * 						ximag[index2];timag = wreal[t] * ximag[index2] + wimag[t] * 						xreal[index2];ureal = xreal[index1];uimag = ximag[index1];xreal[index1] = ureal + treal;ximag[index1] = uimag + timag;xreal[index2] = ureal - treal;ximag[index2] = uimag - timag;}}}
}void FFT_test()
{float xreal[N] = {}, ximag[N] = {};int n = 8;int i = 0;printf("请输入数据,格式(实部 虚部) : \n");for (i = 0; i < 8; i++){scanf("%f%f",xreal+i,ximag+i);}n = i;    // 要求 n 为 2 的整数幂while (i > 1){if (i % 2){printf("%d is not a power of 2! ", n);}i /= 2;}FFT(xreal, ximag, n);printf("=====================================\n");printf("FFT:    i	      实部       虚部 \n");for (i = 0; i < n; i++){printf("     %4d       %8.4f    %8.4f ", i+1, xreal[i], 			ximag[i]);printf("\n");}printf("===================================== \n");printf("DHT:    i	    结果 \n");for (i = 0; i < n; i++){printf("     %4d        %8.4f ", i+1, xreal[i]-ximag[i]);printf("\n");}printf("=====================================\n ");
}int main()
{FFT_test();system("pause");return 0;
}

END

这篇关于离散哈特莱变换(DHT)用C语言实现的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/767690

相关文章

MySQL中查找重复值的实现

《MySQL中查找重复值的实现》查找重复值是一项常见需求,比如在数据清理、数据分析、数据质量检查等场景下,我们常常需要找出表中某列或多列的重复值,具有一定的参考价值,感兴趣的可以了解一下... 目录技术背景实现步骤方法一:使用GROUP BY和HAVING子句方法二:仅返回重复值方法三:返回完整记录方法四:

IDEA中新建/切换Git分支的实现步骤

《IDEA中新建/切换Git分支的实现步骤》本文主要介绍了IDEA中新建/切换Git分支的实现步骤,通过菜单创建新分支并选择是否切换,创建后在Git详情或右键Checkout中切换分支,感兴趣的可以了... 前提:项目已被Git托管1、点击上方栏Git->NewBrancjsh...2、输入新的分支的

Python实现对阿里云OSS对象存储的操作详解

《Python实现对阿里云OSS对象存储的操作详解》这篇文章主要为大家详细介绍了Python实现对阿里云OSS对象存储的操作相关知识,包括连接,上传,下载,列举等功能,感兴趣的小伙伴可以了解下... 目录一、直接使用代码二、详细使用1. 环境准备2. 初始化配置3. bucket配置创建4. 文件上传到os

关于集合与数组转换实现方法

《关于集合与数组转换实现方法》:本文主要介绍关于集合与数组转换实现方法,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、Arrays.asList()1.1、方法作用1.2、内部实现1.3、修改元素的影响1.4、注意事项2、list.toArray()2.1、方

使用Python实现可恢复式多线程下载器

《使用Python实现可恢复式多线程下载器》在数字时代,大文件下载已成为日常操作,本文将手把手教你用Python打造专业级下载器,实现断点续传,多线程加速,速度限制等功能,感兴趣的小伙伴可以了解下... 目录一、智能续传:从崩溃边缘抢救进度二、多线程加速:榨干网络带宽三、速度控制:做网络的好邻居四、终端交互

java实现docker镜像上传到harbor仓库的方式

《java实现docker镜像上传到harbor仓库的方式》:本文主要介绍java实现docker镜像上传到harbor仓库的方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地... 目录1. 前 言2. 编写工具类2.1 引入依赖包2.2 使用当前服务器的docker环境推送镜像2.2

Go语言中nil判断的注意事项(最新推荐)

《Go语言中nil判断的注意事项(最新推荐)》本文给大家介绍Go语言中nil判断的注意事项,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录1.接口变量的特殊行为2.nil的合法类型3.nil值的实用行为4.自定义类型与nil5.反射判断nil6.函数返回的

C++20管道运算符的实现示例

《C++20管道运算符的实现示例》本文简要介绍C++20管道运算符的使用与实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录标准库的管道运算符使用自己实现类似的管道运算符我们不打算介绍太多,因为它实际属于c++20最为重要的

Java easyExcel实现导入多sheet的Excel

《JavaeasyExcel实现导入多sheet的Excel》这篇文章主要为大家详细介绍了如何使用JavaeasyExcel实现导入多sheet的Excel,文中的示例代码讲解详细,感兴趣的小伙伴可... 目录1.官网2.Excel样式3.代码1.官网easyExcel官网2.Excel样式3.代码

Go语言数据库编程GORM 的基本使用详解

《Go语言数据库编程GORM的基本使用详解》GORM是Go语言流行的ORM框架,封装database/sql,支持自动迁移、关联、事务等,提供CRUD、条件查询、钩子函数、日志等功能,简化数据库操作... 目录一、安装与初始化1. 安装 GORM 及数据库驱动2. 建立数据库连接二、定义模型结构体三、自动迁