对分法求解目标函数最优解函数MATLAB实现.docxVIP

%%对分法求解目标函数最优解函数function [t,y]=equation_extremum(x,~,a,b,epsional)%%=======================================%函数默认参数：%x为目标函数，输入字符串型，默认最优解求其最小值%a，b为其初始区间，默认闭区间，ab%对初始区间的要求：初始区间是单谷区间，函数求极小值%则a处的导数=0,b处的导数=0;%%=======================================syms t;x_inline=inline(sym(x));func_diff=diff(sym(x));func_diff_inline=inline(func_diff);t=b;high_value=func_diff_inline(t);t=a;low_value=func_diff_inline(t);ifhigh_value*low_value0 t=本函数不可解; y=给定初始区间非单谷区间;else if high_value*low_value==0ifhigh_value==0 t=b; y=x_inline(t);else t=a; y=x_inline(t);endelsewhile abs(high_value-low_value)=epsional; c=1/2*(b+a); t=c;middle_value=func_diff_inline(t);ifmiddle_value0 high_value=middle_value; b=c;else if middle_value0 low_value=middle_value; a=c;else if middle_value==0break;endendendend c=1/2*(b+a); t=c; y=x_inline(t);endend运行代码：%%========================closeall;clearall;clc;%%目标函数syms t;[t,y]=equation_extremum(t*(t-3),t,-3,5,0.1);disp([t,y]);运行结果：t=1.5y=-2.25%%抛物线法求函数最优解function [t,y]=parabola_extremum(x,~,t1,t2,epsional)%%=======================================%函数默认参数：%x为目标函数，输入字符串型，例如：x^2-sin(x)默认最优解求其最小值%a，b为其初始区间，默认闭区间，t1t2%对初始区间的要求：初始区间是单谷区间，函数求极小值%%=======================================syms t;x_inline=inline(sym(x));%%生成随机数求出t0；high_value=x_inline(t2);low_value=x_inline(t1);t0=parabola_overture(x,t,t1,t2);%自定义函数，试探法求t0middle_value=x_inline(t0);%%拟合曲线，待定系数法%p(t0)=a0+a1*t0+a2*t0^2=f(t0)%p(t1)=a0+a1*t1+a2*t2^2=f(t1)%p(t2)=a0+a1*t2+a2*t2^2=f(t2)%推出：%a1=((t0^2-t2^2)*f(t1)+(t2^2-t1^2)*f(t0)+(t1^2-t0^2)*f(t2))/((t1-t0)*(t0-t2)*(t2-t1))%a2=((t0-t2)*f(t1)+(t2-t1)*f(t0)+(t1-t0)*f(t2))/((t1-t0)*(t0-t2)*(t2-t1))while 1t_ba=0.5*((t0^2-t2^2)*low_value+(t2^2-t1^2)*middle_value+(t1^2-t0^2)*high_value)/((t0-t2)*low_value+(t2-t1)*middle_value+(t1-t0)*high_value);t_ba_value=x_inline(t_ba);if abs(t_ba-t0)epsional t=t_ba; y=t_ba_value;break;else if t_bat0ift_ba_value=middle_va