趣味三角——第8章——高斯的一个求和法

Python 趣味百题 趣味整数 1 不重复的3位数 -易 2 水仙花数 -易 3 完全数 -中 4 相亲数 -中 5 黑洞数 -中 6 勾股数 -易 7 自守数 -易 8 3位反序数 -中 趣味素数 1 素数 阅读详情

第8章 高斯的一个求和法

(Variations on a Theme by Gauss)

The solving of an astronomical problem [proposed by

the French Academy of Sciences in 1735], for which

several eminent mathematicians had demanded several

months’ time. . . was solved by the illustrious Gauss in

one hour.

(解决一个天文问题[由法国科学院于1735年提出],几位著名的数学家为此花费了几个月的时间……被杰出的高斯在一小时内给解决了。)

——Florian Cajori引自R. E. Moritz的<<On Mathematics and Mathematicians>>(论数学和数学家)第155页

有一个关于伟大的德国数学家卡尔·弗里德里希·高斯(Carl Friedrich Gauss,1777-1855年)的故事,他在上学时被老师要求对1到100的数字求和,他几乎立即就给出了正确答案 5,050。 Gauss向惊讶的老师解释说,他只是注意到通过将数字写两次,首先是从1到100,然后是从100到1,并将两个数字垂直相加,每一对加起来就是101。因为有一百个这样的数字对,我们得到100×101 = 10,100,因为这是所需总和的两倍,答案是这个的一半,即 5,050。

就像许多关于名人的故事一样,这个故事可能真的发生过,也可能没有发生过;然而,重要的是要从中吸取教训——寻找模式的重要性。这种情况下的模式是楼梯模式,我们在一端添加的内容在另一端减去:

S = 1 +  2    +   3   + ... + n

S = n + (n-1) + (n-2) + ... + 1

------------------------------------------------------------------------

2 S = (n + 1) + (n + 1) + (n + 1) + ... + (n + 1) = n(n + 1)

---------------------------- n项----------------------------------------

S = n(n + 1)/2 。--------------------------------------------(1)

有一天,我在浏览一本序列和级数手册时想起了Gauss的故事,在那里我找到了以下求和公式:[1]

sin\alpha+sin2\alpha+sin3\alpha+...+sin(n\alpha)=\frac{sin\frac{n\alpha}{2}.sin\frac{(n+1)\alpha}{2}}{sin\frac{\alpha}{2}}---------(2)

起初对于如何证明这个公式没有任何线索,我开始寻找一个模式。令我印象深刻的是等式 (1)和(2)之间形式上的相似性;事实上,等式(1)的两边分别乘以α得到Sα = n(n + 1) α/2,即,

α + 2α + 3α + ... + nα = n(n + 1) α/2

将最后一个等式“乘以(Multiplying)”(两边取正弦)“sin”并继续化简,就好像“sin”是一个普通的代数量一样,我们得到

sin (α + 2α + 3α + ... + nα) = sin [n(n + 1) α/2]

如果我们在左侧打开括号,在右侧乘以和除以第二个“sin”(将其压缩在n和n + 1之间),然后再乘以α/2,我们得到等式(2)!

当然,我们犯了所有可以想象的数学“sin”错误(此词没有双关语),但我们确实得到了一个正确的公式。 那么,我们能否以类似于高斯对等式(1)中和的处理方式来证明等式(2)呢? 令

S = sin α + sin 2α + ... + sin(n-1)α + sin nα

S = sin nα + sin(n-1)α + ... + sin 2α + sin α

按上下对应成对的方式求和并利用和积公式 sin\alpha+sin\beta=2sin\frac{\alpha+\beta}{2}cos\frac{\alpha-\beta}{2},我们得到

