求经过不共线三点Pi(xi,yi) (i=1,2,3)的圆方程的矩阵形式.用mathematica或matlab求解
来源:学生作业帮 编辑:神马作文网作业帮 分类:数学作业 时间:2024/11/10 19:24:09
求经过不共线三点Pi(xi,yi) (i=1,2,3)的圆方程的矩阵形式.用mathematica或matlab求解
我自己用mathematica做的时候发现这三点必须要具体数值才能得出结果,不知道怎么弄才可以直接用参数xi,yi表示,请赐教~
我自己用mathematica做的时候发现这三点必须要具体数值才能得出结果,不知道怎么弄才可以直接用参数xi,yi表示,请赐教~
先用解析几何推导计算圆心的公式,然后用坐标平移和旋转的方法得到吧,基本推导过程可以这样实现:
(1)求圆心
设O(a,b,c)是圆心坐标,M,N分别是P1P2,P2P3量线段的中点.建立向量MO和NO,根据三个条件写方程:MO点乘P1P2=0;NO点乘P2P3=0;(MO叉乘NO)叉乘(P1P2叉乘P2P3)=0.
三个方程中只有第三个方程是2次方程,用matlab符号函数求解,得到用xi,yi表示的圆心坐标.
(2)求旋转角度
利用Pi平面法向量n=P1P2叉乘P2P3可以求出两个旋转矩阵,使法向量n与某个坐标轴平行.这两个矩阵也可以用xi,yi表示出来
(3)平移和旋转
旋转矩阵有了,平移坐标(a,b,c)也有了,且都是用xi,yi表示的,所以可以自己写出圆的方程了.即设圆的方程为
(X-a)^2+(Y-b)^2+(Z-c)^2=R^2 (1)
其中
[X;Y;Z]=AB[x,y,z] (2) A,B为旋转矩阵
由(2)式可用x,y,z表示出X,Y,Z带入到(1).
ok,圆方程用xi,yi表示出来了.
(1)求圆心
设O(a,b,c)是圆心坐标,M,N分别是P1P2,P2P3量线段的中点.建立向量MO和NO,根据三个条件写方程:MO点乘P1P2=0;NO点乘P2P3=0;(MO叉乘NO)叉乘(P1P2叉乘P2P3)=0.
三个方程中只有第三个方程是2次方程,用matlab符号函数求解,得到用xi,yi表示的圆心坐标.
(2)求旋转角度
利用Pi平面法向量n=P1P2叉乘P2P3可以求出两个旋转矩阵,使法向量n与某个坐标轴平行.这两个矩阵也可以用xi,yi表示出来
(3)平移和旋转
旋转矩阵有了,平移坐标(a,b,c)也有了,且都是用xi,yi表示的,所以可以自己写出圆的方程了.即设圆的方程为
(X-a)^2+(Y-b)^2+(Z-c)^2=R^2 (1)
其中
[X;Y;Z]=AB[x,y,z] (2) A,B为旋转矩阵
由(2)式可用x,y,z表示出X,Y,Z带入到(1).
ok,圆方程用xi,yi表示出来了.
关于偏微分方程组mathematica或matlab编程求解,有没有简单点的办法,方程见图
如何利用Matlab或Mathematica做4y+(xy/siny)=Pi^2的二维x-y函数图.x=0~1之间.
matlab 或 mathematica 接近奇异的矩阵的逆
i,j是两个不共线的向量,已知向量AB=3i+2j,向量CB=i+入j,向量CD=-2i+j,若A,B,D三点共线,试求
matlab求解带符号的矩阵方程
方程为(d^2x)/(dt^2)+a*sin(x)=0的mathematica 和matlab求解命令
已知由样本数据点集{(xi,yi)|i=1,2,……,n},求得的回归直线方程为为^y=1.23x+0.08,且 .
用matlab求解矩阵方程AX=B-2X
i.j是两个不共线向量,已知AB向量=3i+2j,CB向量=i+λj,CD=-2i+j,若ABD三点共线,求实数λ的值.
用matlab求解一个方程:arcsin(1/a)+1/a*sqrt(1-(1/a)^2)=pi/4
用matlab求下面的方程,y=a*exp[-pi*(x-b)^2/c^2]+ d*exp[-pi*(x-e)2/f^2
matlab 矩阵方程求解