Question

我实现了一个定点方法来进行根查找。以下是代码：

clear;
clc;

f = inline('y^4 - 2', 'y');
phii = inline('y - (y^4 -2)/(4*(y^3))', 'y');
x0 = 4;
toler = 10^-10;
nmax = 1000;

a = fixpoint(x0,toler,nmax,f,phii)

function[alpha,res,iters] = fixpoint(x_0,tol,max_iters,fun,phi)

    error = tol +1;
    iters = 0;
    xvect = x_0;
    x = x_0;
    res = fun(x);
    xdif = [];

    while iters < max_iters & error > tol
        iters = iters +1;
        x_n = phi(x);
        error = abs(x_n - x);
        xdif = [xdif; error];
        x = x_n;
        xvect = [xvect;x];
        res = [res;fun(x)];
    end
    alpha = xvect(end);
    residual = res(end);
    return
end

该方法可以正常工作。我唯一的问题是我希望它返回以下所有三个：alpha，res（最后一个条目是重要性的残差，但需要所有值的向量）和迭代。

目前它只返回alpha。

最好的方法是什么？

Answer 1

[alpha,res,iters] = fixpoint(x0,toler,nmax,f,phii)

此外，此处不需要return。

Matlab - 函数：返回变量向量

1 个答案: