再论图像变化和频率的关系(使用数学工具)。

2024-04-22 13:44

本文主要是介绍再论图像变化和频率的关系(使用数学工具)。,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

我之前是做了一些探讨,但是没说清楚,现在再看这个问题。

我先提出这个问题。

以以为点列为例,先写成傅里叶级数的形式,不过这里不是三角函数形式,而是指数形式,是一样的。

对f(n)求导,就可以观察变化率了。但是我暂且不这样做,因为我先从直观感受出发。如果f(n+1)-f(n)较大说明了这个位置的像素变化快,那么在三角函数中该如何显示呢?把上图的指数函数看成是三角函数,所以差值f'(n)是跟频率有关,也跟三角级数的幅值F(k)有关。在连续函数的傅里叶级数求导中,如下图所示:

n跟频率有关,确实也说明了这一点。但是问题是在连续函数中,n是无数多个,而f'(x)的值是确定的,到底是多少个频率nk影响了f'(k)的值呢?连续的不好解决,现在看离散的级数的情况。

由于这里是有限个点,所以问题变简单了。

N个频率,只有一个频率k使得Fk*e^{i*2\pi *n*k/N}的绝对值最接近f'(n),然后其余N-1个F(k)作向量加法等于f'(n)。但是这样想无助于问题,那就利用方程组的思想吧。

已知有N个不同频率的正弦函数ck(n)和余弦函数sk(n)建立方程组

\sum_{k=0}^{N-1}Fk*[ck(n)+i*sk(n)]=fn

这样看不方便,用指数函数代替,设为en(k)=ck(n)+i*sk(n),即是

\sum_{k=0}^{N-1}Fk*en,k=fn,写成矩阵形式: E*(F0,F1,...,F(N-1))'=(f0,f1,...,f(N-1))'。

Fk的下标k表示频率,fn的下标n表示位置。

情况本身是这样的,首先是取了fn的N个点,然后由于N确定了矩阵E。所以可以求出来唯一解Fk。

但是问题是若f(m+1)-f(m)较大,则可能只是存在极少数比如l个的kl,跟这个差值非常接近。

E中的(n,k)元是第n行第k列元素,代表在复平面上x轴上的单位向量逆时转旋转2\pikn/N个角度,

可以看出来这是个对称矩阵。我不分析了,网上有傅里叶变换的矩阵分析,是范德蒙矩阵,还是个正交矩阵,也是对称矩阵。

设En表示矩阵的行向量, n是空间域的位置。

现在计算f(n+1)-f(n)=[E(n+1)-E(n)]*F

En看不出来是什么,但是矩阵E具有对称性,所以En(k)=E(k,n), 设Ek=E(k,n),实际上Ek代表的是在空间域上频率为k的一位置n为定义域的正弦函数和余弦函数对。

所以f(n+1)-f(n)=[E(n+1)-E(n)]*F表示: 当f(n)分解为三角函数的时候,变化率为两个相邻的正弦函数和余弦函数对的差值和F的内积。而三角函数早就已经固定了。

反正正弦函数与余弦函数都是在一个周期内,等分成了N个点。

f(n+1)-f(n)=[E(n+1)-E(n)]*E^(-1)*f', 由于E是对称矩阵,且是正定矩阵,则E^(-1)=E。

所以f(n+1)-f(n)=f(n+1)-f(n)。这只能说明我没有推到错误。但是k呢?我希望看到的是关于k的函数。由于对称性En已可以理解为在频率为n的时候,不同位置的三角函数值,Fk的k也可以理解为位置上的权重。所以也可以理解为不同频率上的差值的向量的内积。

所以到底是理解为在位置n和n+1处的两个相同频率的三角函数的差值,还是理解为在频率为n+1和n处的两个相同位置的三角函数的差值。为了不累加频率,那只能理解为Fk是关于位置的权重。但是依然没用。我需要的是在频率为n的时候,减少该频率的F(n)的值或者F(n)附近的值,就能改变位置n的或者附近的变化率。(利用数学工具都无法说明,那我只能从算法看了。)

这篇关于再论图像变化和频率的关系(使用数学工具)。的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



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

