How to get positive and negative area under the curve separately in simulink ?

I got the following output in simulink. I want seperate the positive and negative area under the curve separately and also want to colour it with diffrent colours.

Answer 1

[Mathieu NOE]

hello

see on FEX :

climanomaly - File Exchange - MATLAB Central (mathworks.com)

save your simulink data as array to the workspace and use the climanomaly function

exemple :

x = 1:.1:20;
y = sin(x);
ref = 0;
figure
climanomaly(x,y,ref);

[Swaraj Shevgan]

but i still cant figure out about how to do for x - axis as it is the simulation time. When i try to do it , it says dimensions aren't matching

[Mathieu NOE]

can you share your code and data from simulink ?

[Swaraj Shevgan]

This link is my reference model . Here i want to check for used battery power and regenerated power.I want to know the values seperate and also want to colour the positive and negative part seperate.Thank you .

https://www.mathworks.com/matlabcentral/fileexchange/82250-battery-electric-vehicle-model-in-simscape

[Mathieu NOE]

hello

I am not a Simscape user

could you provide the data from the simulation ?

[Swaraj Shevgan]

I have uploaded the excel file of all the values of power obtained from the simulation. Its the plot of power vs time , where the time ranges from 0 to 1180 seconds.

I want seperate the positive and negative area under the curve separately and also want to get the separate values and also colour it with diffrent colours. Following is the curve for your reference. Thank you alot for helping.You can also contact me on Email: [e-mail]

[Mathieu NOE]

hello again

this is what you can achieve with this code.

Others options are possible - edit the climanomaly function to see other coding alternatives .

clc
clearvars
y = readmatrix('New Microsoft Excel Worksheet.xlsx');
dt = 1180/(length(y)-1); % time ranges from 0 to 1180 seconds

% resampling (decimation) - optionnal , use it at your own convenience
% benefit : makes the below code runs faster without compromising the
% plot quality
R = 20;
y = decimate(y,R);

samples = length(y);
dt = dt*R;
x = (0:samples-1)*dt;
ref = 0;
figure
climanomaly(x,y,ref);
xlabel('Time (s)');
ylabel('Power (W)');

[Swaraj Shevgan]

Hello,

Wow , Thank you very much. Can you also help me with getting the area values of both positive and negative part separate?

[Mathieu NOE]

hello

this is basically the data splitted in positive and negative components plotted separately

clc
clearvars
y = readmatrix('New Microsoft Excel Worksheet.xlsx');
dt = 1180/(length(y)-1); % time ranges from 0 to 1180 seconds

% resampling (decimation) - optionnal , use it at your own convenience
% benefit : makes the below code runs faster without compromising the
% plot quality
R = 20;
y = decimate(y,R);

samples = length(y);
dt = dt*R;
x = (0:samples-1)*dt;

y_pos = zeros(samples,1) ;
ind = find(y>=0);
y_pos(ind) = y(ind);

y_neg = zeros(samples,1) ;
ind = find(y<0);
y_neg(ind) = y(ind);

ref = 0;
figure(1)
subplot(211),climanomaly(x,y_pos,ref);
xlabel('Time (s)');
ylabel('Power (W)');
subplot(212),climanomaly(x,y_neg,ref);
xlabel('Time (s)');
ylabel('Power (W)');

[Swaraj Shevgan]

Hello,

Thank you very much . This will help me alot in my simulation.

[Mathieu NOE]

My pleasure