myline.m

MDPSAS工具箱是马里兰大学开发的

function A = myline(x,y)
% myline.m Constructor method for myline class. Called as
%
%          A = myline(x,y)
%
% INPUT PARAMETERS
%         x : x values or a (slope)
%         y : y values or b (intercept)

if nargin == 0; x = 0; y = 0; end % default values
if prod(size(x)) == 1
    A.a = x; A.b = y;
else 
    B(:,1) = x(:); B(:,2) = 1;
% Note how the line is computed by least squares, equivalent
% to coeff = (B'*B) \ (B'*y)
    coeff = B\y(:);
    A.a = coeff(1); A.b = coeff(2);
end
A = class(A,'myline');