下面是我的思路,盡量用Matlab語(yǔ)言敘述的,方便你作圖.
假設(shè)：(x1,y1,z1),(x2,y2,z2),(x3,y3,z3),(x0,y0,z0),R,(A,B,C,D)均已知.
法向量(A,B,C)歸一化后,設(shè)單位向量K=[A B C]'/sqrt(A^2+B^2+C^2).設(shè)單位向量I=[x1-x0 y1-y0 z1-z0]'/sqrt((x1-x0)^2+(y1-y0)^2+(z1-z0)^2),設(shè)J=K × I（Matlab里的向量叉乘表示為J=cross(K,I)）,則亦為單位向量.且I,J和K滿(mǎn)足右手定則,可立為新坐標(biāo)系.
則T=[I J K] 為正交矩陣,而且是坐標(biāo)旋轉(zhuǎn)的矩陣,把[1 0 0]'變?yōu)镮,[0 1 0]'變?yōu)镴,[0 0 1]'變?yōu)镵.
給定一個(gè)角度Theta,如下計(jì)算式
[x0 y0 z0]'+T*[R*cos(Theta) R*sin(Theta) 0]'
即為該圓的參數(shù)方程,把Theta取個(gè)0到2*pi的循環(huán)即可畫(huà)出此圓.謝謝，計(jì)算式不太懂，這個(gè)計(jì)算式和參數(shù)方程是什么關(guān)系？參數(shù)方程不是應(yīng)該是x=f(t),y=f(t),z=f(t),t=?是的，計(jì)算式子里的Theta（角度），就是你說(shuō)的t?？赡銢](méi)有給參數(shù)方程啊設(shè)空間旋轉(zhuǎn)矩陣T=[T11 T12 T13 T21 T22 T23 T31 T32 T33]為已知。圓心坐標(biāo)(x0,y0,z0)為已知，圓半徑R為已知，則x(t)=x0+(T11*R*cos(t)+T12*R*sin(t))y(t)=y0+(T21*R*cos(t)+T22*R*sin(t))z(t)=z0+(T31*R*cos(t)+T32*R*sin(t))即為參數(shù)方程。