运行时错误：未加括号的根

Question

class ComputeIV
{
    public:    
        typedef std::pair<SimpleQuote,SimpleQuote> BidAsk;    

        static Volatility ComputeImpliedVol(const Date evalDate, const Date expiration, ptime quoteTime, const Option::Type optionType, 
                  const Real underlyingPrice, const Real optionPrice, const Real strike, const Rate riskFree) 
        {
            ActualActual actualActual;
            Settings::instance().evaluationDate() = evalDate;    

            Time timeToMaturity = actualActual.yearFraction(Settings::instance().evaluationDate(), expiration);

            time_duration timeOfDayDuration = quoteTime.time_of_day();
            timeToMaturity += (timeOfDayDuration.hours() + timeOfDayDuration.minutes()/60.0)/(24.0 * 365.0);               

            DiscountFactor discount = std::exp(-riskFree * timeToMaturity);

            Bisection bisection;
            Real accuracy = 0.000001, guess = .20;
            Real min = .05, max = .40;
            Volatility sigma = bisection.solve([&](const Volatility & sigma) {
                Real stdDev = sigma * std::sqrt(timeToMaturity);
                BlackCalculator blackCalculator(optionType, strike, underlyingPrice, stdDev, discount);
                    return blackCalculator.value() - optionPrice;
            }, accuracy, guess, min, max);  

            return sigma;    

        } 

    static const double& pTest(const std::string evalDateStr, const std::string expirationStr, const std::string quoteTimeStr, 
                               const int optType, const Real forwardBid, const Real forwardAsk, const Rate riskFree, const Real strike, 
                               const Real oBid, const Real oAsk)
    {
        std::cout << "Computing IV" << std::endl << std::flush;
        ActualActual actualActual;

        std::cout << evalDateStr << " " << expirationStr << " " << quoteTimeStr << std::endl << std::flush;
        std::cout << optType << " "  << riskFree << " " << forwardBid << " " << forwardAsk << " " << strike 
                  << " " <<  oBid << " " << oAsk << std::endl << std::flush;

        Date evalDate = DateParser::parseFormatted(evalDateStr.c_str(), "%d/%m/%Y");
        Settings::instance().evaluationDate() = evalDate;

        Date expiration = DateParser::parseFormatted(expirationStr.c_str(), "%d/%m/%Y");

        ptime quoteTime(from_iso_string(quoteTimeStr));
        time_duration timeOfDayDuration = quoteTime.time_of_day();

        Real price = (oBid + oAsk) / 2.0;    

        Option::Type oType = (optType > 0 ? Option::Call : Option::Put);
        Volatility *sigma = new Volatility();
        *sigma = ComputeIV::ComputeImpliedVol(evalDate, expiration, quoteTime, oType, 
                                forwardAsk, price, strike, riskFree);

        return *sigma;                        
    }

    private:        

};

当我用这些参数打电话时

^{pr2}$

我得到一个运行时错误：

Traceback (most recent call last):
  File "cboelivedata.py", line 575, in <module>
    main()
  File "cboelivedata.py", line 400, in main
    sigma = ComputeIV.pTest(evalDate, expiration, quoteTime, oType, forwardBid, forwardAsk, riskFree, strike, oBid, oAsk)
RuntimeError: root not bracketed: f[0.05,0.4] -> [-1.350000e-01,-2.434993e-02]
[idf@node3 python]$

我有什么地方做得不对吗？

Answer 1

你用的是黑色公式。在Black-Scholes欧式期权框架中，期权价格是波动性的递增函数。在

平分波动率的最小和最大范围是0.05和{}。它们够好吗？你的看跌期权已严重缺钱，因此需要很高的波动性。在

让我们检查一下你的射程。转到http://www.erieri.com/blackscholes，然后键入如下信息：

你的解算器能给你的最大看跌期权价格约为0.0797，但你的报价是0.135。因此，没有解，你的二等分根解算器告诉你是正确的。在

你需要增加波动范围。尝试：

Real min = .05, max = 1.00;

你可能不需要1.00（这太大了），但是你明白了——你需要调整你的根解算范围。在

试试看，你就会得到隐含波动率。在