具有形状参数的新型广义混合三角贝齐尔样曲线（GHT_Bezier_curves） matlab实现

2023-10-20 18:40
文章标签 具有 参数 matlab 实现 新型 曲线 混合 形状 三角 广义 curves bezier ght 贝齐尔样

本文主要是介绍具有形状参数的新型广义混合三角贝齐尔样曲线（GHT_Bezier_curves） matlab实现，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

具有形状参数的新型广义混合三角贝齐尔样曲线

    • 一、n维 GHT-Bernstein 基函数
      • 1、定义
      • 2、性质：
    • 二、基函数的matlab程序实现
      • 1、程序：
      • 2.2 基函数图像
    • 三、n维GHT-Bézier曲线（GHT-Bézier Curves of Degree n）
      • 相应程序：
      • 相应结果：
      • 注：得到的曲线应该满足凸包性质，但是有几个参数得到的曲线不在控制多边形里面（如b图中红色曲线），程序检查了几遍没查出来原因，因为这篇文章就是借鉴一下想法，就不继续研究下去了。

       
       

参考文献：Geometric Modeling of Novel Generalized Hybrid Trigonometric Bézier-Like Curve with Shape Parameters and Its Applications
       
本博客是对上述参考文献的部分程序实现

一、n维 GHT-Bernstein 基函数

1、定义

在这里插入图片描述
       论文中给出的 q 2 , 2 ( θ ) q_{ 2,2}(θ) q2,2​(θ)定义不对，不满足端点性质。需要修改一点。
       下面是我修改的定义，以及写出的三阶公式：

在这里插入图片描述

2、性质：

  1. 单位分解性。
  2. 正性。
  3. 对称性。
  4. 端点插值性。
    在这里插入图片描述
    在这里插入图片描述

二、基函数的matlab程序实现

1、程序：

function GHT_Bernstein
%  三阶GHT_Bernstein基函数图像
clear;clc;
% Gma = 1;        %Gma = γ
% v = 0.2;
% beta = 0;         %β
Gma = 0.5;                                   %  Gma = γ
beta = 0;                                    %  β
v = -0.5;
x = (0:0.01:1)';
y = F0_3(x,v,Gma);
y1 = F2_3(x,v,beta,Gma);
y2 = F1_3(x,v,beta,Gma);
y3 = F3_3(x,v,Gma);
plot(x,y,'y')
title('三阶GHT-Bernstein基函数')
hold on
plot(x,y1,'b')
hold on
plot(x,y2,'g')
hold on
plot(x,y3,'r')% 二阶GHT_Bernstein基函数图像  （可以直接在命令行使用）
Gma = 0.5;                                   %  Gma = γ
beta = 0;                                    %  β
v = -0.5;x = (0:0.01:1)';
y = F1_2(x,v,beta,Gma);
y1 = F2_2(x,beta,Gma);
y2 = F0_2(x,v,Gma);
figure(2)
plot(x,y,'y')
title('二阶GHT-Bernstein基函数')
hold on
plot(x,y1,'b')
hold on
plot(x,y2,'r')
endfunction g = F0_2(x,v,Gma)
g = (1-sin(0.5*pi*x)).*(1-v*sin(0.5*pi*x)).*exp(Gma*x);
endfunction g = F1_2(x,v,beta,Gma)
g = 1 - F0_2(x,v,Gma) - F2_2(x,beta,Gma);
endfunction g = F2_2(x,beta,Gma)
g = (1-cos(0.5*pi*x)).*(1-beta*cos(0.5*pi*x)).*exp((1-x)*Gma);
endfunction g = F0_3(x,v,Gma)
g = (1-x).*F0_2(x,v,Gma);
endfunction g = F1_3(x,v,beta,Gma)
g = (1-x).*F1_2(x,v,beta,Gma) + x.*F0_2(x,v,Gma);
endfunction g = F2_3(x,v,beta,Gma)
g = (1-x).*F2_2(x,v,Gma)+x.*F1_2(x,v,beta,Gma);
endfunction g = F3_3(x,v,Gma)
g =  x.*F2_2(x,v,Gma);
end

2.2 基函数图像

在这里插入图片描述
在这里插入图片描述

三、n维GHT-Bézier曲线（GHT-Bézier Curves of Degree n）

