Click here to Skip to main content
15,895,709 members
Articles / Programming Languages / C#

Grandiose Projects 3. Compatibility with Simulink

Rate me:
Please Sign up or sign in to vote.
4.27/5 (11 votes)
8 Feb 2010CPOL23 min read 47.8K   5.9K   38  
Import of Simulink files
using System;
using System.Collections.Generic;
using System.Text;
using System.Runtime.Serialization;

using Motion6D.Intrefaces;

namespace Motion6D
{
    public class Motion6DFrame : RotatedFrame, IVelocity
    {

        #region Fields
        public static readonly Motion6DFrame Base = new Motion6DAcceleratedFrame();

        double[] velocity = new double[] { 0, 0, 0 };
        double[] hv = new double[3];

        //protected double[] relativeVelocity = new double[] { 0, 0, 0 };


        protected double[] der = new double[4];
        protected double[] qd = new double[4];

        #endregion

        #region Ctor

        #endregion

        #region IVelocity Members

        double[] IVelocity.Velocity
        {
            get { return velocity; }
        }

       /* double[] IVelocity.RevativeVelocity
        {
            get
            {
                return relativeVelocity;
            }
        }*/


        #endregion

        #region Overriden Members

        

        public override void Set(ReferenceFrame baseFrame, ReferenceFrame relative)
        {
            base.Set(baseFrame, relative);
            IOrientation bo = baseFrame as IOrientation;
            IOrientation ro = relative as IOrientation;
            IVelocity bv = baseFrame as IVelocity;
            IVelocity rv = relative as IVelocity;
            IAngularVelocity oa = baseFrame as IAngularVelocity;
            IAngularVelocity ra = relative as IAngularVelocity;
            double[] vb = bv.Velocity;
            double[] vr = rv.Velocity;
            double[,] mb = bo.Matrix;
            double[,] m = ro.Matrix;
            double[] om = oa.Omega;
            double[] pos = relative.Position;
            Vector3D.V3DOperations.VectorPoduct(om, pos, hv);
            for (int i = 0; i < 3; i++)
            {
                velocity[i] = vb[i];
                for (int j = 0; j < 3; j++)
                {
                    velocity[i] += mb[i, j] * (vr[j] + hv[j]);
                }
            }
        }

        #endregion

    }
}

By viewing downloads associated with this article you agree to the Terms of Service and the article's licence.

If a file you wish to view isn't highlighted, and is a text file (not binary), please let us know and we'll add colourisation support for it.

License

This article, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)


Written By
Architect
Russian Federation Russian Federation
Ph. D. Petr Ivankov worked as scientific researcher at Russian Mission Control Centre since 1978 up to 2000. Now he is engaged by Aviation training simulators http://dinamika-avia.com/ . His additional interests are:

1) Noncommutative geometry

http://front.math.ucdavis.edu/author/P.Ivankov

2) Literary work (Russian only)

http://zhurnal.lib.ru/editors/3/3d_m/

3) Scientific articles
http://arxiv.org/find/all/1/au:+Ivankov_Petr/0/1/0/all/0/1

Comments and Discussions