Question

您好我正在尝试使用scipy.optimize.fmin来最小化一个函数。但事情进展不顺利，因为我的计算似乎有分歧而不是收敛，我得到了一个错误。我试图修复一个公差，但它不起作用。这是我的代码（主程序）：

import sys,os
import numpy as np
from math import exp
import scipy
from scipy.optimize import fmin

from carlo import *

A=real()

x_r=0.11245
x_i=0.14587

#C=A.minim
part_real=0.532
part_imag=1.2
R_0 = fmin(A.minim,[part_real,part_imag],xtol=0.0001)

上课：

import sys,os
import numpy as np
import random, math
import matplotlib.pyplot as plt
import cmath
#import pdb
#pdb.set_trace()

class real:
    def __init__(self):
        self.nmodes = 4
        self.L_ch = 1
        self.w = 2

    def minim(self,p):

        x_r=p[0]
        x_i=p[1]
        x=complex(x_r,x_i)
        self.a=complex(3,4)*(3*np.exp(1j*self.L_ch))
        self.T=np.array([[0.0,2.0*self.a],[(0.00645+(x)**2), 4.3*x**2]])
        self.Id=np.array([[1,0],[0,1]])
        self.disp=np.linalg.det(self.T-self.Id)
        print self.disp
        return self.disp

错误是：

(-2.16124712985-8.13819476595j)
/usr/local/lib/python2.7/site-packages/scipy/optimize/optimize.py:438: ComplexWarning: Casting complex values to real discards the imaginary part
  fsim[0] = func(x0)
(-1.85751684826-8.95377303768j)
/usr/local/lib/python2.7/site-packages/scipy/optimize/optimize.py:450: ComplexWarning: Casting complex values to real discards the imaginary part
  fsim[k + 1] = f
(-2.79592712985-8.13819476595j)
(-3.08484130014-7.36240080015j)
(-3.68788935914-6.62639114029j)
/usr/local/lib/python2.7/site-packages/scipy/optimize/optimize.py:475: ComplexWarning: Casting complex values to real discards the imaginary part
  fsim[-1] = fxe
(-2.62046851255e+87-1.45013007728e+88j)
(-4.037931857e+87-2.2345341712e+88j)
(-7.45017628087e+87-4.12282179854e+88j)
(-1.14801242605e+88-6.35293780534e+88j)
(-2.11813751435e+88-1.17214723347e+89j)
Warning: Maximum number of function evaluations has been exceeded.

实际上我并没有理解为什么计算会发生分歧，也许我不得不使用其他东西而不是使用fmin进行最小化？有人有个主意吗？非常感谢你。

Answer 1

尝试优化绝对值而不是复数值。这为我带来了不错的结果。

f = lambda x: abs(A.minim(x))
R_0 = fmin(f,[part_real,part_imag],xtol=0.0001)

我猜fmin不能很好地处理复杂的值。

如何使用scipy.optimize.fmin

1 个答案: