抛物线法程序的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

function [x,minf]=minpw(f,a,b,eps)
format long;
if nargin==3
    eps=1.0e-4;
end
t0=(a+b)/2;
while 1
    fa=subs(f,findsym(f),a);
    fb=subs(f,findsym(f),b);
    ft0=subs(f,findsym(f),t0);
    u=fa*(b^2-t0^2)+fb*(t0^2-a^2)+ft0*(a^2-b^2);
    v=fa*(b-t0)+fb*(t0-a)+ft0*(a-b);
    t1=0.5*u/v;
    ft1=subs(f,findsym(f),t1);
    tol=abs(t1-t0);
    if ft1<=ft0
        if tol<eps
           x=t1;
            break;
        else
            if t1<=t0
                b=t0;
               t0=t1;

            else a=t0;
               t0=t1;
            end

        end
    else
       if tol<eps
           x=t0;
           break;
       else
           if t1<=t0
              a=t1;

          else b=t1;

           end
      end
    end
end

minf=subs(f,findsym(f),x);
format short;

一直在运行，出不来结果。麻烦大神看一下哪里有问题。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：抛物线 function format nargin FORMA 抛物线 程序

加几个断点看看吧，这种情况大多是break条件没有满足造成的。