polyshape 

函数可创建由二维顶点定义的多边形,并返回具有描述其顶点、实心区域和孔的各种属性的 polyshape 对象。例如,pgon = polyshape([0 0 1 1],[1 0 0 1]) 将创建由四个点 (0,1)、(0,0)、(1,0) 和 (1,1) 定义的实心正方形。

pgon = polyshape(x,y) 从由 x 坐标向量和对应的 y 坐标向量定义的二维顶点创建 polyshapex 和 y 的长度必须相同,且至少要有三个元素。

pgon = polyshape([0 0 1 1],[1 0 0 1]);
figure(4)
plot(pgon)

polyshape 二维多边形构建

nsidedpoly 正多边形

pgon = nsidedpoly(n,Name,Value)

pgon1 = nsidedpoly(6);
pgon2 = nsidedpoly(6,'Center',[5 0],'SideLength',3);
figure(4)
plot([pgon1 pgon2])
axis equal

nsidedpoly(正多边形)

polyarea 计算多边形的面积

L = linspace(0,2.*pi,100);
xv = 1.2*cos(L)';
yv = 1.2*sin(L)';
xx=1.1*cos(L)';
yy=1.1*sin(L)';
AA=polyshape(xv,yv);
AAA=polyshape(xx,yy); plot(AA);
hold on;
plot(AAA);
title(['Area = ' num2str(A)])
axis image A = polyarea(xv,yv) %计算面积

polyarea 计算多边形面积

inpolygon 位于多边形区域边缘内部或边缘上的点

in = inpolygon(xq,yq,xv,yv)  返回 in,以指明 xq 和 yq 所指定的查询点是在 xv 和 yv 定义的多边形区域的边缘内部还是在边缘上

[in,on] = inpolygon(xq,yq,xv,yv) 还返回 on 以指明查询点是否位于多边形区域的边缘。

clc,clear all;
L = linspace(0,2.*pi,100);
xv1 = 1.2*cos(L)';
yv1 = 1.2*sin(L)';
xv2=0.4*cos(L)';
yv2=0.4*sin(L)'; Round1=polyshape(xv1,yv1);
Round2=polyshape(xv2,yv2);
Round1Area = polyarea(xv1,yv1)
Round2Area = polyarea(xv2,yv2) rng default %恢复matlab启动时默认的全局随机流
xq = randn(250,1); %随机生成250个正态分布点
yq = randn(250,1); %随机生成250个正态分布点 in = inpolygon(xq,yq,xv1,yv1); % numel(xq(in))查看在几何体内部的数量
figure(4)
plot(Round1); %画第一个圆
hold on;
plot(Round2); %画第二个圆
plot(xq(in),yq(in),'r+') % 标记在内部的点
plot(xq(~in),yq(~in),'bo') % 标记在外部的点
hold off
title(['Area1,Area2 = ' num2str(Round1Area),num2str(Round2Area)])
axis equal
axis image

凸多边形内的点

axes('xlim',[0 10],'ylim',[0 10])
a=[1 1 2 2];
b=[3 4 5 6];
c=[1 1 3 3];
figure
rectangle('position',a,'EdgeColor', 'r');
rectangle('position',b,'EdgeColor', 'b');
rectangle('position',c,'EdgeColor', 'k'); area_ab = rectint(a,b)
area_ac = rectint(a,c)

rectint 矩形交叉区域

area = rectint(A,B) 返回位置向量 A 和 B 指定的矩形的交叉区域。

如果 A 和 B 分别指定一个矩形,则输出 area 为标量。

A 和 B 还可以是矩阵,其中每行是一个位置向量。area 是一个矩阵,为 B 指定的所有矩形提供 A 指定的所有矩阵的交叉区域。即,如果 A 是 n×4 并且 B 是 m×4,则 area 是 n×m 矩阵,其中 area(i,j) 是 A 的第 i 行和 B 的第 j 行指定的矩形的交叉区域。

a=1;b=1.5;d=0:360;
x=a*cosd(d);
y=b*sind(d);
figure(1),cla
patch(x+0.6,y+0.2,'b','edgecolor','none','facealpha',0.15);%不透明度0.15
patch(x-0.6,y-0.2,'b','edgecolor','none','facealpha',0.15);%不透明度0.15
patch(x-0.2,y+0.9,'b','edgecolor','none','facealpha',0.15);%不透明度0.15
patch(x+0.2,y-0.9,'b','edgecolor','none','facealpha',0.15);%不透明度0.15
axis equal
box on;
set(gcf,'color','w');

