Matematické Fórum

↑ xm111:Třeba sem?

↑ xm111:Myslíš, že ti tu všichni budou psát a prosit tě, abys jim to poslal a oni ti to mohli vyřešit? Nenapadá tě nic, jak to případným pomocníkům nekomplikovat?

Můžeš to nahrát na nějakou úschovnu a hodit sem odkaz...

↑ xm111:No vidíš, že jsi to zvládl, jsi šikulka. Teď ještě pořádně popsat o co ti vlastně v té tabulce jde a bude to perfektní.

Zdravím, ve vzorci máte psáno (pouze výběr):

EXP(1)^-E9^2*H9

ale mělo by být:

EXP(-(E9^2*H9))

a tato chyba se několikrát opakuje, dál jsem se nedíval ...

Mimochodem výsledkem výrazu -E9^2 v Excelu je kladné číslo a to je špatně, správně by mělo být -(E9^2).

Vzorec je poměrně nepřehledný, snadno se v něm udělá chyba a ještě nesnadněji se kontroluje, dalo by se to napsat jednodušeji a přehledněji jako makro ve smyčce, s využitím názvů proměnných (případně i s indexy)

↑ xm111: Tak to nebude jediná chyba. místo 22,84 my vyšlo něco kolem 55.  Zkusím to ještě nekolikrát přepočítat a zkusit to udělat třeba jinak.

Rozdělil bych výpočet do dvou částí, stejně jak to tam máš.
První část bych počítal vzorec:

[mathjax]2\,\sum_{n=1}^{10}{{{e^ {- F_{0}\,q_{n}^2 }\,\sin q_{n}\,\cos  \left({{q_{n}\,x}\over{b}}\right)}\over{\cos q_{n}\,\sin q_{n}+q_{n} }}}[/mathjax]

Ten se mi fakt kontrolovat nechce. Takže jestli výsledek sedí, tak tuhle část můžeme přeskočit.
Napíšu jenom makro jak jsem to počítal já:

Function Vzorec(x#, b#, F0#, q As Variant) As Variant
    Dim T#, qn#, Z As Integer 
    On Error GoTo Chyba
    Z=UBound(q) 
    T=0
    For i=1 To Z
        qn = q(i,1)    
        T = T + (Sin(qn)*Cos(qn*x/b)*exp(-((qn)^2*F0))) / (qn+Sin(qn)*Cos(qn))
    Next     
    Vzorec = 2*T
    Exit Function
Chyba:
    Vzorec = "Err" 
End Function

Nevím jak si na tom s makry, ale už od pohledu vidíš, že je to mnohem přehlednější.
Makro je napsáno pro LibreOffice, tak nezaručuji, že bude pracovat správně v Excelu.
Podle názvů proměnných asi poznáš jak pracuje.
Zádrhel asi bude v poslední předávané proměnné [mathjax]q[/mathjax], které se předává pole dat (tj. buňky $E$9:$E$18) - program si je načítá přímo z listu. (možná pomůže odstranit "As Variant", nebo nahradit za "As Object", ale to fakt nevím, možná poradí někdo jiný)

Druhá část je jednodušší tu počítáš správně dle vzorce (24):
[mathjax]t = t^{*} (t_{p}-t_{o}) - t_{o}[/mathjax]

Přičemž [mathjax]t_{o}[/mathjax] je počáteční teplota a [mathjax]t_{p}[/mathjax] je konečná, takže změna ohřevu za chlazení se provede pouze výměnou těchto teplot.

Ještě připomenu, že aby to vyšlo podle nich, musí mít proměnná [mathjax]b[/mathjax] hodnotu 0.006 (ty tam máš 0.004)

Odkaz na soubor