How can I device a time series A by a Time series B if they are the same size?

I have 2 time series, they are both the same size but I don't know how to devide time series A by time series B.

Also, Once is A/B. How can I get the max value regardless of time?

Thanks

Answer 1

[thoughtGarden]

Time series division is the same as general variable division

>> timeS1 = timeseries(linspace(1,2));
>> timeS2 = timeseries(linspace(1,5));
>> newSeries = timeS1/timeS2; %this works...

to find the max value, use the max() function

timeSeriesMax = max(newSeries);

[Diego Dranuta]

I am a little conffuesed when looking ar this example.

Do you mind putting it with the code I am using?

In my case I have 2 time series out.Vehicle_Speed_mph and out.Driver_Speed_Desired_mph

i have to find the percent difference between them. What I did first was to synchronize the time since they were off. Then to find the difference at each second between them. After that, I Devided the difference by one of the time series. and then I took the max value of that and multi[plied by 100 to get the percentage.

However I am doing something wrong because i get inf for the answer.

[out.Vehicle_Speed_mph,out.Driver_Speed_Desired_mph]=synchronize(out.Vehicle_Speed_mph,out.Driver_Speed_Desired_mph,'intersection')
Speed_Difference=(out.Vehicle_Speed_mph-out.Driver_Speed_Desired_mph)
a=Speed_Difference/out.Driver_Speed_Desired_mph;
b=max(a);
Maximum_Percent_Driver_Error_relative_to_Desired_Driving_Cycle=100*b

Thanks

[thoughtGarden]

Do you want the average difference or the difference at each sample time?

[Diego Dranuta]

the difference at each sample time, so I can take the max difference possible

[thoughtGarden]

Below is an example using your variable names.

% Build sample data using OP variable names
t = 1:.001:10;

out.Vehicle_Speed_mph = timeseries(abs(sin(t))+.1);
out.Driver_Speed_Desired_mph = timeseries(abs(cos(t))+.1);

% Compute Percent of Desired that has been achieved
percentOfDesired = 100.*(out.Vehicle_Speed_mph - out.Driver_Speed_Desired_mph)/out.Driver_Speed_Desired_mph;

% Plot Results
figure(1);clf;
subplot(2,1,1)
plot(out.Vehicle_Speed_mph,'displayName','Speed Achieved');
hold on;
plot(out.Driver_Speed_Desired_mph,'displayName','Speed Desired');
legend('show');
grid on;
title('Both Data Sets')

subplot(2,1,2)
plot(percentOfDesired,'displayName','Percent of Desired')
legend('show')
grid on;
title('The Percentage of Desired')

[Diego Dranuta]

in this case it would be if I am looking fot viualization of the difference. I just want the maximum percent difference in a scalar

[thoughtGarden]

Then simply use the max() function.

% Build sample data using OP variable names
t = 1:.001:10;

out.Vehicle_Speed_mph = timeseries(abs(sin(t))+.1);
out.Driver_Speed_Desired_mph = timeseries(abs(cos(t))+.1);

% Compute Percent of Desired that has been achieved
percentOfDesired = 100.*(out.Vehicle_Speed_mph - out.Driver_Speed_Desired_mph)/out.Driver_Speed_Desired_mph;

% Plot Results
figure(1);clf;
subplot(2,1,1)
plot(out.Vehicle_Speed_mph,'displayName','Speed Achieved');
hold on;
plot(out.Driver_Speed_Desired_mph,'displayName','Speed Desired');
legend('show');
grid on;
title('Both Data Sets')

subplot(2,1,2)
plot(percentOfDesired,'displayName','Percent of Desired')
legend('show')
grid on;
title('The Percentage of Desired')

disp(['The Maximum percentage difference: ' num2str(max(percentOfDesired)) '%']);

[thoughtGarden]

If you found this answer to be correct, please accept it. Otherwise, please add aditional comments so that it may be improved.