如何用Python 求函数 y = sinx 在区间[0, pi/2]上的弧长

最新推荐文章于 2022-11-25 20:51:21 发布
最新推荐文章于 2022-11-25 20:51:21 发布
阅读量3.7k 收藏 6
点赞数 1
分类专栏： 挑战Python 编程练习 文章标签： python 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Alex_mercer_boy/article/details/107577921
版权
本文讲述了如何使用Python求解函数y=sinx在区间[0, π/2]上的弧长近似值。通过将区间等分，利用微积分和弧长公式进行计算，发现当n足够大时，可以达到较高精度。作者还对比了使用勾股定理直接求小三角形斜边的方法，虽然精度相近，但速度较慢。强调了算法设计在提高效率中的重要性。" 107641014,9478665,Java Spring实现MySQL主从分离的代码配置,"['数据库', 'Java', 'Spring框架', 'MyBatis', '数据复制']
摘要由CSDN通过智能技术生成

今天在高数课上老师讲到求光滑曲线弧长问题，老师自己想了一个例子：

求函数 y = sinx 在区间[0, pi/2]上的弧长

但是经过微分和积分发现没法求出确定值，百度后发现不可积，是椭圆积分，只能求近似值。

于是我就心血来潮用Python 进行求近似

 

基本原理：（弧长公式）

