随机产生行和列为固定值的矩阵(附MATLAB源码)

        今天突然需要用一个算法来产生随机矩阵,且矩阵的行和列的和都是固定值。然后百度了很久发现没什么有价值的信息,于是就准备自己动手,说干就干,经过两个小时的思考,想到了一个不是很好的方法,可以产生符合条件的矩阵,也能保证随机性。在这里做个记录,以便给有需要的同学提供写思路,要是有更好的方法和建议,麻烦也共享出来起学习。

        话不多说,直接来介绍我的思路吧:

(1)假设矩阵维度是mxn,行和列的和分别要求为M,N。

(2)产生一个mxn随机矩阵,这个在MATLAB上太容易了,各种概率分布自己选。

(3)对矩阵每一行规划化处理,保证每行的和为1

(4)基于上一步,将矩阵的每一行的每个元素乘以M,则可以得到所有行的和为M的矩阵了,注意,这一步有个需要说明的是,如果要求产生的矩阵所有数都是整数,那么这里就需要对每行的元素再进行修正。修正的方式很简单,判断每行的和,如果大于M,则从这行元素中最大的那个开始递减,直至这样的和等于M。

(5)接下来就是对列进行操作了,每列一次判断是否等于N,,如果不是,再区分大于还是小于N,

如果大于N,这将这列从大到小递减,且将减掉的数再加到所在行的其他列。一个while解决列的修正处理。跳出循环后,整个矩阵就满足要求了。

(6)基于上述策略,即可保证最初产生的随机矩阵与后面修正的矩阵是一一对应的,这样间接保证了最终符合要求的矩阵的随机性

        代码来了:

% sum_row每行的和,可以是一个数,也可以是一个向量
% sum_col每列的和,可以是一个数,也可以是一个向量
%本算例约束条件:
%(1)每一行的和都满足一个固定的值
%(2)每一列的和都满足一个固定的值
%(3)每所有行的和与所有列的和相等
n=length(sum_col);%列数
m=length(sum_row);%行数
x=rand(m,n);%产生随机矩阵
xout=zeros(m,n);%最终要输出的满足要求的矩阵
%行修正
for i=1:m
    %得到行的和一定
    xout(i,:)=round(sum_row(i).*x(i,:)./sum(x(i,:)));
    %修正行的和大于指定值的情况
    while sum(xout(i,:))>sum_row(i)
        [~,ind]=max(xout(i,:));
        xout(i,ind)=xout(i,ind)-1;
    end
    %修正行的和小于指定值的情况
    while sum(xout(i,:))<sum_row(i)
        [~,ind]=min(xout(i,:));
        xout(i,ind)=xout(i,ind)+1;
    end
end
%修正列
while any(sum(xout)~=sum_col)
    suminfo=sum(xout);
    %不满足约束的列
    black=find((suminfo>sum_col)==1);
    %满足约束的列
    white=find((suminfo>sum_col)==0);
    while ~isempty(black)
        col=black(1);
        de=suminfo(col)-sum_col(col);
        list=xout(:,col)';
        [~,ord]=sort(list,'descend');
        for j=1:length(ord)
            if list(ord(j))>=de
                list(ord(j))=list(ord(j))-de;
                for k=1:length(white)
                    de1=abs(suminfo(white(k))-sum_col(white(k)));
                    if de1>=de
                        xout(ord(j),white(k))=xout(ord(j),white(k))+de;
                        de=0;
                        break;
                    else
                        xout(ord(j),white(k))=xout(ord(j),white(k))+de1;
                        de=de-de1;
                    end
                end
                break;
            else
                dec=list(ord(j));
                for k=1:length(white)
                    de1=abs(suminfo(white(k))-sum_col(white(k)));
                    if de1>=dec
                        xout(ord(j),white(k))=xout(ord(j),white(k))+dec;
                        dec=0;
                        break;
                    else
                        xout(ord(j),white(k))=xout(ord(j),white(k))+de1;
                        dec=dec-de1;
                    end
                end
                de=de-(list(ord(j))-dec);
                list(ord(j))=dec;
            end
        end
        xout(:,col)=list';
        suminfo=sum(xout);
        black=find((suminfo>sum_col)==1);
        white=find((suminfo>sum_col)==0);
    end
end
if any(sum(xout,2)~=sum_row) || any(sum(xout,1)~=sum_col)
    disp('行或者列的和不满足要求')
else
    disp('nice^_^')
    xout
end

部分结果截图:

 PS:这里其实还有漏洞,发现了的评论区走起^_^!!!文章来源地址https://www.uudwc.com/A/b1o1R/

原文地址:https://blog.csdn.net/qq_23887115/article/details/128279527

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请联系站长进行投诉反馈,一经查实,立即删除!

h
上一篇 2023年09月26日 15:31
Android Key/Trust Store研究+ssl证书密钥
下一篇 2023年09月26日 15:32