本影演算小记

news/2024/9/19 16:33:46/文章来源:https://www.cnblogs.com/Rolling-star/p/18412639

\( \DeclareMathOperator{\i}{\mathbb{i}} \DeclareMathOperator{\e}{\mathbb{e}} \DeclareMathOperator{\eps}{\varepsilon} \DeclareMathOperator{\P}{\mathbb P} \DeclareMathOperator{\ord}{\operatorname{ord}} \)

一些前置知识
  • 域上的代数定义为:

    \(A\) 是域 \(\mathbb K\) 上的向量空间,再在 \(A\) 上定义乘法 \(\times\) 使得 \((A,+,\times)\) 是一个环,并且 \(\forall a\in\mathbb K,\mathbf {u,v}\in A\) 有:$$a(\mathbf u\times\mathbf v)=(a\mathbf u)\times\mathbf v=\mathbf u\times(a\mathbf v)$$则称 \(A\)\(\mathbb K\) 上的(结合)代数 .

  • 形式幂级数相关

    • 设形式幂级数 \(X=\sum_{k=0}^{\infty}a_kt^k\),则它的序(order)定义为:

      \[\ord (X)=\begin{cases}\infty&X=0\\\min\{n\mid a_n\neq 0\}&\text{else}\end{cases} \]

      请注意不要与多项式的度 \(\deg\) 相混淆 .

    • 如果形式幂级数 \(X\) 满足 \(\ord(X)=0\),则我们称 \(X\)可逆的(invertible) .

    • 如果形式幂级数 \(X\) 满足 \(\ord (X)=1\),则我们称 \(X\)delta 级数(delta series) .

1.本影代数

定义

\(\P\)\(\C\) 上关于单变量多项式的代数,设 \(\P^*\) 为所有 \(\P\) 上线性泛函的向量空间,我们用一个来自物理的记号:

\[\langle L\mid p(x)\rangle \]

来表示一个线性泛函 \(L\) 作用在多项式 \(p(x)\) 上,回过头来,\(\P^*\) 上的操作定义如下:

\[\langle L+M\mid p(x)\rangle=\langle L\mid p(x)\rangle+\langle M\mid p(x)\rangle \]

\[\langle cL\mid p(x)\rangle=c\langle L\mid p(x)\rangle \]

因为线性泛函由它在基上的值唯一确定(因为它本身可以表示成一个列向量),所以 \(L\) 由数列 \(\langle L\mid x^n\rangle\) 唯一确认 .


\(\mathcal F\)\(\C\) 上所有关于 \(t\) 的形式幂级数的集合,则形式幂级数:

\[f(t)=\sum_{k=0}^{\infty}\frac{a_k}{k!}t^k \]

通过定义:

\[\langle f(t)\mid x^n\rangle=a_n \]

定义了 \(\P\) 上的线性泛函,特别的,\(\langle t^k\mid x^n\rangle=[n=k]n!\),实际上,所有 \(\P^*\) 上的线性泛函 \(L\) 都有形式幂级数的表示,因为如果:

\[f_L(t)=\sum_{k=0}^{\infty}\frac{\langle L\mid x^k\rangle}{k!}t^k \]

则:

\[\langle f_L(t)\mid x^n\rangle=\langle L\mid x^n\rangle \]

所以 \(L=f_L(t)\),再加上一些简单的分析,可以得到 \(L\cong f_L\) .


仔细思考,发现对于 \(\mathcal F\) 中的每个元素,代表两个含义—— \(\P\) 上的线性泛函、与关于 \(t\) 的形式幂级数,因此,非常自然地定义了一个 \(\P^*\) 上的代数结构,我们称 \(\mathcal F\)本影代数(umbral algebra).

举个例子,对于 \(\mathcal F\) 上的形式幂级数 \(\e^{yt}\),有:

\[\langle \e^{yt}\mid p(x)\rangle=p(y) \]

因其性质,它被称作赋值泛函(evaluation functional).

性质

注意到 \(\forall f(t)\in \mathcal F\),有:

\[f(t)=\sum_{k=0}^{\infty}\frac{\langle f(t)\mid x^k\rangle}{k!}t^k\tag{1} \]

以及对于所有多项式:

\[p(x)=\sum_{k\ge 0}\frac{\langle t^k\mid p(x)\rangle}{k!}x^k\tag{2} \]

性质 1.1

\[\langle f(t)g(t)\mid x^n\rangle=\sum_{k=0}^n\binom{n}k \langle f(t)\mid x^k\rangle \langle g(t)\mid x^{n-k}\rangle \]

证明由 \((1)\) 展开即得 .

