给定 $k$ 个互不相同的复数 $x_0,\cdots,x_{k-1}$,以及 $k$ 个复数$y_0,\cdots,y_{k-1}$.我们知道存在唯一的复系数 $k-1$ 次多项式
$$
\mathcal{P}_{k-1}(x)=\xi_0+\xi_1x+\cdots+\xi_{k-1}x^{k-1}
$$
使得
$$
\mathcal{P}_{k-1}(x_0)=y_0,\cdots,\mathcal{P}_{k-1}(x_{k-1})=y_{k-1}.
$$
其中 $\xi_0,\cdots,\xi_{k-1}\in\mathbf{C}$.这个结论是范德蒙行列式不为0的一个简单推论.特别的,我们令 $x_i=\omega^i$,其中 $\omega=e^{\frac{2\pi i}{k}}$,我们就得到了三角插值多项式.为了确定三角插值多项式的系数,我们使用 Cramer 法则.我们知道
$$
\begin{cases}
\xi_0+\xi_1\omega^0+\cdots+\xi_{k-1}\omega^0=y_0\\
  \xi_0+\xi_1\omega^1+\cdots+\xi_{k-1}\omega^{k-1}=y_1\\
\xi_0+\xi_1\omega^2+\cdots+\xi_{k-1}\omega^{2(k-1)}=y_2\\
\vdots\\
\xi_0+\xi_1\omega^{k-1}+\cdots+\xi_{k-1}\omega^{(k-1)(k-1)}=y_{k-1}.
\end{cases}
$$
因此
$$
\xi_i=\frac{\begin{vmatrix}
    \omega^{0}&\omega^0&\cdots&y_{0}&\cdots&\omega^{0}\\
    \omega^0&\omega^1&\cdots&y_{1}&\cdots&\omega^{k-1}\\
    \vdots&\vdots&\cdots&\vdots&\cdots&\vdots\\
\omega^0&\omega^{k-1}&\cdots&y_{k-1}&\cdots&\omega^{(k-1)(k-1)}\\
  \end{vmatrix}}{\begin{vmatrix}
    \omega^{0}&\omega^0&\cdots&\omega^{0}\\
\omega^0&\omega^1&\cdots&\omega^{k-1}\\
\vdots&\vdots&\cdots&\vdots\\
\omega^0&\omega^{k-1}&\cdots&\omega^{(k-1)(k-1)}\\
  \end{vmatrix}}.
$$

$$
\begin{pmatrix}
  y_0\\
y_1\\
\vdots\\
y_{k-1}\\
\end{pmatrix}=\alpha_0 \begin{pmatrix}
  \omega^{0}\\
\omega^{0}\\
\vdots\\
\omega^0\\
\end{pmatrix}+\cdots+\alpha_i \begin{pmatrix}
  \omega^{0}\\
\omega^{i}\\
\vdots\\
\omega^{i(k-1)}\\  
\end{pmatrix}+\cdots+\alpha_{k-1}\begin{pmatrix}
  \omega^{0}\\
\omega^{k-1}\\
\vdots\\
\omega^{(k-1)(k-1)}
\end{pmatrix},
$$
则 $\xi_i=\alpha_i$.于是我们只用求 $\alpha_i$ 即可.易得
$$
\alpha_i=\frac{1}{k}\begin{pmatrix}
  y_0\\
y_1\\
\vdots\\
y_{k-1}\\
\end{pmatrix}\cdot \begin{pmatrix}
  \omega^0\\
\omega^{-i}\\
\vdots\\
\omega^{-i(k-1)}\\
\end{pmatrix}.
$$

转载于:https://www.cnblogs.com/yeluqing/p/3827412.html