L = \int_{a}^{b}\sqrt{1 + [f'(x)]^2} dx

基本步骤可以如下：

将区间分割为n份（可以等分）

取每一小区间 Δx （\Delta x = x_{i}-x_{i -1}）以及在x1处求得导函数y' = cox(x)

于是我们可以写出以下代码：

import math
times_n = [1, 10, 100, 1000, 10000, 100000, 1000
最低0.47元/天 解锁文章
Alex_mercer_boy
关注
  • 1
    点赞
  • 踩
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
计算椭圆弧长源码
03-18
椭圆弧长只能求取近似值，计算过程较为复杂。本人参照这篇论文进行了实现。
Python练习题——第六题：编写函数计算弧长的计算公式。弧长计算公式是一个数学公式，为L=n（圆心角度数）× π×2 r（半径）/360（角度制）。其中n是圆心角度数，r是半径，L是圆心角弧长。
阿彻的博客
12-26 2396
Python练习题第六题 题目： 编写函数计算弧长的计算公式。弧长计算公式是一个数学公式，为L=n（圆心角度数）× π×2 r（半径）/360（角度制）。其中n是圆心角度数，r是半径，L是圆心角弧长。 #2019/12/26/18:19 import math math.pi n = int(input("please enter the n:")) r = int(input("please ...
Python 根据弧度和半径求弧长
JiangTao2333的博客
01-20 1924
已知半径R和弧度A，求弧长l l = A * np.pi * R / 180 将直线斜率k转换成角度 karc = math.atan(k) kangle = math.degrees(karc) 根据两个顶点求线段斜率k和b def solveKB(x1, y1, x2, y2): k = (y2 - y1) / (x2 - x1) b = y1 - k * x1 return k, b 求直线的平行 y=ax+b+根号((d*a)^2+d^2) ..
python计算sinx在0-2π_定积分[0,2π]|sinx|
weixin_39661405的博客
12-10 5108
展开全部定积分[0,2π]|sinx|等于4。解：因为|sinx|≥0，而当0≤x≤πe5a48de588b63231313335323631343130323136353331333431363035时，sinx≥0，则|sinx|=sinx，而当π≤x≤2π时，sinx≤0，则|sinx|=-sinx。所以∫(2π，0)|sinx|dx=∫(π，0)sinxdx+∫(2π，π)(-sinx)d...
辛普森法则（Simposon rule）求参数曲线弧长 Python （NURBS Butterfly 为例）
weixin_42532587的博客
11-11 2482
NURBS 非均匀有理样条曲线 （由控制点Pi， 权重wi，节点向量u，阶数k定义） 辛普森法则 已知参数曲线某段弧长的两个端点[a, b]，求这段弧的弧长 在下面的python代码中c=(a+b)/2 Python 代码实现 from geomdl import linalg, NURBS, compatibility import numpy as np import math # 求得曲线一阶导 def first_derivative(a, crv): ders = crv.
正弦函数y=sinx图象与性质.ppt
09-16
为了画出正弦函数 y=sinx 的图象，可以使用“五点法”，即取五个关键点：(0, 0), (π/2, 1), (π, 0), (3π/2, -1), (2π, 0)，然后用光滑的曲线连接这些点。 正弦函数 y=sinx 的图象具有以下特征： * 周期为 2π...
深度学习入门-Pytorch框架前馈神经网络拟合y=sinx+exp(-x)
10-23
本项目将带你入门深度学习，利用PyTorch构建一个前馈神经网络（Feedforward Neural Network，FNN）来拟合一个特定的数学函数：y=sinx+exp(-x)。 首先，理解函数y=sinx+exp(-x)的性质很重要。这是一个复合函数，由...
求根方程（求给定的函数在给定区间内的所有零点。）
10-25
求根方程是数学中的一个重要概念，特别是在数值分析领域，其目标是找到使得函数值为零的点，也就是函数的零点。在给定的区间内寻找一个自定义函数的所有零点，通常需要借助计算机算法来实现，因为有些函数可能有多个...
15函数y=Asin（ωx+φ）的图象课件.ppt
10-23
- 例子3探讨了函数y=sin2x和y=sinx的图象，y=sin2x的图象是y=sinx图象上所有点的横坐标缩短到原来的1/2倍，即周期减半，这是因为ω从1变为2，引起了周期性的变化。 总之，函数y=Asin（ωx+φ）的图象可以通过振幅、...
函数y=Asin(ωχ＋φ)的图象习题精选.doc
08-19
在[0,π]上，y=sinx是奇函数且单调递增。 二、填空题解析： 13. 当正弦函数y=Asin(ωχ＋φ)的值域为[-1,3]时，A=2，因为正弦函数的正常值域是[-1,1]，这里扩大了两倍。 14. 对于函数y=3sin(2x)，当x=π/2时，sin(2...
由A，B经纬度计算A，B两点的距离和弧长_CGW
05-07
由A，B经纬度计算A，B两点的距离和弧长_CGW；简单的excel表批量处理经纬度
求弧长公式
m0_69576880的博客
10-22 326
针对于积分的定义，可以使用牛顿.莱布尼茨公式对其进行求解。那么求该函数在该区间上的弧的长度。
python计算sin程序_Python sin() 函数
weixin_39525313的博客
11-21 3566
描述sin() 返回的x弧度的正弦值。高佣联盟 www.cgewang.com语法以下是 sin() 方法的语法:importmathmath.sin(x)注意：sin()是不能直接访问的，需要导入 math 模块，然后通过 math 静态对象调用该方法。参数x -- 一个数值。返回值返回的x弧度的正弦值，数值在 -1 到 1 之间。实例以下展示了使用 sin() 方法的实例：#!/usr/bin...
求y=sinx在[0,π]上的反函数
weixin_46498268的博客
08-21 4155
为什么在[π/2，π]上y=sinx 的反函数不是x=arcsiny
Python数学基础二、利用正弦sin求曲边图形的面积
红目香薰
11-25 1086
这个小题还是有些麻烦的，理解起来还是没那么容易的。我们看第一个解决方案还是相对理解起来比较容易的，第二个推导式就不是太容易理解的。我们需要对Python的语法非常了解，再加上对题目的理解来搞定这个题目。古代的勾三股四弦五中说的弦就是我们要说的正弦，也就是直角三角形中的斜边，叫做弦，股就是人的大腿，古人称直角三角形长的那个直角边就叫做股。如果把勾股弦放在一个圆里面，弦就是圆上两个点的连线，最大的正弦值就是圆的直径。这个图也能明显的看出来，勾股弦之间的关系。正弦是∠α的对边/斜边的比。∠α的正弦=对边/斜边。
77---Python 计算Sin(x)的积分
u011699626的博客
12-22 2182
计算sin⁡x\sin{x}sinx在区间[0,π]\left[0, \pi\right][0,π]和区间[π,0]\left[\pi, 0\right][π,0]上的积分。 sin⁡x\sin{x}sinx在区间(0,π)\left(0, \pi\right)(0,π)上的积分可写为： ∫0πsin⁡x\int_0^{\pi}\sin{x}∫0π​sinx 令x=−xx=-xx=−x，可得： x∈[0,π]→−x∈[π,0]x\in\left[0, \pi\right]\rightarrow-x\in\l
python绘制正弦函数和余弦函数
m0_66279156的博客
06-08 8719
python绘制正弦函数和余弦函数图像
python坐标轴间隔_使用Python玩转高等数学(5)：三角函数
weixin_39716088的博客
01-08 867
从三角函数开始，我们使用matplotlib绘制函数图像，matplotlib是图形绘制库，使用matplotlib可以方便的绘制函数图形，以及直方图、条形图、散点图等统计图形。matplotlib的绘图命令和matlab的绘图命令基本相似，在使用上可以互相借鉴。安装matplotlib最简单的方法就是在操作系统的命令行窗口分别输入并执行下面的命令：python -m pip insta...
python微元法计算函数曲线长度
ouening的博客
09-14 6295
计算曲线长度，根据线积分公式： ∫BAf(x,y,z)dl∫ABf(x,y,z)dl\int_A^Bf(x,y,z)dl，令积分函数f(x,y,z)f(x,y,z)f(x,y,z)为1，即计算曲线的长度，将其微元化： limn→∞Δl→0∑i=0nf(xi,yi,zi)Δlilimn→∞Δl→0∑i=0nf(xi,yi,zi)Δli \underset{\begin{array}{c} n\...
求函数y=sinx/x在[0,1]上的积分
最新发布
06-09
这个积分是一个比较特殊的积分，它没有原函数，需要通过数值积分方法来求解。一种比较常用的数值积分方法是梯形法则，其基本思路是将积分区间划分成若干个小区间，然后在每个小区间上使用梯形面积来近似计算积分值。 具体来说，在区间[0,1]上均匀取n个点，将区间划分成n-1个小区间，每个小区间的宽度为h=1/n。然后，我们可以将这些小区间上的积分近似地表示为梯形面积之和，即： ∫[0,1] sinx/x dx ≈ h/2 * (y0 + yn + 2(y1+y2+...+yn-1)) 其中，y0 = sin(0)/0 = 1，yn = sin(1)/1 ≈ 0.8415，yi = sin(xi)/xi，xi = i*h。 将上式代入，得到： ∫[0,1] sinx/x dx ≈ h/2 * (1 + 0.8415 + 2(y1+y2+...+yn-1)) 因此，只需要计算出y1+y2+...+yn-1的值即可求得积分的近似值。这可以通过数值计算来实现，例如使用循环来累加每个小区间上的积分值即可。 代码实现如下： ``` #include <stdio.h> #include <math.h> int main() { int n = 100; // 区间划分数 double h = 1.0/n; // 小区间宽度 double sum = 0.0; // 积分值的累加器 // 计算每个小区间上的积分值并累加 for (int i = 1; i < n; i++) { double x = i*h; double y = sin(x)/x; sum += y; } // 计算积分的近似值 double integral = h/2 * (1 + sin(1)/1 + 2*sum); printf("∫[0,1] sinx/x dx ≈ %.6f\n", integral); return 0; } ``` 运行结果为：∫[0,1] sinx/x dx ≈ 0.946083。 注意：由于这是一个数值积分方法，所以计算结果只是一个近似值，并不是精确值。积分的精确值可能需要使用其他更加高级的数值积分方法来求解。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值