这定义了线性泛函的乘积,归纳一下可以得到多项式系数版本,在此不作赘述 .

性质 1.2 \(\quad\) 如果 \(\ord(f(t))>\deg p(x)\),则 \(\langle f(t)\mid p(x)\rangle=0\) .

证明显然 .

有两个意义不明的性质,先放在这里,证明略去:

性质 1.3 \(\quad\)

如果对于所有 \(k\ge 0\)\(\ord(f_k(t))=k\),且

\[\langle f_k(t)\mid p(x) \rangle=\langle f_k(t)\mid q(x) \rangle \]

\(p(x)=q(x)\) .

性质 1.4 \(\quad\)

如果对于所有 \(k\ge 0\)\(\deg p_k(x)=k\),且

\[\langle f(t)\mid p_k(x) \rangle=\langle g(t)\mid q_k(x) \rangle \]

\(f(t)=g(t)\) .

性质 1.5 \(\quad\)

如果 \(f(t)\in\mathcal F\),则

\[\langle f(t)\mid xp(x) \rangle=\langle f'(t)\mid p(x) \rangle \]

对任意多项式 \(p(x)\) 成立 .

证明:

由线性性,只需要考虑 \(p(x)=x^n\) 的情况,设

\[f(t)=\sum_{k=0}^{\infty}\frac{a_k}{k!} t^k \]

则:

\[\begin{aligned} \langle f'(t)\mid x^n \rangle &= \left\langle \sum_{k=0}^{\infty}\frac{a_{k+1}}{k!} t^k\mid x^{n} \right\rangle=a_{n+1}=\langle f(t)\mid x^{n+1}\rangle \end{aligned} \]

得证 .

性质 1.6 \(\quad\)

对于 \(f(t)\in\mathcal F\),与 \(p(x)\in\P\)

\[\langle f(at)\mid p(x) \rangle=\langle f(t)\mid p(ax) \rangle \]

对任意 \(a\in\C\) 成立 .

证明:

\(f(t)=t^n,p(x)=x^k\),我们有:

\[\langle (at)^{n}\mid x^k \rangle=a^nn![n=k]=a^kn![n=k]=\langle t^{n}\mid (ax)^k \rangle \]


下面给出一些 \(\mathcal F\) 上 functional 的例子 .

对于 delta series \(f(t)\in\mathcal F\),我们在下文中称其代表的线性泛函为 delta 泛函(delta functional),类似的,如果 \(f(t)\) 是 invertible series,我们称其代表的线性泛函为可逆泛函(invertible functional) .

  • 赋值泛函(evaluation functional) \(\e^{yt}\)

    前边提到过,其满足:

    \[\langle \e^{yt}\mid p(x)\rangle=p(y) \]

  • 前向差分泛函(forward difference functional) \(\e^{yt}-1\)

    \[\langle \e^{yt}-1\mid p(x)\rangle=p(y)-p(0) \]

  • Abel functional \(t\e^{yt}\)

    \[\langle t\e^{yt}\mid p(x)\rangle=p'(y) \]

  • \(\displaystyle\frac{1}{1-t}\)

    \[\left\langle \frac{1}{1-t}\mid p(x)\right\rangle=\int_0^{\infty} p(u)\e^{-u}\mathbb d u \]

  • \(\displaystyle\frac{\e^{yt}-1}{t}\)

    \[\left\langle \frac{\e^{yt}-1}{t}\mid p(x)\right\rangle=\int_0^{y} p(u)\mathbb d u \]

    顺带一提,其逆 \(\displaystyle\frac{t}{\e^{yt}-1}\) 和伯努利多项式有所关联,实际上,对于伯努利数 \(B_n\) 有:

    \[\left\langle \frac{t}{\e^{t}-1}\mid x^n\right\rangle=B_n \]

    到这里可能有人看出来了,\(\langle f(t)\mid x^n\rangle\Leftrightarrow [t^n]f(t)\),有人好像说过——“发明这玩意的估计是幂级数单推人,什么数列都想 cos 成幂函数。”,的确,从比较常用的本影演算,即用形式幂替换数列中就能看出来,但是,通过我们对 Umbral Algebra 的结构剖析,本影演算本身又是提取形式幂级数的工具,这下闭环了 .

  • \(\displaystyle\frac{1+\e^{yt}}{2}\)

    \[\left\langle \frac{1+\e^{yt}}{2} \mid p(x)\right\rangle=\frac{p(0)+p(y)}2 \]


有时间和精力的情况下我会尽量更新 .

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/799124.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

OpenDiary 24.9

《复活》很久没更日寄了这张是 9.16 下午出门的时候看见拍的 感觉很干净,很漂亮9.17今天突然开始更日寄总之写日记是一鸽再鸽,再加上博客园快完蛋了而新博客还没有眉目,暑假到开学期间写日寄的动力下降了很多 发生了很多事,实在是太多事了,以至于很难理清头绪还有就是,从…

基于FPGA的2ASK调制解调系统,包含testbench,高斯信道模块,误码率统计模块,可以设置不同SNR

1.算法仿真效果 vivado2019.2仿真结果如下(完整代码运行后无水印):本系统在以前写过的ASK调制解调系统的基础上,增加了高斯信道模块,误码率统计模块,可以验证不同SNR情况下的ASK误码情况。设置SNR=20db设置SNR=12db设置SNR=8db设置SNR=4db设置SNR=0dbRTL结构如下:2.算法…

基于心电信号时空特征的QRS波检测算法matlab仿真

1.课题概述通过提取ECG信号的时空特征,并使用QRS波检测算法提取ECG信号的峰值,并在峰值点标记峰值信息。2.系统仿真结果3.核心程序与模型 版本:MATLAB2022aclc; clear; close all; warning off; addpath(genpath(pwd)); rng(default)load ECG.mat Fs = 360; [loc,t…

基于PSO粒子群优化的GroupCNN分组卷积网络时间序列预测算法matlab仿真

1.算法运行效果图预览 (完整程序运行后无水印) 2.算法运行软件版本 matlab2022a3.部分核心程序 (完整版代码包含详细中文注释和操作步骤视频)LR = g1(1);numHiddenUnits1 = floor(g1(2))+1;% numHiddenUnits2 = floor(g1(3))+1;% layers = func_model2(Dim,num…

Qt MetaTypeInterface

Qt中的Meta Type 包含两种类型,每个类型都以一个整型ID来表示, 内建类型和用户自定义类型(自定义类型对应的整型ID 从65536 即 QMetaType::User 开始):1)内建类型: corelib\kernel\qmetatype.h// F is a tuple: (QMetaType::TypeName, QMetaType::TypeNameID, RealType) #d…

Marvelous Designer基础操作3 - 贴图导出

Marvelous Designer基础操作,包括加材质、贴图,导出衣服,导入模特做适配。--本篇导航--加材质 外部贴图 导出衣服、导入到Blender 导入自己的虚拟模特 给虚拟模特加安排板 自动适配板片配合视频看材质和导出的部分 视频中还包含了做双层的方法,也就是【克隆层(外部)】【克…

自动化测试 RobotFramework自定义静态测试类库总结

实践环境 win11 家庭中文版 Python 3.9.13 robotframework6.1.1 说明:为了方便的使用robot命令,安装好robotframwork后,修改系统环境,添加robot.exe(PYTHON_HOME/Scripts/robot.exe)所在路径到系统环境变量path 安装参考连接:https://github.com/robotframework/robotfr…

Marvelous Designer基础操作2 - 制板编辑

Marvelous Designer基础操作,本篇涉及到具体的板片编辑、合并、分割、扩展等等,也包含添加洞、做省、加纽扣、拉链、明线等。​--本篇导航--板片绘制、编辑 缝线 省 内部线、打洞、做褶皱 折叠安排、假缝 扣子 拉链 明线操作可以配合这个视频看一下,注意快捷键。基础操作 双…

图解软件开发的八大模型

图解软件开发的八大模型

[rCore学习笔记 027]地址空间

写在前面 本随笔是非常菜的菜鸡写的。如有问题请及时提出。 可以联系:1160712160@qq.com GitHhub:https://github.com/WindDevil (目前啥也没有 引言 兜兜转转又是新的一章的开始,还是首先要看官方手册里的理论介绍和内容. 这里主要还是提纲挈领地摘抄里面的部分内容,在下面…

如何查找 18 个月之前的 Apple 订单 All In One

如何查找 18 个月之前的 Apple 订单 All In One 为什么 Apple Store 不支持查找 18 个月之前的订单?如何查找 18 个月之前的 Apple 订单 All In One为什么 Apple Store 不支持查找 18 个月之前的订单?errorssolutionsemail history ✅demos(🐞 反爬虫测试!打击盗版⚠️)…

编程日记 更改redis存储默认序列化器

编程日记 更改redis存储默认序列化器 package com.haole.usercenter.service;import com.haole.usercenter.model.domain.User; import jakarta.annotation.Resource; import org.junit.jupiter.api.Assertions; import org.junit.jupiter.api.Test; import org.springframewor…