📄 funmatinmmdpel.m
字号:
0, 0, 0, 0, 0, 0, -d, c*x1+d*y1, -c;
0, -x1, 1, 0, 0, 0, 0, x1*x2, -x2;
1, -y1, 0, 0, 0, 0, -x2, y1*x2, 0;
0, 0, 0, 0, -x1, 1, 0, x1*y2, -y2;
0, 0, 0, 1, -y1, 0, -y2, y1*y2, 0];
case 24,
Jact= [ 0, -c, 0, 0, -a, 0, 0, a*x2+c*y2, 0;
0, -d, 0, 0, -b, 0, 0, b*x2+d*y2, 0;
0, 0, 0, 0, 0, 0, -b, a*x1+b*y1, -a;
0, 0, 0, 0, 0, 0, -d, c*x1+d*y1, -c;
0, 0, 0, 0, -x1, 1, 0, x1*x2, -x2;
0, 0, 0, 1, -y1, 0, -x2, y1*x2, 0;
0, -x1, 1, 0, 0, 0, 0, x1*y2, -y2;
1, -y1, 0, 0, 0, 0, -y2, y1*y2, 0];
case 25,
Jact= [ 0, 0, a*x2+c*y2, 0, 0, -a, 0, 0, -c;
0, 0, b*x2+d*y2, 0, 0, -b, 0, 0, -d;
-a, -b, a*x1+b*y1, 0, 0, 0, 0, 0, 0;
-c, -d, c*x1+d*y1, 0, 0, 0, 0, 0, 0;
-x2, 0, x1*x2, 1, 0, -x1, 0, 0, 0;
0, -x2, y1*x2, 0, 1, -y1, 0, 0, 0;
-y2, 0, x1*y2, 0, 0, 0, 1, 0, -x1;
0, -y2, y1*y2, 0, 0, 0, 0, 1, -y1];
case 26,
Jact= [ 0, 0, a*x2+c*y2, 0, 0, -c, 0, 0, -a;
0, 0, b*x2+d*y2, 0, 0, -d, 0, 0, -b;
-a, -b, a*x1+b*y1, 0, 0, 0, 0, 0, 0;
-c, -d, c*x1+d*y1, 0, 0, 0, 0, 0, 0;
-x2, 0, x1*x2, 0, 0, 0, 1, 0, -x1;
0, -x2, y1*x2, 0, 0, 0, 0, 1, -y1;
-y2, 0, x1*y2, 1, 0, -x1, 0, 0, 0;
0, -y2, y1*y2, 0, 1, -y1, 0, 0, 0];
case 27,
Jact= [ 0, 0, -a, 0, 0, a*x2+c*y2, 0, 0, -c;
0, 0, -b, 0, 0, b*x2+d*y2, 0, 0, -d;
0, 0, 0, -a, -b, a*x1+b*y1, 0, 0, 0;
0, 0, 0, -c, -d, c*x1+d*y1, 0, 0, 0;
1, 0, -x1, -x2, 0, x1*x2, 0, 0, 0;
0, 1, -y1, 0, -x2, y1*x2, 0, 0, 0;
0, 0, 0, -y2, 0, x1*y2, 1, 0, -x1;
0, 0, 0, 0, -y2, y1*y2, 0, 1, -y1];
case 28,
Jact= [ 0, 0, -c, 0, 0, a*x2+c*y2, 0, 0, -a;
0, 0, -d, 0, 0, b*x2+d*y2, 0, 0, -b;
0, 0, 0, -a, -b, a*x1+b*y1, 0, 0, 0;
0, 0, 0, -c, -d, c*x1+d*y1, 0, 0, 0;
0, 0, 0, -x2, 0, x1*x2, 1, 0, -x1;
0, 0, 0, 0, -x2, y1*x2, 0, 1, -y1;
1, 0, -x1, -y2, 0, x1*y2, 0, 0, 0;
0, 1, -y1, 0, -y2, y1*y2, 0, 0, 0];
case 29,
Jact= [ 0, 0, -a, 0, 0, -c, 0, 0, a*x2+c*y2;
0, 0, -b, 0, 0, -d, 0, 0, b*x2+d*y2;
0, 0, 0, 0, 0, 0, -a, -b, a*x1+b*y1;
0, 0, 0, 0, 0, 0, -c, -d, c*x1+d*y1;
1, 0, -x1, 0, 0, 0, -x2, 0, x1*x2;
0, 1, -y1, 0, 0, 0, 0, -x2, y1*x2;
0, 0, 0, 1, 0, -x1, -y2, 0, x1*y2;
0, 0, 0, 0, 1, -y1, 0, -y2, y1*y2];
case 30,
Jact= [ 0, 0, -c, 0, 0, -a, 0, 0, a*x2+c*y2;
0, 0, -d, 0, 0, -b, 0, 0, b*x2+d*y2;
0, 0, 0, 0, 0, 0, -a, -b, a*x1+b*y1;
0, 0, 0, 0, 0, 0, -c, -d, c*x1+d*y1;
0, 0, 0, 1, 0, -x1, -x2, 0, x1*x2;
0, 0, 0, 0, 1, -y1, 0, -x2, y1*x2;
1, 0, -x1, 0, 0, 0, -y2, 0, x1*y2;
0, 1, -y1, 0, 0, 0, 0, -y2, y1*y2];
case 31,
Jact= [ 0, 0, a*x2+c*y2, 0, 0, -a, 0, 0, -c;
0, 0, b*x2+d*y2, 0, 0, -b, 0, 0, -d;
-b, -a, a*x1+b*y1, 0, 0, 0, 0, 0, 0;
-d, -c, c*x1+d*y1, 0, 0, 0, 0, 0, 0;
0, -x2, x1*x2, 0, 1, -x1, 0, 0, 0;
-x2, 0, y1*x2, 1, 0, -y1, 0, 0, 0;
0, -y2, x1*y2, 0, 0, 0, 0, 1, -x1;
-y2, 0, y1*y2, 0, 0, 0, 1, 0, -y1];
case 32,
Jact= [ 0, 0, a*x2+c*y2, 0, 0, -c, 0, 0, -a;
0, 0, b*x2+d*y2, 0, 0, -d, 0, 0, -b;
-b, -a, a*x1+b*y1, 0, 0, 0, 0, 0, 0;
-d, -c, c*x1+d*y1, 0, 0, 0, 0, 0, 0;
0, -x2, x1*x2, 0, 0, 0, 0, 1, -x1;
-x2, 0, y1*x2, 0, 0, 0, 1, 0, -y1;
0, -y2, x1*y2, 0, 1, -x1, 0, 0, 0;
-y2, 0, y1*y2, 1, 0, -y1, 0, 0, 0];
case 33,
Jact= [ 0, 0, -a, 0, 0, a*x2+c*y2, 0, 0, -c;
0, 0, -b, 0, 0, b*x2+d*y2, 0, 0, -d;
0, 0, 0, -b, -a, a*x1+b*y1, 0, 0, 0;
0, 0, 0, -d, -c, c*x1+d*y1, 0, 0, 0;
0, 1, -x1, 0, -x2, x1*x2, 0, 0, 0;
1, 0, -y1, -x2, 0, y1*x2, 0, 0, 0;
0, 0, 0, 0, -y2, x1*y2, 0, 1, -x1;
0, 0, 0, -y2, 0, y1*y2, 1, 0, -y1];
case 34,
Jact= [ 0, 0, -c, 0, 0, a*x2+c*y2, 0, 0, -a;
0, 0, -d, 0, 0, b*x2+d*y2, 0, 0, -b;
0, 0, 0, -b, -a, a*x1+b*y1, 0, 0, 0;
0, 0, 0, -d, -c, c*x1+d*y1, 0, 0, 0;
0, 0, 0, 0, -x2, x1*x2, 0, 1, -x1;
0, 0, 0, -x2, 0, y1*x2, 1, 0, -y1;
0, 1, -x1, 0, -y2, x1*y2, 0, 0, 0;
1, 0, -y1, -y2, 0, y1*y2, 0, 0, 0];
case 35,
Jact= [ 0, 0, -a, 0, 0, -c, 0, 0, a*x2+c*y2;
0, 0, -b, 0, 0, -d, 0, 0, b*x2+d*y2;
0, 0, 0, 0, 0, 0, -b, -a, a*x1+b*y1;
0, 0, 0, 0, 0, 0, -d, -c, c*x1+d*y1;
0, 1, -x1, 0, 0, 0, 0, -x2, x1*x2;
1, 0, -y1, 0, 0, 0, -x2, 0, y1*x2;
0, 0, 0, 0, 1, -x1, 0, -y2, x1*y2;
0, 0, 0, 1, 0, -y1, -y2, 0, y1*y2];
case 36,
Jact= [ 0, 0, -c, 0, 0, -a, 0, 0, a*x2+c*y2;
0, 0, -d, 0, 0, -b, 0, 0, b*x2+d*y2;
0, 0, 0, 0, 0, 0, -b, -a, a*x1+b*y1;
0, 0, 0, 0, 0, 0, -d, -c, c*x1+d*y1;
0, 0, 0, 0, 1, -x1, 0, -x2, x1*x2;
0, 0, 0, 1, 0, -y1, -x2, 0, y1*x2;
0, 1, -x1, 0, 0, 0, 0, -y2, x1*y2;
1, 0, -y1, 0, 0, 0, -y2, 0, y1*y2];
case 37,
disp('Error')
end
vectJ=0;
for jc=1:9,
vectJ(jc,comb)=det([Jact(:,1:(jc-1)) Jact(:,(jc+1):9)]);
end
nor=[nor norm(vectJ(:,comb))];
% nor2=[nor2 (a*d-b*c)^2*sqrt(x1^2+y1^2+1)*sqrt(x2^2+y2^2+1)];
comb=comb+1;
end;
end;
[vm,im]=max(nor);
if vmini<vm,
in1=fix((im-1)/6)+1;
in2=mod(im,6);
if in2==0,
in2=6;
end
j0= indexs(in1,1);
j1= indexs(in1,2);
j2= indexs(in1,3);
i0= indexs(in2,1);
i1= indexs(in2,2);
i2= indexs(in2,3);
if imag(sum(val1(:,1)))==0,
vec1=vect1(:,1)./vect1(j0,1);
else
vec1=vect1(:,3)./vect1(j0,3);
end
x1=vec1(j1);
y1=vec1(j2);
if imag(sum(val2(:,1)))==0,
vec2=vect2(:,1)./vect2(i0,1);
else
vec2=vect2(:,3)./vect2(i0,3);
end
x2=vec2(i1);
y2=vec2(i2);
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -