那位高手用MATLAB帮我解个六元四次方程组eq1='((l-k)*(n*l-p*k)+(h-l)*(p*h-m*l)+(k-h)*(m*k-n*h))/(sqrt((l-k)^2+(h-l)^2+(k-h)^2)*sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2))=a'eq2='(u*(l-k)+v*(h-l)+w*(k-h))/sqrt((l-k)^2+(h-l)^2+(k-h)^2)=b'eq

那位高手用MATLAB帮我解个六元四次方程组eq1='((l-k)*(n*l-p*k)+(h-l)*(p*h-m*l)+(k-h)*(m*k-n*h))/(sqrt((l-k)^2+(h-l)^2+(k-h)^2)*sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2))=a'eq2='(u*(l-k)+v*(h-l)+w*(k-h))/sqrt((l-k)^2+(h-l)^2+(k-h)^2)=b'eq
那位高手用MATLAB帮我解个六元四次方程组
eq1='((l-k)*(n*l-p*k)+(h-l)*(p*h-m*l)+(k-h)*(m*k-n*h))/(sqrt((l-k)^2+(h-l)^2+(k-h)^2)*sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2))=a'
eq2='(u*(l-k)+v*(h-l)+w*(k-h))/sqrt((l-k)^2+(h-l)^2+(k-h)^2)=b'
eq3='(u*(n*l-p*k)+v*(p*h-m*l)+w*(m*k-n*h))/sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2)=c'
eq4='u*h+v*k+w*l=0'
eq5='u^2+v^2+w^2=1'
eq6='h^2+k^2+l^2=1'
eq7='m=-1'
eq8='n=1'
eq9='p=1'
eq10='a=cos(50*3.1415926/180)'
eq11='b=cos(55*3.1415926/180)'
eq12='c=cos(4.9*3.1415926/180)'
其中eq7-eq12是定义a,b,c和m,n,p的值的,是我想自己待定的数值,此方程的变量是u,v,w和h,k,l六个.
最好能把编的程序和解都贴出来,这样我就可以直接用别的a,b,c和m,n,p了,

那位高手用MATLAB帮我解个六元四次方程组eq1='((l-k)*(n*l-p*k)+(h-l)*(p*h-m*l)+(k-h)*(m*k-n*h))/(sqrt((l-k)^2+(h-l)^2+(k-h)^2)*sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2))=a'eq2='(u*(l-k)+v*(h-l)+w*(k-h))/sqrt((l-k)^2+(h-l)^2+(k-h)^2)=b'eq
你的问题可以化为下面向量的问题
已知a=(1,1,1),b=(-1,1,1),c=a×u,d=c×u,
c和d的夹角是50°,c和v的夹角是55°,d和v的夹角是4.9°,
u⊥v,|u|=1,|v|=1
求u,v
题中的a,b,c,d,u,v均为三维向量,×表示向量内积,|u|表示向量u的模
其中,向量b对应你以前的（m,n,p）,向量u对应你以前的（h,k,l）,向量v对应你以前的(u,v,w)
由上题,c=a×u,d=c×u可得c⊥u,d⊥u又u⊥v,且c,d,v有相同的起点即坐标原点,从而c,d,v在同一平面上且有相同的起点,且均与u垂直
所以c,d,v之间的夹角必定满足某个等式,回到题上也就是说,55°=50°+5°,
进一步说,你给的条件是矛盾的,所以matlab找不到解
就算你给出的条件是对的,由于你给出的前三个方程并非完全独立的,也不足以确定你想要的结果

clc
clear all
m=-1;
n=1;
p=1;
a=cos(50*pi/180);
b=cos(55*pi/180);
c=cos(4.9*pi/180);
eq1=sym('((l-k)*(n*l-p*k)+(h-l)*(p*h-m*l)+(k-h)*(m*k-n*h))/(sqrt((l-k)^2+(h-l)^2+(k...

全部展开

clc
clear all
m=-1;
n=1;
p=1;
a=cos(50*pi/180);
b=cos(55*pi/180);
c=cos(4.9*pi/180);
eq1=sym('((l-k)*(n*l-p*k)+(h-l)*(p*h-m*l)+(k-h)*(m*k-n*h))/(sqrt((l-k)^2+(h-l)^2+(k-h)^2)*sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2))-a');
eq2=sym('(u*(l-k)+v*(h-l)+w*(k-h))/sqrt((l-k)^2+(h-l)^2+(k-h)^2)-b');
eq3=sym('(u*(n*l-p*k)+v*(p*h-m*l)+w*(m*k-n*h))/sqrt((n*l-p*k)^2+(p*h-m*l)^2+(m*k-n*h)^2)-c');
eq4=sym('u*h+v*k+w*l');
eq5=sym('u^2+v^2+w^2-1');
eq6=sym('h^2+k^2+l^2-1' );
[u v w h k l]=solve(eq1 ,eq2, eq3, eq4, eq5, eq6)
结果
Warning: Explicit solution could not be found.
> In solve at 140
In sym.solve at 49
In j at 15

u =

[ empty sym ]


v =
[]
w =
[]
h =
[]
k =
[]
l =
[]

收起