https://ja.xen.wiki/index.php?action=history&feed=atom&title=%E3%83%A2%E3%82%B8%E3%83%A5%E3%83%BC%E3%83%AB%3AHarmonic_entropy
モジュール:Harmonic entropy - 版の履歴
2026-06-05T23:38:57Z
このウィキのこのページに関する変更履歴
MediaWiki 1.43.6
https://ja.xen.wiki/index.php?title=%E3%83%A2%E3%82%B8%E3%83%A5%E3%83%BC%E3%83%AB:Harmonic_entropy&diff=1283&oldid=prev
Dummy index: ページの作成:「local limits = require("Module:Limits") local rat = require("Module:Rational") local p = {} -- Compute Harmonic Shannon entropy for an interval of `c` cents -- `c`, `deviation`: in cents -- `ratios`: an array of rational numbers -- `norm`: a function of rational numbers function p.harmonic_entropy(c, ratios, deviation, norm) norm = norm or function(ratio) return math.sqrt(rat.benedetti_height(ratio)) end deviation = deviation or 1200 * math.log(1.01, 2) ratio…」
2025-11-16T03:55:02Z
<p>ページの作成:「local limits = require("Module:Limits") local rat = require("Module:Rational") local p = {} -- Compute Harmonic Shannon entropy for an interval of `c` cents -- `c`, `deviation`: in cents -- `ratios`: an array of rational numbers -- `norm`: a function of rational numbers function p.harmonic_entropy(c, ratios, deviation, norm) norm = norm or function(ratio) return math.sqrt(rat.benedetti_height(ratio)) end deviation = deviation or 1200 * math.log(1.01, 2) ratio…」</p>
<p><b>新規ページ</b></p><div>local limits = require("Module:Limits")<br />
local rat = require("Module:Rational")<br />
local p = {}<br />
<br />
-- Compute Harmonic Shannon entropy for an interval of `c` cents<br />
-- `c`, `deviation`: in cents<br />
-- `ratios`: an array of rational numbers<br />
-- `norm`: a function of rational numbers<br />
function p.harmonic_entropy(c, ratios, deviation, norm)<br />
norm = norm or function(ratio)<br />
return math.sqrt(rat.benedetti_height(ratio))<br />
end<br />
deviation = deviation or 1200 * math.log(1.01, 2)<br />
ratios = ratios<br />
or limits.integer_limit(200, function(ratio)<br />
if math.abs(rat.cents(ratio) - c) > 3 * deviation then<br />
return 1 / 0<br />
end<br />
return norm(ratio)<br />
end, 100)<br />
<br />
local function gaussian(x)<br />
return math.exp(-x * x / (2 * deviation * deviation)) / (deviation * math.sqrt(2 * math.pi))<br />
end<br />
<br />
local function weighted_gaussian(ratio)<br />
return gaussian(rat.cents(ratio) - c) / norm(ratio)<br />
end<br />
<br />
local q_norm = 0<br />
for _, ratio in pairs(ratios) do<br />
q_norm = q_norm + weighted_gaussian(ratio)<br />
end<br />
<br />
local function probability(ratio)<br />
return weighted_gaussian(ratio) / q_norm<br />
end<br />
<br />
local entropy = 0<br />
for _, ratio in pairs(ratios) do<br />
local p_i = probability(ratio)<br />
if p_i > 1e-5 then<br />
entropy = entropy - p_i * math.log(p_i)<br />
end<br />
end<br />
return entropy<br />
end<br />
<br />
return p</div>
Dummy index