蒙特卡罗最优化.pptxVIP

主要内容;1.Numerical optimization methods in R;1.1.1 Bisection method(二分法）; ; 可以看出，二分法不会失效，达到指定精度所需要的迭代次数也是事先可以得到的。如果在区间[a,b]里方程有多个根，则二分; 时停止迭代。此准则可以不考虑x的单位情况下达到指定的精度。 ; 下面我们使用二分法求此方程的一个数值解。我们首先要找到一个区间，比如(0,5n)，使得函数 在区间两端有着不同的符号。然后即可使用二分法。 ;程序： a - 0.5 n - 20 cat(true roots,-a/(n-1)-sqrt(n-2-a^2+(a/(n-1))^2), + -a/(n-1)+sqrt(n-2-a^2+(a/(n-1))^2),\n) bisec-function(b0,b1){ f - function(y, a, n) { a^2 + y^2 + 2*a*y/(n-1) - (n-2) } it - 0 eps - .Machine$double.eps^0.25 r - seq(b0, b1, length=3) y - c(f(r[1], a, n), f(r[2], a, n), f(r[3], a, n)) if (y[