How do I fit a resonance curve?How to guess the correct fitting function to some data?Fitting of exponential data gives me a constant functionCurve Fitting and Multiple ExperimentsPower fit to some experimental dataShould a Gaussian Curve Always Be Drawn Symmetrically?Arrhenius Fit: Linear or exponential form?Reduced chi-squared value for noiseless spectraDamped Harmonic Curve fit and ForceConstant wind drag while falling?Removing zero-counts in exponential decay measurement

Why do these characters still seem to be the same age after the events of Endgame?

Is any special diet a treatment of autism?

Nested loops to process groups of pictures

As a GM, is it bad form to ask for a moment to think when improvising?

Removing racism on a multi raced world

Is it normal for gliders not to have attitude indicators?

Is Soreness in Middle Knuckle of Fretting Hand Index Finger Normal for Beginners?

What was Bran's plan to kill the Night King?

Are pressure-treated posts that have been submerged for a few days ruined?

Would you use "llamarse" for an animal's name?

Should I decline this job offer that requires relocating to an area with high cost of living?

Agena docking and RCS Brakes in First Man

Will 700 more planes a day fly because of the Heathrow expansion?

What is a common way to tell if an academic is "above average," or outstanding in their field? Is their h-index (Hirsh index) one of them?

Why is my arithmetic with a long long int behaving this way?

Why do people keep telling me that I am a bad photographer?

Can there be a single technologically advanced nation, in a continent full of non-technologically advanced nations?

Is disk brake effectiveness mitigated by tyres losing traction under strong braking?

Extra space in cells when using token lists to build tabular content

Snap victim memorial reference in Avengers: Endgame

How does the reduce() method work in Java 8?

Is there an age requirement to play in Adventurers League?

How do I, as a DM, handle a party that decides to set up an ambush in a dungeon?

Voltage Balun 1:1



How do I fit a resonance curve?


How to guess the correct fitting function to some data?Fitting of exponential data gives me a constant functionCurve Fitting and Multiple ExperimentsPower fit to some experimental dataShould a Gaussian Curve Always Be Drawn Symmetrically?Arrhenius Fit: Linear or exponential form?Reduced chi-squared value for noiseless spectraDamped Harmonic Curve fit and ForceConstant wind drag while falling?Removing zero-counts in exponential decay measurement













7












$begingroup$


In an experiment, I collected data points $ (ω,υ(ω))$ that are modeled by the equation:



$$ υ(ω)=fracωCsqrt(ω^2-ω_0^2)^2+γ^2ω^2 ,.$$



How can do I fit the data to the above correlation? And how can I extract $γ$ through this process?










share|cite|improve this question











$endgroup$











  • $begingroup$
    almost looks like the magnitude response of a second-order bandpass filter. other than least-squares fit (or some other $L_p$ metric of error), i dunno how else to get $gamma$. seems to me that the least-squares fit also needs to find $omega_0$ and $C$. but i think $C$ can come out in the wash.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:24











  • $begingroup$
    $C$ can be absorbed into $gamma$.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:28










  • $begingroup$
    You said nothing about how you estimate the errors of your measured points, and what your criterion for a good fit would be. In general, you would want to maximize some likelihood function, in practice (with Gaussian 1D errors on the data points) a least-squares fit may be good enough. Your case probably involves bins ($omega$ being a continuous variable) and the Poisson distribution (few entries in bins far away from $omega_0$), so a more complex Log-Likelihood approach could be called for.
    $endgroup$
    – tobi_s
    Apr 1 at 10:56






  • 3




    $begingroup$
    I strongly disagree that this should be posted on Mathematics, as suggested by the close votes on it. I think such questions are on topic here, but were they not then Cross Validated would be the most obvious choice.
    $endgroup$
    – Kyle Kanos
    Apr 1 at 11:35















7












$begingroup$


In an experiment, I collected data points $ (ω,υ(ω))$ that are modeled by the equation:



$$ υ(ω)=fracωCsqrt(ω^2-ω_0^2)^2+γ^2ω^2 ,.$$



How can do I fit the data to the above correlation? And how can I extract $γ$ through this process?










share|cite|improve this question











$endgroup$











  • $begingroup$
    almost looks like the magnitude response of a second-order bandpass filter. other than least-squares fit (or some other $L_p$ metric of error), i dunno how else to get $gamma$. seems to me that the least-squares fit also needs to find $omega_0$ and $C$. but i think $C$ can come out in the wash.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:24











  • $begingroup$
    $C$ can be absorbed into $gamma$.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:28










  • $begingroup$
    You said nothing about how you estimate the errors of your measured points, and what your criterion for a good fit would be. In general, you would want to maximize some likelihood function, in practice (with Gaussian 1D errors on the data points) a least-squares fit may be good enough. Your case probably involves bins ($omega$ being a continuous variable) and the Poisson distribution (few entries in bins far away from $omega_0$), so a more complex Log-Likelihood approach could be called for.
    $endgroup$
    – tobi_s
    Apr 1 at 10:56






  • 3




    $begingroup$
    I strongly disagree that this should be posted on Mathematics, as suggested by the close votes on it. I think such questions are on topic here, but were they not then Cross Validated would be the most obvious choice.
    $endgroup$
    – Kyle Kanos
    Apr 1 at 11:35













7












7








7


1



$begingroup$


In an experiment, I collected data points $ (ω,υ(ω))$ that are modeled by the equation:



$$ υ(ω)=fracωCsqrt(ω^2-ω_0^2)^2+γ^2ω^2 ,.$$



How can do I fit the data to the above correlation? And how can I extract $γ$ through this process?










share|cite|improve this question











$endgroup$




In an experiment, I collected data points $ (ω,υ(ω))$ that are modeled by the equation:



$$ υ(ω)=fracωCsqrt(ω^2-ω_0^2)^2+γ^2ω^2 ,.$$



How can do I fit the data to the above correlation? And how can I extract $γ$ through this process?







experimental-physics correlation-functions data-analysis






share|cite|improve this question















share|cite|improve this question













share|cite|improve this question




share|cite|improve this question








edited Apr 1 at 9:31









knzhou

48.4k12132234




