Developing simulations for the Missile Defense Agency - Advanced Technology group

30+ years experience in missiles, smart munitions and other aerodynamic systems. Specializing in systems engineering, aerospace systems design, simulation and modeling, aerodynamics, and numerical analysis.

Answered

Kalman Filtering for smoothing GPS speed data

The problem is with this line: P = (1-KG)*P; This causes the gain P to continuously decrease. After a while, the gain P is ve...

Kalman Filtering for smoothing GPS speed data

The problem is with this line: P = (1-KG)*P; This causes the gain P to continuously decrease. After a while, the gain P is ve...

17 días ago | 1

| accepted

Answered

numerical derivate and central difference

Here are some Numerical difference formulas using the central difference method (substitute dt = h). You can find more formulas...

numerical derivate and central difference

Here are some Numerical difference formulas using the central difference method (substitute dt = h). You can find more formulas...

alrededor de 2 meses ago | 1

Answered

Meaning of the following: function [y1,...yn] = input()

It simply means that there are no user inputs to this function. The long list of variable names are returned from the function....

Meaning of the following: function [y1,...yn] = input()

It simply means that there are no user inputs to this function. The long list of variable names are returned from the function....

3 meses ago | 0

| accepted

Answered

Angle Between two vectors.

Given two vectors A and B, the dot product of the two vectors (A dot B) gives the product ABcos(ang), so to get just the angle, ...

Angle Between two vectors.

Given two vectors A and B, the dot product of the two vectors (A dot B) gives the product ABcos(ang), so to get just the angle, ...

3 meses ago | 0

| accepted

Answered

Is it possible to summarize this code?

You could use a switch case structure, that would be a bit more compact. it would look like this: switch davar1.CC_CHECK ...

Is it possible to summarize this code?

You could use a switch case structure, that would be a bit more compact. it would look like this: switch davar1.CC_CHECK ...

3 meses ago | 2

| accepted

Answered

help fixing the Heuns method

Heun's method is otherwise known as "explicit trapezoidal method", "improved Euler's method" or "modified Euler's method". Th...

help fixing the Heuns method

Heun's method is otherwise known as "explicit trapezoidal method", "improved Euler's method" or "modified Euler's method". Th...

3 meses ago | 0

Answered

I need help coding a human cannonball ignoring air resistance

Assuming zero drag and relatively low muzzle velocity ( acceleration due to gravity = constant), the point-mass trajectory can b...

I need help coding a human cannonball ignoring air resistance

Assuming zero drag and relatively low muzzle velocity ( acceleration due to gravity = constant), the point-mass trajectory can b...

4 meses ago | 1

Answered

Runge Kutta 4th Order Method

I think the problem is the T argument. in k1 you have a capital T, which is a 3-vector in k2, k3, and k4 it is a lower case t,...

Runge Kutta 4th Order Method

I think the problem is the T argument. in k1 you have a capital T, which is a 3-vector in k2, k3, and k4 it is a lower case t,...

4 meses ago | 0

Answered

Matlab Homework that im struggling with

To prompt a user to enter a value, try: x = input('Enter the number of equations (from 2 to 5): \n')

Matlab Homework that im struggling with

To prompt a user to enter a value, try: x = input('Enter the number of equations (from 2 to 5): \n')

4 meses ago | 0

Answered

I need help with this script

Try this: T = input('Enter your temperature') if T<=96 disp('You are either hypothermic or a zombie') elseif T<99 d...

I need help with this script

Try this: T = input('Enter your temperature') if T<=96 disp('You are either hypothermic or a zombie') elseif T<99 d...

4 meses ago | 0

Answered

Ballstic Impact of Bullets on armor

I don't know af any specific ballistics tools. You should browse the a list of Matlab products at this link: Matlab Products

Ballstic Impact of Bullets on armor

I don't know af any specific ballistics tools. You should browse the a list of Matlab products at this link: Matlab Products

5 meses ago | 0

| accepted

Answered

How to find theta and theta dot in a simple pendulum simulation in MATLAB.

I'm not sure at all what your system equations are. The first rule of dynamics is where you have translational modes of motion ...

How to find theta and theta dot in a simple pendulum simulation in MATLAB.

I'm not sure at all what your system equations are. The first rule of dynamics is where you have translational modes of motion ...

5 meses ago | 0

Answered

What is a "system variable" in a Simulink model?

In this context, the "system" is the device that you are modelling. Therefore, a "system variable" is one of the variables in y...

What is a "system variable" in a Simulink model?

In this context, the "system" is the device that you are modelling. Therefore, a "system variable" is one of the variables in y...

5 meses ago | 1

| accepted

Answered

Undefined function 'simset' for input arguments of type 'char'.

simset is an obsolete function. See documentation here. It says to use the Configuration Parameters dialog box.

Undefined function 'simset' for input arguments of type 'char'.

simset is an obsolete function. See documentation here. It says to use the Configuration Parameters dialog box.

5 meses ago | 0

Answered

Newtons method for finding minimum of a function.

Here is another variation; function [f,g,h] = myfunction(x) f1 = @(z) exp(2*sin(z))-z; %function dx ...

Newtons method for finding minimum of a function.

