计算几何03_三次参数样条曲线与Cardinal曲线

最新推荐文章于 2024-05-13 09:44:47 发布
最新推荐文章于 2024-05-13 09:44:47 发布
阅读量4.4k 收藏 23
点赞数 3
分类专栏： 计算几何 文章标签： 计算几何
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_35312463/article/details/112388624
版权

在这里插入图片描述

一、三次参数样条曲线

三次样条曲线的唯一缺点就是缺乏几何不变形。即当型值点发生几何变换时不能保证参数递增。因此提出了以弦长为参数的三次参数样条曲线。

1.1 定义

已知n个型值点Pi(xi, yi), i = 1, 2,…, n且相邻型值点不重合；若 p(t) 满足下列条件：
（1）型值点Pi在函数 p(t) 上。
（2）p(t) 在整个区间[P1, Pn]上二阶连续可导。
（3）在每个子区间[Pi, Pi+1], i = 1, 2,…, n-1上，分段函数 p(t) 都是参数t的三次多项式。则称函数 是过型值点的三次参数样条函数，由三次参数样条函数构成的曲线称为三次参数样条曲线。

其中，
在这里插入图片描述
子区间的弦长Li为：
在这里插入图片描述

1.2 表达式

第i段的 pi(t) 表示为：
在这里插入图片描述

总的来说，三次参数样条曲线就是在三次样条曲线基础上对xyz各分量分别进行参数化的曲线表达，如下：

  • 3个坐标分量x，y，z分别是参数 t 的三次样条函数
  • 对型值点做参数化
  • 对3个坐标分量分别处理

1.3 算法

（1）读入n个型值点坐标。
（2）根据实际情况，确定三次参数样条曲线的边界条件。
（3）计算曲线的系数，将其表达为型值点二阶导数的函数。
（4）用追赶法分别沿x方向与y方向求解三弯矩方程。
（5）将求解出的系数代入三次参数样条函数的x分量与y分量表达式中，构造三次参数样条曲线。
（6）循环访问每个节点。在每个子区间内，参数t按照弦长计算增量，根据每组（x，y）坐标，绘制三次参数样条曲线。

1.4 代码实现

1.4.1 读取型值点

首先我们需要来扩充以下二维点类p2.h与p2.cpp(主要如下)：

class CP2  
{
public:
	CP2();
	virtual ~CP2();
	CP2(double x, double y);
	friend CP2 operator +(const CP2 &p0, const CP2 &p1);//运算符重载
	friend CP2 operator -(const CP2 &p0, const CP2 &p1);
	friend CP2 operator -(double scalar, const CP2 &p);
	friend CP2 operator -(const CP2 &p, double scalar);
	friend CP2 operator *(const CP2 &p, double scalar);
	friend CP2 operator *(double scalar, const CP2 &p);
	friend CP2 operator /(const CP2 &p0, CP2 &p1);
	friend CP2 operator /(const CP2 &p, double scalar);
public:
	double x;//直线段端点x坐标
	double y;//直线段端点y坐标
};

CP2::CP2()
{

}

CP2::~CP2()
{

}

CP2::CP2(double x, double y)
{
	this->x = x;
	this->y = y;
}

CP2 operator +(const CP2 &p0, const CP2 &p1)//对象的和
{
	CP2 result;
	result.x = p0.x + p1.x;
	result.y = p0.y + p1.y;
	return result;
}

CP2 operator -(const CP2 &p0, const CP2 &p1)//对象的差
{
	CP2 result;
	result.x = p0.x - p1.x;
	result.y = p0.y - p1.y;
	return result;
}
CP2 operator -(double scalar, const CP2 &p)//常量和对象的差
{
	CP2 result;
	result.x = scalar - p.x;
	result.y = scalar - p.y;
	return result;
}
CP2 operator -(const CP2 &p, double scalar)//对象和常量的差
{
	CP2 result;
	result.x = p.x - scalar;
	result.y = p.y - scalar;
	return result;
}

CP2 operator *(const CP2 &p, double scalar)//对象和常量的积
{
	return CP2(p.x * scalar, p.y * scalar);
}

CP2 operator *(double scalar, const CP2 &p)//常量和对象的积
{	
	return CP2(p.x * scalar, p.y * scalar);
}

