![](https://img-blog.csdnimg.cn/20201014180756738.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
最优化计算
与月遥远
这个作者很懒,什么都没留下…
展开
-
最优化作业||拟牛顿法
拟牛顿法代码function res_x=fxninewton(f,x,n,x0,e)%拟牛顿法(DPF) syms t;%声明变量t grad=gradient(f,x);%一维求导矩阵 xk=x0;%当前点 fk=subs(f,x,x0);%当前值 gk=subs(grad,x,x0);%一维矩阵 ak=eye(n); k=0;%迭代数 while(norm(gk)>e)%开始迭代 dk=-ak*gk;原创 2021-11-02 16:16:10 · 249 阅读 · 0 评论 -
最优化作业||阻尼牛顿法
阻尼牛顿法代码function res_x=fxnewton(f,x,n,x0,e)%阻尼牛顿法 syms t; h=hessian(f); h_1=inv(h);%海塞矩阵的逆 h_1=double(h_1); grad=gradient(f,x);%一维求导矩阵 g1=subs(grad,x,x0); d=-g1 index=0; while(norm(g1)>e)%开始迭代 d=double(d);原创 2021-11-02 16:15:02 · 1072 阅读 · 0 评论 -
最优化作业||最速下降法
最速下降法代码function [res_x,index]=fxsteep(f,x,n,x0,e)%最速下降法 Q=hessian(f); q=double(Q); grad=gradient(f,x);%一维求导矩阵 g1=subs(grad,x,x0); d=-g1; index=0; while(norm(g1)>e)%开始迭代 lamda = (-d)'*d/((-d)'*q*d);%最优步长d d=do原创 2021-11-02 16:12:13 · 346 阅读 · 0 评论 -
最优化作业||抛物线逼近法
抛物线逼近法代码function [x,f_x]=parabola(f,x1,x2,x3)eps = 0.0001;i=0;f=@(x) x.^ 3 - x.^2 - 30*x; %创建题目要求匿名函数,方便使用while 1i=i+1;f1=f(x1); f2=f(x2); f3=f(x3); a=2*(f1*(x2-x3)+f2*(x3-x1)+f3*(x1-x2));b=(f1*(x2^2-x3^2)+f2*(x3^2-x1^2)+f3*(x1^2-x2^2));xp=b/原创 2021-11-02 16:08:54 · 458 阅读 · 0 评论 -
最优化作业||黄金分割法
黄金分割法代码function [x,f_x]=goldfun(f,a0,b0)%黄金分割法a = a0;b = b0; %初始区间index=0;%迭代数eps = 0.01; %区间精度while((b - a) >= eps) a; b; index=index+1; x1 = a + 0.382 * (b - a); x2 = a + 0.618 * (b - a); %黄金分割法主要步骤 if f(x1) < f(x2) %两原创 2021-11-02 16:07:43 · 945 阅读 · 0 评论 -
最优化作业||外推内插法
打工中||最优化作业外推内插法代码// 外推内插法function [a,b]=JinTui(f,x0,h0)%f为函数%x0为初始点%h0为初始步长x1=x0;x2=x1+h0;i=1;if f(x2) < f(x1) %下降沿,前进运算 h=2*h0; x3=x2+h; while f(x3) < f(x2) %迭代过程 h=2*h;%步长翻倍 x1=x2; x2=x3; x3=x3原创 2021-10-25 21:21:28 · 495 阅读 · 0 评论