Re: Optimizing using Fsolve (more equations than u
- From: Pascal <mastersofwar@xxxxxxxxx>
- Date: Tue, 28 Feb 2006 12:19:59 -0500
Function equations got cut off. Sorry for long thread.
The loop inside the optimization functions should be:
F(:,ctr)=
[m(1)*X0 + m(2)*Y0 + m(3)*Z0 + m(4) - x0*(m(9)*X0 + m(10)*Y0 +
m(11)*Z0 + m(12));
m(1)*X1 + m(2)*Y1 + m(3)*Z1 + m(4) - x1*(m(9)*X1 + m(10)*Y1 +
m(11)*Z1 + m(12));
m(1)*X2 + m(2)*Y2 + m(3)*Z2 + m(4) - x2*(m(9)*X2 + m(10)*Y2 +
m(11)*Z2 + m(12));
m(1)*X3 + m(2)*Y3 + m(3)*Z3 + m(4) - x3*(m(9)*X3 + m(10)*Y3 +
m(11)*Z3 + m(12));
m(1)*X4 + m(2)*Y4 + m(3)*Z4 + m(4) - x4*(m(9)*X4 + m(10)*Y4 +
m(11)*Z4 + m(12));
m(1)*X5 + m(2)*Y5 + m(3)*Z5 + m(4) - x5*(m(9)*X5 + m(10)*Y5 +
m(11)*Z5 + m(12));
m(1)*X6 + m(2)*Y6 + m(3)*Z6 + m(4) - x6*(m(9)*X6 + m(10)*Y6 +
m(11)*Z6 + m(12));
m(1)*X7 + m(2)*Y7 + m(3)*Z7 + m(4) - x7*(m(9)*X7 + m(10)*Y7 +
m(11)*Z7 + m(12));
m(1)*X8 + m(2)*Y8 + m(3)*Z8 + m(4) - x8*(m(9)*X8 + m(10)*Y8 +
m(11)*Z8 + m(12));
m(5)*X0 + m(6)*Y0 + m(7)*Z0 + m(8) - y0*(m(9)*X0 + m(10)*Y0 +
m(11)*Z0 + m(12));
m(5)*X1 + m(6)*Y1 + m(7)*Z1 + m(8) - y1*(m(9)*X1 + m(10)*Y1 +
m(11)*Z1 + m(12));
m(5)*X2 + m(6)*Y2 + m(7)*Z2 + m(8) - y2*(m(9)*X2 + m(10)*Y2 +
m(11)*Z2 + m(12));
m(5)*X3 + m(6)*Y3 + m(7)*Z3 + m(8) - y3*(m(9)*X3 + m(10)*Y3 +
m(11)*Z3 + m(12));
m(5)*X4 + m(6)*Y4 + m(7)*Z4 + m(8) - y4*(m(9)*X4 + m(10)*Y4 +
m(11)*Z4 + m(12));
m(5)*X5 + m(6)*Y5 + m(7)*Z5 + m(8) - y5*(m(9)*X5 + m(10)*Y5 +
m(11)*Z5 + m(12));
m(5)*X6 + m(6)*Y6 + m(7)*Z6 + m(8) - y6*(m(9)*X6 + m(10)*Y6 +
m(11)*Z6 + m(12));
m(5)*X7 + m(6)*Y7 + m(7)*Z7 + m(8) - y7*(m(9)*X7 + m(10)*Y7 +
m(11)*Z7 + m(12));
m(5)*X8 + m(6)*Y8 + m(7)*Z8 + m(8) - y8*(m(9)*X8 + m(10)*Y8 +
m(11)*Z8 + m(12));
((m(1)*X1 + m(2)*Y1 + m(3)*Z1 + m(4))/((m(9)*X1 + m(10)*Y1 + m(11)*Z1
+ m(12)))) - ((m(1)*X0 + m(2)*Y0 + m(3)*Z0 + m(4))/((m(9)*X0 +
m(10)*Y0 + m(11)*Z0 + m(12)))) - (x1-x0);
((m(1)*X2 + m(2)*Y2 + m(3)*Z2 + m(4))/((m(9)*X2 + m(10)*Y2 + m(11)*Z2
+ m(12)))) - ((m(1)*X1 + m(2)*Y1 + m(3)*Z1 + m(4))/((m(9)*X1 +
m(10)*Y1 + m(11)*Z1 + m(12)))) - (x2-x1);
((m(1)*X3 + m(2)*Y3 + m(3)*Z3 + m(4))/((m(9)*X3 + m(10)*Y3 + m(11)*Z3
+ m(12)))) - ((m(1)*X2 + m(2)*Y2 + m(3)*Z2 + m(4))/((m(9)*X2 +
m(10)*Y2 + m(11)*Z2 + m(12)))) - (x3-x2);
((m(1)*X4 + m(2)*Y4 + m(3)*Z4 + m(4))/((m(9)*X4 + m(10)*Y4 + m(11)*Z4
+ m(12)))) - ((m(1)*X3 + m(2)*Y3 + m(3)*Z3 + m(4))/((m(9)*X3 +
m(10)*Y3 + m(11)*Z3 + m(12)))) - (x4-x3);
((m(1)*X5 + m(2)*Y5 + m(3)*Z5 + m(4))/((m(9)*X5 + m(10)*Y5 + m(11)*Z5
+ m(12)))) - ((m(1)*X4 + m(2)*Y4 + m(3)*Z4 + m(4))/((m(9)*X4 +
m(10)*Y4 + m(11)*Z4 + m(12)))) - (x5-x4);
((m(1)*X6 + m(2)*Y6 + m(3)*Z6 + m(4))/((m(9)*X6 + m(10)*Y6 + m(11)*Z6
+ m(12)))) - ((m(1)*X5 + m(2)*Y5 + m(3)*Z5 + m(4))/((m(9)*X5 +
m(10)*Y5 + m(11)*Z5 + m(12)))) - (x6-x5);
((m(1)*X7 + m(2)*Y7 + m(3)*Z7 + m(4))/((m(9)*X7 + m(10)*Y7 + m(11)*Z7
+ m(12)))) - ((m(1)*X6 + m(2)*Y6 + m(3)*Z6 + m(4))/((m(9)*X6 +
m(10)*Y6 + m(11)*Z6 + m(12)))) - (x7-x6);
((m(1)*X8 + m(2)*Y8 + m(3)*Z8 + m(4))/((m(9)*X8 + m(10)*Y8 + m(11)*Z8
+ m(12)))) - ((m(1)*X7 + m(2)*Y7 + m(3)*Z7 + m(4))/((m(9)*X7 +
m(10)*Y7 + m(11)*Z7 + m(12)))) - (x8-x7);
((m(5)*X1 + m(6)*Y1 + m(7)*Z1 + m(8))/((m(9)*X1 + m(10)*Y1 + m(11)*Z1
+ m(12)))) - ((m(5)*X0 + m(6)*Y0 + m(7)*Z0 + m(8))/((m(9)*X0 +
m(10)*Y0 + m(11)*Z0 + m(12)))) - (y1-y0);
((m(5)*X2 + m(6)*Y2 + m(7)*Z2 + m(8))/((m(9)*X2 + m(10)*Y2 + m(11)*Z2
+ m(12)))) - ((m(5)*X1 + m(6)*Y1 + m(7)*Z1 + m(8))/((m(9)*X1 +
m(10)*Y1 + m(11)*Z1 + m(12)))) - (y2-y1);
((m(5)*X3 + m(6)*Y3 + m(7)*Z3 + m(8))/((m(9)*X3 + m(10)*Y3 + m(11)*Z3
+ m(12)))) - ((m(5)*X2 + m(6)*Y2 + m(7)*Z2 + m(8))/((m(9)*X2 +
m(10)*Y2 + m(11)*Z2 + m(12)))) - (y3-y2);
((m(5)*X4 + m(6)*Y4 + m(7)*Z4 + m(8))/((m(9)*X4 + m(10)*Y4 + m(11)*Z4
+ m(12)))) - ((m(5)*X3 + m(6)*Y3 + m(7)*Z3 + m(8))/((m(9)*X3 +
m(10)*Y3 + m(11)*Z3 + m(12)))) - (y4-y3);
((m(5)*X5 + m(6)*Y5 + m(7)*Z5 + m(8))/((m(9)*X5 + m(10)*Y5 + m(11)*Z5
+ m(12)))) - ((m(5)*X4 + m(6)*Y4 + m(7)*Z4 + m(8))/((m(9)*X4 +
m(10)*Y4 + m(11)*Z4 + m(12)))) - (y5-y4);
((m(5)*X6 + m(6)*Y6 + m(7)*Z6 + m(8))/((m(9)*X6 + m(10)*Y6 + m(11)*Z6
+ m(12)))) - ((m(5)*X5 + m(6)*Y5 + m(7)*Z5 + m(8))/((m(9)*X5 +
m(10)*Y5 + m(11)*Z5 + m(12)))) - (y6-y5);
((m(5)*X7 + m(6)*Y7 + m(7)*Z7 + m(8))/((m(9)*X7 + m(10)*Y7 + m(11)*Z7
+ m(12)))) - ((m(5)*X6 + m(6)*Y6 + m(7)*Z6 + m(8))/((m(9)*X6 +
m(10)*Y6 + m(11)*Z6 + m(12)))) - (y7-y6);
((m(5)*X8 + m(6)*Y8 + m(7)*Z8 + m(8))/((m(9)*X8 + m(10)*Y8 + m(11)*Z8
+ m(12)))) - ((m(5)*X7 + m(6)*Y7 + m(7)*Z7 + m(8))/((m(9)*X7 +
m(10)*Y7 + m(11)*Z7 + m(12)))) - (y8-y7)];
end
Is this written correctly for the N points analysis?
.
- References:
- Prev by Date: Optimizing using Fsolve (more equations than unknowns)
- Next by Date: Re: How can I remove the padding??
- Previous by thread: Optimizing using Fsolve (more equations than unknowns)
- Next by thread: For loops
- Index(es):