1
在这里插入图片描述
在这里插入图片描述

相应程序：

function [Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma)
px=[0,2,3,5];
py=[0,1,0,1];                              %  控制顶点(0,0),(2,1),(3,0),(5,1)x = (0:0.005:1)';
y0 = F0_3(x,v,Gma);
y1 = F1_3(x,v,beta,Gma);
y2 = F2_3(x,v,beta,Gma);
y3 = F3_3(x,v,Gma);Px = y0*px(1) + y1.*px(2) + y2*px(3) + y3.*px(4);
Py = y0*py(1) + y1.*py(2) + y2*py(3) + y3.*py(4);endfunction g = F0_2(x,v,Gma)
g = (1-sin(0.5*pi*x)).*(1-v*sin(0.5*pi*x)).*exp(Gma*x);
endfunction g = F1_2(x,v,beta,Gma)
g = 1 - F0_2(x,v,Gma) - F2_2(x,beta,Gma);
endfunction g = F2_2(x,beta,Gma)
g = (1-cos(0.5*pi*x)).*(1-beta*cos(0.5*pi*x)).*exp((1-x)*Gma);
endfunction g = F0_3(x,v,Gma)
g = (1-x).*F0_2(x,v,Gma);
endfunction g = F1_3(x,v,beta,Gma)
g = (1-x).*F1_2(x,v,beta,Gma) + x.*F0_2(x,v,Gma);
endfunction g = F2_3(x,v,beta,Gma)
g = (1-x).*F2_2(x,v,Gma)+x.*F1_2(x,v,beta,Gma);
endfunction g = F3_3(x,v,Gma)
g =  x.*F2_2(x,v,Gma);
end

将上述代码保存为函数，用下面代码调用实现：

clear;clc;
%调用GHT_Bezier_curves(v,beta,Gma)函数
%  论文Figure2 中的(a)图像的实现
Gma = 1;                                   %  Gma = γ
beta = 0;                                    %  β
v = 0.95;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
figure(1)
plot(Px,Py,'r',px,py,'b',px,py,'mo')
xlabel('x-axis'), ylabel('y-axes');
title('figure2 (a)')
hold on
v = 0.5;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'y',px,py,'b',px,py,'m*')
hold on
v = -0.1;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'m',px,py,'b',px,py,'m*')
hold on
v = -0.6;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'g',px,py,'b',px,py,'m*')
hold on%  论文Figure2 中的(b)图像的实现
Gma = 0.5;                                   %  Gma = γ
beta = 0;                                    %  β
v = -0.5;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
figure(2)
plot(Px,Py,'r',px,py,'b',px,py,'mo')
xlabel('x-axis'), ylabel('y-axes');
title('figure2 (b)')
hold on
v = 0;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'y',px,py,'b',px,py,'m*')
hold on
v = 0.5;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'m',px,py,'b',px,py,'m*')
hold on
v = 1;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'g',px,py,'b',px,py,'m*')
hold on%  论文Figure2 中的(d)图像的实现
Gma = 0;                                   %  Gma = γ
beta = 1;                                    %  β
v = 1;
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
figure(4)
plot(Px,Py,'r',px,py,'b',px,py,'mo')
xlabel('x-axis'), ylabel('y-axes');
title('figure2 (d)')
hold on
Gma = 0.5;   
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'y',px,py,'b',px,py,'m*')
hold on
Gma = -0.5;   
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'m',px,py,'b',px,py,'m*')
hold on
Gma = -1;   
[Px,Py,px,py] = GHT_Bezier_curves(v,beta,Gma);
plot(Px,Py,'g',px,py,'b',px,py,'m*')
hold on

相应结果：

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

注：得到的曲线应该满足凸包性质，但是有几个参数得到的曲线不在控制多边形里面（如b图中红色曲线），程序检查了几遍没查出来原因，因为这篇文章就是借鉴一下想法，就不继续研究下去了。

      
       
       

这篇关于具有形状参数的新型广义混合三角贝齐尔样曲线（GHT_Bezier_curves） matlab实现的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！


http://www.chinasem.cn/article/248979。

相关文章

Nginx实现高并发的项目实践

Nginx实现高并发的项目实践

