算术几何平均matlab,算术-几何平均数——高斯的发现

“算术-几何平均数”既不是算术平均数,也不是几何平均数,由素有“数学王子”之称的德国数学家高斯首先发现和研究。算术-几何平均数,当然与“算术平均数”和“几何平均数”这两个概念有很深的关系。我们知道,但凡一个数学概念或定理,哪怕再简单不过,只要和高斯扯上关系,那就一定不简单了。带着耐心,我们来看看高斯关于算术-几何平均数的研究。

预备知识

对于两个正实数a和b(不妨设0

我们有基本不等式,

153131015_1_20190129122724477

等号当且仅当a=b时成立。

证明也不难:

从数的角度

153131015_2_20190129122724524

从形的角度

153131015_3_20190129122724602

一目了然。

正文

算术平均数和几何平均数的概念相当简单,绝大部分人认识到基本不等式这一步,可以说是功德圆满了。继续研究的话,无非两个方向:

第一,由两个数向三个、四个乃至任意n个正数的推广:

153131015_4_20190129122724696

第二,研究其他类型的平均,比如立方平均,平方平均,调和平均(倒数平均)以及它们之间的大小关系,得到更高级的基本不等式:

153131015_5_20190129122724758

也就是“立方平均数≥平方平均数≥算术平均数≥几何平均数≥调和平均数”。

上面的不等式同样可以推广到任意n个正数的情形。

绝大部分数学家走到这一步,也可以说是功德圆满了。

高斯,却另辟蹊径。

平均,平均,既然叫做”平均数“,自然介于两者之间,缓和了最大与最小。完整的基本不等式应该是:

153131015_6_20190129122724821

153131015_7_20190129122724883

由a和b,得到(a+b)/2和√ab,显然

153131015_8_20190129122724993

距离不到原来的一半。

令a1=√ab,b1=(a+b)/2,再计算它们的算术平均数和几何平均数,又有

153131015_9_2019012912272555

153131015_10_20190129122725133

同样地,它们之间的距离为

153131015_11_20190129122725243

这个过程可以无限进行下去,也就是

153131015_12_20190129122725305

那么数列{an}单调递增有上界,数列{bn}单调递减有下界,且当n趋于无穷时,

153131015_13_20190129122725368

于是数列{an}和{bn}有同一个极限。

高斯就把这个极限叫做a和b的算术-几何平均数(Arithmetic-Geometric Mean)。记为AGM(a,b)。

高斯当时只研究了算术-几何平均数。但顺着他的这个思路,我们当然还可以发明“算术-平方平均数”,“算术-调和平均数”,“平方-调和平均数”等概念。只需要在上面的迭代过程中,an和bn分别取an-1和bn-1不同的平均数即可。

这些平均数的数值都很容易计算,编个程序,迭代几次就能得到精度相当高的结果,收敛很快。

比如对1和2,小编用MATLAB编程,得到它们的算术-几何平均数约等于1.456791031046907,算术-平方平均数约等于1.540836469462489,平方-调和平均数约等于1.45458688740267。有兴趣的话可以试着计算其他组合的平均数。在计算的过程中,小编发现了一个很有意思的结论。限于篇幅,暂且不表。

本来两个数的平均,算数平均也好,几何平均也好,都很简单,计算简单,结果也简单。对1和2,它们的算术平均是1.5,几何平均是√2,平方平均是√(5/2),调和平均是4/3。然而对如此简单的1和2,它们的算术-几何平均数的卖相却如此“丑陋”!1.540836469462489.......看起来似乎还是个超越数!!!

高斯并不仅仅满足于数值运算。很快,他就找到算术-几何平均数AGM(a,b)的解析表达:

153131015_14_20190129122725430

圆周率π,三角函数,微积分......等等,算术-几何平均数怎么会和这些概念扯到一起???

153131015_15_20190129122725508

当年,高斯22岁。

后续

研究这些平均数,有什么用呢?

对我们来说,可以作为一种数学游戏,具有启发思维的作用。也许,可以应用在某个我们暂时还不知道的领域。

但高斯,他研究算术-几何平均数绝非一时的游戏之作。

作为一个“能从九霄云外的高度按照某种观点掌握星空和深奥数学的天才“,高斯发现,算术-几何平均数跟椭圆积分有很深的联系。

举个例子,有不少人对双纽线比较熟悉,双纽线是平面上到两个定点的距离之积为常数的动点轨迹(类比一下椭圆),长得像一个无穷符号。方程如下:

153131015_16_20190129122725571

学过高数的人应该知道,双纽线的面积是2a^2。但我们这里来看双纽线的周长。

为了简单起见,在上图中取a=1,它的极坐标方程是

153131015_17_20190129122725649

根据对称性,其周长

153131015_18_20190129122725712

利用高斯计算AGM(a,b)的公式,我们很容易得到该双纽线的周长

153131015_19_20190129122725821

为了纪念高斯,称

153131015_20_20190129122725883

为高斯常数(Gauss's Constant)。

双纽线的周长计算其实是一种椭圆积分,而椭圆积分的反演就是椭圆函数。椭圆函数可以说是19世纪的数学界在复变函数论方面取得的最为辉煌壮观的成果,没有之一。

人类历史上第一个被研究的椭圆函数,就是双纽线周长的积分反演。而研究它的,正是高斯。

椭圆函数在数论方面的应用发展出了模函数、模曲线、自守形式等理论。上世纪末,怀尔斯证明了费尔马大定理,应用的基本工具之一正是椭圆函数。

思考题

高斯22岁发现的定理

153131015_21_20190129122725946

有人对证明感兴趣吗?证明仅仅用到了高等数学的基础知识,没有任何知识盲点。如果感兴趣的话,私信或留言告诉我,分享你的思考证明过程。视情况,我将在下一篇贴出。

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

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

相关文章

把nc v6的源码看懂

看懂nc v6的源码! 碧桂园全部的正式环境的补丁都在我手里。 2015-11-18 2:33 谢谢两位一起努力的兄弟 谢谢超哥,谢谢祈冰哥,谢谢连老师,陈明大哥,谢谢龙哥,珍玉,谢谢廖生哥&#x…

.NET Core TLS 协议指定被我钻了空子~~~

【导读】此前,测试小伙伴通过工具扫描,平台TLS SSL协议支持TLS v1.1,这不安全,TLS SSL协议至少是v1.2以上才行,想到我们早已将其协议仅支持v1.3,那应该非我们平台问题。近日,第三方合作伙伴再次…

Android之CheckBox复选框控件使用inelayout.xml Xml代码

linelayout.xml Xml代码<?xml version"1.0" encoding"utf-8"?> <LinearLayout xmlns:android"http://schemas.android.com/apk/res/android" android:layout_width"match_parent" android:layout_height&quo…

linux如何挂载U盘

挂u盘之前&#xff0c;运行命令cat /proc/partitions,看看现在系统中有哪些分区。插上u盘以后&#xff0c;再次运行上述命令&#xff0c;看看多出来什么分区&#xff08;通常是sda1&#xff0c;PS由于偶是在虚拟机下装的&#xff0c;所以是sdb1). 1、插入U盘 2、输入 fdisk -…

光通信调制方式MATLAB仿真,基于LED的紫外光通信调制方式研究

赘固中文核心期刊基于&#xff2c;&#xff25;&#xff24;的紫外光通信调制方式研究何攀&#xff0c;李晓毅&#xff0c;侯倩&#xff0c;肖文林(重庆通信学院&#xff0c;重庆&#xff14;&#xff10;&#xff10;&#xff10;&#xff13;&#xff15;)摘要&#xff1a;紫…

frame buffer编程--画点功能和新增字符串代替RGBT

&#xff08;一&#xff09;&#xff1a;写在前面 这一次进行了一个比较大的迭代&#xff0c;在这里我先把相对于上一个版本的改进做一个解释&#xff1a; 1:进一步完善画点&#xff0c;增加了使用字符串代替RGBT值2:实现字符串解析出RGBT值3:增加了画圆的算法4:进一步增强了代…

这篇博士毕业论文致谢句句诛心......

全世界只有3.14 % 的人关注了爆炸吧知识本文来源&#xff1a;学术志综合自小木虫、网络等&#xff0c;作者yechm AUT博士毕业论文感谢信&#xff1a;我出生的时候是一个下雪的深夜&#xff0c;我的父亲在床边生了一堆柴火&#xff0c;我的母亲躺在四周漏风的那间破屋角落的床上…

30分钟通过Kong实现.NET网关

什么是KongOpenrestry是一个基于Nginx与Lua的高性能平台&#xff0c;内部有大量的Lua库。其中ngx_lua_moudule使开发人员能使用Lua脚本调用Nginx模块。Kong是一个Openrestry程序&#xff0c;而Openrestry运行在Nginx上&#xff0c;用Lua扩展了nginx。所以可以认为Kong Openres…

Proxy模式简介和用例

在软件系统中&#xff0c;有些对象有时候由于某些原因&#xff08;比如对象创建开销很大&#xff0c;或者某些操作需要安全控制&#xff0c;或者需要进程外的访问&#xff09;&#xff0c;如果直接访问会给使用者或者系统结构带来很多麻烦&#xff0c;这时可以在客户程序和目标…

struts2中jsp页面上验证码的生成

如何在jsp页面中显示验证码&#xff1a; 验证码&#xff1a;<img class"yzm_img" align"middle" id"validateImage" src"p_w_picpathCode.action" width"150" height"40" οnclick"this.srcp_w_picpathCo…

挑战程序设计竞赛部分题目选择排序、二分查找找出下标、深度优先算法、求三角形max周长、蚂蚁算法、操作excel

首先提供帮助类 可能这个类还需要更新 、求最大值、求最小值、选择排序、二分查找找出下标、二分查找排好序列里面是否有那个值 public class Util {//求最大值public static int maxValue(int a,int b) {int max =0;if (a > b) {max = a;} else {max = b;}return max;}//求…

nginx php oracle,第8天 ORACLE安装及NGINX整合PHP环境

oracle 11g安装在centos5.9x64下过程1.先检查是否安装开发工具、JAVA、ksh等&#xff0c;及其它如下yum search检查&#xff0c;若没有安装&#xff0c;则用yum groupinstall 开发工具 开发库 Java开发yum install pdkshyum install sysstatyum install libaio-develyum insta…

这个小姐姐真的很火辣......

1 穿着十分火辣的美女...▼2 你还有理了&#xff1f;▼3 学生&#xff1a;我们要把老师照片挂墙上&#xff01;▼4 喜悦中带着阵阵悲伤...▼5 你瞅啥&#xff1f;瞅你咋地&#xff01;▼6 一时竟没发现哪里不对▼7 老板&#xff0c;我们KTV最漂亮的都在这了▼你点的每个赞…

RPC调用框架比较分析

什么是RPC&#xff1a; RPC&#xff08;Remote Procedure Call Protocol&#xff09;——远程过程调用协议&#xff0c;它是一种通过网络从远程计算机程序上请求服务&#xff0c;而不需要了解底层网络技术的协议。 简言之&#xff0c;RPC使得程序能够像访问本地系统资源一样&am…

VS扩展CodeMaid代码整理插件

CodeMaid是Visual Studio的开源扩展&#xff0c;用于清理和简化我们的C&#xff03;&#xff0c;C &#xff0c;F&#xff03;&#xff0c;VB&#xff0c;PHP&#xff0c;PowerShell&#xff0c;JSON&#xff0c;XAML&#xff0c;XML&#xff0c;ASP&#xff0c;HTML&#xff0…

hibernate中的id特殊属性hilo剖解(多用于继承关系)

hilo是hibernate中最长用的一种生成方式&#xff0c;hibernate给出了hilo 和 seqhilo两种生成器&#xff0c;他们的分别向下面一样配置 Hilo配置代码 <id name"id" type"int" column"id"> <generator class"hilo"> <…

Centos5 install vnc

2019独角兽企业重金招聘Python工程师标准>>> 很详细的install steps: http://wiki.centos.org/HowTos/VNC-Server#head-76401321dae4d80916a7fd7e710272a9b85c9485 ---------------------------------------------------------------------------------------------…

搜索引擎(lucene)

专辑&#xff1a;Lucene (全文检索工具)从入门到精通简介:Lucene是apache软件基金会4 jakarta项目组的一个子项目&#xff0c;是一个开放源代码的全文检索引擎工具包&#xff0c;即它不是一个完整的全文检索引擎&#xff0c;而是一个全文检索引擎的架构&#xff0c;提供了完整的…

自动摘要php,修改DEDECMS文章自动摘要长度或者取掉文章摘要

首先&#xff0c;应该在系统→系统基本参数→其它选项“自动摘要长度” 中可以设置&#xff0c;但最高只能实现250&#xff0c;根本没用&#xff0c;所以现在总结一下几种方法方法一&#xff1a;第一步进数据库 找到表 dede_archives 下的字段 description 点更改,修改长度 或 …

解决PKIX:unable to find valid certification path to requested target 的问题

问题的根本是&#xff1a; 缺少安全证书时出现的异常。 解决问题方法&#xff1a; 将你要访问的webservice/url....的安全认证证书导入到客户端即可。 以下是获取安全证书的一种方法&#xff0c;通过以下程序获取安全证书&#xff1a; /** Copyright 2006 Sun Microsystems, In…