MATLAB: How to obtain the value of a vector when a condition is met for multiple times in a row

multiple conditions

Hello,

How can I get the position of averror when the condition of <10^-4 is met 4 times in a row for the first time?

The code below shows the problem. In other words I want to obtain the value of e when the next 3 value of e are also meeting the condition above.

[e]=find(averror<(10^-4))
for k=1:size(e)
    if e(k)==e(k+1)-1 && e(k)==e(k+2)-2 && e(k)==e(k+3)-3
    m=e(k);
    end
end

Does anybody know this?

Best Answer

e = averror < (10^-4);
m = strfind([false e(:).'], [0 1 1 1 1]);
m(1)

Related Solutions

MATLAB: How to manage a .csv imported file

mask = YourSignal(:).' > 0;
starts = strfind([false mask], [0 1]);
stops = strfind([mask false], [1 0]);

Now starts and stops are vectors such that you have non-zero signal between YourSignal(starts(K):stops(K))

This assumes that none of your signal is negative.

MATLAB: Fast Fourier Transform with discrete data

I suspect the first column of ‘data’ is a time vector and the second column is the signal vector.

If that is correct, try this:

t = data(:,1);                                      % Time Vector
s = data(:,2);                                      % Signal Vector
Ts = mean(diff(t));                                 % Sampling Interval
Fs = 1/Ts;                                          % Sampling Frequency
Fn = Fs/2;                                          % Nyquist Frequency
L = numel(t);                                       % Signal Length (Number Of Samples)
FTs = fft(s - mean(s))/L;                           % Fourier Transform
Fv = linspace(0, 1, fix(L/2)+1)*Fn;                 % Frequency Vector
Iv = 1:numel(Fv);                                   % Index Vector
figure
plot(Fv, abs(FTs(Iv))*2)
grid
xlabel('Frequency')
ylabel('Amplitude')

This calculates the Fourier transform, and displays a one-sided Fourier transform (0 Hz to the Nyquist frequency).

This code subtracts the mean of ‘s’ from ‘s’ to remove a constant (d-c) offset (that could obscure other peaks), before taking the Fourier transform.

The sampling interval must be constant, or this code will not give correct results. You can determine that by taking the standard deviation of the time intervals:

Ts_sd = std(mean(diff(t)));

It should be close to 0 (usually less than 1E-12).

Best Answer

Related Solutions

MATLAB: How to manage a .csv imported file

MATLAB: Fast Fourier Transform with discrete data

Related Question