program to create matrix obeying constraints

1 次查看(过去 30 天)
Hello all,
i have following code to create n*n matrix(say n=5 this time) which is following certain constraints.
function selectingset
% n = input ('Enter Num_Nodes :');
% inputData=input ('Enter inputdata :');
n=5;
inputData=[5 1 2];
for z = 1:n
f = inputData(1,1); % f= Demand
S = inputData(1,2); % S= Source node
T = inputData(1,3); % T= Termination node
[idx,p] = sort(rand(n-1,f));
M = accumarray([reshape(p+(p>=T),[],1),repmat([1:S-1,S+1:n]',f,1)],1,[n,n]);
M(1:n+1:n^2) = 0;
end
disp(M);
end
So, my o/p from following code is---
0 0 2 3 0
0 0 0 0 0
0 0 0 0 3
0 4 1 0 0
0 1 0 2 0
Note that as in given code
input S=1, so column 1 is zero and
As T=2, so row 2 is zero.
As f= 5, so sum of row S(i.e. S=1)=Sum of column T(i.e. T=2)= 5
and i am making diagonal elements zero.
all the remaining rows and columns, sum of row(i)=sum of column(i)
Now, How can i modify code such that row(1,3) will always come zero in output but output still follows previous mentioned constraints.
  4 个评论
Matt J
Matt J 2014-9-20
Do all the matrix entries have to be integers?
reshdev
reshdev 2014-9-20
Hello,
S will never be equal to T.
row(1,3) means third element in first row.
all the matrix elements should be integers between 0 to f.

请先登录,再进行评论。

回答(1 个)

Matt J
Matt J 2014-9-20
编辑:Matt J 2014-9-20
If you have the latext version of the Optimization Toolbox, you can use intlinprog to solve for the unknown elements of the matrix. You can choose any objective function that you want. See this recent thread for a very similar problem related to constraining rows and columns of matrices:
This assumes that you require the entries to be integer-valued and that a solution exists, see also My Comments. If they don't have to be integer-valued, you could just use linprog .
  2 个评论
reshdev
reshdev 2014-9-20
Hello,
S will never be equal to T.
row(1,3) means third element in first row.
all the matrix elements should be integers between 0 to f.
Matt J
Matt J 2014-9-20
OK. Then I guess I've answered your question.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Genetic Algorithm 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by