Dotaz: Omezení rozsahu proměnné

[Nemem9aa]

Zdravím,

narazil jsem na asi triviální záležitost, kterou ovšem nejsem schopen vyřešit čistě.

Potřeboval bych omezit vstupní proměnnou do zvoleného intervalu tak, že uprostřed daného intervalu bude proměnná beze změny ale k hranicím intervalu se proměnná bude limitně blížit.

Kupříkladu si zvolím interval 20 až 50, potom potřebuji aby došlo *přibližně* k této transformaci:

-100,0,30,170 --> 20.05, 20.1, 30, 49.95

Děkuji aspoň za nasměřování.

---

---

Odpovědi

[Aleš Kapica]

Kde? V čem?… Zkrátka málo informací.

[Nemem9aa]

Aplikovat to budu v bashi. Neuvedl jsem to, protože jsem měl za to, že implemetace bude stejná ve všech jazycích. Nebo je na toto potřeba nějaká externí knihovna?

[John]

nejdriv si ujasni/ujasnete intervaly:

1) kde budou hodnoty nezmenene(melo by byt jasne)?

2) Kde se hodnoty budou asymptoticky menit? Jakou matematickou funkci presne chcete?

Pak muzete zkoumat, jestli ji Bash umi, nebo ne.

Kazdopadne, neco jineho bude "syrova"/raw hodnota, neco jineho bude zpracovana hodnota...

[Nemem9aa]

ad 1) Na tom příliš nezáleží, může to být omezené a nebo ke změnám může docházet v celém rozsahu, pokud posun nebude příliš veliký (kupříkladu změna maximálně 10% hodnoty v rozsahu 95% intervalu).

ad 2) Toto bude podstata mého dotazu, protože já nevím. Jaká matematická funkce by se na toto dala použít/byla vhodná?

[Andrej]

Aplikovat to budu v bashi.

To asi ne, když tam jsou floating-point čísla.

[Gréta]

ty se jakoby snažíš v bashi namapovat proměnou z jednoho rosahu do druhýho?? :O :O

něco takovýho jak maj jakoby tady v pythonu hele akorátže v bashi nebo to máš nějak víc složitější že to jako nejde takle jednoduše udělat?? :O :O

víc nám to popiš cožeto jakoby chceš :O :O

[Nemem9aa]

Když jsem se snažil googlit řešení, tak jsem na obdobné příklady narazil, ale nelze je na můj případ použít. Lineární interpolace není k řešení vhodná, protože dochází uprostřed intervalu k přílišným chybám:

>>> from scipy.interpolate import interp1d
>>> m = interp1d([-1000,1000],[20,50])
>>> m(30)
array(35.45)

Moje představa je, že cca 95% intervalu by bylo s minimální změnou a k deformaci by docházelo v krajních 5% hodnot.

Omlouvám se za způsob vyjadřování, evidentně masivní absence základních matematických znalostí :(

[OldFrog {Ondra Nemecek}]

Je to primárně matematický problém. Mapování hodnoty z rozsahu na hodnotu v jiném rozsahu se v matematice nazývá funkce. Nakreslete si graf x vs. y, kde x budou vstupní hodnoty a y výstupní hodnoty. Podle toho pak navrhneme odpovídající mat fci. Tu pak použijete v bashi pomocí výpočetních příkazů (nemusí jít o čistý bash - můžete použít příkaz bc apod.)

Jinak nehledáte náhodou narovnávák na ohejbák? Chtělo by to znát detaily, možná se dá najít lepší řešení s ohledem na věcný problém, který řešíte.

[Nemem9aa]

Příloha:

• graph.png (64391 bytů)

Ano, přesně tak jsem to myslel, akorát jsem doufal, že nebudu muset malovat graf a zkoušel jsem to popsat.

Jinak konkrétní problém, kterým tímto řeším není podstatný, můžu ho vyřešit jinak, ale jsem přesvědčen že to lze řešit jednoduše nějakou triviální funkcí. Ale nejsem schopen přijít jakou a především to nejsem schopen ani vygooglit.

Náčet grafu v příloze.

[OldFrog {Ondra Nemecek}]

sigmoida hele :O ;D

[OldFrog {Ondra Nemecek}]

Ta sigmoida podle mě sedí. Detaily konkrétního problému mohou být důležité pro zhodnocení, zda je křivka pro daný problém ospravedlnitelná (pokud máte jiný dost dobrý důvod, můžete to ignorovat).

[Gréta]

vošklivý řešení je přidělat zuby takle třeba

m = interp1d([-1000,-100,100,1000],[20,23,47,50])

jinak si to asi jako musíš nějak vážit tu proměnou funckí nějakou matematickou :O :O

x = np.arange(-1000,1000,1)
y = m(x)

def todle_neni_sigmoida(x):
    return ((15+1)**(abs(x)/1000) - 1) * np.sign(x) + 35

z = todle_neni_sigmoida(x)


for i in range(-1000,1100,100):
    print(f'x: {i} y: {todle_neni_sigmoida(i)}')

x: -1000 y: 20.0
x: -900 y: 23.874267467916816
x: -800 y: 26.81041316002372
x: -700 y: 29.035595493631007
x: -600 y: 30.721968356908423
x: -500 y: 32.0
x: -400 y: 32.9685668669792
x: -300 y: 33.70260329000593
x: -200 y: 34.25889887340775
x: -100 y: 34.68049208922711
x: 0 y: 35.0
x: 100 y: 35.31950791077289
x: 200 y: 35.74110112659225
x: 300 y: 36.29739670999407
x: 400 y: 37.0314331330208
x: 500 y: 38.0
x: 600 y: 39.27803164309158
x: 700 y: 40.96440450636899
x: 800 y: 43.18958683997628
x: 900 y: 46.125732532083184
x: 1000 y: 50.0

[Jendа]

tanh

[Gréta]

taky sigmoida hele :O ;D

[Kit]

https://cs.wikipedia.org/wiki/Logistická funkce

[Gréta]

taky sigmoida hele :O ;D

[wamba]

Rozdělit si to na tři části je čisté řešení?

raku -e '
sub in-range ($_ ) { 
    when * < 21 {20 + 1/(22-$_) }
    when 21 ≤ * ≤ 49 { $_ } 
    when 49 < * { 50 - 1/($_-48) }
};
say $_ => in-range $_ for -100,0,30,35,170
'

-100 => 20.008197
0 => 20.045455
30 => 30
35 => 35
170 => 49.991803

[Gréta]

v raku boužel jako neexistuje žádný čistý řešení :O :/ :/ ;D