python - Scipiy.optimize.minimaze unconstrained error -


i want perform constrained least square methon polynom. before that, decided try unconstrined optimalization. here problem:

my polynom looks

$f(x) = ax^4 + bx^3 + cx^2 + dx + e$

so want find best coefficients a, b, c, d , e. function want minimize looks like

def lsq(args, x, y):     return sum([(y[i] - (args[0]*x[i]**4 + args[1]*x[i]**3 + args[2]*x**2 + args[3]*x + args[4]))**2                 in np.arange(len(x))])

where x , y lists of coordinates of points. thus, code may like:

import numpy np import scipy.optimize ph = np.array([8,8,8,7,7,7,7,7,7,7,7,7,6,3,2,2,2,1]) def rank2(y):     return np.array([(i+1)/len(y) in range(len(y))]) x = rank2(ph) y = ph def lsq(args, x, y):     return sum([(y[i] - (args[0]*x[i]**4 + args[1]*x[i]**3 + args[2]*x**2 + args[3]*x + args[4]))**2                 in np.arange(len(x))])   params = scipy.optimize.minimize(nejmensi_ctverce, [1.0, 1.0, 1.0, 1.0, 1.0], args = (x, y))

but here got error:

file "c:\users\robert\desktop\winpython-64bit-3.6.1.0qt5\python- 3.6.1.amd64\lib\site-packages\scipy\optimize\optimize.py", line 628, in _approx_fprime_helper grad[k] = (f(*((xk + d,) + args)) - f0) / d[k]  valueerror: setting array element sequence

can me? guess not cempletely understand how parse arguments minimize function.

some xs in lsq not indexed i:

args[0]*x[i]**4 + args[1]*x[i]**3 + args[2]*x**2 + args[3]*x + args[4]                                          ---^           ---^

this led lsq returning array of values instead of scalar:

in [9]: lsq([1.0, 1.0, 1.0, 1.0, 1.0], x, y) out[9]:  array([ 468.00714962,  458.38490951,  448.01979096,  436.95911906,         425.25433416,  412.9609918 ,  400.13876277,  386.85143307,         373.16690395,  359.15719185,  344.89842847,  330.47086072,         315.95885075,  301.45087591,  287.0395288 ,  272.82151724,         258.89766428,  245.37290818])

this leads valueerror since scipy.optimize.minimze expects lsq return scalar values minimized. 1 way fix problem replace bare xs x[i]s.

a better way solve problem replace x[i]s xs, remove for in np.arange(len(x)) , use numpy array-based arithmetic:

def lsq(args, x, y):     return ((y - (args[0]*x**4 + args[1]*x**3 + args[2]*x**2                    + args[3]*x + args[4]))**2).sum()

for example,

import numpy np import scipy.optimize optimize  ph = np.array([8,8,8,7,7,7,7,7,7,7,7,7,6,3,2,2,2,1]) def rank2(y):     return np.array([(i+1)/len(y) in range(len(y))]) x = rank2(ph) y = ph  def lsq(args, x, y):     a, b, c, d, e = args     return ((y - (a*x**4 + b*x**3 + c*x**2 + d*x + e))**2).sum()  params = optimize.minimize(lsq, [1.0, 1.0, 1.0, 1.0, 1.0], args = (x, y)) print(params.x)

yields

[  94.48618936 -211.42358992  144.93063545  -37.24078798   10.23934514]

with minimum lsq value of:

print(lsq(params.x, x, y)) # 6.91284752049

Comments

Post a Comment

Popular posts from this blog

Is there a better way to structure post methods in Class Based Views -

i'm using class based views in django , never specify url in ajax calls , omit action parameter in form because know these requests go through post method of class based view. notice view - profile page make many requests depending on user updating on page - post method become convoluted , if/else mess example: def post(self, request, *args, **kwargs): # handle user changing available date if request.post.get('availabledate') != none: self.updateavailabledate(request) return jsonresponse({'result':'success'}) if request.post.get('newprojectname') != none: creative_user = creativeuserprofile.objects.get(id = request.user.id) project = project.create_project(creative_user, request.post['newprojectname']) project.save() return jsonresponse({'projectid': project.id }) if request.post.get('projectimage') != none , request.files['file']: proj...
Read more

performance - Why is XCHG reg, reg a 3 micro-op instruction on modern Intel architectures? -

i'm doing micro-optimization on performance critical part of code , came across sequence of instructions (in at&t syntax): add %rax, %rbx mov %rdx, %rax mov %rbx, %rdx i thought had use case xchg allow me shave instruction , write: add %rbx, %rax xchg %rax, %rdx however, dimay found agner fog's instruction tables , xchg 3 micro-op instruction 2 cycle latency on sandy bridge, ivy bridge, broadwell, haswell , skylake. 3 whole micro-ops , 2 cycles of latency! 3 micro-ops throws off 4-1-1-1 cadence , 2 cycle latency makes worse original in best case since last 2 instructions in original might execute in parallel. now... cpu might breaking instruction micro-ops equivalent to: mov %rax, %tmp mov %rdx, %rax mov %tmp, %rdx where tmp anonymous internal register , suppose last 2 micro-ops run in parallel latency 2 cycles. given register renaming occurs on these micro-architectures, though, doesn't make sense me done way. why wouldn't register renamer ...
Read more

c# - Asp.net web api : redirect unauthorized requst to forbidden page -

im trying redirect unauthorized request forbidden page instead i'm getting forbidden page in response body , how can fix ? here's startup class : app.createperowincontext(storecontext.create); app.createperowincontext<applicationusermanager>(applicationusermanager.create); app.createperowincontext<applicationsigninmanager>(applicationsigninmanager.create); app.usecookieauthentication(new cookieauthenticationoptions { authenticationtype = defaultauthenticationtypes.applicationcookie, expiretimespan = timespan.fromdays(30), }); app.useexternalsignincookie(defaultauthenticationtypes.externalcookie); app.usetwofactorsignincookie(defaultauthenticationtypes.twofactorcookie, timespan.fromminutes(5)); app.usetwofactorrememberbrowsercookie(defaultauthenticationtypes.twofactorrememberbrowsercookie); the method im trying reach : [httpget] [authorize(roles = "admin")] public string getcurrentusername() { return userman...
Read more