Matching columns of arrays to fill up a corresponding empty column

Question

Antonios Kouzelis 2022-5-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1715655-matching-columns-of-arrays-to-fill-up-a-corresponding-empty-column

编辑： Voss 2022-5-17

Hi,

While determining the battery voltage values from battery state of charge values, I'm encountering the following issue:

I have one array of 10.001 rows and 2 colums (10.001x2), giving percentage in steps of 0.01 percent

and one array of 604.800 rows and 2 columns (604.800x2), giving values of the voltage.

The first array describes the relation between state of charge and the voltage potential in the battery. The left column goes from 0 to 100 percent, in steps of 0.01 percent, and the right column gives corresponding voltage values that range from 2.8015 V to 4.1458 V.

0.00 2.8015

0.01 2.8015

0.02 2.8015

... ...

99.98 4.1444

99.99 4.1548

100.00 4.1548

My second array, in the first column, has 604,800 rows of percentage values (one week of state of charge values). The left column starts at 90%, goes up and down a few times, ending up close to 90% again. The second column of this array is empty.

90.00 [empty]

... ...

89.99 [empty]

How can I fill up the second column of the second array ([empty] values), with the values that correspond to those percentages from the first array?

Thanks in advance.

Best, Antonios

5 个评论
显示 3更早的评论隐藏 3更早的评论

Antonios Kouzelis 2022-5-11

SOC.mat

Here it is!

Antonios Kouzelis 2022-5-12

@Image Analyst

Dear ImageAnalyst,

You were able to be of excellent help last time I posted a question on the forum.

Would you be able to help me once more?

Thanks in advance.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Voss 2022-5-13

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1715655-matching-columns-of-arrays-to-fill-up-a-corresponding-empty-column#answer_963410

在 MATLAB Online 中打开

load('OCV_new.mat')
load('SoC_internal_new.mat')
load('SOC.mat')

It looks like SoC_internal_new is not values between 0 and 100 percent but between 0 and 1:

min(SoC_internal_new)
ans = 0
max(SoC_internal_new)
ans = 1

And SOC is in percent:

min(SOC)
ans = 20.7939
max(SOC)
ans = 90

So one of them must be converted in order to compare them. I'll convert SoC_internal_new to percent and construct the first array described in the question:

input_data = [100*SoC_internal_new(:) OCV_new(:)]
input_data = 10001×2
         0    2.8015
0100    2.8015
0200    2.8015
0300    2.8015
0400    2.8015
0500    2.8015
0600    2.8015
0700    2.8015
0800    2.8015
0900    2.8015

Then the second array described in the question has SOC in the first column and interpolated values in the second column:

output_data = [SOC(:) interp1(input_data(:,1),input_data(:,2),SOC)]

output_data = 604800×2

90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878 90.0000 4.0878

% plot both arrays

plot(input_data(:,1),input_data(:,2))

hold on

% SOC (i.e., output_data(:,1)) is unsorted (i.e.,

% it goes up and down, as described in the question),

% so sort output_data according to the values of

% output_data(:,1), for plotting

[~,idx] = sort(output_data(:,1));

plot(output_data(idx,1),output_data(idx,2),'--','LineWidth',2)

2 个评论
显示无隐藏无

Antonios Kouzelis 2022-5-16

Thanks man, with the help of a collegue I ended up doing it as follows:

n = numel(Data.SOC);

soc_voltage_data = zeros(n,2);

for i =1:n

temp1 = Data.SOC(i,1);

temp2 = SoC_internal_new.'*100;

temp3 = temp2 - temp1;

abs_error = abs(temp3);

% ind_min = find(temp3==0);

[~,ind_min] = min(abs_error);

voltage = OCV_new(1,ind_min);

if numel(voltage) == 0

break

end

soc_voltage_data(i,1) = temp1;

soc_voltage_data(i,2) = voltage;

end

Best, Antonis

Voss 2022-5-17

编辑：Voss 2022-5-17

在 MATLAB Online 中打开

You're welcome!

Note that my solution and yours are different, in that I'm doing linear interpolation and you're doing nearest-neighbor (minimum absolute difference) lookup. I don't know if that matters in this case, since the question didn't say which method to use and the data is finely enough spaced that you can't see any difference between the results from the two methods. I just figured I'd point that out in case it matters for your purpose.

However, if you want to use nearest-neighbor, you can do it with the interp1 function much faster than with a for loop (~1000x faster in this case):

Data = load('SOC');
load('OCV_new.mat')
load('SoC_internal_new.mat')
tic
n = numel(Data.SOC);
soc_voltage_data = zeros(n,2);
for i =1:n
    temp1  =  Data.SOC(i,1);
    temp2 = SoC_internal_new.'*100;
    
    temp3 = temp2 - temp1;
    abs_error = abs(temp3);
    
%     ind_min = find(temp3==0);
    [~,ind_min] = min(abs_error);
    voltage = OCV_new(1,ind_min); 
    
    if numel(voltage) == 0
        break
    end
    
    soc_voltage_data(i,1) = temp1;
    soc_voltage_data(i,2) = voltage;
end
toc
Elapsed time is 13.790107 seconds.
tic
soc_voltage_data_test = [Data.SOC ...
    interp1(100*SoC_internal_new(:),OCV_new(:),Data.SOC,'nearest')];
toc
Elapsed time is 0.015937 seconds.
% both ways give the same result
isequal(soc_voltage_data_test,soc_voltage_data)
ans = logical
   1