matlab三维向量,MATLAB三维空间绘制向量和面

http://www.cnblogs.com/shenyuanyuan/p/3524290.html

MATLAB三维空间绘制向量和面

1.绘制带箭头向量图,见下面实例:

clear;clc;

A=[4,5,6];

B=[-10,6,7];

C=[-1 8,0];

scale=1;

quiver3(A(1),A(2),A(3),B(1)-A(1),B(2)-A(2),B(3)-A(3),scale);%表示以A为起点,由A指向B

text(A(1),A(2),A(3),'A');%在A点附件标注字母A

axis([-10 4 5 8 0 7]);xlabel('x');ylabel('y');zlabel('z');%规范x,y,z坐标轴刻度范围,及在各自坐标轴上标注字母x,y,z

grid on;%绘网格

hold on;

quiver3(B(1)-0.8,B(2),B(3),C(1)-B(1),C(2)-B(2),C(3)-B(3),scale);

text(B(1)-0.3,B(2),B(3),'B');

hold on;

quiver3(A(1),A(2),A(3),C(1)-A(1),C(2)-A(2),C(3)-A(3),scale);

text(C(1),C(2),C(3)-0.5,'C');

e5f65329ddf15094f6007df74d6623f7.gif

2.在三维空间中绘制面

可用mesh和surf两个函数,见以下实例:

[x,y]=meshgrid(-5:0.01:5);

z=sin(10.*x)-y.^3+x.^2.*y;

mesh(x,y,z);

fe39431a2d5b7af779e705082c58c716.png

遇到更多的信息会继续补充。

posted on 2014-01-17 16:34 follow your heart 阅读(...) 评论(...) 编辑 收藏

转载于:https://www.cnblogs.com/shenyuanyuan/p/3524290.html

标签:10,scale,和面,quiver3,MATLAB,绘制,向量,三维空间

来源: https://blog.csdn.net/weixin_30254435/article/details/98307765

  • 0
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
假设已知的三维空间点坐标及法向量分别存储在 $P$ 和 $N$ 两个 $n \times 3$ 的矩阵中,其中 $n$ 表示点的数量。 首先,我们需要构造 B 样条基函数。这里我们使用经典的 Cox-de Boor 递推公式进行计算。具体实现如下: ```matlab function [N, U] = bspline_basis(u, k, t) % u: 插值点的位置 % k: B样条的次数 % t: 节点向量 % N: B样条基函数值 % U: 扩展过的节点向量 n = numel(t) - k - 1; % B样条的数量 U = [t(1:k) u t(n+1:n+k)]; % 扩展节点向量 N = zeros(n, numel(u)); % 存储 B样条基函数值 for i = 1:numel(u) for j = 1:n if U(j) <= u(i) && u(i) < U(j+1) N(j,i) = 1; end end end for d = 2:k for j = 1:n if N(j,i) ~= 0 alpha = (u(i) - U(j)) / (U(j+d-1) - U(j)); N(j,i) = alpha * N(j,i-1) + (1-alpha) * N(j+1,i-1); end end end end ``` 接下来,我们可以使用 B 样条插值方法对点和法向量进行插值。具体实现如下: ```matlab function [P_interpolated, N_interpolated] = bspline_interpolation(P, N, num_points, degree) % P: 原始点的坐标 (n x 3) % N: 原始点的法向量 (n x 3) % num_points: 插值后的点数 % degree: B样条的次数 % 构造节点向量 n = size(P, 1); t = linspace(0, 1, n+degree+1); t = [zeros(1,degree) t(1:end-degree)]; % 构造插值点的位置 u = linspace(0, 1, num_points); % 计算 B样条基函数 [Nu, U] = bspline_basis(u, degree, t); % 对点坐标插值 [Nn, ~] = bspline_basis(u, degree, t); % 对法向量插值 % 计算插值后的点坐标和法向量 P_interpolated = zeros(num_points, 3); N_interpolated = zeros(num_points, 3); for i = 1:3 P_interpolated(:,i) = Nu * P(:,i); N_interpolated(:,i) = Nn * N(:,i); end end ``` 这里我们使用了与点坐标插值相同的节点向量和 B 样条基函数对法向量进行插值。实际上,我们也可以使用不同的节点向量和 B 样条基函数对法向量进行插值,以达到更好的效果。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值