> For the complete documentation index, see [llms.txt](https://gtae.gitbook.io/ae4610/llms.txt). Markdown versions of documentation pages are available by appending `.md` to page URLs; this page is available as [Markdown](https://gtae.gitbook.io/ae4610/archive/2-dof-aero/2-dof-helicopter.md).

# A. System Identification

## Introduction

The Quanser Aero Experiment can be configured as a conventional dual-rotor helicopter, as shown in Figure 1a. The front rotor that is horizontal to the ground predominantly affects the motion about the pitch axis while the back or tail rotor mainly affects the motion about the yaw axis (about the vertical shaft).

![Fig. 1a: Quanser Aero Experiment](/files/-MNoecZ-J58ZgtNeU11w)

The tail rotor in helicopters is also known as the anti-torque rotor because it is used to reduce the torque that the main rotor generates about the yaw. Without this, the helicopter would be difficult to stabilize about the yaw axis. The rotors on the Quanser Aero Experiment are the same size and equidistant from the vertical shaft, the tail rotor also generates a torque about the pitch axis. As a result, both the front and back/tail rotors generate torques on each other.

## Background

### Pitch Stiffness

The Quanser 2D AERO is designed with its mass distribution well balanced from the front and back rotors such that its center of gravity is at the mid-point between the two rotors. However, a vertical offset of the cg location is included to result in pitch stiffness from the pendulum effect as depicted in Figure 1b.

<figure><img src="/files/0KDsC7RB0wSVRDiBlV4u" alt=""><figcaption><p>Figure 1b: Pendulum effect creating pitch stiffness due to vertical offset of cg from pitch pivot.</p></figcaption></figure>

The system is in static equilibrium when it is horizontal and parallel to the ground with zero pitch angle. Any change in pitch angle $$\theta\_b$$ gives rise to a restoring moment from the pendulum effect given by restoring pitch moment due to pitch angle change = $$-M\_bgD\_m \sin(\theta\_b)\approx -M\_bgD\_m \theta\_b$$

### Equations of Motion

The free-body diagram of the Quanser Aero Experiment is illustrated in Figure 2.

![Fig. 2: Simple free-body diagram of Quanser Aero Experiment](/files/-MNooya_9MszcrIlGWEO)

The following conventions are used for the modeling:

* The helicopter is horizontal and parallel with the ground when the pitch angle is zero, i.e., $$\theta=0$$ .
* The pitch angle increases positively, $$\theta(t)>0$$ , when the front rotor is moved upwards and the body rotates clockwise (CW) about the Y axis.
* The yaw angle increases positively, $$\psi(t)>0$$ , when the body rotates counter-clockwise (CCW) about  &#x20;the Z axis.
* Pitch increases, $$\theta>0$$ , when the front rotor voltage is positive $$V\_\theta>0$$ .
* Yaw increases, $$\psi>0$$ , when the back (or tail) rotor voltage is positive, $$V\_{\psi}>0$$ .

When voltage is applied to the pitch motor, $$V\_\theta$$ , the speed of rotation results in a force, $$F\_p$$ that acts normal to the body at a distance $$r\_p$$ from the pitch axis as shown in Figure 2. The rotation of the propeller generates a torque about the pitch rotor motor shaft which is in turn seen about the yaw axis. Thus, rotating the pitch propeller does not only cause motion about the pitch axis but also about the yaw axis. As described earlier, that is why conventional helicopters include a tail, or anti-torque, rotor to compensate for the torque generated about the yaw axis by the large, main rotor.

Similarly, the yaw motor causes a force $$F\_y$$ that acts on the body at a distance $$r\_y$$  from the yaw axis as well as a torque about the pitch axis as shown in Figure 2.

The equations of motion can be approximated as:

For pitch axis:

$$
\ddot{\theta} + 2\zeta\_\theta\omega\_{n\_\theta} \displaystyle \dot{\theta} + \omega\_{n\_\theta}^2 \theta = \Tau\_{\theta} \qquad \qquad \tag{1}
$$

For yaw axis:

$$
\ddot{\psi}+ \frac{1}{\tau\_\psi} \dot{\psi}=\Tau\_\psi \qquad \qquad \tag{2}
$$

where the torques acting on the pitch and yaw axes are &#x20;

$$
\Tau\_{\theta} = \overline{K}*{\theta \theta} V*{\theta} + \overline{K}*{\theta \psi} V*{\psi} \qquad \qquad \tag{3}
$$

$$
\Tau\_\psi = \overline{K}*{\psi \theta} V*{\theta} + \overline{K}*{\psi \psi} V*{\psi} \qquad \qquad \tag{4}
$$

The parameters used in the EOMs above are:

* $$\zeta\_\theta$$ : the damping ratio of the pitch dynamics
* $$\omega\_{n\_\theta}$$ : the natural frequency of the pitch dynamics
* $$\tau\_{\psi}$$ : the time constant of the yaw dynamics
* $$\overline{K}*{\theta\theta}$$ : normalized torque thrust gain from the pitch rotor (normalized by $$J*{\theta}$$)
* $$\overline{K}*{\psi\psi}$$ : normalized torque thrust gain from the yaw rotor (normalized by $$J*{\psi}$$)
* $$\overline{K}*{\theta\psi}$$ : normalized cross-torque thrust gain acting on the pitch from the yaw rotor (normalized by $$J*{\theta}$$)
* $$\overline{K}*{\psi\theta}$$ : normalized cross-torque thrust gain acting on the yaw from the pitch rotor (normalized by $$J*{\psi}$$)
* $$V\_\theta$$ : voltage applied to the pitch rotor motor
* $$V\_\psi$$ : voltage applied to the yaw rotor motor
* $$J\_\theta$$ : the total moment of inertia about the pitch axis
* $$J\_{\psi}$$ : the total moment of inertia about the yaw axis

### **First-Order Response**

The step response of a first-order transfer function

$$
Y(s) = \displaystyle \frac{\tau y(0)+KU(s)}{\tau s+1} \qquad \qquad \tag{7}
$$

where y(0) is the initial condition, K is the DC or steady-state gain, and τ is the time constant as illustrated in Figure 3. Figure 3 is for a&#x20;system with y(0) = 0, $$K = 1$$ rad/V and $$\tau = 0.05$$ sec.

To obtain the time constant from the response, find the time it takes to reach $$1-e^{-1}$$or 63.2% of its&#x20;final steady-state value:

$$
y(t\_1)=y\_1=(1-e^{-1})(y\_{\rm ss}-y\_0) +y\_0 \qquad \qquad \tag{8}
$$

The time constant is $$\tau=t\_1-t\_{0 }$$, where $$t\_0$$ is the start time of the step and $$t\_1$$ is the time it takes to reach&#x20;63.2% of the final value, as illustrated in Figure 3.

![Fig. 3 : First-order step response](/files/-MO8ZM_nxDt9dWCWWd9-)

### Second-Order Response

The general equation of motion of a second-order system with input u(t) is described by

$$
\ddot{\alpha} + 2\zeta\omega\_n\dot{\alpha} + \omega\_n^2\alpha = u(t) \qquad \qquad \tag{9}
$$

Assuming the initial conditions $$\alpha(0-)=\alpha\_0$$ and $$\dot\alpha(0-)=0$$, the Laplace transform of Equation (9)&#x20;is

$$
\alpha(s)=\displaystyle{\frac{(s+2\zeta\omega\_n)\alpha\_0 +U(s)}{s^2+2\zeta\omega\_ns+\omega\_n^2}}   \qquad \qquad \tag{10}
$$

The characteristic equation from Equation (10) is obtained as,&#x20;

$$
s^2+2\zeta\omega\_ns+\omega\_n^2 = 0 \qquad \qquad \tag{11}
$$

where $$\zeta$$ is the damping ratio and $$\omega\_n$$ is the natural frequency.

The typical response of an underdamped second order system to a step input with zero initial conditions is shown in Figure 4. The natural frequency and damping ratio can be obtained from transient response using Figure 4 as follows:

#### Finding the Natural Frequency

The period of the oscillations in a system response can be found using the equation

$$
T\_{\rm osc}=\displaystyle\frac{t\_n-t\_1}{n-1} \qquad \qquad \tag{12}
$$

where $$t\_n$$ is the time of the $$n^{th}$$ peak, $$t\_1$$ is the time of the first peak, and $$n$$ is the number of&#x20;oscillations considered. From this, the damped natural frequency (in radians per second) is

$$
\omega\_d=\displaystyle \frac{2\pi}{T\_{\rm osc}} \qquad \qquad \tag{13}
$$

and the undamped natural frequency is

$$
\omega\_n=\displaystyle\frac{\omega\_d}{\sqrt{1-\zeta^2}} \qquad \qquad \tag{14}
$$

#### Finding the Damping Ratio&#xD;

The damping ratio of a second-order system can be found from the transient response to a step input. For a typical second-order&#x20;underdamped system, the subsidence ratio (i.e. decrement ratio) is defined as

$$
\delta = \displaystyle\frac{1}{n-1}\ln\frac{O\_1}{O\_n} \qquad \qquad \tag{15}
$$

![Fig. 4: Step Response of an underdamped second order system.](/files/-MO8i_-_o1XPQpRvSpBR)

where $$O\_1$$ is the peak magnitude of the first oscillation and $$O\_n$$ is the peak magnitude of the $$n^{th}$$ oscillation. Note that $$O\_1$$ > $$O\_n$$ ,&#x20;as this is a decaying response. The damping ratio can then be found using

$$
\zeta= \displaystyle\frac{\delta}{\sqrt{4\pi^2+\delta^2}} \qquad \qquad \tag{16}
$$

### Estimating the Damping Ratio, Natural Frequency and Time Constant

The damping ratio and the natural frequency of the pitch axis, $$\zeta\_\theta$$ and $$\omega\_{n\_\theta}$$ in Equation (1)&#x20;can be found from the oscillatory portion of the step response.&#x20;

**Pitch Axis:** Apply a voltage only to the pitch rotor motor and get the system&#x20;response in the pitch axis. The&#x20;resulting second-order equation of motion is

$$
\ddot{\theta} + 2\zeta\_\theta\omega\_{n\_\theta} \displaystyle \dot{\theta} + \omega\_{n\_\theta}^2 \theta = \overline{K}*{\theta\theta}V*{\theta}(t) \qquad \qquad \tag{17}
$$

Taking its Laplace transform gives

$$
\theta\left(\Theta(s)s^2-s\theta(0^-)-\dot{\theta}(0^-)\right)+2\zeta\_\theta\omega\_{n\_\theta}\left(\Theta(s)-\theta(0^-)\right)+\omega\_{n\_\theta}^2\Theta(s)=\overline{K}*{\theta\theta}V*{\theta}(s)  \qquad \tag{18}
$$

Assuming the initial velocity is zero, $$\dot{\theta}(0^-)=0$$, and solving for position we get

$$
\Theta(s)=\displaystyle\frac{(s+2\zeta\_\theta \omega\_{n\_\theta})\theta(0^-) +\overline{K}*{\theta\theta}V*{\theta}(s)}{s^2+2\zeta\_\theta \omega\_{n\_\theta}s+\omega^2\_{n\_\theta}}\qquad \qquad \tag{19}
$$

The pitch step response transfer function matches the standard second-order transfer function in Equation (10). Thus, the damping ratio and natural frequency of the pitch axis can be determined from the pitch step response or alternatively from free response to an initial pitch disturbance.&#x20;

**Yaw Axis:** The yaw-only equation of motion by applying a voltage input only to the yaw rotor motor is

$$
\ddot{\psi}+ \frac{1}{\tau\_\psi} \dot{\psi}= \overline{K}*{\psi\psi}V*\psi(t) \qquad \qquad \tag{20}
$$

In terms of angular rate, the equation becomes

$$
\dot\omega\_\psi(t)+\frac{1}{\tau\_\psi}\omega\_\psi(t)=\overline{K}*{\psi\psi}V*\psi(t) \qquad \qquad \tag{21}
$$

where $$\omega\_\psi(t)=\dot{\psi}(t)$$ . Taking its Laplace transform

$$
(\Omega\_\psi(s)s-\omega\_\psi(0^-))+\frac{1}{\tau\_\psi}\Omega\_\psi(s)= \overline{K}*{\psi\psi}V*\psi(s)\qquad \qquad \tag{22}
$$

and solving for the speed we get

$$
\Omega\_\psi(s)=\displaystyle \frac{\omega\_\psi(0^-)+\overline{K}*{\psi\psi}V*\psi(s)}{ s+\displaystyle \frac{1}{\tau\_\psi}} = \displaystyle \frac{(\omega\_\psi(0^-)+\overline{K}*{\psi\psi}V*\psi(s))\tau\_\psi}{ \tau\_\psi s+1}   \qquad \qquad \tag{23}
$$

The yaw transfer function matches the first-order transfer function in Equation (7). Thus, the time constant of the yaw axis can be determined from the yaw step response or alternatively from the yaw response to a yaw rate disturbance.     &#x20;

### Estimating the Thrust Parameters&#xD;

#### A) Estimating $$\overline{K}\_{\theta\theta}$$

Providing step input only in the pitch axis (applying a voltage to the pitch rotor motor) allows us to focus on the pitch-only system. Plugging $$V\_\theta \neq 0$$ and $$V\_\psi=0$$ in Equation (1), the equation of motion&#x20;for the pitch axis is

$$
\ddot{\theta} + 2\zeta\_\theta\omega\_{n\_\theta} \dot{\theta} + \omega\_{n\_\theta}^2 \theta = \overline{K}*{\theta \theta} V*{\theta}(t) \qquad \qquad \tag{24}
$$

Solving for the normalized thrust gain we get

$$
\overline{K}*{\theta\theta}=\displaystyle \frac{\ddot{\theta} + 2\zeta*\theta\omega\_{n\_\theta} \dot{\theta} + \omega\_{n\_\theta}^2 \theta}{V\_\theta(t)} \qquad \qquad \tag{25}
$$

Remark that this is the normalized thrust torque gain parameter. The normalized force thrust gain would be $$\overline{K}\_{\theta\theta}/r\_p$$, where $$r\_p$$is&#x20;the distance between the helicopter pivot and the center of the pitch rotor. \
\
**Steady State Method:** By focusing on the steady-state portion of the system response to a step voltage input to the pitch rotor motor, and neglecting the $$\ddot{\theta}$$ and $$\dot{\theta}$$ (zero at steady-state), the normalized thrust gain is obtained as

$$
\overline{K}*{\theta\theta}=\displaystyle \frac{ \omega*{n\_\theta}^2 \theta\_{ss}}{V\_\theta} \qquad \qquad \tag{26}
$$

where $$V\_\theta$$ is the magnitude of the voltage input to the pitch rotor motor.

**Discrete Derivative Method:** The thrust gain parameter can be obtained as

$$
\overline{K}*{\theta\theta} =\frac{\frac{\Delta \Omega*\theta}{\Delta t} + 2 \zeta\_\theta \omega\_{n\_\theta} \Delta \Omega\_\theta}{V\_\theta} \qquad \qquad \tag{27}
$$

where $$\Omega\_\theta = \dot{\theta}$$ and $$\Delta t$$ is the difference between the time instance when the step input is given $$t\_{\text{step}}$$ and a time instance slightly above $$t\_{\text{step}}$$, i.e., $$t = t\_{\text{step}}^+$$. Thus we can find the normalized pitch thrust gain from the measured pitch rate and derived acceleration (discrete-derivative of pitch rate).

#### B) Estimating $$\overline{K}\_{\psi\psi}$$

Similarly, to find the normalized thrust gain acting on the yaw axis system, apply a&#x20;voltage to the tail rotor motor. Plugging $$V\_\theta = 0$$ and $$V\_\psi \neq 0$$ in Equation (1), the equation of motion&#x20;for the yaw axis is

$$
\ddot{\psi}+ \frac{1}{\tau\_\psi} \dot{\psi}=\overline{K}*{\psi\psi}V*\psi \qquad \qquad \tag{28}
$$

or,

$$
\dot\omega\_\psi+\frac{1}{\tau\_\psi} \omega\_\psi=\overline{K}*{\psi\psi}V*\psi \qquad \qquad \tag{29}
$$

where $$\omega\_\psi=\dot\psi$$ is the angular rate of the yaw axis. The normalized yaw torque thrust gain is

$$
\overline{K}*{\psi\psi}=\displaystyle\frac{\dot\omega*\psi+\displaystyle \frac{1}{\tau\_\psi} \omega\_\psi}{V\_\psi} \qquad \qquad \tag{30}
$$

**Steady State Method:** At steady-state condition, and neglecting the $$\dot{\omega}\_\psi$$ (zero at steady-state), the normalized thrust gain is obtained as

$$
\overline{K}*{\psi\psi}=\displaystyle\frac{\displaystyle \frac{1}{\tau*\psi} \omega\_{\psi\_{ss}}}{V\_\psi} \qquad \qquad \tag{31}
$$

where $$V\_\psi$$is the magnitude of the step voltage input to the yaw rotor motor.

**Discrete Derivative Method:** The thrust gain parameter can be obtained as

$$
\overline{K}*{\psi\psi} =\frac{\frac{\Delta \Omega*\psi}{\Delta t} + \frac{1}{\tau\_{\psi}}\Delta \Omega\_\psi}{V\_\psi}\qquad \qquad \tag{32}
$$

where $$\Omega\_\psi = \dot{\psi}$$ and $$\Delta t$$ is the difference between the time instance when the step input is given $$t\_{\text{step}}$$ and a time instance slightly above $$t\_{\text{step}}$$, i.e., $$t = t\_{\text{step}}^+$$. Thus we can find the normalized yaw thrust gain from the measured yaw rate and derived acceleration (discrete-derivative of yaw rate).

#### C) Estimating  $$\overline{K}*{\theta\psi}$$ and $$\overline{K}*{\psi\theta}$$

The normalized cross-torque thrust parameters, $$\overline{K}*{\theta\psi}$$and $$\overline{K}*{\psi\theta}$$ in Equation (3) and (4), represent the coupling between the&#x20;axes. The cross-torque acting on the pitch axis from a torque applied to the tail rotor motor, can be found by appling a voltage to the tail rotor motor, and examining&#x20;the response of the pitch. The equations representing these dynamics when $$V\_\theta=0$$ and $$V\_\psi \neq 0$$, are

$$
\ddot{\theta} + 2\zeta\_\theta\omega\_{n\_\theta} \dot{\theta} + \omega\_{n\_\theta}^2 \theta = \overline{K}*{\theta \psi} V*{\psi} \qquad \qquad \tag{33}
$$

Similarly to identify the cross-torque acting on the yaw axis from a torque applied to the pitch rotor motor, apply a voltage to the pitch rotor motor, and examine&#x20;the response of the yaw. The equations representing these dynamics when $$V\_\theta \neq 0$$ and $$V\_\psi=0$$ , are

$$
\ddot{\psi}+ \frac{1}{\tau\_\psi} \dot{\psi}=\overline{K}*{\psi\theta}V*\theta \qquad \qquad \tag{34}
$$

**Estimating** $$\overline{K}*{\theta\psi}$$ **and** $$\overline{K}*{\psi\theta}$$ **are out of scope of this lab. The values of these parameters are provided in the** [**Analysis**](#analysis) **section.**

## System Identification

### Experimental Steps for Finding System Parameters

First download the zip file below and extract in your group folder. The experiment data is automatically generated in the folder MATLAB is opened to.

{% file src="/files/SZL9B2mPMOxQDBh6k5jD" %}

#### 1) Pitch Motor Voltage Impulse&#xD;

1. Unlock the pitch axis. Lock the yaw axis.
2. Open the *pitch\_impulse* SIMULINK fil&#x65;*.*
3. Set the amplitude gain block to 18. This applies a pulse input of -18 V for 1.5 sec in the pitch axis.&#x20;
4. Select simulation time 30 sec.
5. To build the model, click the down arrow on **Monitor & Tune** under the Hardware tab and then click **Build** **for monitoring** ![](/files/F2dOEp0PbFGdcI8vs0cM). This generates the controller code.
6. Click **Connect** <img src="/files/uORNyY7T0XmLTYaSZyVj" alt="" data-size="line"> button under **Monitor & Tune** and then run SIMULINK by clicking **Start** <img src="/files/OTdtk1LFdBfwxlcgsJX8" alt="" data-size="line">.
7. Copy *`aero_pitch_impulse.mat`* to your folder.&#x20;
8. If the data is not smooth/clean, repeat steps 3-7 with a different pulse voltage.
9. Close the SIMULINK model. DO NOT SAVE THE CHANGE.

Data is saved in following order: \
1: Time \
2: Pitch motor input (V) \
3: Pitch Angle (rad)\
4: Pitch Speed (rad/s)

<figure><img src="/files/8G1a8MMaNdDtjwLs0XKD" alt=""><figcaption><p>Pitch impulse response</p></figcaption></figure>

#### 2) Pitch Motor Voltage Step&#xD;

1. Lock the yaw axis.
2. Open the *pitch\_step* SIMULINK fil&#x65;*.*
3. Apply a step input in the pitch axis. If your setup is Aero1, apply 18 V. If your setup is Aero2, apply 14 V. &#x20;
4. Select simulation time 100 sec.
5. To build the model, click the down arrow on **Monitor & Tune** under the Hardware tab and then click **Build** **for monitoring** ![](/files/F2dOEp0PbFGdcI8vs0cM). This generates the controller code.
6. Click **Connect** <img src="/files/uORNyY7T0XmLTYaSZyVj" alt="" data-size="line"> button under **Monitor & Tune** and then run SIMULINK by clicking **Start** <img src="/files/OTdtk1LFdBfwxlcgsJX8" alt="" data-size="line">.

   <mark style="background-color:red;">**Note:**</mark> <mark style="background-color:red;">**Pitch angle**</mark> <mark style="background-color:red;"></mark><mark style="background-color:red;">has to</mark> <mark style="background-color:red;"></mark><mark style="background-color:red;">**reach steady state**</mark> <mark style="background-color:red;"></mark><mark style="background-color:red;">i.e. constant pitch angle, or you may have to increase simulation time. If steady state has reached, you can end simulation earlier.</mark>
7. Copy *`aero_pitch_step.mat`* to your folder.&#x20;
8. If the data is not smooth/clean, repeat steps 3-7 with a different step voltage.
9. Close the SIMULINK model. DO NOT SAVE THE CHANGE.

Data is saved in following order: \
1: Time \
2: Pitch motor input (V) \
3: Pitch Angle (rad)\
4: Pitch Speed (rad/s)

<figure><img src="/files/k3Vzybwc3nHriDQzSVie" alt=""><figcaption><p>Pitch step response</p></figcaption></figure>

#### 3) Yaw Motor Voltage Step

1. Unlock the yaw axis. Lock the pitch axis.
2. Open the *yaw\_step* SIMULINK fil&#x65;*.*
3. Apply a step input in the yaw axis. If your setup is Aero1, pick a voltage in the range 15-20 V. If your setup is Aero2, pick a voltage in the range 10-11.5 V.&#x20;
4. Select simulation time 100 sec.
5. To build the model, click the down arrow on **Monitor & Tune** under the Hardware tab and then click **Build** **for monitoring** ![](/files/F2dOEp0PbFGdcI8vs0cM). This generates the controller code.
6. Click **Connect** <img src="/files/uORNyY7T0XmLTYaSZyVj" alt="" data-size="line"> button under **Monitor & Tune** and then run SIMULINK by clicking **Start** <img src="/files/OTdtk1LFdBfwxlcgsJX8" alt="" data-size="line">.\ <mark style="background-color:red;">**Note:**</mark> <mark style="background-color:red;">**Pitch angle**</mark> <mark style="background-color:red;"></mark><mark style="background-color:red;">has to</mark> <mark style="background-color:red;"></mark><mark style="background-color:red;">**reach steady state**</mark> <mark style="background-color:red;"></mark><mark style="background-color:red;">i.e. constant pitch angle, or you may have to increase simulation time. If steady state has reached, you can end simulation earlier.</mark>
7. Copy *`aero_yaw_step.mat`* to your folder.&#x20;
8. If the data is not smooth/clean, repeat steps 3-7 with a different step voltage.
9. Close the SIMULINK model. DO NOT SAVE THE CHANGE.

Data is saved in following order: \
1: Time \
2: Yaw motor input (V) \
3: Yaw Angle (rad)\
4: Yaw Speed (rad/s)

<figure><img src="/files/UnA8BgzwlLrpGLQbGCBB" alt=""><figcaption><p>Yaw step response for Aero1 setup.</p></figcaption></figure>

### Analysis

{% file src="/files/dghCDHRHQZQjbLNrEgJs" %} <mark style="color:red;">**MATLAB script to find**</mark> $$\omega\_n$$<mark style="color:red;">**,**</mark>$$\zeta$$ <mark style="color:red;">**and**</mark> $$\overline{K}\_{\theta\theta}$$
{% endfile %}

{% file src="/files/uSJLDl775XmF4eboyEGE" %} <mark style="color:red;">**MATLAB script to find**</mark> $$\tau$$ <mark style="color:red;">**and**</mark> $$\overline{K}\_{\psi\psi}$$
{% endfile %}

1. The cross-torque thrust parameters are $$\overline{K}*{\theta\psi} = 0.0959$$ Nm/V, $$\overline{K}*{\psi\theta} = -0.1227$$ Nm/V. These will be required for [Controller Design](/ae4610/archive/2-dof-aero/b.-control-design.md).
2. From the pitch response due to pitch motor impulse input, find the **natural frequency** and **damping ratio**&#x20;
   * **Hint:** Use equations [12-16](#finding-the-natural-frequency)
   * You can also use the above-attached `Aero_Pitch_SystemID.mlx` script to find $$\omega\_n$$ and $$\zeta$$ . Complete the lines marked with **TODO**.
   * Fill the corresponding values in [Table B.1](#table-b.1)\
     &#x20;
3. From the pitch response due to pitch motor step input find the normalized **thrust** gain $$\overline{K}\_{\theta\theta}$$, using both the steady state method and discrete derivative method.&#x20;
   * **Hint:** For the steady-state method, use equation [26](#a-estimating) and for the discrete derivative method, use equation [27](#a-estimating).
   * You can also use the above-attached `Aero_Pitch_SystemID.mlx` script to find $$\overline{K}\_{\theta\theta}$$. Complete the lines marked with **TODO**.
   * Fill in the values in [Table B.2](#table-b.2).<br>
4. From the yaw **speed** response due to yaw motor step input find the time constant and normalized **thrust** gain $$\overline{K}\_{\psi\psi}$$
   * **Hint:** For time constant. use equation [8](#first-order-response). For $$\overline{K}\_{\psi\psi}$$, use equation [31](#b-estimating) for the steady-state method and use equation [32](#b-estimating) for the discrete derivative method.
   * You can use the above-attached `Aero_Yaw_SystemID.mlx` script to find $$\tau$$ and $$\overline{K}\_{\psi\psi}$$. Complete the lines marked with **TODO**.
   * Fill the respective values in [Table B.](#table-b.1)[1](#table-b.1) and [Table B.2](#table-b.2).

#### Table B.1

<table><thead><tr><th width="365.5">System ID Parameters</th><th>Value</th></tr></thead><tbody><tr><td><span class="math">\omega_{n_\theta}</span></td><td></td></tr><tr><td><span class="math">\zeta_\theta</span></td><td></td></tr><tr><td><span class="math">\tau_\psi</span></td><td></td></tr></tbody></table>

#### Table B.2

<table><thead><tr><th width="234">Thrust Gain Parameters</th><th width="227">  Steady State Method</th><th>   Discrete Derivative Method</th></tr></thead><tbody><tr><td><span class="math">\overline{K}_{\theta\theta}</span></td><td></td><td></td></tr><tr><td><span class="math">\overline{K}_{\psi\psi}</span></td><td></td><td></td></tr></tbody></table>

## Results for Report

1. Equations used for calculating $$\omega\_n$$, $$\zeta$$, $$\tau$$, $$\overline{K}*{\theta\theta},\overline{K}*{\psi\psi}$$.&#x20;
2. Table B.1 and Table B.2.

## Questions for Report

1. The given value of $$\overline{K}*{\theta\psi}$$ is positive and $$\overline{K}*{\psi\theta}$$ is negative for our setup (see Step 1 of [Analysis](#analysis) section). Why?
2. Compare the normalized thrust gains obtained from the steady-state method and discrete derivative method and explain any differences.


---

# Agent Instructions
This documentation is published with GitBook. GitBook is the documentation platform designed so that both humans and AI agents can read, navigate, and reason over technical content effectively. Learn more at gitbook.com.

## Querying This Documentation
If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question.

Perform an HTTP GET request on the current page URL with the `ask` query parameter:

```
GET https://gtae.gitbook.io/ae4610/archive/2-dof-aero/2-dof-helicopter.md?ask=<question>
```

The question should be specific, self-contained, and written in natural language.
The response will contain a direct answer to the question and relevant excerpts and sources from the documentation.

Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.
