Download Reference Manual

121
⎛x ⎞
⎛ cos θ (α ) − sen θ (α )⎞ ⎡
⎟ ⎢(1 + α )⎜ 1 ⎟ + x12 + x 22
⎝ sen θ (α ) cos θ (α ) ⎠ ⎣
⎝ x2 ⎠
(
f(x, α) = ⎜
J=
)⎛⎜⎝ ab ((αα))
− b (α )⎞ ⎛ x1 ⎞ ⎤
⎟⎜ ⎟
a (α ) ⎠ ⎝ x 2 ⎠ ⎥⎦
d
f(x, α)
dx
Y se presenta del siguiente modo en el archivo SACKERJD.M,
function [J,ValPr]=SackerJD(y)
x1=y(1); x2=y(2); alfa=y(3);
teta=pi/2; a=-alfa; b=0;
dAdx1=(1+alfa)+
(x1^2+x2^2)*(a)+ (2*x1)*(a*x1-b*x2);
dAdx2=
0
+
(x1^2+x2^2)*(-b)+ (2*x2)*(a*x1-b*x2);
dBdx1=
0
+
(x1^2+x2^2)*(b)+ (2*x1)*(b*x1+a*x2);
dBdx2=(1+alfa)+
(x1^2+x2^2)*(a)+ (2*x1)*(b*x1+a*x2);
J=[cos(teta)*dAdx1-sin(teta)*dAdx1, cos(teta)*dAdx2-sin(teta)*dAdx2;
sin(teta)*dBdx1+cos(teta)*dBdx1,cos(teta)*dBdx2-sin(teta)*dBdx2];
ValPr=eig(J);
7.6.5 Ver también
FUN, GRADFUN, FUNP, GRADFUNP, FUND
7.6.6 Referencias
Capítulo 4