您所在位置：主题：[讨论]求教 解二元三次方程
apcz
[专家分：0] 发布于 2007-12-11 23:04:00
我写的计算步骤如下:(一开如可以很快运行,但试了几次就运行不了了,不知为什么,大家帮帮我!)
>>syms x y x1 y1 x2 y2 x3 y3 x4 y4 x5 y5 x6 y6 x7 y7 x8 y8 a b
>>N1=1/4*(1-a)*(1-b)*(-a-b-1);
>>N2=1/4*(1+a)*(1-b)*(a-b-1);
>>N3=1/4*(1+a)*(1+b)*(a+b-1);
>>N4=1/4*(1-a)*(1+b)*(-a+b-1);
>>N5=1/2*(1-a^2)*(1-b);
>>N6=1/2*(1-b^2)*(1+a);
>>N7=1/2*(1-a^2)*(1+b);
>>N8=1/2*(1-b^2)*(1-a);
>>eq1=N1*x1+N2*x2+N3*x3+N4*x4+N5*x5+N6*x6+N7*x7+N8*x8-x;
>>eq2=N1*y1+N2*y2+N3*y3+N4*y4+N5*y5+N6*y6+N7*y7+N8*y8-y;
>>s=solve(eq1,eq2,a,b)
>>syms x y x1 y1 x2 y2 x3 y3 x4 y4 x5 y5 x6 y6 x7 y7 x8 y8 a b
>>N1=1/4*(1-a)*(1-b)*(-a-b-1);
>>N2=1/4*(1+a)*(1-b)*(a-b-1);
>>N3=1/4*(1+a)*(1+b)*(a+b-1);
>>N4=1/4*(1-a)*(1+b)*(-a+b-1);
>>N5=1/2*(1-a^2)*(1-b);
>>N6=1/2*(1-b^2)*(1+a);
>>N7=1/2*(1-a^2)*(1+b);
>>N8=1/2*(1-b^2)*(1-a);
>>eq1=N1*x1+N2*x2+N3*x3+N4*x4+N5*x5+N6*x6+N7*x7+N8*x8-x;
>>eq2=N1*y1+N2*y2+N3*y3+N4*y4+N5*y5+N6*y6+N7*y7+N8*y8-y;
>>s=solve(eq1,eq2,a,b)