草庐IT

python - 在非常具体的约束下生成随机数

coder 2023-08-15 原文

我面临以下编程问题。我需要生成 n (a, b) 元组，所有 a 的总和是给定的 A 和所有 b 的总和是给定的 B 并且对于每个元组，a/b 的比率在 范围内>(c_min，c_max)。 A/B 也在同一范围内。我还试图确保除了约束引入的结果之外没有任何偏差，并且 a/b 值或多或少均匀分布在给定范围内。

一些说明和元约束:

    给出了
  • A、B、c_min 和c_max。
  • 比率 A/B 在 (c_min, c_max) 范围内。如果问题是在给定其他约束的情况下找到解决方案，则必须如此。
  • a 和 b 是 >0 且非整数。

我正在尝试用 Python 实现它，但非常感谢任何语言(包括英语)的想法。

最佳答案

我们寻找满足 a_i 和 b_i 的元组

  • (a_1, ... a_n) 和 (b_1, ... b_n) 具有在指数排列下不变的分布(您称之为“无偏”)
  • 比率 a_i/b_i 在 [cmin, cmax] 上均匀分布
  • sum(a_i) = A, sum(b_i) = B

如果 c_min 和 c_max 条件不是太差(即它们彼此不是很接近)，并且 n 不是很大，以下作品:

  • “均匀地”生成 a_i 使得 sum a_i = A:
    • 从某些分布(例如均匀分布)中抽取 n 个样本 aa_i (i = 1..n)
    • 将它们除以它们的总和并乘以 A:a_i = A * aa_i/sum(aa_i) 具有所需的属性。
  • 生成 b_i 使得 sum b_i = B 用同样的方法。
  • 如果存在 i 使得 a_i/b_i 不在区间 [cmin, cmax] 中，丢弃所有 a_i 和 b_i 并从头开始重试。

它不能与n 很好地扩展，因为a_i 和b_i 的集合满足约束随着 n 的增加，变得越来越窄(因此您拒绝了更多候选人)。

老实说，我没有看到任何其他简单的解决方案。如果 n 变大并且 cmin ~ cmax，那么你将不得不使用大锤(例如 MCMC)从你的分布中生成样本，除非我们做了一些技巧没看到。

如果你真的想使用 MCMC 算法，请注意你可以将 cmin 更改为 cmin * B/A(对于 cmax 也是如此)并假设 A == B == 1。然后问题是统一绘制两个单位 n-单纯形 (u_1...u_n, v_1...v_n) 的乘积，使得

u_i / v_i \in [cmin, cmax].

所以你必须对两个单位 n-单纯形与密度的乘积使用 MCMC 算法(Metropolis-Hastings 似乎更适合)

f(u_1, ..., u_n, v_1, ..., v_n) = \prod indicator_{u_i/v_i \in [cmin, cmax]}

这绝对可行(尽管涉及)。

关于python - 在非常具体的约束下生成随机数，我们在Stack Overflow上找到一个类似的问题： https://stackoverflow.com/questions/7908800/

python在codesectionlialgorithmrandom

