Deriving a value for Full Width Half Maximum

Click For Summary
SUMMARY

The discussion focuses on deriving the Full Width Half Maximum (FWHM) for the response function |R(ω)|² of a driven damped harmonic oscillator, where |R(ω)|² = γ at FWHM. The equation of motion is given by $$\ddot x + \gamma \dot x + \omega_0^2x = \frac{F_{ext}(t)}{m}$$, leading to the response function $$R(ω) = \frac{f_0}{\omega_0^2 - \omega^2 + (i\omega) \gamma}$$. The FWHM is determined by finding the maximum of |R(ω)|² and setting the derivative to zero, ultimately showing that the FWHM equals the damping factor γ.

PREREQUISITES
  • Understanding of damped harmonic oscillators
  • Familiarity with complex numbers in physics
  • Knowledge of calculus, specifically derivatives
  • Ability to manipulate and solve differential equations
NEXT STEPS
  • Study the derivation of the response function for driven damped harmonic oscillators
  • Learn about the significance of Full Width Half Maximum in signal processing
  • Explore the effects of varying the damping factor γ on system behavior
  • Investigate the application of complex analysis in solving differential equations
USEFUL FOR

Students and professionals in physics, particularly those studying oscillatory systems, engineers working with damped systems, and anyone interested in the mathematical modeling of physical phenomena.

Emspak
Messages
240
Reaction score
1

Homework Statement



In a damped harmonic oscillator, show that [itex]|R(\omega)|^2 = \gamma[/itex] at full width half maximum.

2. My attempt at a solution

This is another one of those questions that I feel like I am almost there, but not quite, and it's the math that gets me. But here goes:

For a driven damped harmonic oscillator, show that the full width at half maximum of the response function [itex]| R(\omega)|^2[/itex] is [itex]\gamma[/itex]. Where [itex]\gamma[/itex] is the damping factor. So we start off with: $$\ddot x + \gamma \dot x + \omega_0x = \frac{F_{ext}(t)}{m} = f(t)$$

And I move to solve this. I know that the solution will be [itex]x_0 e^{i\omega t}[/itex] and that I should get the following:
$$-\omega^2 x_0 e^{i\omega t} + (i\omega)x_0 e^{i\omega t}\gamma+ \omega_0^2x_0e^{i\omega t} = \frac{F_0 e^{i\omega t}}{m}$$
which turns into $$x_0 (\omega_0^2 -\omega^2 + (i\omega) \gamma ) = \frac{F_0}{m}= f_0$$ which then becomes $$x_0 = \frac{f_0}{\omega_0^2 -\omega^2 + (i\omega) \gamma }$$
The denominator should give me the response function $R(\omega)$. So I take the absolute value and square it and I get: $$|R(\omega^2)| = \frac{f_0}{\omega_0^2 -\omega^2 t + (i\omega) \gamma } \frac{f_0}{\omega_0^2 -\omega^2 t - (i\omega) \gamma } = \frac{f_0}{(\omega_0^2 -\omega^2)^2 + \gamma^2 \omega^2 }$$

So far so good. I want to find the maximum of this function and that means I want to know where the derivative of the denominator is zero. Taking that derivative: $$\frac{d}{d\omega}(\omega_0^2 -\omega^2)^2 + \gamma^2 \omega^2) = 2((\omega_0^2 -\omega^2)(2\omega) + 2\gamma^2 \omega = 0$$ and since it = 0 we can divide thru by [itex]2\omega[/itex] and we are left with $$-2(\omega_0^2 - \omega^2) + \gamma^2 = 0$$ and solving for [itex]\omega[/itex]: [itex]\omega = \sqrt{\omega_0^2 - \frac{\gamma^2}{2}}[/itex]. We will replace [itex]\gamma[/itex] with [itex]2\beta[/itex], leaving us [itex]\omega = \sqrt{\omega_0^2 - 2 \beta^2}[/itex]

we know where the maximum is, but now we want to find the half-maximum. Plugging my [itex]\omega[/itex] back into my [itex]|R(\omega)|^2[/itex] expression, I have $$\frac{1}{(\omega_0^2 - \omega_0^2 + 2 \beta^2)^2 + 4\beta^4 (\omega_0^2 - 2\beta^2)}=\frac{1}{(2 \beta^2)^2 + 4\beta^4 (\omega_0^2 - 2\beta^2)}=\frac{1}{4 \beta^4 + 4\beta^2 \omega_0^2 - 8\beta^6}$$

which shows me the maximum of [itex]|R(\omega)|^2[/itex]. At 1/2 that is the half maximum and I want to know what [itex]\omega[/itex] is at that point. So going back to my original equation I posit: $$|R(\omega^2)| = \frac{1}{(\omega_0^2 -\omega^2)^2 + 4\beta^2 \omega^2 }=\frac{1}{8 \beta^4 + 8\beta^2 \omega_0^2 - 16\beta^6}$$

But at this point I feel I have lost the plot. Going through this the whole thing struck me as more complicated than it needs to be. My text says showing this should be "an easy exercise." So I am turning to people here to see where I messed up.

Best to you all, and thanks.
 
Last edited:
Physics news on Phys.org
The attempt at a solutionFor a driven damped harmonic oscillator, the equation of motion is given by:$$\ddot x + \gamma \dot x + \omega_0^2x = \frac{F_{ext}(t)}{m} = f(t)$$The solution of this equation can be written as $x(t) = x_0 e^{i\omega t}$ where $\omega$ is the angular frequency. Substituting this into the equation of motion, we get:$$-\omega^2 x_0 e^{i\omega t} + (i\omega)x_0 e^{i\omega t}\gamma+ \omega_0^2x_0e^{i\omega t} = \frac{F_0 e^{i\omega t}}{m}$$ Rearranging this, we get:$$x_0 (\omega_0^2 -\omega^2 + (i\omega) \gamma ) = \frac{F_0}{m}= f_0$$ which yields the response function:$$R(\omega) = \frac{f_0}{\omega_0^2 -\omega^2 + (i\omega) \gamma }$$ Taking the absolute value and squaring this equation gives us:$$|R(\omega)|^2 = \frac{f_0^2}{(\omega_0^2 -\omega^2)^2 + \gamma^2 \omega^2 }$$Now we need to find the full width at half maximum (FWHM) of the response function. To do this, we set the derivative of the above equation to zero in order to find the maximum. Taking the derivative with respect to $\omega$ yields:$$\frac{d}{d\omega}(\omega_0^2 -\omega^2)^2 + \gamma^2 \omega^2) = 2((\omega_0^2 -\omega^2)(2\omega) + 2\gamma^2 \omega
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
1K
Replies
9
Views
2K
Replies
17
Views
3K
Replies
18
Views
2K
  • · Replies 6 ·
Replies
6
Views
2K
  • · Replies 17 ·
Replies
17
Views
1K
Replies
4
Views
3K
Replies
7
Views
1K
Replies
3
Views
1K
  • · Replies 12 ·
Replies
12
Views
2K