Module:Harmonic entropy - Revision history

ArrowHead294: Remove already finished Todo

2024-06-13T13:53:51Z

Remove already finished Todo

← Older revision		Revision as of 13:53, 13 June 2024
Line 1:		Line 1:
	~~-- TODO: this module should just be called harmonic entropy!~~
	local limits = require("Module:Limits")		local limits = require("Module:Limits")
	local rat = require("Module:Rational")		local rat = require("Module:Rational")
	local p = {}		local p = {}

	-- ~~compute harmonic~~ Shannon entropy for an interval of `c` cents		-- Compute Harmonic Shannon entropy for an interval of `c` cents
	-- `c`, `deviation`: in cents		-- `c`, `deviation`: in cents
	-- `ratios`: an array of rational numbers		-- `ratios`: an array of rational numbers

2024-06-12T11:10:18Z

← Older revision		Revision as of 11:10, 12 June 2024
Line 12:		Line 12:
	return math.sqrt(rat.benedetti_height(ratio))		return math.sqrt(rat.benedetti_height(ratio))
	end		end
	deviation = deviation or 17.~~22635157248405~~		deviation = deviation or 1200 * math.log(1.01, 2)
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)

2024-06-12T07:25:08Z

FloraC moved page Module:Interval to Module:Harmonic entropy without leaving a redirect: See talk

← Older revision	Revision as of 07:25, 12 June 2024
(No difference)

2024-05-10T13:11:12Z

i give up

← Older revision		Revision as of 13:11, 10 May 2024
Line 12:		Line 12:
	return math.sqrt(rat.benedetti_height(ratio))		return math.sqrt(rat.benedetti_height(ratio))
	end		end
	deviation = deviation or 0.01		deviation = deviation or 17.22635157248405
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(rat.~~as_float~~(ratio) * 2^(rat.as_float(c~~) / 1200)~~) > 3 * deviation then		if math.abs(rat.cents(ratio) - c) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(rat.~~as_float~~(ratio) * 2^(rat.as_float(c~~) / 1200)~~) / norm(ratio)		return gaussian(rat.cents(ratio) - c) / norm(ratio)
	end		end

2024-05-10T13:09:48Z

← Older revision		Revision as of 13:09, 10 May 2024
Line 15:		Line 15:
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(rat.as_float(ratio) * 2^(c / 1200)) > 3 * deviation then		if math.abs(rat.as_float(ratio) * 2^(rat.as_float(c) / 1200)) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(rat.as_float(ratio) * 2^(c / 1200)) / norm(ratio)		return gaussian(rat.as_float(ratio) * 2^(rat.as_float(c) / 1200)) / norm(ratio)
	end		end

2024-05-10T13:08:38Z

← Older revision		Revision as of 13:08, 10 May 2024
Line 15:		Line 15:
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(rat.~~abs~~(ratio) * 2^(c / 1200)) > 3 * deviation then		if math.abs(rat.as_float(ratio) * 2^(c / 1200)) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(rat.~~abs~~(ratio) * 2^(c / 1200)) / norm(ratio)		return gaussian(rat.as_float(ratio) * 2^(c / 1200)) / norm(ratio)
	end		end

2024-05-10T13:06:54Z

← Older revision		Revision as of 13:06, 10 May 2024
Line 15:		Line 15:
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(ratio~~[1]~~ * 2^(c / 1200)) > 3 * deviation then		if math.abs(rat.abs(ratio) * 2^(c / 1200)) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(ratio~~[1]~~ * 2^(c / 1200)) / norm(ratio)		return gaussian(rat.abs(ratio) * 2^(c / 1200)) / norm(ratio)
	end		end

2024-05-10T13:04:49Z

← Older revision		Revision as of 13:04, 10 May 2024
Line 15:		Line 15:
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(ratio[0] * 2^(c / 1200)) > 3 * deviation then		if math.abs(ratio[1] * 2^(c / 1200)) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(ratio[0] * 2^(c / 1200)) / norm(ratio)		return gaussian(ratio[1] * 2^(c / 1200)) / norm(ratio)
	end		end

2024-05-10T13:04:24Z

← Older revision		Revision as of 13:04, 10 May 2024
Line 15:		Line 15:
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(ratio * 2^(c / 1200)) > 3 * deviation then		if math.abs(ratio[0] * 2^(c / 1200)) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(ratio * 2^(c / 1200)) / norm(ratio)		return gaussian(ratio[0] * 2^(c / 1200)) / norm(ratio)
	end		end

2024-05-10T12:59:43Z

Try this out

← Older revision		Revision as of 12:59, 10 May 2024
Line 12:		Line 12:
	return math.sqrt(rat.benedetti_height(ratio))		return math.sqrt(rat.benedetti_height(ratio))
	end		end
	deviation = deviation or 17.~~22635157248405~~		deviation = deviation or 0.01
	ratios = ratios		ratios = ratios
	or limits.integer_limit(200, function(ratio)		or limits.integer_limit(200, function(ratio)
	if math.abs(~~rat.cents~~(~~ratio~~) ~~- c~~) > 3 * deviation then		if math.abs(ratio * 2^(c / 1200)) > 3 * deviation then
	return 1 / 0		return 1 / 0
	end		end
Line 26:		Line 26:

	local function weighted_gaussian(ratio)		local function weighted_gaussian(ratio)
	return gaussian(~~rat.cents~~(~~ratio~~) ~~- c~~) / norm(ratio)		return gaussian(ratio * 2^(c / 1200)) / norm(ratio)
	end		end