我有一个矩阵,可以说
value = [ 1,0,0;0,1,0]
,并且需要在这些向量生成的空间中绘制平面,如下所示:
加上下一个代码
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);
但是我所得到的只是下一个数字
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之类的方法。
答案 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