《Nginx实现高并发的项目实践》本文主要介绍了Nginx实现高并发的项目实践,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录使用最新稳定版本的Nginx合理配置工作进程(workers)配置工作进程连接数(worker_co
阅读更多...
python中列表list切分的实现

python中列表list切分的实现

《python中列表list切分的实现》列表是Python中最常用的数据结构之一,经常需要对列表进行切分操作,本文主要介绍了python中列表list切分的实现,文中通过示例代码介绍的非常详细,对大家... 目录一、列表切片的基本用法1.1 基本切片操作1.2 切片的负索引1.3 切片的省略二、列表切分的高
阅读更多...
基于Python实现一个PDF特殊字体提取工具

基于Python实现一个PDF特殊字体提取工具

《基于Python实现一个PDF特殊字体提取工具》在PDF文档处理场景中,我们常常需要针对特定格式的文本内容进行提取分析,本文介绍的PDF特殊字体提取器是一款基于Python开发的桌面应用程序感兴趣的... 目录一、应用背景与功能概述二、技术架构与核心组件2.1 技术选型2.2 系统架构三、核心功能实现解析
阅读更多...
使用Python实现表格字段智能去重

使用Python实现表格字段智能去重

《使用Python实现表格字段智能去重》在数据分析和处理过程中,数据清洗是一个至关重要的步骤,其中字段去重是一个常见且关键的任务,下面我们看看如何使用Python进行表格字段智能去重吧... 目录一、引言二、数据重复问题的常见场景与影响三、python在数据清洗中的优势四、基于Python的表格字段智能去重
阅读更多...
Spring AI集成DeepSeek实现流式输出的操作方法

Spring AI集成DeepSeek实现流式输出的操作方法

《SpringAI集成DeepSeek实现流式输出的操作方法》本文介绍了如何在SpringBoot中使用Sse（Server-SentEvents）技术实现流式输出,后端使用SpringMVC中的S... 目录一、后端代码二、前端代码三、运行项目小天有话说题外话参考资料前面一篇文章我们实现了《Spring
阅读更多...
Nginx中location实现多条件匹配的方法详解

Nginx中location实现多条件匹配的方法详解

《Nginx中location实现多条件匹配的方法详解》在Nginx中,location指令用于匹配请求的URI,虽然location本身是基于单一匹配规则的,但可以通过多种方式实现多个条件的匹配逻辑... 目录1. 概述2. 实现多条件匹配的方式2.1 使用多个 location 块2.2 使用正则表达式
阅读更多...
使用Apache POI在Java中实现Excel单元格的合并

使用Apache POI在Java中实现Excel单元格的合并

《使用ApachePOI在Java中实现Excel单元格的合并》在日常工作中,Excel是一个不可或缺的工具,尤其是在处理大量数据时,本文将介绍如何使用ApachePOI库在Java中实现Excel... 目录工具类介绍工具类代码调用示例依赖配置总结在日常工作中，Excel 是一个不可或缺的工http://
阅读更多...
SpringBoot实现导出复杂对象到Excel文件

SpringBoot实现导出复杂对象到Excel文件

《SpringBoot实现导出复杂对象到Excel文件》这篇文章主要为大家详细介绍了如何使用Hutool和EasyExcel两种方式来实现在SpringBoot项目中导出复杂对象到Excel文件,需要... 在Spring Boot项目中导出复杂对象到Excel文件，可以利用Hutool或EasyExcel
阅读更多...
Python如何实现读取csv文件时忽略文件的编码格式

Python如何实现读取csv文件时忽略文件的编码格式

《Python如何实现读取csv文件时忽略文件的编码格式》我们再日常读取csv文件的时候经常会发现csv文件的格式有多种,所以这篇文章为大家介绍了Python如何实现读取csv文件时忽略文件的编码格式... 目录1、背景介绍2、库的安装3、核心代码4、完整代码1、背景介绍我们再日常读取csv文件的时候经常
阅读更多...
Golang中map缩容的实现

Golang中map缩容的实现

《Golang中map缩容的实现》本文主要介绍了Go语言中map的扩缩容机制,包括grow和hashGrow方法的处理,具有一定的参考价值,感兴趣的可以了解一下... 目录基本分析带来的隐患为什么不支持缩容基本分析在 Go 底层源码 src/runtime/map.go 中，扩缩容的处理方法是 grow
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2