如何在Matlab中绘制由两个向量跨越的平面?

时间:2019-05-04 20:01:02

标签: matlab user-interface linear-algebra

我有一个矩阵,可以说

value = [ 1,0,0;0,1,0]

,并且需要在这些向量生成的空间中绘制平面,如下所示: enter image description here
加上下一个代码

 tt=-0.5:0.1:1.5;
 ss=-0.5:0.1:1.5;
 [s,t]=meshgrid(ss,tt);
 grid on
 mesh(0+(value(1:1)*t)+(value(2:1)*s),0+(value(1:2)*t)+(value(2:2)*s),(0+value(1:3)*t)+(value(2:3)*s))
 view(45, 45);

但是我所得到的只是下一个数字 enter image description here
Matlab在mesh命令行中向我发送错误。

说明了网格指令
我曾经用过这样的想法:平面需要一个点和两个向量
点(0,0,0)
u = 1,0,0
v = 0,1,0
因此在参数方程式中您将拥有

w =(x =(0 + t + 0s),y =(0 + 0t + s),z =(0 + 0t + 0s)
但是这个想法可以推广到任何大小的矩阵,因此,它不是固定值,而是根据网格部分中的矩阵元素表示的。 我见过几种方法,但是这些方法涉及到平面的极限(取决于矢量,因此并不总是知道),并使用了诸如patch之类的方法。

1 个答案:

答案 0 :(得分:1)

就像@ Dev-il所说的,您只需要使用cross函数

% Generate vectors
v1 = rand(3,1) - 0.5 ;
v1 = v1 / norm(v1) ;
v2 = rand(3,1) - 0.5 ;
v2 = v2 / norm(v2 ) ;

% Generate point
p = 10*(rand(3,1) - 0.5);

% Start actual code
v3 = cross(v1,v2); % get plane's normal vector

% Points on plane
[ x , y ] = meshgrid( p(1)+(-5:5) , p(2)+(-5:5) );

% Equation for a plane
z = p(3) - (v3(1)*(x-p(1)) + v3(2)*(y-p(2)))/v3(3);

% plot plane
surf(x,y,z)

% plot vectors
hold all
v1Plot = num2cell( [ p , p + v1 ] , 2 );
v2Plot = num2cell( [ p , p + v2 ] , 2 );
v3Plot = num2cell( [ p , p + v3 ] , 2 );
plot3( v1Plot{:}, v2Plot{:}, v3Plot{:} )
hold off

axis equal