那位高手用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
那位高手用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了,
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 =
[]
你的问题可以化为下面向量的问题已知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...