清空筛选条件
清空筛选条件

Estimating the period and the damping of the system

3 次查看（过去 30 天）
Obsidian
Obsidian 2021-7-19
评论： Obsidian 2021-7-19
Hello,
I am new with MatLab. I want to Estimate the period and the damping of the system from a set of data using three values of j cycles:
j = 10 cycles, 20 cycles, 30 cycles
For each case c.I, c.II, and c.III I am supposed to use different starting peak (i.e., moving window of peaks) and calculate the median value of damping ratio and natural period for each case C.I, C.II, and C.III. and also tabulate the results.
Does anyone have a suggestion about how to approad this problem?
I appreciate the help!
I used following code to extract and plot data and I am stuck on how to find period and damping:
code:
clc; clear all; close all
load Vibration_Data.txt
acc = Vibration_Data(:,2);
dt = 0.02;
time = (0:numel(acc)-1) * dt;
plot (time,acc,'b-');
findpeaks(acc,time)
peaks=findpeaks(acc,time)
period=time(peaks)-time(peaks(0:size(peaks)-1))
error:
last line: array indices must be positive integers or logical values.
  2 个评论
Sameer Pujari
Sameer Pujari 2021-7-19
编辑：Sameer Pujari 2021-7-19
As per my understanding, you are trying to calculate period(s) between two peaks.
You could use numel(peaks) or length(peaks) to find number of elements in the peaks array.
Then use for loop to generate an array of period.
for i=1:length(peaks)-1
period(i)=time(i+1)-time(i)
end
In matlab, Array indexing starts from index '1'. index '0' is not valid. So, you might be getting the error

请先登录，再进行评论。

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

回答（1 个）

Simon Chan
Simon Chan 2021-7-19
This line will give an error:
period=time(peaks)-time(peaks(0:size(peaks)-1))
So just try the following:
period=time(peaks)-time(peaks-1)

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by