Combine two different matrices to create one new matrix of all matrix combinations

1 次查看(过去 30 天)
Hi, I have the following matrices
WT_energy_supply(:,:,M) and PV_power_output(:,:,K)
which are 365 x 24 matrices and M = 3 and K = 10. Therefore there are 30 possible combinations of WY_energy_supply and PV_power_output. I am looking to create a set of matrices (in the form (:,:,Z)) that add corresponding elements in each matrix to form this new set (again 365x24 with 30 different matrices) Does anyone have any idea where to start coding this?
Thank you

采纳的回答

Sean de Wolski
Sean de Wolski 2012-5-11
Nice brain teaser!
%Sample data:
A = ones(4,4,5);
B = bsxfun(@times,ones(4,4,3),cat(3,1,2,3));
%Nice to know:
sz = size(B);
%Engine:
C = reshape(bsxfun(@plus,A, reshape(B,sz(1),sz(2),1,sz(3))),sz(1),sz(2),[]);
  18 个评论

请先登录,再进行评论。

更多回答(2 个)

Andrei Bobrov
Andrei Bobrov 2012-5-11
eg:
M = 3;
K = 10;
WT_energy_supply = randi(10,10,4,M);
PV_power_output = randi([10 50],10,4,K);
%solution
ij = fullfact([M K]);
out = WT_energy_supply(:,:,ij(:,1)) + PV_power_output(:,:,ij(:,2)) ;
  5 个评论
Andrew Alkiviades
Andrew Alkiviades 2012-5-11
having tried this code i get the following error "Array dimensions must match for binary array op."
Do I need to define M and K as 3 and 10 since I have them previously in the code for working out WT_energy_supply and PV_power_output?
This may be a stupid question but why can't I just add WT_energy_supply(:,:,M) to PV_power_output(:,:,K) since this is really what I am trying to do?
Thank you
Sean de Wolski
Sean de Wolski 2012-5-11
That's what bsxfun does for you!
think about it, how can you add two things that are different sizes?

请先登录,再进行评论。


owr
owr 2012-5-11
To generate all combinations of K and M, check out the file exchange function "allcomb": href=""<http://www.mathworks.com/matlabcentral/fileexchange/10064</a>>
Ive been using it for years, works with more than 2 parameters as well.

类别

Help CenterFile Exchange 中查找有关 Matrices and Arrays 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by