半透明重叠加深 patch函数

v1 = [2 4; 2 8; 8 4];
f1 = [1 2 3];
figure
patch('Faces',f1,'Vertices',v1,'FaceColor','red','FaceAlpha',.3); v2 = [2 4; 2 8; 8 8];
f2 = [1 2 3];
patch('Faces',f2,'Vertices',v2,'FaceColor','blue','FaceAlpha',.5);

通过将 FaceAlpha 属性设置为 0 和 1 之间的值,创建两个半透明的多边形。

boundary 二维或三维空间内的一组点的边界

x = gallery('uniformdata',3000,1,1);
y = gallery('uniformdata',3000,1,10); plot(x,y,'k.')
xlim([-0.3 1.2])
ylim([-0.3 1.2])
k = boundary(x,y); %返回一个表示包围点 (x,y) 的单个相容二维边界的点索引向量。收缩因子为1
hold on;
line(x(k),y(k),'color','b'); j = boundary(x,y,0.01); %收缩因子为0.1 最大为1 越大包裹性越紧
line(x(j),y(j),'color','r'); [~, vol] = boundary(x,y); %包围点形成的形状的体积 或 面积
axis equal;
vol %返回 包围的体积或者面积

boundary 二维或三维空间内的一组点的边界

[~, vol] = boundary(x,y); %包围点形成的形状的体积
axis equal;
vol 可以输出二维平面 各个点包围的面积

P = gallery('uniformdata',30,3,5);
subplot(1,2,1);
plot3(P(:,1),P(:,2),P(:,3),'.','MarkerSize',10)
axis equal;
grid on
k = boundary(P);
hold on
subplot(1,2,2);
trisurf(k,P(:,1),P(:,2),P(:,3),'Facecolor','red','FaceAlpha',0.1) %收缩因子默认0
[~, vol] = boundary(P); %包围点形成的形状的体积
axis equal;
vol %体积

三维

[~, vol] = boundary(P);   vol 返回体积

convhull 凸包

xx = -1:.05:1;
yy = abs(sqrt(xx));
[x,y] = pol2cart(xx,yy);
k = convhull(x,y);
plot(x(k),y(k),'r-',x,y,'bo')

二维凸包

alphaShape 依据二维和三维中的点构建的多边形和多面体

th = (pi/12:pi/12:2*pi)';
x1 = [reshape(cos(th)*(1:5), numel(cos(th)*(1:5)),1); 0];
y1 = [reshape(sin(th)*(1:5), numel(sin(th)*(1:5)),1); 0];
x = [x1; x1+15];
y = [y1; y1];
subplot(1,2,1)
plot(x,y,'.')
axis equal
subplot(1,2,2)
shp = alphaShape(x,y); %默认 alpha 半径可生成带不规则边界的 alpha 形状。
shp.Alpha = 2.5; % 要更好地捕获点集边界,请尝试更大的alpha 半径。
plot(shp)
axis equal

alphaShape alphaShape 创建一个可将一组二维或三维点包围起来的边界面或三维体。

delaunayTriangulation

使用 delaunayTriangulation 对象可以基于一组点创建二维或三维 Delaunay 三角剖分。对于二维数据,您也可以指定边约束。

P = gallery('uniformdata',[30 2],0);
DT = delaunayTriangulation(P);
IC = incenter(DT);
triplot(DT) %绘制 Delaunay 三角剖分
hold on
plot(IC(:,1),IC(:,2),'*r')

绘图

