Bezier曲线插值拟合

最新推荐文章于 2024-04-23 19:01:02 发布

我的女友叫苏苏

最新推荐文章于 2024-04-23 19:01:02 发布

阅读量420

点赞数 1

分类专栏： matlab 插值文章标签：算法 matlab

本文链接：https://blog.csdn.net/qq_58533871/article/details/128218822

版权

matlab 同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

插值

2 篇文章 0 订阅

订阅专栏

基本原理

传统的插值方法是经过所有的给定点 $(x_{i},f(x_{i}))$ 而形成的光滑曲线，但随着给定点的变多，简单曲线的求解过程也变得越来越繁琐和不易计算并且会产生一些误差，例如龙格现象，而Bezier曲线与传统插值方法不一样，他只经过给定点中的控制点 $(x_{i_{j}},f(x_{i_{j}}))$ ,而对于给定点中的非控制点并不要求必须经过所求出的插值曲线

插值公式

根据上述原理介绍，我们要构造一条经过那些控制点且对那些非控制点有良好拟合效果的插值曲线，我们不妨先考虑简单的情况：只有两个控制点 $min(x_{i}),max(x_{i})$ 的情况。

在这种简单情况插值公式为

$\begin{bmatrix} x(t)\\y(t) \end{bmatrix}=\sum_{i=0}^{n}P_{i}B_{i,n}(t),P_{i}= \begin{bmatrix} x_{i}\\f(x_{i}) \end{bmatrix},B_{i,n}(t)=C^{i}_{n}t^{i}(1-t)^{n-i},t\in[0,1]$

$P(t)=\begin{bmatrix} x(t) \\ y(t) \end{bmatrix}$ , 显然 $P(0)=\begin{bmatrix} x(0) \\ y(0) \end{bmatrix}=\begin{bmatrix}x_{0}\\f(x_{0}) \end{matrix}, P(1)=\begin{bmatrix} x(1) \\ y(1) \end{bmatrix}=\begin{bmatrix}x_{n}\\f(x_{n}) \end{matrix}$

显然满足上述Bezier曲线的要求。

下面列举关于 $B_{i,n}(t)$ 的一些性质，

这些性质很容易验证，这里不在阐述。

代码部分

clear;clc;close all;
% Bezier曲线
P=[0,2;2,4;4,3;6,0];
pause(1)
t=linspace(0,1,11);
n=3 
x=[];y=[];
for j=1:11
    b=0;b1=0;
    for i=1:n+1
        a=factorial(n)/(factorial(i-1)*factorial(n-i+1));
        b=b+P(i,1)*a*t(j)^(i-1)*(1-t(j))^(n-i+1);
        b1=b1+P(i,2)*a*t(j)^(i-1)*(1-t(j))^(n-i+1);
    end
    x=[x,b];
    y=[y,b1];
end
plot(x,y,'r*-')