数值分析笔记(三)函数逼近

news2024/10/7 20:27:32

最佳平方逼近

函数逼近是使用一种简单易算的函数来近似表示一个复杂函数。

在这里插入图片描述

该问题可转化为求解线性方程组

G n C = F n ​ G_{n}C=F_{n}​ GnC=Fn

其中,系数 C = ( c 0 , c 1 , ⋯   , c n ) T , F n = ( ( f , φ 0 ) , ( f , φ 1 ) , ⋯   , ( f , φ n ) ) T C=(c_{0},c_{1},\cdots,c_{n})^{\mathrm{T}},F_{n}=((f,\varphi_{0}),(f,\varphi_{1}),\cdots,(f,\varphi_{n}))^{\mathrm{T}} C=(c0,c1,,cn)T,Fn=((f,φ0),(f,φ1),,(f,φn))T

G n G_n Gn是格拉姆矩阵。称该线性方程组为法方程组或正规方程组。

最佳平方逼近的解函数为 φ ∗ = ∑ i = 0 n c i ∗ φ i \varphi^*=\sum_{i=0}^nc_i^*\varphi_i φ=i=0nciφi

最佳平方逼近函数,继承内积,即 ( φ ∗ , φ ∗ ) = ( φ ∗ , f ) (\varphi^*,\varphi^*)=(\varphi^*,f) (φ,φ)=(φ,f)

取逼近区间[a,b]为[0,1]时,其平方误差为:
∥ φ ∗ − f ∥ 2 2 = ( f , f ) − F n T C ∗ = ∫ 0 1 f 2 ( x ) d x − F n T C ∗ . \parallel\varphi^*-f\parallel_2^2=(f,f)-F_n^\mathrm{T}C^*=\int_0^1f^2(x) \mathrm{d}x-F_n^\mathrm{T} C^* . φf22=(f,f)FnTC=01f2(x)dxFnTC.
在这里插入图片描述

正交系