CP2 operator /(const CP2 &p0, CP2 &p1)//对象的商
{
	if(fabs(p1.x)<1e-6)
		p1.x = 1.0;
	if(fabs(p1.y)<1e-6)
		p1.y = 1.0;
	CP2 result;
	result.x = p0.x / p1.x;
	result.y = p0.y / p1.y;
	return result;
}

CP2 operator /(const CP2 &p, double scalar)//对象数除
{
	if(fabs(scalar)<1e-6)
		scalar = 1.0;
	CP2 result;
	result.x = p.x / scalar;
	result.y = p.y / scalar;
	return result;
}

之后在主函数中定义读取函数：

// CGeometricfiguretestViewmessage handlers
void CGeometricfiguretestView::ReadPoint()
{
	P[1].x = -400, P[1].y = -150;
	P[2].x = -100, P[2].y =  0;
	P[3].x = -300, P[3].y =  100;
	P[4].x =  0,   P[4].y =  200;
	P[5].x =  300, P[5].y =  100;
	P[6].x =  100, P[6].y =  0;
	P[7].x =  400, P[7].y =  -150;
}

1.4.2 绘制型值点

//绘制型值点
void CGeometricfiguretestView::DrawDataPoint(CDC* pDC)
{
	CBrush NewBrush, *OldBrush;
	NewBrush.CreateSolidBrush(RGB(0, 0, 0));
	OldBrush = pDC->SelectObject(&NewBrush);
	for(int i = 1; i < 8; i++)
		pDC->Ellipse(ROUND(P[i].x - 5), ROUND(P[i].y - 5), ROUND(P[i].x + 5), ROUND(P[i].y + 5));
	pDC->SelectObject(OldBrush);
}

1.4.3 三次参数样条曲线绘制

三次参数样条曲线计算与绘制，具体过程见注释(计算与推导和三次样条曲线致)：

