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

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

相关文章

Java使用Javassist动态生成HelloWorld类

《Java使用Javassist动态生成HelloWorld类》Javassist是一个非常强大的字节码操作和定义库,它允许开发者在运行时创建新的类或者修改现有的类,本文将简单介绍如何使用Javass... 目录1. Javassist简介2. 环境准备3. 动态生成HelloWorld类3.1 创建CtC

使用Python批量将.ncm格式的音频文件转换为.mp3格式的实战详解

《使用Python批量将.ncm格式的音频文件转换为.mp3格式的实战详解》本文详细介绍了如何使用Python通过ncmdump工具批量将.ncm音频转换为.mp3的步骤,包括安装、配置ffmpeg环... 目录1. 前言2. 安装 ncmdump3. 实现 .ncm 转 .mp34. 执行过程5. 执行结

Java使用jar命令配置服务器端口的完整指南

《Java使用jar命令配置服务器端口的完整指南》本文将详细介绍如何使用java-jar命令启动应用,并重点讲解如何配置服务器端口,同时提供一个实用的Web工具来简化这一过程,希望对大家有所帮助... 目录1. Java Jar文件简介1.1 什么是Jar文件1.2 创建可执行Jar文件2. 使用java

C#使用Spire.Doc for .NET实现HTML转Word的高效方案

《C#使用Spire.Docfor.NET实现HTML转Word的高效方案》在Web开发中,HTML内容的生成与处理是高频需求,然而,当用户需要将HTML页面或动态生成的HTML字符串转换为Wor... 目录引言一、html转Word的典型场景与挑战二、用 Spire.Doc 实现 HTML 转 Word1

Java中的抽象类与abstract 关键字使用详解

《Java中的抽象类与abstract关键字使用详解》:本文主要介绍Java中的抽象类与abstract关键字使用详解,本文通过实例代码给大家介绍的非常详细,感兴趣的朋友跟随小编一起看看吧... 目录一、抽象类的概念二、使用 abstract2.1 修饰类 => 抽象类2.2 修饰方法 => 抽象方法,没有

MyBatis ParameterHandler的具体使用

《MyBatisParameterHandler的具体使用》本文主要介绍了MyBatisParameterHandler的具体使用,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参... 目录一、概述二、源码1 关键属性2.setParameters3.TypeHandler1.TypeHa

Spring 中的切面与事务结合使用完整示例

《Spring中的切面与事务结合使用完整示例》本文给大家介绍Spring中的切面与事务结合使用完整示例,本文通过实例代码给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考... 目录 一、前置知识:Spring AOP 与 事务的关系 事务本质上就是一个“切面”二、核心组件三、完

使用docker搭建嵌入式Linux开发环境

《使用docker搭建嵌入式Linux开发环境》本文主要介绍了使用docker搭建嵌入式Linux开发环境,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面... 目录1、前言2、安装docker3、编写容器管理脚本4、创建容器1、前言在日常开发全志、rk等不同

使用Python实现Word文档的自动化对比方案

《使用Python实现Word文档的自动化对比方案》我们经常需要比较两个Word文档的版本差异,无论是合同修订、论文修改还是代码文档更新,人工比对不仅效率低下,还容易遗漏关键改动,下面通过一个实际案例... 目录引言一、使用python-docx库解析文档结构二、使用difflib进行差异比对三、高级对比方

sky-take-out项目中Redis的使用示例详解

《sky-take-out项目中Redis的使用示例详解》SpringCache是Spring的缓存抽象层,通过注解简化缓存管理,支持Redis等提供者,适用于方法结果缓存、更新和删除操作,但无法实现... 目录Spring Cache主要特性核心注解1.@Cacheable2.@CachePut3.@Ca