48.4k12132234










asked Mar 31 at 19:47









Unstoppable TachyonUnstoppable Tachyon

987




987











  • $begingroup$
    almost looks like the magnitude response of a second-order bandpass filter. other than least-squares fit (or some other $L_p$ metric of error), i dunno how else to get $gamma$. seems to me that the least-squares fit also needs to find $omega_0$ and $C$. but i think $C$ can come out in the wash.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:24











  • $begingroup$
    $C$ can be absorbed into $gamma$.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:28










  • $begingroup$
    You said nothing about how you estimate the errors of your measured points, and what your criterion for a good fit would be. In general, you would want to maximize some likelihood function, in practice (with Gaussian 1D errors on the data points) a least-squares fit may be good enough. Your case probably involves bins ($omega$ being a continuous variable) and the Poisson distribution (few entries in bins far away from $omega_0$), so a more complex Log-Likelihood approach could be called for.
    $endgroup$
    – tobi_s
    Apr 1 at 10:56






  • 3




    $begingroup$
    I strongly disagree that this should be posted on Mathematics, as suggested by the close votes on it. I think such questions are on topic here, but were they not then Cross Validated would be the most obvious choice.
    $endgroup$
    – Kyle Kanos
    Apr 1 at 11:35
















  • $begingroup$
    almost looks like the magnitude response of a second-order bandpass filter. other than least-squares fit (or some other $L_p$ metric of error), i dunno how else to get $gamma$. seems to me that the least-squares fit also needs to find $omega_0$ and $C$. but i think $C$ can come out in the wash.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:24











  • $begingroup$
    $C$ can be absorbed into $gamma$.
    $endgroup$
    – robert bristow-johnson
    Apr 1 at 7:28










  • $begingroup$
    You said nothing about how you estimate the errors of your measured points, and what your criterion for a good fit would be. In general, you would want to maximize some likelihood function, in practice (with Gaussian 1D errors on the data points) a least-squares fit may be good enough. Your case probably involves bins ($omega$ being a continuous variable) and the Poisson distribution (few entries in bins far away from $omega_0$), so a more complex Log-Likelihood approach could be called for.
    $endgroup$
    – tobi_s
    Apr 1 at 10:56






  • 3




    $begingroup$
    I strongly disagree that this should be posted on Mathematics, as suggested by the close votes on it. I think such questions are on topic here, but were they not then Cross Validated would be the most obvious choice.
    $endgroup$
    – Kyle Kanos
    Apr 1 at 11:35















$begingroup$
almost looks like the magnitude response of a second-order bandpass filter. other than least-squares fit (or some other $L_p$ metric of error), i dunno how else to get $gamma$. seems to me that the least-squares fit also needs to find $omega_0$ and $C$. but i think $C$ can come out in the wash.
$endgroup$
– robert bristow-johnson
Apr 1 at 7:24





$begingroup$
almost looks like the magnitude response of a second-order bandpass filter. other than least-squares fit (or some other $L_p$ metric of error), i dunno how else to get $gamma$. seems to me that the least-squares fit also needs to find $omega_0$ and $C$. but i think $C$ can come out in the wash.
$endgroup$
– robert bristow-johnson
Apr 1 at 7:24













$begingroup$
$C$ can be absorbed into $gamma$.
$endgroup$
– robert bristow-johnson
Apr 1 at 7:28




$begingroup$
$C$ can be absorbed into $gamma$.
$endgroup$
– robert bristow-johnson
Apr 1 at 7:28












$begingroup$
You said nothing about how you estimate the errors of your measured points, and what your criterion for a good fit would be. In general, you would want to maximize some likelihood function, in practice (with Gaussian 1D errors on the data points) a least-squares fit may be good enough. Your case probably involves bins ($omega$ being a continuous variable) and the Poisson distribution (few entries in bins far away from $omega_0$), so a more complex Log-Likelihood approach could be called for.
$endgroup$
– tobi_s
Apr 1 at 10:56




$begingroup$
You said nothing about how you estimate the errors of your measured points, and what your criterion for a good fit would be. In general, you would want to maximize some likelihood function, in practice (with Gaussian 1D errors on the data points) a least-squares fit may be good enough. Your case probably involves bins ($omega$ being a continuous variable) and the Poisson distribution (few entries in bins far away from $omega_0$), so a more complex Log-Likelihood approach could be called for.
$endgroup$
– tobi_s
Apr 1 at 10:56




3




3




$begingroup$
I strongly disagree that this should be posted on Mathematics, as suggested by the close votes on it. I think such questions are on topic here, but were they not then Cross Validated would be the most obvious choice.
$endgroup$
– Kyle Kanos
Apr 1 at 11:35




$begingroup$
I strongly disagree that this should be posted on Mathematics, as suggested by the close votes on it. I think such questions are on topic here, but were they not then Cross Validated would be the most obvious choice.
$endgroup$
– Kyle Kanos
Apr 1 at 11:35










4 Answers
4






active

oldest

votes


















10












$begingroup$

What you want to find is the parameters $theta=(C, omega_0, gamma)$ that minimizes the difference between $nu(omega|theta)$ (the curve given the parameters) and the measured $nu_i$ values.



The most popular method is least mean square fitting, which minimizes the sum of the squares of the differences. One can also do it by formulating the normal equations and solve it as a (potentially big) linear equation system. Another approach is the Gauss-Newton algorithm, a simple iterative method to do it. It is a good exercise to implement the solution oneself, but once you have done it once or twice it is best to rely on some software package.



Note that this kind of fitting works well when you know the functional form (your equation for $nu(omega)$), since you can ensure only that the parameters that matter are included. If you try to fit some general polynomial or function you can get overfitting (some complex curve that fits all the data but has nothing to do with your problem) besides the problem of identifying the parameters you care about.






share|cite|improve this answer









