具有形状参数的新型广义混合三角贝齐尔样曲线（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。

相关文章

Qt实现对Word网页的读取功能

Qt实现对Word网页的读取功能

《Qt实现对Word网页的读取功能》文章介绍了几种在Qt中实现Word文档（.docx/.doc）读写功能的方法,包括基于QAxObject的COM接口调用、DOCX模板替换及跨平台解决方案,重点讨论... 目录1. 核心实现方式2. 基于QAxObject的COM接口调用（Windows专用）2.1 环境
阅读更多...
MySQL查看表的历史SQL的几种实现方法

MySQL查看表的历史SQL的几种实现方法

《MySQL查看表的历史SQL的几种实现方法》：本文主要介绍多种查看MySQL表历史SQL的方法,包括通用查询日志、慢查询日志、performance_schema、binlog、第三方工具等,并... 目录mysql 查看某张表的历史SQL1.查看MySQL通用查询日志（需提前开启）2.查看慢查询日志3.
阅读更多...
Java实现字符串大小写转换的常用方法

Java实现字符串大小写转换的常用方法

《Java实现字符串大小写转换的常用方法》在Java中,字符串大小写转换是文本处理的核心操作之一,Java提供了多种灵活的方式来实现大小写转换,适用于不同场景和需求,本文将全面解析大小写转换的各种方法... 目录前言核心转换方法1.String类的基础方法2. 考虑区域设置的转换3. 字符级别的转换高级转换
阅读更多...
使用Python实现局域网远程监控电脑屏幕的方法

使用Python实现局域网远程监控电脑屏幕的方法

《使用Python实现局域网远程监控电脑屏幕的方法》文章介绍了两种使用Python在局域网内实现远程监控电脑屏幕的方法,方法一使用mss和socket,方法二使用PyAutoGUI和Flask,每种方... 目录方法一：使用mss和socket实现屏幕共享服务端（被监控端）客户端（监控端）方法二：使用PyA
阅读更多...
MyBatis-Plus逻辑删除实现过程

MyBatis-Plus逻辑删除实现过程

《MyBatis-Plus逻辑删除实现过程》本文介绍了MyBatis-Plus如何实现逻辑删除功能,包括自动填充字段、配置与实现步骤、常见应用场景,并展示了如何使用remove方法进行逻辑删除,逻辑删... 目录1. 逻辑删除的必要性编程1.1 逻辑删除的定义1.2 逻辑删php除的优点1.3 适用场景2.
阅读更多...
C#借助Spire.XLS for .NET实现在Excel中添加文档属性

C#借助Spire.XLS for .NET实现在Excel中添加文档属性

《C#借助Spire.XLSfor.NET实现在Excel中添加文档属性》在日常的数据处理和项目管理中,Excel文档扮演着举足轻重的角色,本文将深入探讨如何在C#中借助强大的第三方库Spire.... 目录为什么需要程序化添加Excel文档属性使用Spire.XLS for .NET库实现文档属性管理Sp
阅读更多...
Python+FFmpeg实现视频自动化处理的完整指南

Python+FFmpeg实现视频自动化处理的完整指南

《Python+FFmpeg实现视频自动化处理的完整指南》本文总结了一套在Python中使用subprocess.run调用FFmpeg进行视频自动化处理的解决方案,涵盖了跨平台硬件加速、中间素材处理... 目录一、 跨平台硬件加速：统一接口设计1. 核心映射逻辑2. python 实现代码二、 中间素材处
阅读更多...
Java数组动态扩容的实现示例

Java数组动态扩容的实现示例

《Java数组动态扩容的实现示例》本文主要介绍了Java数组动态扩容的实现示例,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 目录1 问题2 方法3 结语1 问题实现动态的给数组添加元素效果，实现对数组扩容，原始数组使用静态分配
阅读更多...
Python实现快速扫描目标主机的开放端口和服务

Python实现快速扫描目标主机的开放端口和服务

《Python实现快速扫描目标主机的开放端口和服务》这篇文章主要为大家详细介绍了如何使用Python编写一个功能强大的端口扫描器脚本,实现快速扫描目标主机的开放端口和服务,感兴趣的小伙伴可以了解下... 目录功能介绍场景应用1. 网络安全审计2. 系统管理维护3. 网络故障排查4. 合规性检查报错处理1.
阅读更多...
Python轻松实现Word到Markdown的转换

Python轻松实现Word到Markdown的转换

《Python轻松实现Word到Markdown的转换》在文档管理、内容发布等场景中,将Word转换为Markdown格式是常见需求,本文将介绍如何使用FreeSpire.DocforPython实现... 目录一、工具简介二、核心转换实现1. 基础单文件转换2. 批量转换Word文件三、工具特性分析优点局
阅读更多...

Copyright China编程 23002807@qq.com

沪ICP备2023033072号-2