IT数码 购物 网址 头条 软件 日历 阅读 图书馆
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
图片批量下载器
↓批量下载图片,美女图库↓
图片自动播放器
↓图片自动播放器↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁
 
   -> 数据结构与算法 -> 数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法 -> 正文阅读

[数据结构与算法]数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法

正交多项式

上一节中我们学会了最佳一致逼近和平方逼近的方式求解函数逼近多项式,但是发现在求解法方程时,多元方程求解难度较大,如果我们能够将对应的法方法的矩阵直接变为只有对角线元素不为零的矩阵,将极大降低计算的运算量,因此我们引入了正交函数族的概念,从而得出正交多项式。下面简单介绍两种正交多项式勒让德多项式切比雪夫多项式,其他的正交多项式如:
L a g u e r r e 多 项 式 Laguerre多项式 Laguerre H e r m i t e 多 项 式 Hermite多项式 Hermite有兴趣的可以自行了解。

勒让德多项式(Legendre)

在这里插入图片描述
在这里插入图片描述

勒让德多项式的性质

在这里插入图片描述

在这里插入图片描述

其中最为重要的就是性质一到三和五,第四条基本上不咋用。然后最重要的表格来了,其实也可以用性质三来推出来:
在这里插入图片描述
整个题目来理解一下其中的妙处
在这里插入图片描述

性质

在这里插入图片描述

切比雪夫多项式(Shebyshev)

和上面学勒让德多项式一样来学习切比雪夫多项式;
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
同样的,咱搞个实战试一试
在这里插入图片描述

最小二乘法

定义;
在这里插入图片描述
求函数 S ? ( x ) = ∑ j = 0 n α j ? φ j ( x ) S^*(x)=\sum\limits_{j=0}^n\alpha_j^*\varphi_j(x) S?(x)=j=0n?αj??φj?(x)的问题可归结为求它的系数 α 0 ? , . . . , α n ? \alpha_0^*,...,\alpha_n^* α0??,...,αn??使得多元函数
F ( α 0 , α 1 , . . . , α n ) ? ∑ i = 1 m [ f ( x i ) ? ∑ j = 0 n α j φ j ( x j ) ] 2 \qquad\qquad\qquad\qquad\qquad\Large F(\alpha_0,\alpha_1,...,\alpha_n)\triangleq\sum\limits_{i=1}^m[f(x_i)-\sum\limits_{j=0}^n\alpha_j\varphi_j(x_j)]^2 F(α0?,α1?,...,αn?)?i=1m?[f(xi?)?j=0n?αj?φj?(xj?)]2取得最小值。由于 F ( α 0 , α 1 , . . . , α n ) F(\alpha_0,\alpha_1,...,\alpha_n) F(α0?,α1?,...,αn?)是关于 α 0 , α 1 , . . . , α n \large \alpha_0,\alpha_1,...,\alpha_n α0?,α1?,...,αn?的二次函数,所以可以使用多元函数取得极值的必要条件得:
? F ? α k = 0 , ( k = 0 , 1 , 2 , . . . , n ) \qquad\qquad\qquad\qquad\qquad\frac{\large\partial F}{\large\partial \alpha_k}=0 ,(k=0,1,2,...,n) ?αk??F?=0,(k=0,1,2,...,n)则通过
F ( α 0 , α 1 , . . . , α n ) ? ∑ i = 1 m [ f ( x i ) ? ∑ j = 0 n α j φ j ( x j ) ] 2 \qquad\qquad\qquad\qquad\large F(\alpha_0,\alpha_1,...,\alpha_n)\triangleq\sum\limits_{i=1}^m[f(x_i)-\sum\limits_{j=0}^n\alpha_j\varphi_j(x_j)]^2 F(α0?,α1?,...,αn?)?i=1m?[f(xi?)?j=0n?αj?φj?(xj?)]2有:
? F ? α k = ? 2 ∑ i = 1 n { [ f ( x i ) ? ∑ j = 0 n α j φ j ( x i ) ] φ k ( x i ) ] } = 0 , ( k = 0 ~ n ) \qquad\qquad\qquad\qquad\frac{\large\partial F}{\large\partial \alpha_k}=-2\sum\limits_{i=1}^n \lbrace [f(x_i)-\sum\limits_{j=0}^n\alpha_j\varphi_j(x_i)]\varphi_k(x_i)]\rbrace=0,(k=0 \sim n) ?αk??F?=?2i=1n?{[f(xi?)?j=0n?αj?φj?(xi?)]φk?(xi?)]}=0,(k=0n)
得到方程组:
∑ i = 1 m ( ∑ j = 0 n φ k ( x i ) φ j ( x i ) ) α j = ∑ i = 1 m f ( x i ) φ k ( x k ) \qquad\qquad\qquad\sum\limits_{i=1}^m(\sum\limits_{j=0}^n\varphi_k(x_i)\varphi_j(x_i))\alpha_j=\sum\limits_{i=1}^mf(x_i)\varphi_k(x_k) i=1m?(j=0n?φk?(xi?)φj?(xi?))αj?=i=1m?f(xi?)φk?(xk?)在离散情况下, f ( x ) , g ( x ) f(x),g(x) f(x),g(x)关于点 x 1 , x 2 , x 3 , . . . , x m x_1,x_2,x_3,...,x_m x1?,x2?,x3?,...,xm?的内积为 ( f , g ) = ∑ i = 1 m f ( x i ) g ( x i ) (f,g)=\sum\limits_{i=1}^mf(x_i)g(x_i) (f,g)=i=1m?f(xi?)g(xi?)也就是说可以将方程组改写为如下形式: ∑ j = 0 n ( φ k , φ j ) φ j = ( f , φ k ) , ( k = 0 ~ n ) \sum\limits_{j=0}^n(\varphi_k,\varphi_j)\varphi_j=(f,\varphi_k),(k=0\sim n) j=0n?(φk?,φj?)φj?=(f,φk?),(k=0n),再改写为矩阵形式:
在这里插入图片描述
也就得到了对应的线性拟合模型。