//三次参数样条曲线
void CGeometricfiguretestView::DrawCubicSpline(CDC* pDC)
{
	int n = 7;
	const int dim = 8;//二维数组维数
	double b1 = 0, bn = 0;//起点和终点的一阶导数
	double L[dim];//参数样条曲线的弦长
	double lambda[dim], mu[dim];
	double l[dim], m[dim], u[dim];
	CP2 c1, cn;//起点和终点的方向余弦
	CP2	D[dim];
    CP2 M[dim], K[dim];//追赶法过渡矩阵
    CP2 B1[dim], B2[dim], B3[dim], B4[dim];//函数的系数
	CP2 delt[dim];
	for(int i=1; i<n;i++)//计算弦长
    {
		delt[i]=P[i+1]-P[i];
		L[i]=sqrt(delt[i].x*delt[i].x+delt[i].y*delt[i].y);
	}
	//边界条件的投影
    c1.x = b1*cos(delt[1].x/L[1]);//起点
	c1.y = b1*cos(delt[1].y/L[1]);
   	cn.x = bn*cos(delt[n-1].x/L[n-1]);//终点
	cn.y = bn*cos(delt[n-1].y/L[n-1]);
	for(int i = 2; i < n; i++)
	{
		lambda[i] = L[i-1]/(L[i-1]+L[i]);//计算λ
		mu[i]=L[i]/(L[i-1]+L[i]);//计算μ
		D[i]=6/(L[i-1]+L[i])*((P[i+1]-P[i])/L[i]-(P[i]-P[i-1])/L[i-1]);//计算D
    }
	D[1]=6*((P[2]-P[1])/L[1]-b1)/L[1];//夹持端的D[1]
    D[n]=6*(bn-(P[n]-P[n-1])/L[n-1])/L[n-1];//夹持端的D[n]
    mu[1]=1;//夹持端的μ[1],
	lambda[n]=1;//夹持端的λ[n]
	//追赶法求解三弯矩方程
    l[1]=2;
	u[1]=mu[1]/l[1];
    for(int i=2; i <= n; i++)
    {
		m[i]=lambda[i];
        l[i]=2-m[i]*u[i-1];
        u[i]=mu[i]/l[i];
    }
    K[1] = D[1]/l[1];//解LK=D
    for(int i = 2; i <= n;i++)
    {
		K[i]=(D[i]-m[i]*K[i-1])/l[i];
	}
	M[n] = K[n];//解UM=K
	for(int i = n-1; i >= 1;i--)
	{
		M[i]=K[i]-u[i]*M[i+1];
	}
	//计算三次样条函数的系数
	for(int i = 1; i < n; i++)
    {
		B1[i]=P[i];
        B2[i]=(P[i+1]-P[i])/L[i]-L[i]*(M[i]/3+M[i+1]/6);
        B3[i]=M[i]/2;
        B4[i]=(M[i+1]-M[i])/(6*L[i]);
    }
	CPen pen(PS_SOLID,3,RGB(0,0,0));
	pDC->SelectObject(&pen);
	pDC->MoveTo(ROUND(P[1].x), ROUND(P[1].y));
	double tStep = 0.5;//步长
	CP2 p;
	for(int i = 1; i < n; i++)//循环访问每个节点
    {
		for(double t=0; t <= L[i]; t += tStep)
		{
			p =B1[i]+B2[i]*t+B3[i]*t*t+B4[i]*t*t*t;
			pDC->LineTo(ROUND(p.x), ROUND(p.y));//绘制参数样条曲线
        }
    }

1.4.4 主函数调用

void CGeometricfiguretestView::OnDraw(CDC* pDC)
{
	CTestDoc* pDoc = GetDocument();
	ASSERT_VALID(pDoc);
	if (!pDoc)
		return;
	// TODO: add draw code for native data here
	CRect rect;//定义客户区矩形
	GetClientRect(&rect);//获得客户区矩形的信息
	pDC->SetMapMode(MM_ANISOTROPIC);//自定义二维坐标系
	pDC->SetWindowExt(rect.Width(), rect.Height());//设置窗口范围
	pDC->SetViewportExt(rect.Width(), -rect.Height());//设置视区范围，且x轴水平向右为正，y轴垂直向上为正
	pDC->SetViewportOrg(rect.Width() / 2, rect.Height() / 2);//设置客户区中心为二维坐标系原点
	rect.OffsetRect(-rect.Width() / 2, -rect.Height() / 2);//rect矩形与客户区重合
	ReadPoint();
	DrawDataPoint(pDC);
	DrawCubicParameterSpline(pDC);
}

编译运行，可见如下：

在这里插入图片描述
我们也可以改变一下边界约束条件：

double b1 = 1, bn = 1;//起点和终点的一阶导数

在这里插入图片描述

二、Cardinal曲线

2.1 Hermite基矩阵

Hermite的每个曲线两个端点的坐标和导数来定义。
在这里插入图片描述
假设Hermite样条曲线的方程可以写成矩阵形式，如下：
在这里插入图片描述
将边界条件p(0)=yi和p(1)=yi+1代入上式，得：

yi ﹦d，yi+1﹦a﹢b﹢c﹢d

且上式的一阶导数为：

在这里插入图片描述
将边界条件p’(0)=Ri和p’(1)=Ri+1带入上式，得：

Ri ﹦c，Ri+1﹦3a﹢2b﹢c

Hermite边界条件的矩阵表示为：
在这里插入图片描述

该方程对多项式系数求解，有：
在这里插入图片描述
式中，Mh称为Hermite基矩阵，是边界约束矩阵的逆矩阵。需要说明的是，弗格森于1963年用于飞机设计的均匀参数插值三次样条方程就是Mh。
在这里插入图片描述
因此，可得：
在这里插入图片描述
最终可得到：
在这里插入图片描述

2.2 Cardinal曲线

Cardinal样条是Hermite样条的变形，仅由相邻型值点的坐标就可以计算出导数。

设相邻的四个型值点分别记为Pi-1, Pi, Pi+1,Pi+2。Cardinal样条插值方法规定Pi, Pi+1两型值点间插值多项式的边界条件为：
在这里插入图片描述

式中, u为可调参数,称为张力参数,可以控制Cardinal样条曲线型值点间的松紧程度。
记s=(1-u)/2,用类似Hermite曲线样条中的方法,将Cardinal边界条件代入下边的参数方程,
在这里插入图片描述

可以得到矩阵表达式：
在这里插入图片描述

式中，Mc称为Cardinal矩阵。

一条Cardinal样条曲线完全由四个连续的型值点给出，中间两个型值点是曲线段端点，另外两个型值点用来辅助计算端点导数。只要给出一组型值点的坐标值，就可以分段绘制出Cardinal样条曲线。

2.3 Cardinal 算法

（1）读入n个型值点坐标Pi，i=1,2，…n。
（2）在两端各加入一个虚拟点P0和Pn+1。
（3）设置张力参数u值，计算Cardinal基矩阵Mc。
（4）计算Mc矩阵与边界条件矩阵的算法，每4个顶点构成一个条件矩阵。
（5）根据型值点个数计算每段曲线中插值参数t的值，使用直线段连接对应的点p(x(t)，y(t))绘制曲线。

2.4 代码实现

1.4.1 读取型值点

void CGeometricfiguretestView::ReadPoint()
{
	P[0].x =-500, P[0].y =0;
	P[1].x =-450, P[1].y =-100;
	P[2].x =-350, P[2].y =-250;
	P[3].x =-250, P[3].y = 200;
	P[4].x =-150, P[4].y = 100;
	P[5].x =-50,  P[5].y = 150;
	P[6].x = 40, P[6].y = 10;
	P[7].x = 100, P[7].y = -60;
	P[8].x = 150, P[8].y = 80;
	P[9].x = 200, P[9].y = 70;
	P[10].x =230,P[10].y =-10;
	P[11].x =300,P[11].y =-200;
	P[12].x =400,P[12].y = 150;
	P[13].x =520,P[13].y = 250;
	n = 13;
}

1.4.2 绘制型值点

void CGeometricfiguretestView::DrawDataPoint(CDC* pDC)//绘制型值点
{
	CBrush NewBrush, *pOldBrush;
	NewBrush.CreateSolidBrush(RGB(0, 0, 0));
	pOldBrush = pDC->SelectObject(&NewBrush);
	for( int i = 1; i < 13; i++)
		pDC->Ellipse(ROUND(P[i].x - 5), ROUND(P[i].y - 5), ROUND(P[i].x + 5), ROUND(P[i].y + 5));
	pDC->SelectObject(pOldBrush);
	pDC->Ellipse(ROUND(P[0].x - 5), ROUND(P[0].y - 5), ROUND(P[0].x + 5), ROUND(P[0].y + 5));
	pDC->Ellipse(ROUND(P[13].x - 5), ROUND(P[13].y - 5), ROUND(P[13].x + 5), ROUND(P[13].y + 5));
}

1.4.3 计算与绘制

void CGeometricfiguretestView::DrawCardinalCurve(CDC* pDC)//Cardinal曲线
{
	double M[4][4];
	double u=0.2;
	double s=(1-u)/2;
	CP2 p;
	M[0][0]=-s ,M[0][1]=2-s,M[0][2]= s-2 ,M[0][3]= s;//Mc矩阵
	M[1][0]=2*s,M[1][1]=s-3,M[1][2]=3-2*s,M[1][3]=-s;
	M[2][0]=-s ,M[2][1]= 0 ,M[2][2]=  s  ,M[2][3]= 0;
	M[3][0]= 0 ,M[3][1]= 1 ,M[3][2]=  0  ,M[3][3]= 0;
	double t3,t2,t1,t0;
	CPen pen(PS_SOLID,2,RGB(0,0,0));
	CPen* pOldPen=pDC->SelectObject(&pen);
	pDC->MoveTo(ROUND(P[1].x),ROUND(P[1].y));
	for(int i=0; i< n-2; i++)
	{
		Point[0]=P[i],Point[1]=P[i+1],Point[2]=P[i+2],Point[3]=P[i+3];
		MultiplyMatrix(M, P, i);
		double tStep = 0.001;
		for(double t=0.0; t<1; t += tStep)
		{
			t3 = t*t*t; t2 = t*t; t1 = t; t0 = 1;
			p = t3*Point[0] + t2*Point[1] + t1*Point[2] + t0*Point[3];
			pDC->LineTo(ROUND(p.x),ROUND(p.y));
		}
	}
	pDC->SelectObject(pOldPen);
}

void CGeometricfiguretestView::MultiplyMatrix(double M0[][4], CP2 P0[4], int n)
{
	for(int i=0;i<4;i++)
		Point[i] = M0[i][0]*P0[n] + M0[i][1]*P0[n+1] + M0[i][2]*P0[n+2] + M0[i][3]*P0[n+3];
}

其中定义MultiplyMatrix函数为了实现4x4的Mc矩阵与1x4的边界约束矩阵的乘法计算。

编译运行，如下：
在这里插入图片描述

我们也可以改变一下张力参数来试一下：

double u=1.0;

在这里插入图片描述

沉默的舞台剧
关注
  • 3
    点赞
  • 踩
  • 23
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
专栏目录
B-spline曲线详解
03-09
给资料为B-spline曲线详细解释.B-spline基本特性， B-spline曲線數學模式，B-spline模型特性。
机器人轨迹规划：三次样条曲线
热门推荐
xuuyann
05-05 1万+
机器人轨迹规划：三次样条曲线 写在前面 在一些避障的应用场景下，一般都是先在任务空间中对多轴机械臂的末端进行路径规划，得到的是末端的运动路径点数据。这条轨迹只包含位置关系，并没有告诉机器人应该以怎样的速度、加速度运动，这就需要进行带时间参数的轨迹规划处理，也就是对这条空间轨迹进行速度、加速度约束，并且计算运动到每个路点的时间，高级的算法有TOPP等，一般的呢就是贝塞尔、三次准/非/均匀B、五次及三...
曲线生成 图解三次样条曲线生成原理(附ROS C++ Python Matlab仿真)
最新发布
HUAXIAL的博客
05-13 976
🔥附C++/Python/Matlab全套代码🔥课程设计、毕业设计、创新竞赛必备！详细介绍全局规划(图搜索、采样法、智能算法等)；局部规划(DWA、APF等)；曲线优化(贝塞尔曲线、B样条曲线等)。一、Python所有方向的学习路线Python所有方向的技术点做的整理，形成各个领域的知识点汇总，它的用处就在于，你可以按照下面的知识点去找对应的学习资源，保证自己学得较为全面。二、Python必备开发工具工具都帮大家整理好了，安装就可直接上手！三、最新Python学习笔记。
三次参数样条曲线与Cardinal曲线
danshiming的博客
10-16 551
已知n个型值点Pi(xi, yi), i = 1, 2,…, n且相邻型值点不重合；型值点Pi在函数 p(t) 上。p(t) 在整个区间[P1, Pn]上二阶连续可导。在每个子区间[Pi, Pi+1], i = 1, 2,…, n-1上，分段函数 p(t) 都是参数t的三次多项式。则称函数是过型值点的三次参数样条函数，由三次参数样条函数构成的曲线称为三次参数样条曲线。其中，
三次参数样条曲线拟合（任意控制点）
jiangjjp2812的专栏
08-29 1万+
三次参数样条曲线拟合，主要是为解决三次样条函数不能解决的问题而提出的。三次样条函数要求x满足单调递增，即x序列满足x0<x1<x2<...<xn。物理上的意义是，曲线不可以出现绕回或打圈。三次参数样条采用参数方程表示曲线，较为方便解决此问题。 三次样条函数原理百度里比较多，这里就不讲了。三次参数样条实现的原理是在三次样条函数的基础上进行参数化求解的...
[计算机动画] 路径曲线与运动物体控制（Cardinal样条曲线）
(ง •̀_•́)ง
10-05 8330
•路径曲线与运动物体控制 •掌握Cardinal样条曲线的表示和算法，了解控制参数对曲线形状的影响。 •对照Cardinal样条曲线的数学表示和程序代码的对应关系。 •在路径曲线上放置一小汽车，使其在路径上运动起来，汽车运动速度可调。
三次样条曲线CubicSpline
weixin_62705892的博客
05-04 2878
三次样条曲线CubicSpline
绘制三次参数样条插值曲线 (3KB).ZIP_三次 参数 样条 曲线_三次样条_三次样条插值_插值_样条插值
09-19
三次参数样条插值曲线是一种在数学和计算机图形学中广泛应用的技术，用于平滑地连接一系列离散数据点，形成连续且光滑的曲线。这种插值方法特别适用于数据拟合和曲线绘制，尤其是在工程、科学计算以及动画制作等领域...
statistics.zip_三次样条拟合_曲线拟合_样条曲线
07-14
用此程序，可将一条曲线离散成散点，并用三次样条曲线重新拟合，使其具有良好而形态。
三次样条曲线.rar_三次样条_三次样条代码_三次样条曲线_曲线平滑_样条曲线
07-14
3种源码，三次样条函数算法C++源代码,用于工程上计算平滑曲线的插值计算.
huatu.zip_B样条曲线_三次B样条曲线_三次bezier曲线_三次样条_样条曲线
09-14
本篇将深入探讨"三次Bezier样条曲线"和"三次B样条曲线"，它们是样条曲线家族中的两个关键成员。 首先，我们来理解什么是Bezier曲线。Bezier曲线是由法国工程师Pierre Bezier在1962年提出的，它是一种参数化的曲线，...
cardinal曲线JS算法和demo
07-12
cardinal曲线JS算法和demo演示。 JavaScript版的。 Java版的请看 http://blog.csdn.net/aofavx/article/details/75012648
Cardinal曲线（计算机图形学作业）
12-07
c++ QT实现可交互的Cardinal曲线（计算机图形学作业），鼠标点击屏幕自动生成点，点击绘制自动生成cardinal曲线，拖动滑动条修改张力参数t，曲线随之变化，同时点也可以用鼠标点击拖动
三次样条法求函数的导数
10-19
建模基础算法，数值微分，用三次样条法求函数的导数。
7种插值算法的c++代码实现
09-13
7种插值算法的c++代码实现，1 拉格朗日插值(POLINT) 2 有理函数插值(RATINT) 3 三次样条插值(SPLINE(二阶导数值)->SPLINT(函数值)) 4 有序表的检索法(LOCATE(二分法), HUNT(关联法)) 5 插值多项式(POLCOE(n2), POLCOF(n3)) 6 二元拉格朗日插值(POLIN2) 7 双三次样条插值(SPLIE2)
Cardinal样条曲线Qt代码
10-15
掌握Cardinal样条曲线的表示和算法，了解控制参数对曲线形状的影响。 对照Cardinal样条曲线的数学表示和程序代码的对应关系。 在路径曲线上放置一小汽车，使其在路径上运动起来，汽车运动速度可调。
ThreeCurve.rar_三次样条 曲线_曲线 插值_样条平滑_样条插值
07-14
用MFC做的可以实现画三种平滑曲线的程序(三次样条插值，贝塞尔曲线，GDI+曲线平滑)
MFC绘制三次Cardinal曲线
ha0ha0xuexi的博客
05-21 231
MFC绘制三次Cardinal曲线 参考《计算几何算法与实现》–孔令德 给定12个型值点，绘制Cardinal曲线： 1、新建MFC单文档程序 2、添加二维点类： #pragma once //为了避免按照x和y方向进行重复运算，重载运算对象 class CP2 { public: CP2(void); ~CP2(void); CP2(double x,double y); friend CP2 operator+(const CP2&p0,const CP2&p1);//运算符
计算几何02_三次样条曲线
qq_35312463的博客
01-07 8879
一、样条 样条（Spline）函数是由舍恩伯格于1946年提出的。样条是富有弹性的细木条或有机玻璃条，它的作用相当于“万能”曲线板。早期船舶、汽车、飞机放样时用铅压铁压住样条，使其通过一系列型值点，调整压铁达到设计要求后绘制其曲线，称为样条曲线。这样设计曲线的方法在20世纪六七十年代得到了广泛应用。 二、几何连续性 2.1 连续性条件 通常单一的曲线段或曲面片难以表达复杂的形状，必须将一些曲线段拼接成组合曲线，或将一些曲面片拼接成组合曲面，才能表达复杂的形状。为了保证在结合点处光滑过渡，需要满足连续性条件
三次参数样条曲线拟合代码
05-21
以下是使用三次参数样条曲线进行拟合的 Python 代码示例： ```python import numpy as np from scipy.interpolate import CubicSpline import matplotlib.pyplot as plt # 生成一些随机数据 x = np.array([1, 2, 3, 4, 5, 6]) y = np.array([3, 5, 4, 6, 2, 4]) # 对数据进行三次参数样条曲线拟合 cs = CubicSpline(x, y) # 生成拟合曲线上的点 x_fit = np.linspace(x.min(), x.max(), 100) y_fit = cs(x_fit) # 绘制原始数据和拟合曲线 plt.plot(x, y, 'o', label='Data') plt.plot(x_fit, y_fit, label='Cubic Spline') plt.legend() plt.show() ``` 该代码使用了 `scipy` 库中的 `CubicSpline` 方法进行拟合，生成拟合曲线上的点，最后使用 `matplotlib` 库绘制原始数据和拟合曲线。你可以将其中的 x 和 y 替换为你自己的数据。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值