内积空间 V V V上的两个元素 f f f g g g,如果有内积 ( f , g ) = 0 (f,g)=0 (f,g)=0,则称 f f f g g g关于内积 ( ⋅ , ⋅ ) (\cdot,\cdot) (,)正交。若内积空间上的元素系 { f i } \{f_{i}\} {fi}满足两两正交
{ ( f i , f j ) = 0 ( i ≠ j ) , ( f i , f i ) = γ i > 0 , \begin{cases}(f_i,f_j)=0\quad(i\neq j) ,\\(f_i,f_i)=\gamma_i>0 ,\end{cases} {(fi,fj)=0(i=j),(fi,fi)=γi>0,
则称 { f i } \{f_{i}\} {fi}为正交系,若有 ( f i , f i ) = 1 ( i = 0 , 1 , 2 , 3... ) (f_i,f_i)=1(i=0,1,2,3...) (fi,fi)=1(i=0,1,2,3...),则称 { f i } \{f_{i}\} {fi}为标准正交系。

给定一组正交基,法方程组系数矩阵 G n G_n Gn为对角矩阵,其解向量为:
C ∗ = ( ( φ 0 , f ) ( φ 0 , φ 0 ) , ( φ 1 , f ) ( φ 1 , φ 1 ) , ⋯   , ( φ n , f ) ( φ n , φ n ) ) T . C^* = \left(\frac{(\varphi_0,f)}{(\varphi_0,\varphi_0)},\frac{(\varphi_1,f)}{(\varphi_1,\varphi_1)},\cdots,\frac{(\varphi_n,f)}{(\varphi_n,\varphi_n)}\right)^\mathrm{T}. C=((φ0,φ0)(φ0,f),(φ1,φ1)(φ1,f),,(φn,φn)(φn,f))T.
函数 f ( x ) f(x) f(x)的最佳平方逼近函数为
φ ∗ = ( φ 0 , f ) ( φ 0 , φ 0 ) φ 0 + ( φ 1 , f ) ( φ 1 , φ 1 ) φ 1 + ⋯ + ( φ n , f ) ( φ n , φ n ) φ n . \varphi^*=\frac{(\varphi_0,f)}{(\varphi_0,\varphi_0)}\varphi_0+\frac{(\varphi_1,f)}{(\varphi_1,\varphi_1)}\varphi_1+\cdots+\frac{(\varphi_n,f)}{(\varphi_n,\varphi_n)}\varphi_n. φ=(φ0,φ0)(φ0,f)φ0+(φ1,φ1)(φ1,f)φ1++(φn,φn)(φn,f)φn.
平方误差为
∥ f − φ ∗ ∥ 2 2 = ( f , f ) − ∑ i = 0 n ( f , φ i ) 2 ( φ i , φ i ) . \parallel f-\varphi^*\parallel_2^2=(f,f)-\sum_{i=0}^n\frac{(f,\varphi_i)^2}{(\varphi_i,\varphi_i)}. fφ22=(f,f)i=0n(φi,φi)(f,φi)2.

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

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

私域电商的新篇章:构建深度连接与高效生态

大家好,我是电商领域的探索者,今天我想和大家分享关于私域电商的一些心得与洞见。在这个数字化飞速发展的时代,如何构建与用户之间更为紧密、深入的连接,以及如何通过私域生态来挖掘用户的更大价值,成为了我们关注的焦…

Studio One 6.6.2中文破解版安装图文激活教程

Studio One 6.6.2中文破解版做为新生代音乐工作站,凭借更低的价格和完备的功能,获得了音乐人和直播行业工作者的青睐,尤其是对硬件声卡的适配支持更好,特别适合用来配合线上教学和电商带货。 最近网上出现不少关于StudioOne不能用…

springboot民宿信息管理系统-计算机毕业设计源码08818

摘 要 信息化社会内需要与之针对性的信息获取途径,但是途径的扩展基本上为人们所努力的方向,由于站在的角度存在偏差,人们经常能够获得不同类型信息,这也是技术最为难以攻克的课题。针对民宿信息管理系统等问题,对民宿…

Nuxt3 的生命周期和钩子函数(一)

title: Nuxt3 的生命周期和钩子函数(一) date: 2024/6/25 updated: 2024/6/25 author: cmdragon excerpt: 摘要:本文是关于Nuxt3的系列文章之一,主要探讨Nuxt3的生命周期和钩子函数,引导读者深入了解其在前端开发中…

Windows server 由于没有远程桌面授权服务器可以提供许可证,远程会话连接已断开。

问题现象: 解决办法 临时远程方式1: 打开 mstsc 时带上 /admin 等参数,如下图所示: 使用“mstsc /admin /v:目标ip”来强制登录服务器,但只能是管理员身份。 远程方式2: 通过VM远程登陆系统后,运行输入R…

如何将本地的Django项目部署到阿里云服务器上?

场景:在本地的pycharm上已经写好了一个Django架构的网站,现在要把它放到公网上 一、阿里云服务器 选择云服务器ECS,新用户可以免费使用三个月 购买时选择预装宝塔面板 买好后,进入云服务器控制台 重置实例密码 远程连接至服务…

【ACM出版】第13届亚洲膜计算会议(ACMC2024)暨 2024年机器学习、模式识别与自动化工程国际学术会议(MLPRAE 2024,8月7日-9)

第13届亚洲膜计算会议(ACMC2024)暨2024年机器学习、模式识别与自动化工程国际学术会议(MLPRAE 2024) 将于2024年8月7日-9日在新加坡举行。它致力于为机器学习、模式识别与自动化工程领域的专家和学者之间的学术交流创造一个平台。 会议的理念是让来自世…

JVM专题十:JVM中的垃圾回收机制

在JVM专题九:JVM分代知识点梳理中,我们主要介绍了JVM为什么采用分代算法,以及相关的概念,本篇我们将详细拆分各个算法。 垃圾回收的概念 垃圾回收(Garbage Collection,GC)确实是计算机编程中的…

视频录制软件哪个好用?5款简单好用软件推荐

在我们的日常生活中,都有哪些好用的视频录制软件?在很多场合中我们都会用电脑记录下重要的时刻。比如,在电脑上听老师讲解一道难题的方法时,怕自己会忘记,想要录制下来进行重复的观看。这时,选择一款好用的…

震惊!CURRENT_TIMESTAMP不能乱用

事情发生在签到和查询签到记录. 设置mysql时间默认值为CURRENT_TIMESTAMP可以随系统生成默认时间戳,即生成该数据的时间戳, 但是有些特殊场景要避免由mysql给我们生成默认时间: 1、首先签到成功之后返回给前端, 2、前端收到执行成功之后立马去查询签到记录, 3、发现并没有…

【软件工具】Xshell安装教程

1、安装软件:Xshell-5.0.1337p.exe,双击安装即可,可以选择安装到D盘; 2、在D盘安装完成后,将文件nslicense.dll拷贝到对应的安装目录下; 3、打开快捷方式即可打开应用软件。

“代码规范”这样做【高级前端必备软技能之一】

✨在前端开发领域,良好的代码规范是团队协作、项目可维护性的基石。最近在梳理我们团队关于代码规范相关的一些文档,顺便给大家分享一下我们是从哪些方面入手来从一个草台班子,到开发规范的高效团队。 ✨本文旨在提供一份详尽的前端代码规范指…

Linux的免交互

交互:我们发出指令控制程序的运行,程序在接收到指令之后按照指令的效果做出对应的反应。 免交互:间接的通过第三方的方式把指令传送给程序,不用直接的下达指令。 1、here document免交互 ere document免交互:是命令…

Studying-代码随想录训练营day20| 235.二叉搜索树的最近公共祖先、701.二叉搜索树中的插入操作、450.删除二叉搜索树中的节点

第二十天,二叉树part07,二叉树搜索树加油加油💪 目录 235.二叉搜索树的最近公共祖先 701.二叉搜索树中的插入操作 450.删除二叉搜索树中的节点 拓展:普通二叉树的删除方式 总结 235.二叉搜索树的最近公共祖先 文档讲解&…

在FlowUs息流,让知识库为你所用|如何打造个人知识库|如何打造企业知识库

🎉 在 FlowUs 的世界中,知识绽放出无限的可能!🚀 在当今信息爆炸的时代,知识的更新换代速度极快,我们每天都面临着海量的信息冲击。拥有一个属于自己的知识库变得至关重要。 首先,打造自己的知…

Linux_应用篇(27) CMake 入门与进阶

在前面章节内容中,我们编写了很多示例程序,但这些示例程序都只有一个.c 源文件,非常简单。 所以,编译这些示例代码其实都非常简单,直接使用 GCC 编译器编译即可,连 Makefile 都不需要。但是,在实…

解决IMX6ULL GPIO扩展板PWM7/8中的pwm0/period后卡死问题

前言 本篇文章主要是记录解决百问网论坛上面设置 IMX6ULL GPIO扩展板PWM7/8中的pwm0/period后卡死问题,如下图: 一、查看原理图,找出对应引脚 在这里我们如何确定哪个扩展口中的引脚输出PWM波呢?我们可以通过查看原理图。 查看…

【Java】解决Java报错:IllegalMonitorStateException in Synchronization

文章目录 引言一、IllegalMonitorStateException的定义与概述1. 什么是IllegalMonitorStateException?2. IllegalMonitorStateException的常见触发场景3. 示例代码 二、解决方案1. 确保在同步代码块或方法中调用wait()、notify()和notifyAll()2. 使用同步方法3. 使用…

办公人导航-上网导航,找网站,下软件,找资源!

办公人导航是一个专门为办公人员设计的实用导航网站,旨在帮助用户高效地找到各种优质的办公资源和工具。无论是需要查找办公软件、学习资源还是娱乐工具,在办公人导航上都能找到你需要的内容。 办公人导航-实用的办公生活导航网站!https://ww…

Linux内核测试技术

Linux 内核是Linux操作系统的核心部分,负责管理硬件资源和提供系统调用接口。随着 Linux 内核的不断发展和更新,其复杂性和代码规模也在不断增加。因此,确保内核的稳定性和可靠性变得尤为重要。内核测试技术是实现这一目标的关键手段。本文将…