多项式拟合

在离散数据 ( x i , y i ) i = 1 ? m (x_i,y_i)_{i=1-m} (xi?,yi?)i=1?m?的最小二乘拟合中,最简单,最常用的数学模型就是多项式,取 φ 0 ( x ) = 1 , φ 1 ( x ) = x , φ 2 ( x ) = x 2 , . . . , φ n ( x ) = x n \varphi_0(x)=1,\varphi_1(x)=x,\varphi_2(x)=x^2,...,\varphi_n(x)=x^n φ0?(x)=1,φ1?(x)=x,φ2?(x)=x2,...,φn?(x)=xn得到的拟合函数 f ( x ) = c 0 + c 1 x + c 2 x 2 + . . . + c n x n f(x)=c_0+c_1x+c_2x^2+...+c_nx^n f(x)=c0?+c1?x+c2?x2+...+cn?xn
称为多项式拟合。
加权的最小二乘法一般是在通常拟合函数前加一个加权函数乘以拟合部分,具体的可以自己百度一下,嘻嘻。

最小二乘法实战

在这里插入图片描述
函数逼近我们将告一段落,后面我们将开启数值积分的版块,船新版本哦,等你来学!!!

欢乐的时光总是短暂的,让我们下一次再见!!!

good good study,day day up! (study hard, improve every day)

预知后事,请听下回分解!!!!

  数据结构与算法 最新文章
【力扣106】 从中序与后续遍历序列构造二叉
leetcode 322 零钱兑换
哈希的应用:海量数据处理
动态规划|最短Hamilton路径
华为机试_HJ41 称砝码【中等】【menset】【
【C与数据结构】——寒假提高每日练习Day1
基础算法——堆排序
2023王道数据结构线性表--单链表课后习题部
LeetCode 之 反转链表的一部分
【题解】lintcode必刷50题<有效的括号序列
上一篇文章      下一篇文章      查看所有文章
加:2021-12-23 15:58:58  更:2021-12-23 16:01:02 
 
开发: C++知识库 Java知识库 JavaScript Python PHP知识库 人工智能 区块链 大数据 移动开发 嵌入式 开发工具 数据结构与算法 开发测试 游戏开发 网络协议 系统运维
教程: HTML教程 CSS教程 JavaScript教程 Go语言教程 JQuery教程 VUE教程 VUE3教程 Bootstrap教程 SQL数据库教程 C语言教程 C++教程 Java教程 Python教程 Python3教程 C#教程
数码: 电脑 笔记本 显卡 显示器 固态硬盘 硬盘 耳机 手机 iphone vivo oppo 小米 华为 单反 装机 图拉丁

360图书馆 购物 三丰科技 阅读网 日历 万年历 2024年11日历 -2024/11/26 16:55:47-

图片自动播放器
↓图片自动播放器↓
TxT小说阅读器
↓语音阅读,小说下载,古典文学↓
一键清除垃圾
↓轻轻一点,清除系统垃圾↓
图片批量下载器
↓批量下载图片,美女图库↓
  网站联系: qq:121756557 email:121756557@qq.com  IT数码