Fitting a damping function to a sinusoidal response

7 次查看（过去 30 天）
Diogo Carvalho
Diogo Carvalho 2020-1-31
回答： Star Strider 2020-1-31
I'm attempting to fit an exponential damping function to a sinusoidal wave by iterating zeta values.
However, it seems I can never get it to properlly 'ride' along the peaks, meaning it always crosses the response function (pic)
As I understand, the damping function looks like exp(-zeta*wn*t). I have estimated the value of wn and I have all values of time.
The function code looks like this:
xi1 = 0.01;
y = exp(-xi1*wn1*(time1-4.3))-0.6568;
plot(time1-4.3,Response1)
hold on
plot(time1-4.3,y)
hold off
grid on
I subtract 0.6568 to have the damping function begin on the first peak.

请先登录，再回答此问题。

回答（2 个）

Image Analyst
Image Analyst 2020-1-31
Try using findpeaks(), and then just pass the peak locations into fitnlm(), like in my attached demo.
Star Strider
Star Strider 2020-1-31

类别

在 Help Center 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by