Question

所以我四处寻找导致这个错误的原因，并且人们继续列出罪魁祸首是matlab不知道函数的目录路径。我检查过，这似乎不是问题，因为matlab知道路径。我这里有一些代码，假设要进行高斯消除

classdef gaussfunctions
    %UNTITLED Summary of this class goes here
    %   Detailed explanation goes here

properties
end

methods (Static)
    function solution_vec = gauss_reduce(param_mat, const_vec)
        %check for conistent sixe of param_mat and const_vec
        [n,m] = size(param_mat);

        if n ~= m | n ~= size(const_vec)
            disp('Non square matrix or constant vector was incorrect size');
            solution_vec = zeros(size(const_vec));
            return

        end

        %reduce coeffeicient matrix to upper triangular
        [ut_mat, new_const_vec] = ut_reduce(param_mat, const_vec);

        %compute the solution vector
        solution_vec = back_subst(ut_mat, new_const_vec);
    end

    function [ut_mat, new_const_vec] = ut_reduce(param_mat, const_vec)
        %get the size of the parameter matrix
        [m,n] = size(param_mat);

        %get the values of the arguments to change and the return
        cur_mat = param_mat;
        cur_vec = const_vec;

        %loop through all of the columns to reduce them
        for cur_col = 1:m,
            [cur_mat, cur_vec] = reduce_column(cur_mat, cur_vec, cur_col);
        end

        %return the values
        ut_mat = cur_mat;
        new_const_vec = cur_vec;

    end

    function [new_mat, new_const_vec] = reduce_column(param_mat, const_vec, column)
        %get the size 
        [m,n] = size(param_mat);

        %copy the argument values to change them and the return
        cur_mat = param_mat;
        cur_vec = const_vec;

        %check to see if the pivot is zero or not
        if param_mat(column,column) == 0

            disp('Zero pivot encourtered');
            new_mat = param_mat;
            new_const_vec = zeros(size(const_vec));

            return
        end

        %loops down calling reduce_row_at_col on all row below the
        %current column
        for i = column + 1:n,
            [cur_mat, cur_vec] = reduce_row_at_col(cur_mat, cur_vec, column, column, i);
        end

        %return the end values
        new_mat = cur_mat;
        new_const_vec = cur_vec;

    end


    function [new_mat, new_const_vec] = reduce_row_at_col(param_mat, const_vec, col, row_added, row_reduced)
        %transfer over the values
        new_mat = param_mat;
        new_const_vec = const_vec;

        %get the multiple that will be multiplied by the row that will
        %be added
        m = (-param_mat(row_reduced,col))/param_mat(row_added, col);

        %change the value of the new_mat at the row that's going to be reduced 
        new_mat(row_reduced,col) = param_mat(row_reduced,col) + m*param_mat(row_added,col);

        %change the value of the new_const_vec at the row that's going
        %to be reduced
        new_const_vec(row_reduced) = const_vec(row_reduced) + m*const_vec(row_reduced);

    end

    function solution_vec = back_subst(ut_mat, new_const_vec)
        %get the size of the matrix
        [n,m] = size(ut_mat);
        %set the partial soltuions to zeroes
        part_solution_vec = zeros(new_const_vec);

        %start from the last column and work backwards
        for i = n: -1: 1,
            part_solution_vec = back_subst_for_col(ut_mat, new_const_vec, i, part_solution_vec);
        end

        %return the final answer
        solution_vec = part_solution_vec;

    end

    function new_part_solution_vec = back_subst_for_col(ut_mat, new_const_vec, column, part_solution_vec)
        %get size
        [n,m] = size(ut_mat);
        %value used for calculation
        subst_val;
        %copy over the partial solution
        new_part_solution_vec = part_solution_vec;

        %if we are at the last element then we want to do a simplified
        %calculation
        if column == n
            new_part_solution_vec(column) = new_const_vec(column)/ut_mat(column,column);
            return
        end

        %otherwise calculate the soltuion through back substituion
        for i = column + 1: n,
            subst_val = subst_val + ut_mat(column,i)*part_solution_vec(i);
        end

        %do the final calculation
        new_part_solution_vec(column) = (new_const_vec(column) - subst_val)/ut_mat(column,column);

    end

end

端

所以问题是可以调用一个函数，但是当它调用一个辅助函数时它会通过一个错误。这是我在命令窗口中调用gauss_reduce()的示例

gaussfunctions.gauss_reduce(A,B)
??? Undefined function or method 'ut_reduce' for input arguments of type 'double'.

Error in ==> gaussfunctions>gaussfunctions.gauss_reduce at 21
            [ut_mat, new_const_vec] = ut_reduce(param_mat, const_vec);

与我联系是否将所有内容都放在一个班级？我以前从未在matlab中编程，所以请原谅我，如果这是一个相当明显的错误，但我无法弄清楚发生了什么。

Answer 1

MATLAB与C ++和Java®等语言的区别在于没有传递给所有方法的特殊隐藏类实例。你必须通过该类的对象显式地为该方法。最左边的论点不需要是类实例，参数列表可以有多个对象。

来自http://www.mathworks.com/help/matlab/matlab_oop/specifying-methods-and-functions.html

此外：

classdef文件中的本地函数对于仅在该文件中使用的实用程序函数很有用。这些功能可以采取或返回作为类的实例的参数，但是，没有必要，作为在普通方法的情况下。例如，以下代码在classdef块之外定义myUtilityFcn：

这说的是，如果你在效用函数之前放置一个end，那么它将在classdef之外_并且一切都会好的。它将在文件中，但在课外。因此它被“隐藏”在世界之外（除了类定义之外没有看到）但是你的类可以使用它：

  % ... earlier code ...
        %compute the solution vector
        solution_vec = back_subst(ut_mat, new_const_vec);
    end
end  % <<<< add this one

function [ut_mat, new_const_vec] = ut_reduce(param_mat, const_vec)

%....
end
end % <<<<<< remove the last end statement

未定义的函数或方法xxx用于类型＆＃39; double＆＃39;的输入参数。

1 个答案: