数据拟合

2744557306

这段代码通过一次函数对给定的数据进行拟和，并绘制原始数据与拟和曲线的图像。下面是对代码的解释：
x = 1:8;
) G. z, v# @4 X! x5 y  uy = [15.3,20.5,27.4,36.6,49.1,65.6,87.8,117.6];
y1 = log(y);
n = 1; % 表示用一次函数拟和
! F- r' n9 z9 P& w) j0 wm = length(x);
! l5 Z. W  t! j! CA = zeros(m, n+1);

8 A. v0 [) @' o1 m% 构建矩阵 A
% w1 B( t+ K9 r+ s. R8 Efor i = 1:n+1
    A(:,i) = x.^(i-1)';
end

2 b; P9 d8 \" S, K. u1 S% 计算矩阵 B 和向量 C
B = A' * A;
: r; t( Z7 B1 q4 O- v2 MC = A' * y1';

% 求解参数向量 re
re = B \ C;

% 计算拟和曲线
/ e0 D" k& _) Z# U( M! k% Y4 N8 n: va = exp(re(1));
# Y$ _1 E; K) k0 ~& l" `4 Hb = re(2);
Y = a * exp(b * x);
0 m6 t2 U6 f4 {( R0 f
% 计算拟合误差
err = norm(Y - y, 2);

% 绘图
/ M/ G' n. [+ e) y+ `plot(x, y, 'r+', x, Y, 'g');
) G4 z: S: I& ^xlabel('x');
ylabel('y');
legend('原始数据', '拟和曲线');

% 输出拟合参数和拟合误差
7 t5 V5 c1 O! i' K6 Pdisp('拟合参数:');
: c" |" ?* G4 cdisp(['a = ', num2str(a)]);
6 y1 b' @) {: t1 E( Tdisp(['b = ', num2str(b)]);
disp(['拟合误差 (L2 范数): ', num2str(err)]);
- J) [2 g2 P- |$ i4 F( h7 G+ V) }
这段代码使用一次函数 (Y = a \cdot e^{b \cdot x}) 对给定数据进行拟和。通过最小二乘法，计算了参数 (a) 和 (b) 的值，然后计算拟和曲线 (Y)。最后，绘制了原始数据和拟和曲线的图像，并输出了拟合参数和拟合误差。