Multiply and integrate two plots for single output value

Below is illustrated the code used to compute functions, with a comment at the bottom summarising the aim. To obtain a single value for annual energy output from two plots.

%% Clear Workspace 
clc; clear all; clear;

% Define Parameters 
K_0 = 3149/1080;
K_1 = 469/1080;
T_1 = 12.4/24; %hours 
T_0 = 353/24;  %hours 

% Define Timestep
ts_min = 2; %mins
ts_hours = ts_min/60;
ts_days = ts_hours/24;
ts_years = ts_days/365;

start_val = 0; % days 
end_val = 365; % days  
nElements = 1/ts_years; % intervals in a year
ts = linspace(start_val, end_val, nElements); % create timestep array

% Define Harmonic Analysis
V = (K_0 + K_1.*cos((2.*pi.*ts)/T_1)).*cos((2.*pi.*ts)/T_0);
V2 = V.^2;
Vnorm = sqrt(V2);
plot(ts,Vnorm);
grid on;

% Plot Histogram of Wind Velocity Frequencies 
hist(Vnorm);
grid on;

%% Tidal Power Curves for 4m and 6.5m from Excel 
PowerTable = readtable('TidalData.xlsx','Sheet','Sheet1');
x4 = PowerTable.Header1; % Flow Speed [m/s] at 4m turbine
y4 = PowerTable.Header2; % Electrical Power [kW] at 4m turbine
x63 = PowerTable.Header3; %: Flow Speed [m/s] at 6.3m turbine
y63 = PowerTable.Header4; %: Electrical Power [kW] at 6.3m turbine

figure;
plot(x4,y4,'DisplayName','4m');
hold on;
grid on;
plot(x63,y63,'DisplayName','6.3m');
hold off;
legend

%% Multiplication Calculation power curve by yearly quantity
% aim to multiply plot(x4,y4)- flow speed (m/s) vs. electrical power (kW),
% by hist(Vnorm)- flow speed (m/s) vs. annual number of hours, to get the total annual
% electrical power kWh in a year..

[Dyuman Joshi]

What is the question? What is the problem you are facing?

[Ben Bawtree]

Calculate the energy captured of a tidal turbine per annum (kWh/yr)?

I have plotted the distrubution of velocities over a year, and then formed a histogram plot, for a distrubution of the frequencies of occurence for each velocity. Then plotted a power curve of the output power against flow speed.

Then the issue i am facing it what function/plot to use to combine said plots (histogram - frequency distrubution of tidal velocities & power curve of flow speed vs. tidal power), to obtain a singular value for the total annual power. Hope that makes sense.

[Mathieu NOE]

hello

first comment, you did not provide the excel file needed to run your code

second comment, why would we need the histogram distribution if you know the velocity all year long and you know the turbine power output vs velocity.

Combining both data will give you the power output every time increment (hour ?) of the year. and the end you just have to sum all those numbers to have the total kWh in one year.

[Ben Bawtree]

Understood thank you, will attempt this.

Answer 1

[Peter Perkins]

Any time the number 365 appears in code, thatr's probably a bug.

You might try to use datetimes and durations in your code, including the plot.