2S = 2[sin\frac{(1+n)\alpha}{2}.cos\frac{(1-n)\alpha}{2}+sin\frac{(1+n)\alpha}{2}.cos\frac{(3-n)\alpha}{2}+...+sin\frac{(n+1)\alpha}{2}.cos\frac{(n-3)\alpha}{2}+sin\frac{(n+1)\alpha}{2}.cos\frac{(n-1)\alpha}{2}=2[sin\frac{(1+n)\alpha}{2}+cos\frac{(1-n)\alpha}{2}+cos\frac{(3-n)\alpha}{2}+...+cos\frac{(n-3)\alpha}{2}+cos\frac{(n-1)\alpha}{2}]

为了消除掉余弦项中烦人的1/2,3/2,...,让我们将最后一个等式乘以sin\frac{\alpha}{2}并利用积和公式(product-to-sum) sin\alpha.cos\beta = \frac{1}{2}[sin(\alpha-\beta)+sin(\alpha+\beta)];我们得到

2Ssin\frac{\alpha}{2} = sin\frac{(n+1)\alpha}{2}[ sin\frac{n\alpha}{2}+sin(1-\frac{n}{2})\alpha+sin(-1+\frac{n}{2})\alpha]+sin(2-\frac{n}{2})\alpha+...+sin(2-\frac{n}{2})\alpha+sin(-1+\frac{n}{2})\alpha+sin(1-\frac{n}{2})\alpha+sin\frac{n\alpha}{2}

但是,sin(-1+\frac{n}{2})\alpha=-sin(1-\frac{n}{2})\alpha, 以及类似的其它项;因此,括号中的表达式是“伸缩总和(telescopic sum)”,除了第一项和最后一项之外,其它项都抵消掉了(cancel out)。因此,我们有

2Ssin\frac{\alpha}{2}=2sin\frac{(n+1)\alpha}{2}sin\frac{n\alpha}{2}

或者

S=\frac{sin\frac{(n+1)\alpha}{2}sin\frac{n\alpha}{2}}{sin\frac{\alpha}{2}},

这就是我们想要证明的公式。

可以用类似的方式证明余弦求和的类似公式:[2]

cos \alpha + cos 2\alpha + cos 3\alpha + ...+ cos n\alpha = \frac{cos\frac{(n+1)\alpha}{2}sin\frac{n\alpha}{2}}{sin\frac{\alpha}{2}}--------(3)

如果我们现在将等式(2)除以等式(3),我们将得到漂亮的公式

tan\frac{(n+1)\alpha}{2}=\frac{sin\alpha+sin2\alpha+sin3\alpha+...+sinn\alpha}{cos\alpha+cos2\alpha+cos3\alpha+...+cosn\alpha}-----(4)

但这还没有结束。遵循每个三角函数公式最终都来自几何的想法,我们转向图51。从原点开始(基于我们的目的,我们标为P_{0})我们画了一条单位长度的线段P_{0}P_{1},使其与正向x轴构成角α。在P_{1}点我们画了第二条单位长的线段,使其与第一条线段构成角α,因此,与x轴构成的角为2α。持续n次这样的操作,我们来到点P_{n},我们将其坐标记为X和Y。很显然,X是n条线段的水平投影之和,而Y是n条线段的垂直投影之和,因此,我们有

X = cos α + cos 2α + cos 3α + ...+ cos nα

Y = sin α + sin 2α + sin 3α + ... + sin nα 。 ------------------------(5)

-------图51 S = sin α + sin 2α + ... + sin(n-1)α + sin nα 的几何构成-------------

现在,P_{i} 点位于一个规则的多边形上,且这个规则的多边形(polygon)内接于(inscribed)以O为圆心半径r为的圆中。每条线段P_{i-1}P_{i}正对向着一个以O点为顶点大小为α的角,因此,线段P_{0}P_{n}正对向着一个以O点为顶点大小为nα的角。但是这条线段是连接对角线(diagonal) P_{0} 和 P_{n};我们将其记为d。在这个等腰三角形(isosceles) P_{0}OP_{n}中,我们有

d=2rsin\frac{n\alpha}{2} ,

而在等腰三角形P_{0}OP_{1}中,我们有

1=2rsin\frac{\alpha}{2} 。

在这些等式中消掉r,我们得到

d = \frac{sin\frac{n\alpha}{2}}{sin\frac{\alpha}{2}}   。

为了求得线段P_{0}P_{n}的垂直和水平投影,我们必须求得它与x轴构成的这个角。这个角为α + β,其中\angle P_1P_0P_n。现在,角β正对向的弦P_1Pn内接于圆中,因此,它等于同一弦所正对向的圆心角的一半,即,等于(n-1)α/2。因此,α + β  = α + (n-1)α/2 = (n+1)α/2 。

因此

X=dcos[\frac{(n+1)\alpha}{2}]=\frac{cos\frac{(n+1)\alpha}{2}sin\frac{n\alpha}{2}}{sin\frac{\alpha}{2}}

和-----------------------------------------------------------(6)

Y=dsin[\frac{(n+1)\alpha}{2}]=\frac{sin\frac{(n+1)\alpha}{2}sin\frac{n\alpha}{2}}{sin\frac{\alpha}{2}}

如果我们将等式(5)中的 X 和 Y 的表达式代入等式(6),我们得到等式(2)和(3)。

假如我们将每条线段P_{i-1}P_{i}看成是从P_{i-1}到P_{i} 的向量, 则线段P_{0}P_{n}是它们的向量和。那么,等式(2)和(3)表示各个线段的(垂直或水平)投影之和等于它们的向量和的(垂直或水平)投影。这表明投影是一种“线性操作(linear operation)”——即遵循分配律(distributive law)p(u + v) = p(u) + p(v)的一种操作,其中,p表示“……的投影”, u 和 v表示任意两个向量——就像所有的线性操作——行为和普通乘法无异。我们可以用求和的“高斯法”来证明其他的三角求和公式。 这里有一些例子:

sin\alpha+sin3\alpha+sin5\alpha+...+sin (2n - 1)\alpha =\frac {sin^{2}n\alpha}{sin\alpha}------(7)

cos\alpha+cos3\alpha+cos5\alpha+...+cos (2n - 1)\alpha =\frac {sin2n\alpha}{2sin\alpha}------(8)

sin\frac{\pi }{n}+sin\frac{2\pi }{n}+...+sin\frac{n\pi }{n}=cot\frac{\pi }{2n}-------------(9)

cos\frac{\pi }{n}+cos\frac{2\pi }{n}+...+cos\frac{n\pi }{n}=-1-----------------(10)

cos\frac{\pi }{2n+1}+cos\frac{3\pi }{2n+1}+...+cos\frac{(2n-1)\pi }{2n+1}=\frac{1}{2}----------(11)

最后两个分别是等式(3)和(8)的特例;它们非常出色,因为每个中的总和与n无关。

匈牙利数学家Lip´ot Fej´er (1880–1959)结合他在傅立叶级数求和方面的工作对三角求和进行了研究,我们将在第 15 章中回到这个主题。

注释和资料来源:

1. <<Summation of Series>>(序列求和),由L.B.W. Jolley收集整理(1925初版; New York: Dover, 1961重印),序列号147。

2. 两个公式也可以通过取几何级数和 e^{ix}+e^{2ix}+...+e^{nix} 的实部和虚部来证明,其中i=\sqrt{-1};参见Richard Courant,<<Differential and Integral Calculus>>(微分与积分)( 1934初版;London: Blackie & Son, 1956重印),卷1,第436页。

内容来源:

<<Trigonometric Delights>> 作者:Eli Maor

高斯消元快速入门 高斯消元快速入门 一、基本描述 学习一个算法/技能,首先要知道它是干什么的,那么高斯消元是干啥的呢? 高斯消元主要用来求解线性方程组,也可以求解矩阵的秩,矩阵的逆。在ACM中是一个有力的数学武器. 它的时间复杂度是n^3,主要与方程组的个数,未知数的个数有关。 那么什么是线性方程组呢? 简而言之就是有多个未知数,并且每个未知数的次数均为一次,这样多个未知数组成的方程组为线性... 阅读详情

相关推荐

数字图像处理_06_(阈值分割——直方图阈值、三角法阈值,迭代法阈值,大津法阈值和自适应阈值)

本章将学习阈值分割,即找到一个合适的阈值实现数字图像的二值化。阈值分割有很多种方法,本章学习直方图阈值、三角法阈值,迭代法阈值,大律法阈值,和自适应阈值。 1.直方图阈值 图像直方图阈值分割法: 算法描述:根据图像的灰度直方图寻找阈值。 算法特点:适应于于直方图为双峰的图像。 绘制直方图代码: plt.hist(img.ravel(),256,[0,256]) 或者 plt.hist(img.flatten(),np.arange(-0.5,256,1),color='g') 绘制二值图: _,img_b

qq_40331078的博客 1万+

高次Gauss和估计简介

1953年采用错误的算术代数几何结果,还有一个基本技术错误,[N]的引理2的证明(因此[N]的定理1是无效的),因为在p.488的第一个不等式中,...。关于n阶多项式素数模完整三角和估计工作基础上得到任意整数模的估计,连 d(q)因子都无法改进为阶数n的函数,因此。[Lu]陆鸣皋,完整三角和的估计,《中国科学》A辑,1985年第1期pp.1-17。等人的工作全错,原因是最初企图改进的一篇论文,苏联人V.I.陆鸣皋,陈景润,戚鸣皋,丁平。

ATINER的博客 757

高斯求和法

高斯求和公式:和=(首项+末项)x项数÷2 例, 1+2+3+4+5+…+100 =(1+100)x100÷2 =101x100÷2 =101x50 =5050 这公式怎么来的? 1+2+3..+100 =(1+100)+(2+99)..(50+51) =101*50 =5050 图示如下; 下面来用C语言实现一下; void CGsqhView::OnDraw(CDC* pDC) { CGsqhDoc* pDoc = GetDocument(); ASSERT_VALID(p...

bcbobo21cn的专栏 1万+

歪写数学史(现代数学第一人——柯西)

柯西的这个名头有点大,虽然柯西的贡献毋庸置疑,第一人并不是指柯西对现代数学的贡献,而只是他出现的时间。如同从文言文到白话文,即使有新文化运动的推波助澜,改变也并非一夜之间,数学发展绝对是一个连续而非离散的过程,因此我们不可能用历史上的某一天划分数学。柯西与三L都有着或多或少的交集,他的一生的大部分时间与高斯重叠,他的数学成就毫无疑问是建立在众多前辈所打下的基础之上,但仅就今天数学严谨的形态而言,可

Ivan 的专栏 6278

【评价模型】模糊综合评价模型(模糊加权求和法、模糊层次分析法、模糊贝叶斯网络、模糊 TOPSIS 法)

模糊综合评价模型(Fuzzy Comprehensive Evaluation Model)是一种基于模糊数学理论的决策分析方法,用于评估和比较多个因素对目标的综合影响。它能够处理不确定性和模糊性的问题,适用于复杂的决策和评价任务。

Fanjufei的博客 2382

Gaussian Quadrature for Triangles:计算三角形的高斯节点和权重-matlab开发

此脚本计算高斯节点和权重,以对任意三角形域上的函数进行数值积分。 注意:此方法使用折叠的正方形而不是一般的体积情况。

三角形数

三角形数[编辑] 维基百科,自由的百科全书 本条目需要扩充。(2013年2月14日)  请协助改善这篇条目,更进一步的信息可能会在讨论页或扩充请求中找到。请在扩充条目后将此模板移除。 一定数目的点或圆在等距离的排列下可以形成一个等边三角形,这样的数被称为三角形数。比如10个点可以组成一个等边三角形,因此10是一个三角形数:

Domacles 3955

有限元高斯积分点MATLAB

在设计有限元分析程序时,要使用高斯积分代替传统积分。

高斯消去、列主元消去、Lu分解、追赶法(matlab)

基于高斯消去法的改进,当列主元出现数值较小时,在消去的过程列主元会作为分母,导致误差的出现。基于此,列主元消去是将每一次的主元通过行互换,用绝对值最大的列主元进行消去。去未知数的方法把原方程组 Ax = b 化为与其等价的三角形方程组,而求解三角形方程组可用回代的方法.。矩阵为三对角矩阵时 ,用a,b,c 存储。a为下对角线元素,b为对角线,c为上对角线。比如对与上面的这个方程组,用消去法解方程组的基本思想是用逐次消。追赶法只有当矩阵为三对角矩阵时才能使用(形如这样的)大体过程主要分为消去与回代。

weixin_53026957的博客 9207

多维高斯点和权重:给出 1-D 到 3-D 的高斯点和权重以及 6 节点三角形有限元-matlab开发

程序计算 1-D、2-D、3-D 的高斯点及其权重以用于数值积分。 最初是为有限元程序编写的,因此能够为 6 节点三角形单元提供积分点。 如果需要,可以以符号形式输出。 FEM 在http://www.colorado.edu/engineering/CAS/courses.d/IFEM.d/由科罗拉多大学博尔德分校的 Carlos Felippa 教授描述

01-SA8155P 冷启动强制进入EDL模式

冷启动EDL模式

静虚待令,有物混成 3095

高斯消去matlab,高斯消去、追赶法 matlab

1. 分别用Gauss消去法、列主元Gauss消去法、三角分解方法求解方程组程序:(1)Guess消去法:function x=GaussXQByOrder(A,b)%Gauss消去法N = size(A);n = N(1);x = zeros(n,1);for i=1:(n-1)for j=(i+1):nif(A(i,i)==0)disp(‘对角元不能为0‘);return;endm = A(j...

weixin_28737603的博客 642

在 STM32CubeIDE 中针对 STM32F407 移植 ThreadX

在 STM32CubeIDE 中针对 STM32F407 移植 ThreadX,不依赖盗版破解的Keil MDK、IAR等软件。

wangyuewu08的博客 3376

台湾国立大学郭彦甫Matlab教程笔记(21)linear equations(高斯消去法和追赶法)

本文记录了线性方程组的一些解法。 【1】高斯消去法rref()函数。 【2】追赶法:上下三角阵:Ax=b ,A =inv(L)*U 得到 y=U*x 先求解:inv(L)*y=b 再求解:Ux=y

阿正的梦工坊 2835

2025年最新Cesium学习示例:Cesium-Examples

2025年最新Cesium学习示例:Cesium-Examples

2403_88103571的博客 1011

centos7 python3.12.1 报错 No module named _ssl,坑我几个小时

网上的教程很多,但几乎是无效的,跟着我这教程,可以完美解决。

Amio_的博客 1万+

matlab高斯消去,高斯消去、追赶法 matlab

1. 分别用Gauss消去法、列主元Gauss消去法、三角分解方法求解方程组 程序:(1)Guess消去法:function x=GaussXQByOrder(A,b)%Gauss消去法N = size(A);n = N(1);x = zeros(n,1);for i=1:(n-1)for j=(i+1):nif(A(i,i)==0)disp('对角元不能为0');return;endm = A(...

weixin_42395213的博客 1134

程序员实战:DINOv3 权重文件下载与本地化部署全流程指南

Meta开源的DINOv3视觉大模型凭借70亿参数和17亿图像训练,在语义分割、目标检测等任务表现优异。本文详解预训练权重的获取渠道、校验方法和加载实践。开发者可通过HuggingFace或GitHub获取不同版本权重,包括针对服务器端的ViT系列和边缘设备的ConvNeXt系列。文章提供MD5校验代码确保文件完整性,并给出PyTorch和Transformers两种加载方式。特别强调7B模型需16GB GPU内存,建议根据硬件条件选择合适版本,如边缘设备使用ConvNeXt-Tiny。最后指出DINOv3

xiugou3365的博客 5911

高斯消元(Gaussian elimination)

我们将分两个部分来说明高斯消去,第一部分是数学方面,第二个部分是如何用 C++ 来实现。 第一部分 数学知识 什么是高斯消去 高斯消元法(或译:高斯消去法),是线性代数规划中的一个算法,可用来为线性方程组求解。 高斯消元的目标是将矩阵变为上三角矩阵,例如,对应的矩阵为就是一个上三角矩阵(upper triangular matrix)。 对于一个线性方程组,我们可以写出两个矩阵(系数矩阵 coefficient matrix 和增广矩阵 Augmented Matrix),如下图所示。 .

努力中的老周的专栏 4万+

新能源汽车电池及充电原理详解

本文详细介绍了新能源汽车的电池包结构、铅酸电池充放电原理,以及交流和直流充电的工作原理和优缺点。内容涵盖电池串联并联原因、充电设备如车载充电机、交流充电桩和直流充电桩的介绍,以及不同充电模式的连接方式。通过对充电过程的深入探讨,有助于读者理解新能源汽车的充电技术。

让你爱上电路设计 4528

SAP_MM模块非限制库存转冻结库存

项目非限制 → 冻结事务码移动类型343影响库存不可用,不参与MRP常见用途质量问题封存、盘点冻结、待退货品通过本文,我们了解了 SAP 的[SAP_MM模块非限制库存转冻结库存]的基本概念和应用方法。希望这些信息能帮助您在实际工作中更好地使用 SAP 系统,提高工作效率和管理水平。更多SAP文章请点击更多SAP文章更新,大家一起学习进步!!!文章中如有错误。敬请指正!搜索公众号:SAP资料文库微信号:SPRO_PP时间:2025-10-14。

我小时候很黑的博客 588

计算方法——三角分解法与高斯迭代法求解线性方程组

三角分解法求解线性方程组 【问题描述】为求解一个线性方程组,首先采用偏序选主元策略的三角分解法构造矩阵L,U和P,再用前向替换法对方程组LY=PB求解Y,最后用回代法对方程组UX=Y求解X。 【输入形式】在屏幕上依次输入方阵阶数n,系数矩阵A和常数矩阵B。 二、高斯赛德尔迭代法 【问题描述】为求解一个线性方程组,使用高斯赛德尔迭代法,采用欧几里得距离判定是否收敛。精度delta为1E-9,最大迭代次数为20。 【输入形式】在屏幕上依次输入方阵阶数n,系数矩阵A,常数矩阵B和起始点P。

啊噗呲咔的博客 934

什么是华为OD?怎么才能通过华为OD机试?华为OD机试怎么刷题?(2024最新华为OD机试E卷更新中)

华为OD是Outsourcing Dispacth模式,是华为和外企德科联合招聘的简称。目前华为大多数是OD招聘。OD模式也是华为提出的一种新的用工形式,每年都会从OD项目挑优秀员工转为正编。所有OD和正式员工一样对待,没有差别,部门平时开会、评审等活动都是一起参加的,是同薪同酬,都是由华为的HR直接定的薪资,看的是我们技术能力。华为是大企业,能到华为上班也是有能力的,每个部门都有好多职位,是金子,在哪都会发光。

KJ.JK 11万+

高斯消元法与LU分解

1.高斯消去法 有一般方程:Ax=b; 图1 ​​ 正常思路消元到最后,可以求出,回带之前的方程,可以解出. 图2 但是一旦数据规模大了,那么便很麻烦,很难求解。 ...

微凉的博客 1万+

数学建模模型总结

这是本人在两年内准备参加数学建模竞赛搜集到的模型资料,已经经过了本人的筛选和整理,希望能够对准备参加数学建模竞赛的学弟学妹们提供一些帮助!

win10 跳过U8环境检测的IE WebControls.rar

win10 安装u8的时候 提示IE WebControls 不通过,这个工具可以欺骗u8的检查。 跳过U8环境检测的IE WebControls。提示导入失败后,重新检测通过。

高斯-克吕格投影

由于这个投影是由德国数学家、物理学家、天文学家高斯于19 世纪20 年代拟定,后经德国大地测量学家克吕格于1912 年对投影公式加以补充,故称为高斯-克吕格投影。   即等角横切椭圆柱投影。假想用一个圆柱横切于地球椭球体的某一经线上,这条与圆柱面相切的经线,称中央经线。以中央经线为投影的对称轴,将东西各3°或1°30′的两条子午线所夹经差6°或3°的带状地区按数学法则、投影法则投影到圆柱面上

yusongwhu的专栏 3422

IEM模型,地表粗糙度模拟

利用微波反演土壤水分时候,利用IEM模型可以实现对随机地表的模拟

上一篇: 趣味三角——第7章——圆的内外旋轮线
下一篇: 趣味三角——第9章——Zeno错失无穷小
ComputerInBook
博客等级 码龄19年 1168粉丝 464原创
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付:10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值