$endgroup$




















    8












    $begingroup$

    Don't try using any general-purpose curve fitting algorithm for this.



    The form of your function looks like a frequency response function, with the two unknown parameters $omega_0$ and $gamma$ - i.e. the resonant frequency, and the damping parameter. The function you specified omits an important feature if this is measured data, namely the relative phase between the "force" driving the oscillation and the response.



    If you didn't measure the phase at each frequency, repeat the experiment, because that is critical information.



    When you have the amplitude and phase data, there are curve fitting techniques devised specifically for this problem of "system identification" in experimental modal analysis. A simple one is the so-called "circle fitting" method. If you make a Nyquist plot of your measured data (i.e. plot imaginary part of the response against the real part), the section of the curve near the resonance is a circle, and you can fit a circle to the measured data and find the parameters from it.



    In practice, a simplistic approach assuming the system only has one resonance often doesn't work well, because the response of a real system near resonance also includes the off-resonance response to all the other vibration modes. If the resonant frequencies are well separated and lightly damped, it is possible to correct for this while fitting "one mode at a time". If this is not the case, you need methods that can identify several resonances simultaneously from one response function.



    Rather than re-invent the wheel, use existing code. The signal processing toolbox in MATLAB would be a good starting point - for example https://uk.mathworks.com/help/signal/ref/modalfit.html






    share|cite|improve this answer











    $endgroup$








    • 7




      $begingroup$
      That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
      $endgroup$
      – Emilio Pisanty
      Mar 31 at 22:14










    • $begingroup$
      what is a well-known method for identifying several closely spaced resonances at the same time?
      $endgroup$
      – AmIAStudent
      Mar 31 at 22:56






    • 2




      $begingroup$
      What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
      $endgroup$
      – Federico Poloni
      Apr 1 at 13:40



















    6












    $begingroup$

    If we put:



    $$Y = fracomega^2u(omega)^2$$



    and



    $$X = omega^2$$



    the equation becomes:



    $$Y =fracX^2C^2 +frac(gamma^2 - 2 omega_0^2)C^2 X + fracomega_0^4C^2$$



    You can then extract the coefficients using polynomial fitting. To get the least-squares fit right, you have to compute the errors in $Y$ and $X$ for each data point from the measurement errors in $omega$ and $u(omega)$.






    share|cite|improve this answer









    $endgroup$








    • 4




      $begingroup$
      This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
      $endgroup$
      – tobi_s
      Apr 1 at 10:43


















    0












    $begingroup$

    Are you looking for something like polynomial regression? The general idea is, if you have measured pairs of (x, y(x)) and you are looking for find a fit of the form:



    $$y = alpha_0 + alpha_1 x + alpha_2 x^2 ...$$



    You can write this in matrix form as:



    $$beginbmatrix y_1 \ y_2 \ y_3 \ vdots \ y_n endbmatrix = beginbmatrix 1 & x_1 & x_1^2 & cdots \ 1 & x_2 & x_2^2 & cdots \ 1 & x_3 & x_3^2 & cdots \ vdots & vdots & vdots & vdots \ 1 & x_n &x_n^2 & cdots endbmatrix beginbmatrix alpha_0 \ alpha_1 \ alpha_2 \ vdots \ alpha_m endbmatrix$$



    This can now be solved for your coefficients, $alpha_i$. That being said, and as was hinted at in your comments, I've never actually done this, and have instead used non-linear fitting functions provided by libraries.



    More information on polynomial regression on the wikipedia page.



    Edit: As you say in the comments, this method is only applicable if you can write your function that you wish to fit in polynomial form, which I don't think you can do for your example. In which case you are best off referring to the other answers to this question.






    share|cite|improve this answer











    $endgroup$








    • 2




      $begingroup$
      The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
      $endgroup$
      – Unstoppable Tachyon
      Mar 31 at 20:28










    • $begingroup$
      Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
      $endgroup$
      – Anon1759
      Mar 31 at 20:32










    • $begingroup$
      The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
      $endgroup$
      – Vladimir F
      Apr 1 at 13:04










    • $begingroup$
      @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
      $endgroup$
      – Vladimir F
      Apr 1 at 13:06










    • $begingroup$
      @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
      $endgroup$
      – Federico Poloni
      Apr 1 at 13:37












    Your Answer








    StackExchange.ready(function()
    var channelOptions =
    tags: "".split(" "),
    id: "151"
    ;
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function()
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled)
    StackExchange.using("snippets", function()
    createEditor();
    );

    else
    createEditor();

    );

    function createEditor()
    StackExchange.prepareEditor(
    heartbeatType: 'answer',
    autoActivateHeartbeat: false,
    convertImagesToLinks: false,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: null,
    bindNavPrevention: true,
    postfix: "",
    imageUploader:
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    ,
    noCode: true, onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    );



    );













    draft saved

    draft discarded


















    StackExchange.ready(
    function ()
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fphysics.stackexchange.com%2fquestions%2f469754%2fhow-do-i-fit-a-resonance-curve%23new-answer', 'question_page');

    );

    Post as a guest















    Required, but never shown

























    4 Answers
    4






    active

    oldest

    votes








    4 Answers
    4






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    10












    $begingroup$

    What you want to find is the parameters $theta=(C, omega_0, gamma)$ that minimizes the difference between $nu(omega|theta)$ (the curve given the parameters) and the measured $nu_i$ values.



    The most popular method is least mean square fitting, which minimizes the sum of the squares of the differences. One can also do it by formulating the normal equations and solve it as a (potentially big) linear equation system. Another approach is the Gauss-Newton algorithm, a simple iterative method to do it. It is a good exercise to implement the solution oneself, but once you have done it once or twice it is best to rely on some software package.



    Note that this kind of fitting works well when you know the functional form (your equation for $nu(omega)$), since you can ensure only that the parameters that matter are included. If you try to fit some general polynomial or function you can get overfitting (some complex curve that fits all the data but has nothing to do with your problem) besides the problem of identifying the parameters you care about.






    share|cite|improve this answer









    $endgroup$

















      10












      $begingroup$

      What you want to find is the parameters $theta=(C, omega_0, gamma)$ that minimizes the difference between $nu(omega|theta)$ (the curve given the parameters) and the measured $nu_i$ values.



      The most popular method is least mean square fitting, which minimizes the sum of the squares of the differences. One can also do it by formulating the normal equations and solve it as a (potentially big) linear equation system. Another approach is the Gauss-Newton algorithm, a simple iterative method to do it. It is a good exercise to implement the solution oneself, but once you have done it once or twice it is best to rely on some software package.



      Note that this kind of fitting works well when you know the functional form (your equation for $nu(omega)$), since you can ensure only that the parameters that matter are included. If you try to fit some general polynomial or function you can get overfitting (some complex curve that fits all the data but has nothing to do with your problem) besides the problem of identifying the parameters you care about.






      share|cite|improve this answer









      $endgroup$















        10












        10








        10





        $begingroup$

        What you want to find is the parameters $theta=(C, omega_0, gamma)$ that minimizes the difference between $nu(omega|theta)$ (the curve given the parameters) and the measured $nu_i$ values.



        The most popular method is least mean square fitting, which minimizes the sum of the squares of the differences. One can also do it by formulating the normal equations and solve it as a (potentially big) linear equation system. Another approach is the Gauss-Newton algorithm, a simple iterative method to do it. It is a good exercise to implement the solution oneself, but once you have done it once or twice it is best to rely on some software package.



        Note that this kind of fitting works well when you know the functional form (your equation for $nu(omega)$), since you can ensure only that the parameters that matter are included. If you try to fit some general polynomial or function you can get overfitting (some complex curve that fits all the data but has nothing to do with your problem) besides the problem of identifying the parameters you care about.






        share|cite|improve this answer









        $endgroup$



        What you want to find is the parameters $theta=(C, omega_0, gamma)$ that minimizes the difference between $nu(omega|theta)$ (the curve given the parameters) and the measured $nu_i$ values.



        The most popular method is least mean square fitting, which minimizes the sum of the squares of the differences. One can also do it by formulating the normal equations and solve it as a (potentially big) linear equation system. Another approach is the Gauss-Newton algorithm, a simple iterative method to do it. It is a good exercise to implement the solution oneself, but once you have done it once or twice it is best to rely on some software package.



        Note that this kind of fitting works well when you know the functional form (your equation for $nu(omega)$), since you can ensure only that the parameters that matter are included. If you try to fit some general polynomial or function you can get overfitting (some complex curve that fits all the data but has nothing to do with your problem) besides the problem of identifying the parameters you care about.







        share|cite|improve this answer












        share|cite|improve this answer



        share|cite|improve this answer










        answered Mar 31 at 20:28









        Anders SandbergAnders Sandberg

        10.7k21532




        10.7k21532





















            8












            $begingroup$

            Don't try using any general-purpose curve fitting algorithm for this.



            The form of your function looks like a frequency response function, with the two unknown parameters $omega_0$ and $gamma$ - i.e. the resonant frequency, and the damping parameter. The function you specified omits an important feature if this is measured data, namely the relative phase between the "force" driving the oscillation and the response.



            If you didn't measure the phase at each frequency, repeat the experiment, because that is critical information.



            When you have the amplitude and phase data, there are curve fitting techniques devised specifically for this problem of "system identification" in experimental modal analysis. A simple one is the so-called "circle fitting" method. If you make a Nyquist plot of your measured data (i.e. plot imaginary part of the response against the real part), the section of the curve near the resonance is a circle, and you can fit a circle to the measured data and find the parameters from it.



            In practice, a simplistic approach assuming the system only has one resonance often doesn't work well, because the response of a real system near resonance also includes the off-resonance response to all the other vibration modes. If the resonant frequencies are well separated and lightly damped, it is possible to correct for this while fitting "one mode at a time". If this is not the case, you need methods that can identify several resonances simultaneously from one response function.



            Rather than re-invent the wheel, use existing code. The signal processing toolbox in MATLAB would be a good starting point - for example https://uk.mathworks.com/help/signal/ref/modalfit.html






            share|cite|improve this answer











            $endgroup$








            • 7




              $begingroup$
              That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
              $endgroup$
              – Emilio Pisanty
              Mar 31 at 22:14










            • $begingroup$
              what is a well-known method for identifying several closely spaced resonances at the same time?
              $endgroup$
              – AmIAStudent
              Mar 31 at 22:56






            • 2




              $begingroup$
              What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:40
















            8












            $begingroup$

            Don't try using any general-purpose curve fitting algorithm for this.



            The form of your function looks like a frequency response function, with the two unknown parameters $omega_0$ and $gamma$ - i.e. the resonant frequency, and the damping parameter. The function you specified omits an important feature if this is measured data, namely the relative phase between the "force" driving the oscillation and the response.



            If you didn't measure the phase at each frequency, repeat the experiment, because that is critical information.



            When you have the amplitude and phase data, there are curve fitting techniques devised specifically for this problem of "system identification" in experimental modal analysis. A simple one is the so-called "circle fitting" method. If you make a Nyquist plot of your measured data (i.e. plot imaginary part of the response against the real part), the section of the curve near the resonance is a circle, and you can fit a circle to the measured data and find the parameters from it.



            In practice, a simplistic approach assuming the system only has one resonance often doesn't work well, because the response of a real system near resonance also includes the off-resonance response to all the other vibration modes. If the resonant frequencies are well separated and lightly damped, it is possible to correct for this while fitting "one mode at a time". If this is not the case, you need methods that can identify several resonances simultaneously from one response function.



            Rather than re-invent the wheel, use existing code. The signal processing toolbox in MATLAB would be a good starting point - for example https://uk.mathworks.com/help/signal/ref/modalfit.html






            share|cite|improve this answer











            $endgroup$








            • 7




              $begingroup$
              That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
              $endgroup$
              – Emilio Pisanty
              Mar 31 at 22:14










            • $begingroup$
              what is a well-known method for identifying several closely spaced resonances at the same time?
              $endgroup$
              – AmIAStudent
              Mar 31 at 22:56






            • 2




              $begingroup$
              What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:40














            8












            8








            8





            $begingroup$

            Don't try using any general-purpose curve fitting algorithm for this.



            The form of your function looks like a frequency response function, with the two unknown parameters $omega_0$ and $gamma$ - i.e. the resonant frequency, and the damping parameter. The function you specified omits an important feature if this is measured data, namely the relative phase between the "force" driving the oscillation and the response.



            If you didn't measure the phase at each frequency, repeat the experiment, because that is critical information.



            When you have the amplitude and phase data, there are curve fitting techniques devised specifically for this problem of "system identification" in experimental modal analysis. A simple one is the so-called "circle fitting" method. If you make a Nyquist plot of your measured data (i.e. plot imaginary part of the response against the real part), the section of the curve near the resonance is a circle, and you can fit a circle to the measured data and find the parameters from it.



            In practice, a simplistic approach assuming the system only has one resonance often doesn't work well, because the response of a real system near resonance also includes the off-resonance response to all the other vibration modes. If the resonant frequencies are well separated and lightly damped, it is possible to correct for this while fitting "one mode at a time". If this is not the case, you need methods that can identify several resonances simultaneously from one response function.



            Rather than re-invent the wheel, use existing code. The signal processing toolbox in MATLAB would be a good starting point - for example https://uk.mathworks.com/help/signal/ref/modalfit.html






            share|cite|improve this answer











            $endgroup$



            Don't try using any general-purpose curve fitting algorithm for this.



            The form of your function looks like a frequency response function, with the two unknown parameters $omega_0$ and $gamma$ - i.e. the resonant frequency, and the damping parameter. The function you specified omits an important feature if this is measured data, namely the relative phase between the "force" driving the oscillation and the response.



            If you didn't measure the phase at each frequency, repeat the experiment, because that is critical information.



            When you have the amplitude and phase data, there are curve fitting techniques devised specifically for this problem of "system identification" in experimental modal analysis. A simple one is the so-called "circle fitting" method. If you make a Nyquist plot of your measured data (i.e. plot imaginary part of the response against the real part), the section of the curve near the resonance is a circle, and you can fit a circle to the measured data and find the parameters from it.



            In practice, a simplistic approach assuming the system only has one resonance often doesn't work well, because the response of a real system near resonance also includes the off-resonance response to all the other vibration modes. If the resonant frequencies are well separated and lightly damped, it is possible to correct for this while fitting "one mode at a time". If this is not the case, you need methods that can identify several resonances simultaneously from one response function.



            Rather than re-invent the wheel, use existing code. The signal processing toolbox in MATLAB would be a good starting point - for example https://uk.mathworks.com/help/signal/ref/modalfit.html







            share|cite|improve this answer














            share|cite|improve this answer



            share|cite|improve this answer








            edited Mar 31 at 21:03

























            answered Mar 31 at 20:58









            alephzeroalephzero

            6,08521123




            6,08521123







            • 7




              $begingroup$
              That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
              $endgroup$
              – Emilio Pisanty
              Mar 31 at 22:14










            • $begingroup$
              what is a well-known method for identifying several closely spaced resonances at the same time?
              $endgroup$
              – AmIAStudent
              Mar 31 at 22:56






            • 2




              $begingroup$
              What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:40













            • 7




              $begingroup$
              That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
              $endgroup$
              – Emilio Pisanty
              Mar 31 at 22:14










            • $begingroup$
              what is a well-known method for identifying several closely spaced resonances at the same time?
              $endgroup$
              – AmIAStudent
              Mar 31 at 22:56






            • 2




              $begingroup$
              What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:40








            7




            7




            $begingroup$
            That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
            $endgroup$
            – Emilio Pisanty
            Mar 31 at 22:14




            $begingroup$
            That is, of course, if the phase information is experimentally accessible. It's measurable in plenty of systems, but there are also many cases where it is either inaccessible or much more expensive to access.
            $endgroup$
            – Emilio Pisanty
            Mar 31 at 22:14












            $begingroup$
            what is a well-known method for identifying several closely spaced resonances at the same time?
            $endgroup$
            – AmIAStudent
            Mar 31 at 22:56




            $begingroup$
            what is a well-known method for identifying several closely spaced resonances at the same time?
            $endgroup$
            – AmIAStudent
            Mar 31 at 22:56




            2




            2




            $begingroup$
            What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
            $endgroup$
            – Federico Poloni
            Apr 1 at 13:40





            $begingroup$
            What are the advantages of these algorithms with respect to the general-purpose ones? What cost function do they minimize?
            $endgroup$
            – Federico Poloni
            Apr 1 at 13:40












            6












            $begingroup$

            If we put:



            $$Y = fracomega^2u(omega)^2$$



            and



            $$X = omega^2$$



            the equation becomes:



            $$Y =fracX^2C^2 +frac(gamma^2 - 2 omega_0^2)C^2 X + fracomega_0^4C^2$$



            You can then extract the coefficients using polynomial fitting. To get the least-squares fit right, you have to compute the errors in $Y$ and $X$ for each data point from the measurement errors in $omega$ and $u(omega)$.






            share|cite|improve this answer









            $endgroup$








            • 4




              $begingroup$
              This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
              $endgroup$
              – tobi_s
              Apr 1 at 10:43















            6












            $begingroup$

            If we put:



            $$Y = fracomega^2u(omega)^2$$



            and



            $$X = omega^2$$



            the equation becomes:



            $$Y =fracX^2C^2 +frac(gamma^2 - 2 omega_0^2)C^2 X + fracomega_0^4C^2$$



            You can then extract the coefficients using polynomial fitting. To get the least-squares fit right, you have to compute the errors in $Y$ and $X$ for each data point from the measurement errors in $omega$ and $u(omega)$.






            share|cite|improve this answer









            $endgroup$








            • 4




              $begingroup$
              This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
              $endgroup$
              – tobi_s
              Apr 1 at 10:43













            6












            6








            6





            $begingroup$

            If we put:



            $$Y = fracomega^2u(omega)^2$$



            and



            $$X = omega^2$$



            the equation becomes:



            $$Y =fracX^2C^2 +frac(gamma^2 - 2 omega_0^2)C^2 X + fracomega_0^4C^2$$



            You can then extract the coefficients using polynomial fitting. To get the least-squares fit right, you have to compute the errors in $Y$ and $X$ for each data point from the measurement errors in $omega$ and $u(omega)$.






            share|cite|improve this answer









            $endgroup$



            If we put:



            $$Y = fracomega^2u(omega)^2$$



            and



            $$X = omega^2$$



            the equation becomes:



            $$Y =fracX^2C^2 +frac(gamma^2 - 2 omega_0^2)C^2 X + fracomega_0^4C^2$$



            You can then extract the coefficients using polynomial fitting. To get the least-squares fit right, you have to compute the errors in $Y$ and $X$ for each data point from the measurement errors in $omega$ and $u(omega)$.







            share|cite|improve this answer












            share|cite|improve this answer



            share|cite|improve this answer










            answered Apr 1 at 0:07









            Count IblisCount Iblis

            8,55411440




            8,55411440







            • 4




              $begingroup$
              This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
              $endgroup$
              – tobi_s
              Apr 1 at 10:43












            • 4




              $begingroup$
              This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
              $endgroup$
              – tobi_s
              Apr 1 at 10:43







            4




            4




            $begingroup$
            This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
            $endgroup$
            – tobi_s
            Apr 1 at 10:43




            $begingroup$
            This is a beautiful transformation, but it will also distort the error distributions of the data points, rendering the fit much harder as result. The preferred approach to fitting depends on your measurement errors (or error estimates) and on whether your data is binned. For unbinned data, if the $omega$ values are exact (or if the error is negligible compared to $v(omega)$), and if the errors on $v(omega)$ are drawn from a Gaussian distribution, a (non-linear) least-squares fit to your data points is hard to beat, as it will also be a maximum-likelihood fit.
            $endgroup$
            – tobi_s
            Apr 1 at 10:43











            0












            $begingroup$

            Are you looking for something like polynomial regression? The general idea is, if you have measured pairs of (x, y(x)) and you are looking for find a fit of the form:



            $$y = alpha_0 + alpha_1 x + alpha_2 x^2 ...$$



            You can write this in matrix form as:



            $$beginbmatrix y_1 \ y_2 \ y_3 \ vdots \ y_n endbmatrix = beginbmatrix 1 & x_1 & x_1^2 & cdots \ 1 & x_2 & x_2^2 & cdots \ 1 & x_3 & x_3^2 & cdots \ vdots & vdots & vdots & vdots \ 1 & x_n &x_n^2 & cdots endbmatrix beginbmatrix alpha_0 \ alpha_1 \ alpha_2 \ vdots \ alpha_m endbmatrix$$



            This can now be solved for your coefficients, $alpha_i$. That being said, and as was hinted at in your comments, I've never actually done this, and have instead used non-linear fitting functions provided by libraries.



            More information on polynomial regression on the wikipedia page.



            Edit: As you say in the comments, this method is only applicable if you can write your function that you wish to fit in polynomial form, which I don't think you can do for your example. In which case you are best off referring to the other answers to this question.






            share|cite|improve this answer











            $endgroup$








            • 2




              $begingroup$
              The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
              $endgroup$
              – Unstoppable Tachyon
              Mar 31 at 20:28










            • $begingroup$
              Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
              $endgroup$
              – Anon1759
              Mar 31 at 20:32










            • $begingroup$
              The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
              $endgroup$
              – Vladimir F
              Apr 1 at 13:04










            • $begingroup$
              @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
              $endgroup$
              – Vladimir F
              Apr 1 at 13:06










            • $begingroup$
              @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:37
















            0












            $begingroup$

            Are you looking for something like polynomial regression? The general idea is, if you have measured pairs of (x, y(x)) and you are looking for find a fit of the form:



            $$y = alpha_0 + alpha_1 x + alpha_2 x^2 ...$$



            You can write this in matrix form as:



            $$beginbmatrix y_1 \ y_2 \ y_3 \ vdots \ y_n endbmatrix = beginbmatrix 1 & x_1 & x_1^2 & cdots \ 1 & x_2 & x_2^2 & cdots \ 1 & x_3 & x_3^2 & cdots \ vdots & vdots & vdots & vdots \ 1 & x_n &x_n^2 & cdots endbmatrix beginbmatrix alpha_0 \ alpha_1 \ alpha_2 \ vdots \ alpha_m endbmatrix$$



            This can now be solved for your coefficients, $alpha_i$. That being said, and as was hinted at in your comments, I've never actually done this, and have instead used non-linear fitting functions provided by libraries.



            More information on polynomial regression on the wikipedia page.



            Edit: As you say in the comments, this method is only applicable if you can write your function that you wish to fit in polynomial form, which I don't think you can do for your example. In which case you are best off referring to the other answers to this question.






            share|cite|improve this answer











            $endgroup$








            • 2




              $begingroup$
              The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
              $endgroup$
              – Unstoppable Tachyon
              Mar 31 at 20:28










            • $begingroup$
              Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
              $endgroup$
              – Anon1759
              Mar 31 at 20:32










            • $begingroup$
              The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
              $endgroup$
              – Vladimir F
              Apr 1 at 13:04










            • $begingroup$
              @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
              $endgroup$
              – Vladimir F
              Apr 1 at 13:06










            • $begingroup$
              @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:37














            0












            0








            0





            $begingroup$

            Are you looking for something like polynomial regression? The general idea is, if you have measured pairs of (x, y(x)) and you are looking for find a fit of the form:



            $$y = alpha_0 + alpha_1 x + alpha_2 x^2 ...$$



            You can write this in matrix form as:



            $$beginbmatrix y_1 \ y_2 \ y_3 \ vdots \ y_n endbmatrix = beginbmatrix 1 & x_1 & x_1^2 & cdots \ 1 & x_2 & x_2^2 & cdots \ 1 & x_3 & x_3^2 & cdots \ vdots & vdots & vdots & vdots \ 1 & x_n &x_n^2 & cdots endbmatrix beginbmatrix alpha_0 \ alpha_1 \ alpha_2 \ vdots \ alpha_m endbmatrix$$



            This can now be solved for your coefficients, $alpha_i$. That being said, and as was hinted at in your comments, I've never actually done this, and have instead used non-linear fitting functions provided by libraries.



            More information on polynomial regression on the wikipedia page.



            Edit: As you say in the comments, this method is only applicable if you can write your function that you wish to fit in polynomial form, which I don't think you can do for your example. In which case you are best off referring to the other answers to this question.






            share|cite|improve this answer











            $endgroup$



            Are you looking for something like polynomial regression? The general idea is, if you have measured pairs of (x, y(x)) and you are looking for find a fit of the form:



            $$y = alpha_0 + alpha_1 x + alpha_2 x^2 ...$$



            You can write this in matrix form as:



            $$beginbmatrix y_1 \ y_2 \ y_3 \ vdots \ y_n endbmatrix = beginbmatrix 1 & x_1 & x_1^2 & cdots \ 1 & x_2 & x_2^2 & cdots \ 1 & x_3 & x_3^2 & cdots \ vdots & vdots & vdots & vdots \ 1 & x_n &x_n^2 & cdots endbmatrix beginbmatrix alpha_0 \ alpha_1 \ alpha_2 \ vdots \ alpha_m endbmatrix$$



            This can now be solved for your coefficients, $alpha_i$. That being said, and as was hinted at in your comments, I've never actually done this, and have instead used non-linear fitting functions provided by libraries.



            More information on polynomial regression on the wikipedia page.



            Edit: As you say in the comments, this method is only applicable if you can write your function that you wish to fit in polynomial form, which I don't think you can do for your example. In which case you are best off referring to the other answers to this question.







            share|cite|improve this answer














            share|cite|improve this answer



            share|cite|improve this answer








            edited Apr 1 at 6:52

























            answered Mar 31 at 20:24









            Anon1759Anon1759

            513




            513







            • 2




              $begingroup$
              The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
              $endgroup$
              – Unstoppable Tachyon
              Mar 31 at 20:28










            • $begingroup$
              Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
              $endgroup$
              – Anon1759
              Mar 31 at 20:32










            • $begingroup$
              The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
              $endgroup$
              – Vladimir F
              Apr 1 at 13:04










            • $begingroup$
              @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
              $endgroup$
              – Vladimir F
              Apr 1 at 13:06










            • $begingroup$
              @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:37













            • 2




              $begingroup$
              The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
              $endgroup$
              – Unstoppable Tachyon
              Mar 31 at 20:28










            • $begingroup$
              Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
              $endgroup$
              – Anon1759
              Mar 31 at 20:32










            • $begingroup$
              The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
              $endgroup$
              – Vladimir F
              Apr 1 at 13:04










            • $begingroup$
              @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
              $endgroup$
              – Vladimir F
              Apr 1 at 13:06










            • $begingroup$
              @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
              $endgroup$
              – Federico Poloni
              Apr 1 at 13:37








            2




            2




            $begingroup$
            The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
            $endgroup$
            – Unstoppable Tachyon
            Mar 31 at 20:28




            $begingroup$
            The answer is yes if the equation can be reduced to a polynomial one. I don't think it can be though.
            $endgroup$
            – Unstoppable Tachyon
            Mar 31 at 20:28












            $begingroup$
            Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
            $endgroup$
            – Anon1759
            Mar 31 at 20:32




            $begingroup$
            Then I think your only choice is to follow the advice as given in Anders Sandberg's answer and use one of the fitting techniques suggested there.
            $endgroup$
            – Anon1759
            Mar 31 at 20:32












            $begingroup$
            The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
            $endgroup$
            – Vladimir F
            Apr 1 at 13:04




            $begingroup$
            The formula with the Vandermonde matrix is for Linear interpolation, not for linear regeression. Or what am I missing?
            $endgroup$
            – Vladimir F
            Apr 1 at 13:04












            $begingroup$
            @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
            $endgroup$
            – Vladimir F
            Apr 1 at 13:06




            $begingroup$
            @AndreasMastronikolis You can always connect n+1 points with a Lagrange polynomial of degree n. But I doubt it makes much sense here.
            $endgroup$
            – Vladimir F
            Apr 1 at 13:06












            $begingroup$
            @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
            $endgroup$
            – Federico Poloni
            Apr 1 at 13:37





            $begingroup$
            @VladimirF The Vandermonde formula works also for linear regression. You just need to take the pseudoinverse $V^+ = (V^TV)^-1V^T$ of that (rectangular) matrix rather than its classical inverse, i.e., solve the system in the least-squares sense.
            $endgroup$
            – Federico Poloni
            Apr 1 at 13:37


















            draft saved

            draft discarded
















































            Thanks for contributing an answer to Physics Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid


            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.

            Use MathJax to format equations. MathJax reference.


            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fphysics.stackexchange.com%2fquestions%2f469754%2fhow-do-i-fit-a-resonance-curve%23new-answer', 'question_page');

            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown







            Popular posts from this blog

            Færeyskur hestur Heimild | Tengill | Tilvísanir | LeiðsagnarvalRossið - síða um færeyska hrossið á færeyskuGott ár hjá færeyska hestinum

            He _____ here since 1970 . Answer needed [closed]What does “since he was so high” mean?Meaning of “catch birds for”?How do I ensure “since” takes the meaning I want?“Who cares here” meaningWhat does “right round toward” mean?the time tense (had now been detected)What does the phrase “ring around the roses” mean here?Correct usage of “visited upon”Meaning of “foiled rail sabotage bid”It was the third time I had gone to Rome or It is the third time I had been to Rome

            Slayer Innehåll Historia | Stil, komposition och lyrik | Bandets betydelse och framgångar | Sidoprojekt och samarbeten | Kontroverser | Medlemmar | Utmärkelser och nomineringar | Turnéer och festivaler | Diskografi | Referenser | Externa länkar | Navigeringsmenywww.slayer.net”Metal Massacre vol. 1””Metal Massacre vol. 3””Metal Massacre Volume III””Show No Mercy””Haunting the Chapel””Live Undead””Hell Awaits””Reign in Blood””Reign in Blood””Gold & Platinum – Reign in Blood””Golden Gods Awards Winners”originalet”Kerrang! Hall Of Fame””Slayer Looks Back On 37-Year Career In New Video Series: Part Two””South of Heaven””Gold & Platinum – South of Heaven””Seasons in the Abyss””Gold & Platinum - Seasons in the Abyss””Divine Intervention””Divine Intervention - Release group by Slayer””Gold & Platinum - Divine Intervention””Live Intrusion””Undisputed Attitude””Abolish Government/Superficial Love””Release “Slatanic Slaughter: A Tribute to Slayer” by Various Artists””Diabolus in Musica””Soundtrack to the Apocalypse””God Hates Us All””Systematic - Relationships””War at the Warfield””Gold & Platinum - War at the Warfield””Soundtrack to the Apocalypse””Gold & Platinum - Still Reigning””Metallica, Slayer, Iron Mauden Among Winners At Metal Hammer Awards””Eternal Pyre””Eternal Pyre - Slayer release group””Eternal Pyre””Metal Storm Awards 2006””Kerrang! Hall Of Fame””Slayer Wins 'Best Metal' Grammy Award””Slayer Guitarist Jeff Hanneman Dies””Bullet-For My Valentine booed at Metal Hammer Golden Gods Awards””Unholy Aliance””The End Of Slayer?””Slayer: We Could Thrash Out Two More Albums If We're Fast Enough...””'The Unholy Alliance: Chapter III' UK Dates Added”originalet”Megadeth And Slayer To Co-Headline 'Canadian Carnage' Trek”originalet”World Painted Blood””Release “World Painted Blood” by Slayer””Metallica Heading To Cinemas””Slayer, Megadeth To Join Forces For 'European Carnage' Tour - Dec. 18, 2010”originalet”Slayer's Hanneman Contracts Acute Infection; Band To Bring In Guest Guitarist””Cannibal Corpse's Pat O'Brien Will Step In As Slayer's Guest Guitarist”originalet”Slayer’s Jeff Hanneman Dead at 49””Dave Lombardo Says He Made Only $67,000 In 2011 While Touring With Slayer””Slayer: We Do Not Agree With Dave Lombardo's Substance Or Timeline Of Events””Slayer Welcomes Drummer Paul Bostaph Back To The Fold””Slayer Hope to Unveil Never-Before-Heard Jeff Hanneman Material on Next Album””Slayer Debut New Song 'Implode' During Surprise Golden Gods Appearance””Release group Repentless by Slayer””Repentless - Slayer - Credits””Slayer””Metal Storm Awards 2015””Slayer - to release comic book "Repentless #1"””Slayer To Release 'Repentless' 6.66" Vinyl Box Set””BREAKING NEWS: Slayer Announce Farewell Tour””Slayer Recruit Lamb of God, Anthrax, Behemoth + Testament for Final Tour””Slayer lägger ner efter 37 år””Slayer Announces Second North American Leg Of 'Final' Tour””Final World Tour””Slayer Announces Final European Tour With Lamb of God, Anthrax And Obituary””Slayer To Tour Europe With Lamb of God, Anthrax And Obituary””Slayer To Play 'Last French Show Ever' At Next Year's Hellfst””Slayer's Final World Tour Will Extend Into 2019””Death Angel's Rob Cavestany On Slayer's 'Farewell' Tour: 'Some Of Us Could See This Coming'””Testament Has No Plans To Retire Anytime Soon, Says Chuck Billy””Anthrax's Scott Ian On Slayer's 'Farewell' Tour Plans: 'I Was Surprised And I Wasn't Surprised'””Slayer””Slayer's Morbid Schlock””Review/Rock; For Slayer, the Mania Is the Message””Slayer - Biography””Slayer - Reign In Blood”originalet”Dave Lombardo””An exclusive oral history of Slayer”originalet”Exclusive! Interview With Slayer Guitarist Jeff Hanneman”originalet”Thinking Out Loud: Slayer's Kerry King on hair metal, Satan and being polite””Slayer Lyrics””Slayer - Biography””Most influential artists for extreme metal music””Slayer - Reign in Blood””Slayer guitarist Jeff Hanneman dies aged 49””Slatanic Slaughter: A Tribute to Slayer””Gateway to Hell: A Tribute to Slayer””Covered In Blood””Slayer: The Origins of Thrash in San Francisco, CA.””Why They Rule - #6 Slayer”originalet”Guitar World's 100 Greatest Heavy Metal Guitarists Of All Time”originalet”The fans have spoken: Slayer comes out on top in readers' polls”originalet”Tribute to Jeff Hanneman (1964-2013)””Lamb Of God Frontman: We Sound Like A Slayer Rip-Off””BEHEMOTH Frontman Pays Tribute To SLAYER's JEFF HANNEMAN””Slayer, Hatebreed Doing Double Duty On This Year's Ozzfest””System of a Down””Lacuna Coil’s Andrea Ferro Talks Influences, Skateboarding, Band Origins + More””Slayer - Reign in Blood””Into The Lungs of Hell””Slayer rules - en utställning om fans””Slayer and Their Fans Slashed Through a No-Holds-Barred Night at Gas Monkey””Home””Slayer””Gold & Platinum - The Big 4 Live from Sofia, Bulgaria””Exclusive! Interview With Slayer Guitarist Kerry King””2008-02-23: Wiltern, Los Angeles, CA, USA””Slayer's Kerry King To Perform With Megadeth Tonight! - Oct. 21, 2010”originalet”Dave Lombardo - Biography”Slayer Case DismissedArkiveradUltimate Classic Rock: Slayer guitarist Jeff Hanneman dead at 49.”Slayer: "We could never do any thing like Some Kind Of Monster..."””Cannibal Corpse'S Pat O'Brien Will Step In As Slayer'S Guest Guitarist | The Official Slayer Site”originalet”Slayer Wins 'Best Metal' Grammy Award””Slayer Guitarist Jeff Hanneman Dies””Kerrang! Awards 2006 Blog: Kerrang! Hall Of Fame””Kerrang! Awards 2013: Kerrang! Legend”originalet”Metallica, Slayer, Iron Maien Among Winners At Metal Hammer Awards””Metal Hammer Golden Gods Awards””Bullet For My Valentine Booed At Metal Hammer Golden Gods Awards””Metal Storm Awards 2006””Metal Storm Awards 2015””Slayer's Concert History””Slayer - Relationships””Slayer - Releases”Slayers officiella webbplatsSlayer på MusicBrainzOfficiell webbplatsSlayerSlayerr1373445760000 0001 1540 47353068615-5086262726cb13906545x(data)6033143kn20030215029