Fitting a damping function to a sinusoidal response
显示 更早的评论
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
2020-1-31
0 个投票
Try using findpeaks(), and then just pass the peak locations into fitnlm(), like in my attached demo.
Star Strider
2020-1-31
0 个投票
类别
在 帮助中心 和 File Exchange 中查找有关 Measurements and Feature Extraction 的更多信息
2020-1-31
2020-1-31
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