有关python - 在非常具体的约束下生成随机数的更多相关文章

  1. ruby - 使用 RubyZip 生成 ZIP 文件时设置压缩级别 - 2

    我有一个Ruby程序，它使用rubyzip压缩XML文件的目录树。gem。我的问题是文件开始变得很重，我想提高压缩级别，因为压缩时间不是问题。我在rubyzipdocumentation中找不到一种为创建的ZIP文件指定压缩级别的方法。有人知道如何更改此设置吗？是否有另一个允许指定压缩级别的Ruby库？ 最佳答案 这是我通过查看ruby​​zip内部创建的代码。level=Zlib::BEST_COMPRESSIONZip::ZipOutputStream.open(zip_file)do|zip|Dir.glob("**/*")d

  2. python - 如何使用 Ruby 或 Python 创建一系列高音调和低音调的蜂鸣声？ - 2

    关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点？我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等)，但我确实想创建一个输出文件。

  3. ruby - 在 jRuby 中使用 'fork' 生成进程的替代方案？ - 2

    在MRIRuby中我可以这样做:deftransferinternal_server=self.init_serverpid=forkdointernal_server.runend#Maketheserverprocessrunindependently.Process.detach(pid)internal_client=self.init_client#Dootherstuffwithconnectingtointernal_server...internal_client.post('somedata')ensure#KillserverProcess.kill('KILL',

  4. ruby - 如何使用 Ruby aws/s3 Gem 生成安全 URL 以从 s3 下载文件 - 2

    我正在编写一个小脚本来定位aws存储桶中的特定文件，并创建一个临时验证的url以发送给同事。(理想情况下，这将创建类似于在控制台上右键单击存储桶中的文件并复制链接地址的结果)。我研究过回形针，它似乎不符合这个标准，但我可能只是不知道它的全部功能。我尝试了以下方法:defauthenticated_url(file_name,bucket)AWS::S3::S3Object.url_for(file_name,bucket,:secure=>true,:expires=>20*60)end产生这种类型的结果:...-1.amazonaws.com/file_path/file.zip.A

  5. ruby-on-rails - Ruby on Rails - 为文本区域和图片生成列 - 2

    我是Rails的新手，所以请原谅简单的问题。我正在为一家公司创建一个网站。那家公司想在网站上展示它的客户。我想让客户自己管理这个。我正在为“客户”生成一个表格，我想要的三列是:公司名称、公司描述和Logo。对于名称，我使用的是name:string但不确定如何在脚本/生成脚手架终端命令中最好地创建描述列(因为我打算将其设置为文本区域)和图片。我怀疑描述(我想成为一个文本区域)应该仍然是描述:字符串，然后以实际形式进行调整。不确定如何处理图片字段。那么……说来话长:我在脚手架命令中输入什么来生成描述和图片列？ 最佳答案 对于“文本”数

  6. ruby-on-rails - 如何生成传递一些自定义参数的 `link_to` URL？ - 2

    我正在使用RubyonRails3.0.9，我想生成一个传递一些自定义参数的link_toURL。也就是说，有一个articles_path(www.my_web_site_name.com/articles)我想生成如下内容:link_to'Samplelinktitle',...#HereIshouldimplementthecode#=>'http://www.my_web_site_name.com/articles?param1=value1¶m2=value2&...我如何编写link_to语句“alàRubyonRailsWay”以实现该目的？如果我想通过传递一些

  7. Python 相当于 Perl/Ruby ||= - 2

    这个问题在这里已经有了答案:关闭10年前。PossibleDuplicate:Pythonconditionalassignmentoperator对于这样一个简单的问题表示歉意，但是谷歌搜索||=并不是很有帮助；)Python中是否有与Ruby和Perl中的||=语句等效的语句？例如:foo="hey"foo||="what"#assignfooifit'sundefined#fooisstill"hey"bar||="yeah"#baris"yeah"另外，类似这样的东西的通用术语是什么？条件分配是我的第一个猜测，但Wikipediapage跟我想的不太一样。

  8. ruby-on-rails - 如何在 Rails 3 中创建自定义脚手架生成器？ - 2

    有这些railscast。http://railscasts.com/episodes/218-making-generators-in-rails-3有了这个，你就会知道如何创建样式表和脚手架生成器。http://railscasts.com/episodes/216-generators-in-rails-3通过这个，您可以了解如何添加一些文件来修改脚手架View。我想把两者结合起来。我想创建一个生成器，它也可以创建脚手架View。有点像RyanBates漂亮的生成器或web_app_themegem(https://github.com/pilu/web-app-theme)。我

  9. java - 什么相当于 ruby​​ 的 rack 或 python 的 Java wsgi？ - 2

    什么是ruby​​的rack或python的Java的wsgi？还有一个路由库。 最佳答案 来自Python标准PEP333:Bycontrast,althoughJavahasjustasmanywebapplicationframeworksavailable,Java's"servlet"APImakesitpossibleforapplicationswrittenwithanyJavawebapplicationframeworktoruninanywebserverthatsupportstheservletAPI.ht

  10. 报告回顾丨模型进化狂飙，DetectGPT能否识别最新模型生成结果？ - 2

    导读语言模型给我们的生产生活带来了极大便利，但同时不少人也利用他们从事作弊工作。如何规避这些难辨真伪的文字所产生的负面影响也成为一大难题。在3月9日智源Live第33期活动「DetectGPT：判断文本是否为机器生成的工具」中，主讲人Eric为我们讲解了DetectGPT工作背后的思路——一种基于概率曲率检测的用于检测模型生成文本的工具，它可以帮助我们更好地分辨文章的来源和可信度，对保护信息真实、防止欺诈等方面具有重要意义。本次报告主要围绕其功能，实现和效果等展开。（文末点击“阅读原文”，查看活动回放。）Ericmitchell斯坦福大学计算机系四年级博士生，由ChelseaFinn和Chri

随机推荐