quadi.m

来自「数字通信第四版原书的例程」· M 代码 · 共 19 行

function [stepmin]=quadinter(steps,matf)
%QUADINTER Quadraticly interpolates three points to estimate minimum.
%         This function uses QUADRATIC interpolation and
%         the values of three unevenly spaced points 
%         in order estimate the minimum of a 
%         a function along a line.
%    syntax: [stepmin]=quadinter(steps,matf)
%    where steps is a matrix constraining the spacing of the points.
%          matf is a vector which contains the corresponding elements of f 
%          stepmin is the step to the minimum.

%	Copyright (c) 1990 by the MathWorks, Inc.
%	Andy Grace 7-9-90.
steps=steps(:);
% Solve simultaneous equations:
      amat=[0.5*steps.*steps, steps, ones(3,1)];
      abc=amat\matf(:);
      stepmin=-abc(2)/abc(1);