数值分析复习:逼近理论的应用——最小二乘问题、解超定、欠定方程组

本文主要是介绍数值分析复习:逼近理论的应用——最小二乘问题、解超定、欠定方程组,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

文章目录

  • 逼近理论的应用——最小二乘问题、解超定、欠定方程组
    • 离散平方逼近
    • 最小二乘解

本篇文章适合个人复习翻阅,不建议新手入门使用
本专栏:数值分析复习 的前置知识主要有:数学分析、高等代数、泛函分析

逼近理论的应用——最小二乘问题、解超定、欠定方程组

离散平方逼近

设全空间 X = R n X=\mathbb{R}^n X=Rn, 在 R n \mathbb{R}_n Rn 中取 m < n m<n m<n 个线性无关的向量 ( X 1 , … , X m ) (X_1,\dots,X_m) (X1,,Xm),令 M = s p a n { X 1 , … , X m } M=span\{X_1,\dots,X_m\} M=span{X1,,Xm},则对任意 Y ∈ X \ M Y\in X\backslash M YX\M M M M 中存在唯一的最佳逼近元 X ∗ = ∑ i = 1 m c i X i X^*=\sum\limits_{i=1}^mc_iX_i X=i=1mciXi,其满足以下法方程组
∑ i = 1 m < X i , X j > c i = < Y , X j > \sum\limits_{i=1}^m<X_i,X_j>c_i=<Y,X_j> i=1m<Xi,Xj>ci=<Y,Xj>若设 A = [ X 1 , … , X m ] , C = [ c 1 , … , c m ] T A=[X_1,\dots,X_m],C=[c_1,\dots,c_m]^T A=[X1,,Xm],C=[c1,,cm]T,则方程组等效于
A T A C = A T Y A^TAC=A^TY ATAC=ATY

最小二乘解

求如下的最小化问题的解
x ∈ R n , s . t . min ⁡ ∣ ∣ A x − b ∣ ∣ 2 x\in \mathbb{R}^n,s.t.\min||Ax-b||_2 xRn,s.t.min∣∣Axb2由离散平方逼近的理论,其解满足
A T A x = A T b A^TAx=A^Tb ATAx=ATb

应用:求解超定、欠定方程组

我们把线性方程组 A x = b Ax=b Ax=b 中,
未知数多于方程个数的方程组称为欠定方程组
未知数多于方程个数且有矛盾方程的方程组称为超定方程组

欠定方程组一般有多个解,超定方程组一般无解,故在工程上常用1范数或2范数意义下的最佳逼近解来作为解,即上述的最小二乘解
x ∈ R n , s . t . min ⁡ ∣ ∣ A x − b ∣ ∣ 2 x\in \mathbb{R}^n,s.t.\min||Ax-b||_2 xRn,s.t.min∣∣Axb2其解满足
A T A x = A T b A^TAx=A^Tb ATAx=ATb

参考书籍:《数值分析》李庆扬 王能超 易大义 编

这篇关于数值分析复习:逼近理论的应用——最小二乘问题、解超定、欠定方程组的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/860233

相关文章

MybatisGenerator文件生成不出对应文件的问题

《MybatisGenerator文件生成不出对应文件的问题》本文介绍了使用MybatisGenerator生成文件时遇到的问题及解决方法,主要步骤包括检查目标表是否存在、是否能连接到数据库、配置生成... 目录MyBATisGenerator 文件生成不出对应文件先在项目结构里引入“targetProje

C#使用HttpClient进行Post请求出现超时问题的解决及优化

《C#使用HttpClient进行Post请求出现超时问题的解决及优化》最近我的控制台程序发现有时候总是出现请求超时等问题,通常好几分钟最多只有3-4个请求,在使用apipost发现并发10个5分钟也... 目录优化结论单例HttpClient连接池耗尽和并发并发异步最终优化后优化结论我直接上优化结论吧,

Java内存泄漏问题的排查、优化与最佳实践

《Java内存泄漏问题的排查、优化与最佳实践》在Java开发中,内存泄漏是一个常见且令人头疼的问题,内存泄漏指的是程序在运行过程中,已经不再使用的对象没有被及时释放,从而导致内存占用不断增加,最终... 目录引言1. 什么是内存泄漏?常见的内存泄漏情况2. 如何排查 Java 中的内存泄漏?2.1 使用 J

numpy求解线性代数相关问题

《numpy求解线性代数相关问题》本文主要介绍了numpy求解线性代数相关问题,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 在numpy中有numpy.array类型和numpy.mat类型,前者是数组类型,后者是矩阵类型。数组

解决systemctl reload nginx重启Nginx服务报错:Job for nginx.service invalid问题

《解决systemctlreloadnginx重启Nginx服务报错:Jobfornginx.serviceinvalid问题》文章描述了通过`systemctlstatusnginx.se... 目录systemctl reload nginx重启Nginx服务报错:Job for nginx.javas

Redis主从/哨兵机制原理分析

《Redis主从/哨兵机制原理分析》本文介绍了Redis的主从复制和哨兵机制,主从复制实现了数据的热备份和负载均衡,而哨兵机制可以监控Redis集群,实现自动故障转移,哨兵机制通过监控、下线、选举和故... 目录一、主从复制1.1 什么是主从复制1.2 主从复制的作用1.3 主从复制原理1.3.1 全量复制

Redis缓存问题与缓存更新机制详解

《Redis缓存问题与缓存更新机制详解》本文主要介绍了缓存问题及其解决方案,包括缓存穿透、缓存击穿、缓存雪崩等问题的成因以及相应的预防和解决方法,同时,还详细探讨了缓存更新机制,包括不同情况下的缓存更... 目录一、缓存问题1.1 缓存穿透1.1.1 问题来源1.1.2 解决方案1.2 缓存击穿1.2.1

将Python应用部署到生产环境的小技巧分享

《将Python应用部署到生产环境的小技巧分享》文章主要讲述了在将Python应用程序部署到生产环境之前,需要进行的准备工作和最佳实践,包括心态调整、代码审查、测试覆盖率提升、配置文件优化、日志记录完... 目录部署前夜:从开发到生产的心理准备与检查清单环境搭建:打造稳固的应用运行平台自动化流水线:让部署像

vue解决子组件样式覆盖问题scoped deep

《vue解决子组件样式覆盖问题scopeddeep》文章主要介绍了在Vue项目中处理全局样式和局部样式的方法,包括使用scoped属性和深度选择器(/deep/)来覆盖子组件的样式,作者建议所有组件... 目录前言scoped分析deep分析使用总结所有组件必须加scoped父组件覆盖子组件使用deep前言

解决Cron定时任务中Pytest脚本无法发送邮件的问题

《解决Cron定时任务中Pytest脚本无法发送邮件的问题》文章探讨解决在Cron定时任务中运行Pytest脚本时邮件发送失败的问题,先优化环境变量,再检查Pytest邮件配置,接着配置文件确保SMT... 目录引言1. 环境变量优化:确保Cron任务可以正确执行解决方案:1.1. 创建一个脚本1.2. 修