[matlab] 19.matlab 基础几何学的更多相关文章

  1. 第二章 Matlab面向对象编程基础

    DeepLab是一款基于Matlab面向对象编程的深度学习工具箱,所以了解Matlab面向对象编程的特点是必要的.笔者在做Matlab面向对象编程的时候发现无论是互联网上还是书店里卖的各式Matlab ...

  2. Matlab面向对象编程基础

    DeepLab是一款基于Matlab面向对象编程的深度学习工具箱,所以了解Matlab面向对象编程的特点是必要的.笔者在做Matlab面向对象编程的时候发现无论是互联网上还是书店里卖的各式Matlab ...

  3. [.net 面向对象编程基础] (19) LINQ基础

    [.net 面向对象编程基础] (19)  LINQ基础 上两节我们介绍了.net的数组.集合和泛型.我们说到,数组是从以前编程语言延伸过来的一种引用类型,采用事先定义长度分配存储区域的方式.而集合是 ...

  4. 机器学习及其matlab实现—从基础到实践

    第1周 MATLAB入门基础 第2周 MATLAB进阶与提高 第3周 BP神经网络 第4周 RBF.GRNN和PNN神经网络 第5周 竞争神经网络与SOM神经网络 第6周 支持向量机(Support ...

  5. Python matlab octave 矩阵运算基础

    基础总结,分别在三种软件下,计算 求逆矩阵 矩阵转置 等运算,比较异同 例子:正规方程法求多元线性回归的最优解 θ=(XTX)-1XTY octave: pwd()当前目录 ones() zeros( ...

  6. 2014.08.04,读书,读书笔记-《Matlab概率与数理统计分析》-第1章 MATLAB的数据基础

    第1章 MATLAB数据基础 虽然一直间或使用MATLAB,但从来没有系统的学习过,现在开始也不晚.先对几个重点或者平时忽略的要点做下笔记. %后的所有文字为注释,多条命令可以放在一行,但要用逗号或分 ...

  7. 【Matlab】调试基础

    1.matlab 调试子程序 在主程序进入子程序前一句加断点,然后step in,可以进入子程序. 但是直接在子程序里设置断点,运行主程序是不能进入子程序的.

  8. MATLAB的一些基础知识

    1.已知a1=sin(sym(pi/4)+exp(sym(0.7)+sym(pi/3)))产生精准符号数字,请回答:以下产生的各种符号数哪些是精准的?若不精准,误差又是多少?能说出产生误差的原因吗? ...

  9. MATLAB数字图像处理基础

    图像的输入.输出和显示 1.图像的输入    imread('filename'),  实际中写的是 >> f = imread('sky.jpg'); 2.图像的显示    imshow ...

随机推荐

  1. STL中的Set用法(详+转)

    set是STL中一种标准关联容器(vector,list,string,deque都是序列容器,而set,multiset,map,multimap是标准关联容器),它底层使用平衡的搜索树——红黑树实 ...

  2. pullMsg有感

    开发功能过程中,始终会有些东西是确认的,比如美丑.业务是否合理.对错. 如果明知道不合理,却按照已有规定.框架.设计去开发,其实是不够职业. 好的做法是朝对的方向去push,并落地: 次之是去push ...

  3. MongoDB复合索引详解

    摘要: 对于MongoDB的多键查询,创建复合索引可以有效提高性能. 什么是复合索引? 复合索引,即Compound Index,指的是将多个键组合到一起创建索引,这样可以加速匹配多个键的查询.不妨通 ...

  4. git入门 多人协作

    当你从远程仓库克隆时,实际上Git自动把本地的master分支和远程的master分支对应起来了,并且,远程仓库的默认名称是origin. 要查看远程库的信息,用git remote: $ git r ...

  5. Laravel 系列入门教程(一)【最适合中国人的 Laravel 教程】

    热烈庆祝 Laravel 5.5 LTS 发布! 实际上 Laravel 上一个 LTS 选择 5.1 是非常不明智的,因为 5.2 增加了许许多多优秀的特性.现在好了,大家都用最新的长期支持版本 5 ...

  6. canvas学习笔记之2d画布基础的实现

    一. Canvas是啥 < canvas > 是一个可以使用脚本(通常是js)来绘图的HTML元素 < canvas > 最早由Apple引入WebKit,用于Mac OS X ...

  7. P1993 小 K 的农场

    题目描述 小 K 在 Minecraft 里面建立很多很多的农场,总共 n 个,以至于他自己都忘记了每个 农场中种植作物的具体数量了,他只记得一些含糊的信息(共 m 个),以下列三种形式描 述: 农场 ...

  8. django请求和响应

    本文转载自https://blog.csdn.net/xiaogeldx/article/details/88096341 HttpRequest对象 服务器接收到http协议的请求后,会根据报文创建 ...

  9. 18.Odoo产品分析 (二) – 商业板块(10) – 电子商务(2)

    查看Odoo产品分析系列--目录 接上一篇Odoo产品分析 (二) – 商业板块(10) – 电子商务(1) 6. 高级属性 除了我们到目前为止已经覆盖基本选项,Odoo在产品页面还提供了一些高级选项 ...

  10. Android assets文件夹之位置放置和作用

    Android 的assets文件夹的放置位置,Eclipse创建项目时就生成了的,Android Studio则不太一样,AS可以包含几种方式, 1:可以在build.gradle文件下配置,加如下 ...