Here is another variation; function [f,g,h] = myfunction(x) f1 = @(z) exp(2*sin(z))-z; %function dx ...

5 meses ago | 0

Answered

help me in trying to solve the secant method

Try preallocating x, right after you define n: n = 100; x = nan(1,n); ...

help me in trying to solve the secant method

Try preallocating x, right after you define n: n = 100; x = nan(1,n); ...

5 meses ago | 0

Answered

trying to calculate the central diff approximation

There is a Matlab function "diff" which will do this. Otherwise, your subscripts must all match (2:end-1) Adi_1 = diff(Vdi_1...

trying to calculate the central diff approximation

There is a Matlab function "diff" which will do this. Otherwise, your subscripts must all match (2:end-1) Adi_1 = diff(Vdi_1...

5 meses ago | 0

| accepted

Answered

Performing calculations (Newton's method) over values in loop

Yes, your loop is not quite right. Notice your placement of the three lines i=1; converge = false; po = p_init(1); These t...

Performing calculations (Newton's method) over values in loop

Yes, your loop is not quite right. Notice your placement of the three lines i=1; converge = false; po = p_init(1); These t...

5 meses ago | 0

| accepted

Answered

Draw a circle for arbitrary orientation on spherical surface

My approach is to define the circle around the X-axis, then rotate it into the desired position through a Yaw-Pitch rotation. ...

Draw a circle for arbitrary orientation on spherical surface

My approach is to define the circle around the X-axis, then rotate it into the desired position through a Yaw-Pitch rotation. ...

5 meses ago | 1

| accepted

Answered

the meaning of basic matlab words

an argument is a constant or variable that you provide to a function; e.g. sin(x) % here x is the argument to the function sin...

the meaning of basic matlab words

an argument is a constant or variable that you provide to a function; e.g. sin(x) % here x is the argument to the function sin...

5 meses ago | 0

| accepted

Answered

How to alter Comet Speed

I don't see a way to control the speed of the comet function. There is a discussion of controlling the speed of an animated lin...

How to alter Comet Speed

I don't see a way to control the speed of the comet function. There is a discussion of controlling the speed of an animated lin...

5 meses ago | 0

Answered

Incorrect Answer Calculated For Force Calculation

It is assumed that you want to find F3 and y such that the system is in equilibrium (this is not stated) (and I disagree with D...

Incorrect Answer Calculated For Force Calculation

It is assumed that you want to find F3 and y such that the system is in equilibrium (this is not stated) (and I disagree with D...

5 meses ago | 0

Answered

Transformation Matrix Between Two Cartesian Systems

This is a very generic question and there are many ways to answer. It depends on what information you have about the two coordi...

Transformation Matrix Between Two Cartesian Systems

This is a very generic question and there are many ways to answer. It depends on what information you have about the two coordi...

5 meses ago | 0

Answered

Why is ode45 changing the value of the state variable between calls to my dynamics function?

ode45 implements the Dormand-Prince method which is an embedded 4th / 5th order Runge-Kutta method. For every "step" in the sol...

Why is ode45 changing the value of the state variable between calls to my dynamics function?

ode45 implements the Dormand-Prince method which is an embedded 4th / 5th order Runge-Kutta method. For every "step" in the sol...

5 meses ago | 0

| accepted

Answered

Implement friction into a model, related to the velocity. (Missile Launch)

This question is familiar. look here.

Implement friction into a model, related to the velocity. (Missile Launch)

This question is familiar. look here.

5 meses ago | 0

Answered

How to make a marker circle bold?

plot(x, y, 'ro', 'linewidth',2);

How to make a marker circle bold?

plot(x, y, 'ro', 'linewidth',2);

5 meses ago | 0

| accepted

Answered

how can i resolve this equation Runge kutta method

Ths is the equation for a parabola. y'' = -1 y' = -y y = -1/2 y^2 No need for a numerical approximation.

how can i resolve this equation Runge kutta method

Ths is the equation for a parabola. y'' = -1 y' = -y y = -1/2 y^2 No need for a numerical approximation.

5 meses ago | 0

Answered

Missile Lauch friction gamma factor - Output deviating from desired output, please help!

Step 1: Draw a free body diagram that shows the body with all of the forces that act on it and the coordinate frame(s) that are...

Missile Lauch friction gamma factor - Output deviating from desired output, please help!

Step 1: Draw a free body diagram that shows the body with all of the forces that act on it and the coordinate frame(s) that are...

6 meses ago | 0

| accepted

Answered

How to calculate the slope from geopraphical coordinates?

First of all, you should consider the relationship between the WGS84 lat/lon/alt and the ECEF frame. Altitude is perpendicular ...

How to calculate the slope from geopraphical coordinates?

First of all, you should consider the relationship between the WGS84 lat/lon/alt and the ECEF frame. Altitude is perpendicular ...

6 meses ago | 1

| accepted

Answered

calculate rotation matrix in 3D

I'll assume that we make the reference frame right-handed by negating the axis_3 vector such that axis_1 = [2.7239 -0.2014 -...

calculate rotation matrix in 3D

I'll assume that we make the reference frame right-handed by negating the axis_3 vector such that axis_1 = [2.7239 -0.2014 -...

6 meses ago | 1

| accepted