相关文章

Spring Boot整合Redis注解实现增删改查功能(Redis注解使用)

《SpringBoot整合Redis注解实现增删改查功能(Redis注解使用)》文章介绍了如何使用SpringBoot整合Redis注解实现增删改查功能,包括配置、实体类、Repository、Se... 目录配置Redis连接定义实体类创建Repository接口增删改查操作示例插入数据查询数据删除数据更

使用python生成固定格式序号的方法详解

《使用python生成固定格式序号的方法详解》这篇文章主要为大家详细介绍了如何使用python生成固定格式序号,文中的示例代码讲解详细,具有一定的借鉴价值,有需要的小伙伴可以参考一下... 目录生成结果验证完整生成代码扩展说明1. 保存到文本文件2. 转换为jsON格式3. 处理特殊序号格式(如带圈数字)4

Java使用Swing生成一个最大公约数计算器

《Java使用Swing生成一个最大公约数计算器》这篇文章主要为大家详细介绍了Java使用Swing生成一个最大公约数计算器的相关知识,文中的示例代码讲解详细,感兴趣的小伙伴可以了解一下... 目录第一步:利用欧几里得算法计算最大公约数欧几里得算法的证明情形 1:b=0情形 2:b>0完成相关代码第二步:加

Java中流式并行操作parallelStream的原理和使用方法

《Java中流式并行操作parallelStream的原理和使用方法》本文详细介绍了Java中的并行流(parallelStream)的原理、正确使用方法以及在实际业务中的应用案例,并指出在使用并行流... 目录Java中流式并行操作parallelStream0. 问题的产生1. 什么是parallelS

Linux join命令的使用及说明

《Linuxjoin命令的使用及说明》`join`命令用于在Linux中按字段将两个文件进行连接,类似于SQL的JOIN,它需要两个文件按用于匹配的字段排序,并且第一个文件的换行符必须是LF,`jo... 目录一. 基本语法二. 数据准备三. 指定文件的连接key四.-a输出指定文件的所有行五.-o指定输出

Linux jq命令的使用解读

《Linuxjq命令的使用解读》jq是一个强大的命令行工具,用于处理JSON数据,它可以用来查看、过滤、修改、格式化JSON数据,通过使用各种选项和过滤器,可以实现复杂的JSON处理任务... 目录一. 简介二. 选项2.1.2.2-c2.3-r2.4-R三. 字段提取3.1 普通字段3.2 数组字段四.

Linux kill正在执行的后台任务 kill进程组使用详解

《Linuxkill正在执行的后台任务kill进程组使用详解》文章介绍了两个脚本的功能和区别,以及执行这些脚本时遇到的进程管理问题,通过查看进程树、使用`kill`命令和`lsof`命令,分析了子... 目录零. 用到的命令一. 待执行的脚本二. 执行含子进程的脚本,并kill2.1 进程查看2.2 遇到的

详解SpringBoot+Ehcache使用示例

《详解SpringBoot+Ehcache使用示例》本文介绍了SpringBoot中配置Ehcache、自定义get/set方式,并实际使用缓存的过程,文中通过示例代码介绍的非常详细,对大家的学习或者... 目录摘要概念内存与磁盘持久化存储:配置灵活性:编码示例引入依赖:配置ehcache.XML文件:配置

Java 虚拟线程的创建与使用深度解析

《Java虚拟线程的创建与使用深度解析》虚拟线程是Java19中以预览特性形式引入,Java21起正式发布的轻量级线程,本文给大家介绍Java虚拟线程的创建与使用,感兴趣的朋友一起看看吧... 目录一、虚拟线程简介1.1 什么是虚拟线程?1.2 为什么需要虚拟线程?二、虚拟线程与平台线程对比代码对比示例:三

k8s按需创建PV和使用PVC详解

《k8s按需创建PV和使用PVC详解》Kubernetes中,PV和PVC用于管理持久存储,StorageClass实现动态PV分配,PVC声明存储需求并绑定PV,通过kubectl验证状态,注意回收... 目录1.按需创建 PV(使用 StorageClass)创建 StorageClass2.创建 PV