最小二乘拟合平面与直线matlab.docVIP

利用Matlab实现直线和平面的拟合 1、直线拟合的matlab代码% Fitting a best-fit line to data, both noisy and non-noisyx = rand(1,10);n = rand(size(x)); % Noisey = 2*x + 3; % x and y satisfy y = 2*x + 3yn = y + n; % x and yn roughly satisfy yn = 2*x + 3 due to the noise% Determine coefficients for non-noisy line y=m1*x+b1Xcolv = x(:); % Make X a column vectorYcolv = y(:); % Make Y a column vectorConst = ones(size(Xcolv)); % Vector of ones for constant termCoeffs = [Xcolv Const]\Ycolv; % Find the coefficientsm1 = Coeffs(1);b1 = Coeffs(2);% To fit another function to this data, simply change the first % matrix on the line defining Coeffs% For example, this code would fit a quadratic % y = Coeffs(1)*x^2+Coeffs(2)*x+Coeffs(3)% Coeffs = [Xcolv.^2 Xcolv Const]\Ycolv; % Note the .^ before the exponent of the first term% Plot the original points and the fitted curvefigureplot(x,y,ro)hold onx2 = 0:0.01:1;y2 = m1*x2+b1; % Evaluate fitted curve at many pointsplot(x2, y2, g-)title(sprintf(Non-noisy data: y=%f*x+%f,m1,b1))% Determine coefficients for noisy line yn=m2*x+b2Xcolv = x(:); % Make X a column vectorYncolv = yn(:); % Make Yn a column vectorConst = ones(size(Xcolv)); % Vector of ones for constant termNoisyCoeffs = [Xcolv Const]\Yncolv; % Find the coefficientsm2 = NoisyCoeffs(1);b2 = NoisyCoeffs(2);% Plot the original points and the fitted curvefigureplot(x,yn,ro)hold onx2 = 0:0.01:1;yn2 = m2*x2+b2;plot(x2, yn2, g-)title(sprintf(Noisy data: y=%f*x+%f,m2,b2))2、平面拟合matlab代码x = rand(1,10);y = rand(1,10);z = (3-2*x-5*y)/4; % Equation of the plane containing % (x,y,z) points is 2*x+5*y+4*z=3Xcolv = x(:); % Make X a column vectorYcolv = y(:); % Make Y a column vectorZcolv = z(:); % Make Z a column vectorConst = ones(size(Xcolv)); % Vector of ones for constant termCoefficients = [Xcolv Ycolv Const]\Zcolv; % Find the coefficientsXCoeff = Coefficients(1); % X coefficientYCoeff = Coefficients(2); % X coefficientCCoeff = Coefficients(3); % constant term% Using the above variables, z = XCoeff * x + YCoeff *