For an idealized pendulum with massless arm of length L, traveling along a vertical circular track of radius R with constant angular velocity W, this program will integrate the equations of motion given the above parameters and initial conditions, and will then animate the system, tracing the pendulum's path.
Good demo for introductory dynamics or numerical integration.
[alpha,delta] = eclip2equat(lambda,beta) returns the right ascension and declination in the equatorial coordinate system corresponding to the longitudinal and latitudinal coordinates in the ecliptic coordinate system (lambda and beta).
[t,z,r_PO] = wattLinkage(n,npathpoints,lengths,omega) integrates the equations of motion for a 3 bar Watt linkage and animates the motion of the system over n cycles (one cycle is defined as a full sweep top to bottom).
[coordSys,as] = animEuler(h,rotSet,angs) animates the rotation of a dextral coordinate system using three arbitrary Euler angles (angs) about three axes (rotSet) in figure (h). rotSet is an array of 3 numbers in the range of [1, 3] representing...
SAMPLEDIST Sample from an arbitrary distribution
sampleDist(f,M,N,b) retruns an array of size X of random values sampled from the distribution defined by the probability density function referred to by handle f, over the range b = [min,...
vec2orbElem(rs,vs,mus) converts positions (rs) and velocities (vs)
of bodies with gravitational parameters (mus) to Keplerian orbital elements.
rs: 3n x 1 stacked initial position vectors:
This simulation finds the angular velocity time histories required
to raise a glass from ground level to table height in a straight 45 degree trajectory. Build Robot constructs the robot from inputs, and Go! starts the simulation.
findEulerAngs generates Euler angle sets and animates rotations in response to user manipulation of a 3D rigid body (a box).
Use the mouse to change the orientation of the box. A wireframe of the original box position will appear. The...
[U D] = UDFactor(P) returns matrices U and D such that U.'*D*U = P
[U D] = UDFactor(P,uflag) returns matrices U and D such that U*D*U.' = P
when uflag is set to TRUE. Setting uflag to FALSE is equivalent to
running UDFactor with...
Phi = keplerSTM(x0,dt,mu) will return the state transition matrix Phi
for a set of objects described by initial positions and velocities (x0)
with gravitational parameters mu, propogated along Keplerian orbits for
time dt. The...