Das deutsche Python-Forum

Code: Alles auswählen

import numpy as np

def lma(fs, dfs, initial, lmbd, min_gain):
    def error(x):
        return 0.5*sum(f(x)**2 for f in fs)
    
    def evalualte(x):
        return np.matrix([[f(x)] for f in fs])
    
    def jacobian(x):
        return np.matrix([[f(x) for f in r] for r in dfs])
    
    def step(x):
        j = jacobian(x)
        jtj = j.T*j
        li = lmbd * np.identity(jtj.shape[0])
        left = jtj+li
        
        k = evalualte(x)
        right = -j.T*k
        
        return left.I*right
    
    
    x = np.matrix(initial).T
    gain = 2 * min_gain
    new_error = None
    
    while gain >= min_gain:
        delta = step(x)
        old_error = new_error
        new_error = error(x)
        x += delta
        
        if old_error is not None:
            gain = old_error - new_error
    
    return x
        

def main():
    f1 = lambda x: x[0,0]**2
    f2 = lambda x: x[1,0]**2
    f3 = lambda x: x[0,0]*x[1,0]
    
    f1dx = lambda x: 2*x[0,0]
    f1dy = lambda x: 0.0
    f2dx = lambda x: 0.0
    f2dy = lambda x: 2*x[1,0]
    f3dx = lambda x: x[1,0]
    f3dy = lambda x: x[0,0]
    
    fs = [f1, f2, f3]
    dfs = [[f1dx, f1dy], [f2dx, f2dy], [f3dx, f3dy]]
    
    print lma(fs, dfs, (10.0, 20.0), 1.0, 0.00001)


if __name__ == "__main__":
    main()

Code: Alles auswählen

 y_data = columns[1]
IndexError: list index out of range

Code: Alles auswählen

y_data = []
t = []
crs = open("Datei.txt", "r")
for columns in ( raw.strip().split() for raw in crs ):
    y_data = columns[1]
    t = columns[0]
    
data.append(y_data)
time.append(time)

Code: Alles auswählen

for columns in ( raw.strip().split() for raw in crs ):
    ...

Code: Alles auswählen

with open(...) as handler:
    # handler ist hier File-Object
    # es wird nach Durchlaufen des Blocks automatisch geschlossen

Code: Alles auswählen

0.1837
-1.6762
-1.7639
0.3118
4.0473
8.4418
12.4349
15.0443
15.4261
13.2425
9.1291
4.5988
1.2063
-0.3603
-0.4915
-0.0904
0.1993
0.3526
0.7231
1.4638
2.0657
1.4138
-1.5479
-6.9588
-13.6981
-19.8780
-23.7647
-24.4734
-22.1911
-18.0425
-13.6499
-10.4237
-8.9728
-9.1255
-10.4815
-12.8756
-16.3304
-20.6567
-25.1772
-28.9096
-31.0687
-31.4088
-30.2043
Traceback (most recent call last):
  File "/home/CurveFiting/src/Read.py", line 29, in <module>
    data_load = columns[3]
IndexError: list index out of range

Code: Alles auswählen

>>> data = """1 eins
... 2 zwei
... 3
... 4 vier"""
>>>
>>> for line in data.split("\n"):
...     print(line)
...
1 eins
2 zwei
3
4 vier
>>> for line in data.split("\n"):
...     print(line.split())
...
['1', 'eins']
['2', 'zwei']
['3']
['4', 'vier']
>>> for line in data.split("\n"):
...     print(line.split()[1])
...
eins
zwei
Traceback (most recent call last):
  File "<stdin>", line 2, in <module>
IndexError: list index out of range