Skip to content
Snippets Groups Projects
Commit fd812b03 authored by jcandrsn's avatar jcandrsn
Browse files

Corrections to linode boundary conditions to correctly apply boundary...

Corrections to linode boundary conditions to correctly apply boundary condition to segment period (forgot to divide by frequency) and to update Jacobian to correct initial size (ran anyway due to automatically expanding it, but needed to be corrected for clarity).
parent 12a5e7cf
No related branches found
No related tags found
1 merge request!1Development
......@@ -7,7 +7,7 @@ x1 = x(1,:);
p2 = p(2,:);
p3 = p(3,:);
J = zeros(2,2,numel(x1));
J = zeros(2,3,numel(x1));
J(2,1,:) = -x1;
J(2,2,:) = -sin(p3.*t + p2);
J(2,3,:) = -t.*sin(p3.*t + p2);
......
......@@ -5,8 +5,9 @@ function J = linode_het_DFDT(t, x, p)
x1 = x(1,:);
phi = p(2,:);
om = p(3,:);
J = zeros(2,numel(x1));
J(2,:) = -sin(t + phi);
J(2,:) = -sin(om.*t + phi);
end
......@@ -7,7 +7,8 @@ x0 = u(1:2);
x1 = u(3:4);
T0 = u(5);
T = u(6);
p3 = u(7);
y = [x1(1:2)-x0(1:2); x0(2); T0; T-2*pi];
y = [x1(1:2)-x0(1:2); x0(2); T0; T-2*pi/p3];
end
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment