自适应辛普森迭代求曲线的长度

已于 2023-06-14 15:43:09 修改
阅读量1.3k 收藏 6
点赞数 1
分类专栏： 算法 文章标签： 算法
于 2017-08-31 09:49:10 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/rui307/article/details/77737257
版权
本文探讨了如何使用自适应辛普森迭代来求解曲线的长度，主要针对贝塞尔曲线和样条曲线的弧长计算。通过递归方式计算参数曲线的弧长，并利用导矢和辛普森公式。在满足精度要求时，递归终止，确保算法效率。文章强调了算法实现和优化的重要性。
摘要由CSDN通过智能技术生成

自适应辛普森迭代求曲线的长度

最近工作都在研究几何相关的算法：需要计算贝塞尔，样条等曲线的弧长。贴一块最近写的递归计算长度的代码。

递归计算弧长的代码：

a,b=[0,1];  //
der1=ArcLenDer(0);//曲线导矢的模
der2= ArcLenDer(0.5);
der3=ArcLenDer(1);
//自适应辛普森计算公式
double simpsonCalLength(double SP,double EP
最低0.47元/天 解锁文章
超级替补
关注
专栏目录
MATLAB数值计算实战---2019年8月
qq_34965996的博客
01-23 2097
《MATLAB数值计算实战》—2019.8.14 第一章 MATLAB程序设计基础
机器人轨迹规划：求解样条曲线长度——自适应辛普森公式求积分
xuuyann
03-17 3723
写在前面 挖个坑~ 在样条学习过程中遇到了积分求解样条长度的问题，一般曲线长度的求解可以采用对速度积分的方式，但如果直接由速度解析表达式计算路径长度非常困难，尤其当表达式复杂的时候，对于计算机编程而言非常不友好。为此，可通过数值积分的方式对表达式进行近似积分，常见方法有复合梯形积分、复合辛普森积分。下面尝试采用自适应辛普森公式，先简单用matlab测试一下~ 自适应辛普森公式求积分 matlab简单测试代码如下 % 自适应辛普森求积分 % https://blog.csdn.net/frosero/art
Simpson积分方法计算NURBS曲线弧长，详细原理+代码实现
weixin_45650767的博客
04-13 1366
自适应 Simpson 积分方法是一种递归方法，将整个曲线分成若干小段，对每一小段采用 Simpson 公式进行数值积分，若相邻两次积分的结果相差较大，则将这一小段再分成两半，递归计算，直到精度满足要求为止。它的基本思想是将曲线分成若干小段，对每一小段采用 Simpson 公式进行数值积分，然后将各小段的积分结果加起来得到整个曲线的积分结果。函数中，采用递归方法计算曲线弧长，当相邻两次积分的结果相差较小时，直接将两段的积分结果加起来，否则将这一小段再分成两半，递归调用。在计算导数时，可以使用。
辛普森法则（Simposon rule）求参数曲线弧长 Python （NURBS Butterfly 为例）
weixin_42532587的博客
11-11 2584
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.
曲线积分求长度simpson公式
12-04 4115
曲线积分求长度 辛普森公式 基本思想就是把复杂的函数f近似成二次函数 简化得到 simpson公式对于在区间[a,b]越小，则求出的积分长度越接近实际长度。 但若[a,b]区间较大，则需要使用复化的辛普森公式。 复化辛普森公式 顾名思义是指将[a,b]区间划分为n等份，然后分别使用辛普森公式进行求和。 ...
hdu4498 求曲线长度 数值积分 simpson公式
weixin_30329623的博客
03-28 256
http://acm.hdu.edu.cn/showproblem.php?pid=4498 simpson公式可以求连续曲线的定积分。 这里注意解ax^2+bx+c=0时a=0,退化为一次的情况。 #include<iostream> #include<cstdio> #include<cstring> #include<cstdli...
LA 3485 辛普森公式求积分
退役ACMer
04-09 1222
传送门 : LA 3485题解 曲线积分知识（弧长公式）： 若可导曲线函数为f(x)其弧长公式为 L=∫ba1+f′(x)2−−−−−−−−√L = \int_a^b\sqrt{1 + f\prime(x)^2} 所以这题抛物线设为f(x) = a(x - d)(x +d) 弧长为2∗∫d01+4a2x2−−−−−−−−√ 2 * \int _0^d\sqrt{1 + 4a^2x^2}先预
数值积分和微分方程的解法
在数学中，积分是微积分的一个重要内容，用于求曲线下面积、求函数的定积分和求函数的不定积分等。而数值积分则是利用数值方法对积分进行近似计算的过程，常用于求解无法通过解析方法得到精确解的积分问题。数值积分...
Lecture.Matlab数值积分.pdf
02-03
辛普森规则是一种更为精确的数值积分方法，它通过拟合抛物线来近似积分区间内的曲线形状。辛普森规则通常比矩形公式或梯形公式更准确。在MATLAB中，可以编写相应的函数来实现辛普森规则，从而计算积分值。 ### 6. ...
Matlab学习
m0_75082186的博客
10-15 1103
2subs函数可以用指定符号替换符号表达式中的某一特定符号，调用格式：R=subs（s）%用工作空间中的变量值替代符号表达式s中的所有符号变量R=subs（S，New）%用新符号变量New替代原来符号表达式S中的默认变量R=subs（S,Old,New）%用符号变量New替代原来符号表达式S中的变量Old【9 7】替换函数subs应用示例syms a b t;subs(a^2+a*b+8,a,1) %将a^2+a*b+8中的a替换为1magic函数。
UVA - 1356 Bridge （辛普森公式）
lljjccsskk的专栏
11-03 597
Bridge Time Limit: 3000MS     64bit IO Format: %lld & %llu Submit Status uDebug Description 题意： 桥上等距离摆放着塔高为H的若干个塔，忽略宽度，相邻塔的距离不超过D，塔之间绳索形成全等的对称抛物线。桥长度为B，
辛普森法则
jackyvan的专栏
09-20 2万+
辛普森法则是求定积分近似值的一种方法，也叫抛物线法 它将区间[a,b] n 等分 ，n为偶数。设 ，过三点做一条抛物线，抛物线的函数表达式用 表示，此函数为偶函数。 则 根据 牛顿-莱布尼茨公式 可以推算 取 计算积分得到 最后得出 所以 $$ ...
simpson公式求定积分（模板）
weixin_33949359的博客
10-06 761
1 #include&lt;cstdio&gt; 2 #include&lt;cmath&gt; 3 #include &lt;algorithm&gt; 4 using namespace std; 5 6 double r,R1;//公式需要用到的变量 7 8 // simpson公式用到的函数，就是待积分函数 9 double F(double x)...
Simpson’s Rule （辛普森法则）
热门推荐
liangdaojun的博客
09-13 2万+
  In Simpson’s Rule, we will use parabolas to approximate each part of the curve. This proves to be very efficient since it’s generally more accurate than the other numerical methods such as straight ...
辛普森法求定积分
木质的旋律
10-19 7963
#include #include using namespace std; double f(double x); double getResult(double a, double b); int main() { double a = 0; double b = 2; double result = getResult(a, b); cout << "3 * x^2 在 " <<
数值分析-基本算法流程图：二分法、简单迭代法、复化梯形公式、弦截法、辛普森
Iris54Iris的博客
06-28 1万+
1.二分法2.复化Simpson公式3.复化梯形公式4.简单迭代法5.牛顿迭代法6.弦截法
黎曼积分求解可微曲线的弧线长度
Zhang Phil
12-31 5572
黎曼积分求解可微曲线的弧线长度假设曲线y=f(x)在区间[a,b]内光滑、可微且连续。那么可以根据微积分求解y=f(x)在a如图：从微分的思想入手建立数学函数式，假设s为曲线上(x,f(x))到(x+dx,f(x+dx))两点连线。这两点在水平方向的长度为dx，在垂直方向的y坐标轴长度为dy，根据直角三角形的勾股定理可知：其中，由f’(x)=dy/dx，得到dy =f’(x) dx从而：即ds的长
计算几何_求面积_辛普森积分公式
weixin_30319097的博客
08-04 900
附上一道水(模板)题 https://vjudge.net/problem/HDU-1724 和一道我认为卡精度要命的题(当然 是用这个方法被卡得要命的...) https://vjudge.net/problem/HDU-5858 /* %.17f的答案 附上 HDU-5858的答案 :0.29276251905757444 /*----3-----0.291392...
自适应辛普森求积分matlab
最新发布
11-11
以下是使用自适应辛普森公式求解积分的Matlab代码示例： ```matlab function [I, err] = adaptive_simpson(f, a, b, tol) % 自适应辛普森公式求解积分 % f: 被积函数 % a, b: 积分区间 % tol: 允许误差 % I: 积分...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值