三角插值的 Fourier 系数推导相关推荐

  1. 【六更完结!由于字数限制开新文章继续】零基础信号与系统学习笔记:复指数信号、傅里叶级数的系数推导、三角函数正交性、离散傅里叶变换、相位补偿、z变换表、逆变换表、常见序列及其作用

    零基础信号与系统学习笔记:复指数信号.傅里叶变换.三角函数正交性 基础1:复指数信号 复指数信号基础知识 复指数信号推导1 虚指数信号 虚指数信号特性和作用 直流信号 基础2:傅里叶级数 推导傅里叶级 ...

  2. 转贴: 傅立叶级数(Fourier Series) 推导

    转自: http://hi.baidu.com/hinus/blog/item/e9885e0eca9b30c97acbe173.html 傅立叶级数实在是我最痛恨的一种学问之一,来得突兀之至,一点兆 ...

  3. 图说Fourier变换

    傅里叶变换是傅里叶发现的,说起傅里叶这位大佬,他是一位研究热能的科学家,在研究热传导的时候发现了傅里叶变换,但傅里叶变换不仅只应用于热能,在声学,光学,电子学中也有广泛的应用. 傅里叶老先生的原话是: ...

  4. Part 13 (1) Fourier级数基本概念与应用

    文章目录 1. Fourier级数基本概念 1.1. 概念准备 1.1.1. 三角函数系的正交 1.1.2. 一种特殊的函数项级数--三角级数 1.1.3. 分段分析性质 1.2. Fourier系数 ...

  5. 【Fourier变换】傅里叶变换的性质与常用变换对(附注意事项)

    1. 定义 (1)Fourier正变换 F(ω)=F(f(t))=∫−∞+∞f(t)e−jωtdtF\left( \omega \right) =\mathscr{F} \left( f\left( ...

  6. 应用Maple软件推导直升机旋翼的气动力和力矩

    首先,看一下应用Maple推导的直升机旋翼的拉力.侧向力.后向力和反扭矩系数的总体效果: 应用Maple软件推导直升机旋翼的气动力和力矩 如果视频不能观看,请在浏览器地址栏直接复制视频地址:https ...

  7. matlab Fourier变换--方波信号的分解与重建

    N=200; %采用的数据点数 dt=4/N; %采样的间隔 for n=1:N %生成一个方波if (n*dt>=2)x(n)=0.8;elsex(n)=-0.8;end endfigure ...

  8. Fourier变换极其应用(Brad G. Osgood)——第1章——Fourier级数

    目录 第1章    Fourier级数 1.1 选择:"欢迎入局"(Choices: Welcome Aboard) 1.2 周期性现象(Periodic phenomena) 1 ...

  9. 理解神经网络函数高频成分的收敛率界限

    ©作者 | 王志伟.罗涛.许志钦 单位 | 上海交通大学 神经网络的频率原则 深度神经网络(DNN)在监督学习问题上展现出了其广泛的应用前景.近期的一系列的研究表明,神经网络的输出关于频率存在一种隐式 ...

最新文章

  1. no segments* file found in SimpleFSDirectory问题总结
  2. SpringBoot Scheduled Cron表达式范例记录
  3. COGS——T 803. [USACO Hol10] 政党 || 1776: [Usaco2010 Hol]cowpol 奶牛政坛
  4. laravel 安装
  5. MVC学习一:MVC简单流程
  6. 华锋e路航x10升级工具_E路航电子狗升级工具下载
  7. 深度linux操作系统
  8. 作业管理系统设计报告_Struts2_设计报告
  9. firefox浏览器一分钟去广告--去广告插件安装教程(adblock plus)
  10. IOS中impactor报81错误解决方法
  11. 无主之地3重型武器节奏评分榜(9.25) 枪械名 红字效果 元素属性 清图评分 Boss战评分 泛用性评分 特殊性评分 最终评级 掉落点 掉率 图片 瘟疫传播
  12. Windows中使用Java执行shell命令运行检测,通过sonarqube的webapi获取扫描结果
  13. app live photo_iphone live photo app恢复
  14. Polhemus PATRIOT 电磁位置跟踪系统
  15. PIL gray img /灰度图 转OpenCV格式灰度图
  16. MOSFET是什么器件,它的作用是什么
  17. DNS服务器未响应,电脑网页打不开
  18. Python大数据培训班特色优势及工作方向
  19. jzoj9382泡泡龙
  20. 小白易语言post培训接码登录day03

热门文章

  1. 插入顶部_声屏障顶部斜插式安装可获得10-25dB的降噪量
  2. 物流设计大赛优秀作品_独家!亲历全国大学生物流设计大赛,队员、老师、专家这样说...
  3. LetCode: 227. 简单计算器2
  4. python聚类分析如何确定分类个数_Python数据挖掘—聚类—KMeans划分法
  5. 获取crm服务器信息失败,无法连接到 Dynamics CRM 服务器,因为凭据身份验证 - Dynamics 365 Sales | Microsoft Docs...
  6. python求数列的积_python小游戏学习笔记4-2(列表【】,数列矩阵排列)
  7. mysql多实例安装配置演示
  8. 启动vm虚拟机里的系统时,提示此主机支持 AMD-V,但 AMD-V 处于禁用状态。
  9. 微信公众号点击菜单即可打开并登录微站的实现方法
  10. Boost.PropertyTree 属性树的简介