相切于球体上定点的平面

文章目录


1.问题描述

该问题同二维平面上求解过圆上一点的直线问题类似。

2.已知条件

已知球体表面上一点P, 球心C,求球面上P点的切面M。

3.求解

令向量 n = P − C n = P-C n=P−C,则 u = ( n y , − n x , 0 ) , v = u × n u=(n_y,-n_x,0),v=u\times n u=(ny,−nx,0),v=u×n

平面 M = P + t ⋅ u + m ⋅ v M = P+t\cdot u +m\cdot v M=P+t⋅u+m⋅v

4.Code

matlab 复制代码
close all;
clear;


rng default;

C = rand(1,3);
r = rand();
t = rand(1,2)*pi;
P = C + [r*cos(t(1))*sin(t(2)), r*sin(t(1))*sin(t(2)), r*cos(t(2))];

t = linspace(0, 2*pi, 50);
u = linspace(0, 2*pi, 50);
[t,u] = meshgrid(t,u);

x = C(1) + r*cos(t).*sin(u);
y = C(2) + r*sin(t).*sin(u);
z = C(3) + r*cos(u);

n = normalize(P - C,'norm');
m = [0, n(3), -n(2)];
l = cross(m,n);
c1 = linspace(-1,1,20);
c2 = linspace(-1,1,20);
[c1,c2] = meshgrid(c1,c2);
xx = P(1) + c1*m(1) + c2*l(1);
yy = P(2) + c1*m(2) + c2*l(2);
zz = P(3) + c1*m(3) + c2*l(3);


figure;
mesh(x,y,z);
axis equal
hold on
mesh(xx,yy,zz)
hold on
scatter(P(1),P(2), P(3),'filled')
相关推荐
云 无 心 以 出 岫15 分钟前
贪心算法QwQ
数据结构·c++·算法·贪心算法
俏布斯29 分钟前
算法日常记录
java·算法·leetcode
独好紫罗兰32 分钟前
洛谷题单3-P5719 【深基4.例3】分类平均-python-流程图重构
开发语言·python·算法
SheepMeMe1 小时前
蓝桥杯2024省赛PythonB组——日期问题
python·算法·蓝桥杯
随便昵称1 小时前
蓝桥杯专项复习——前缀和和差分
c++·算法·前缀和·蓝桥杯
脑子慢且灵1 小时前
蓝桥杯冲刺:一维前缀和
算法·leetcode·职场和发展·蓝桥杯·动态规划·一维前缀和
姜威鱼1 小时前
蓝桥杯python编程每日刷题 day 21
数据结构·算法·蓝桥杯
CYRUS STUDIO2 小时前
Unidbg Trace 反 OLLVM 控制流平坦化(fla)
android·汇编·算法·网络安全·逆向·ollvm
ゞ 正在缓冲99%…2 小时前
leetcode22.括号生成
java·算法·leetcode·回溯
小卡皮巴拉2 小时前
【力扣刷题实战】矩阵区域和
开发语言·c++·算法·leetcode·前缀和·矩阵