Approximate experimental data from multiple cells

Hello. I have several cell arrays with probe data from .mat-file.
I need to find resonance frequencies and q-factors for each of the 1000 inputs, and then estimate (approximate) them. However, the code below provides only 1 value of fres instead of 1000, and no plot for (i, fres) is shown. (I attached a frequency-amplitude plot)
  • FStart - cell array (1000*1)
  • FEnd - cell array (1000*1)
  • Amp - multiarray (1000*100)
for i = 1:1:1000
f1 = FStart(i):10:FEnd(i);
grid on
y1 = plot(f1,Amp(i,:));
[maxValue, maxIndex] = max(Amp(i,:)); %find maximum value of amplitude for each i
[Q_Value, Q_Index] = max(0.5*AmpA2(i,:)); %also tried 0.5*max() and /2
fres = f1(maxIndex); %by index of max amplitude value find resonance frequency
plot(i,fres) %plot resonance frequency for each i
hold on
end
NB: the last FStart value is less than the first FreqEnd value
Futhermore, I try to esteem the Q factor as: Max(frequency on level = 1/2*MaxAmplitude)-Min(frequency on level = 1/2*MaxAmplitude)
fmin = min(f1(Q_Index))
fmax = max(f1(Q_Index))
But it shows fmin = fmax
Could you please tell, what's the problem here?

3 个评论

YOu have not shown what is f2 ? And the question is based on f2.
Excuse me, that's f1, too, I'll edit it.
You say that FStart and Fend are cell arrays, yet your code has
f1 = FStart(i):10:FEnd(i);
which would cause an error if they were indeed cell arrays. Are they actually just plain matrices (in this case vectors)?

请先登录,再进行评论。

回答(1 个)

But it shows fmin = fmax Well, of course. I don't understand what you're trying to do but you have:
[Q_Value, Q_Index] = max(0.5*AmpA2(i,:))
fmin = min(f1(Q_Index))
fmax = max(f1(Q_Index))
so AmpA2(i, :) is a column vector from which you're taking the maximum and its location Q_index. There is only ever one Q_Index value returned by max. Therefore f1(Q_Index) is a scalar value. The minimum and maximum of a single value is the same: that single value.
Also: see comment to the question, above, about the cell arrays.

1 个评论

Thank you for the answer. I understand that there was 1 value, I pointed that out. fres(i) instead of fres solved my problem. I was trying to get an array for fres.

请先登录,再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Particle & Nuclear Physics 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by