线性代数 习题六.doc

第六章习题解答与问题 一、习题解答 用最小二乘法求解超定方程组 2x + 4y = 11 3x-5y = 3 x + 2y = 6 4x + 2y = 14 解：超定方程组的矩阵形式为 将方程两端同乘以系数矩阵的转置矩阵，可得正规方程组 「30 3_ X ~93~ _3 49_ _y_ 69 解之，得 * = 2.9774, 1.2259。 观测一个作直线运动的物体，测得以下数据: 时间t 0 0.9 1.9 3.0 3.9 5.0 距离S 0 10 30 50 80 110 在表中，时间单位为秒，距离单位为米。假若加速度为常数，求这物体的初速度和加速度。 解：设物体运动的初速度和加速度分别为vo和a,初始时刻距离为0,则距离函数为 1 , S(/) = Vgt + ~Clt~ 用后5个点的数据作曲线拟合 t 0.9 1.9 3.0 3.9 5.0 S 10 30 50 80 110 可得，v0= 10.6576, a= 4.6269 用最小二乘法求一个形如y-的经验公式，使与下列数据相拟合 X 1 2 3 4 y 60 30 20 15 贝!J z = In A + Bxo 数 据变换如下 X 1 2 3 4 z = lny 4.0943 3.4012 2.9957 2.7081 山最小二乘法作线性拟合得，In A = 4.4409, B = -0.4564?所以A =84.8528=故，所求经 难公式为 二 84.25 04564x o 4已知实验观测数据(花，yi) (I = L 2,…，m)o令 ( (p\ (x) = X 取拟合函数为 0（x） = ao0o（x） + ai? （x） 试利用曲线拟合的最小二乘法确定组合系数ao, ai （推导出计算公式）。 解：记 00=［00（兀）00 （*2）??- 0o（x,”）r 01= Wig） eg） ??- ?（x,”）r y = lyi y2 ??- yJ 显然，00是元素全为“1”的列向量。将所有实验数据的X坐标代入拟合函数，并令其分 别等于实验数据的Y坐标值，得超定方程组 将方程组两端同乘以矩阵［00 ?r，得正规方程组 (00,00)(0o，0i) X 「風，y）~ _（01，00）（01，010）_ @, y)_ 1 m 记x=—Yx.,由于系数矩阵中两个非对角元素为 m曰 TOC \o 1-5 \h \z m m （0o，0i）= （0i ,00）=工（兀一元）=工兀—?元=0 Z=1 1=1 所以 _風,刃 _ （?,刃 vtn — ， ⑷— （观，观） （D 5.对某个物体的长度测量”次后，得n个近似值x 1，X2,……x,?,通常取平均值作为所 求长度的值。试用最小二乘法原理说明其理由。 解：利用最小二乘原理，设物体的长度为X,记 6k = x - X/, （ k = 1, 2, , m） 则残差平方和为 m S（x） = ^（x —xj k=l 为了求上面函数极小值，由极值必要条件，令S（x） = O,得 m Y（x-xj = o k=l 由此得 解：利用勒让德多项式作基函数，即P(x) = aopo(x) + aipi(x) + t?2P2(x) + a3P3(x)， 中 Po(x)=l, Pi(x) = x, P-. (x) — —x2 — , p (x) — —X’ — 2 2 2 2 2 利用正交性，得系数为 (P”，/) 2n + l ^pn(x)f(x)dx ( n = 0, 1, 2, 3) ￡p0(x) f(x)dx = Xdx = e-e-1 ^p1(x)f(x)dx = ^xexdx = 2e~r ^p2(x)f(x)dx = ￡弓* - ￡)e*dx = e-7e ^p3(x)f (x)dx = x3 ~^x)exdx = 37e_1 -5e a0 a0 = — x(e-e) a 1.1752, ax = —x2e_1 ? 1.1036, 25 2 ? 7 a2 = —x(e -7eJ)a0.3578, a3 =—x(37e — 5e) a0.0705 2 2 所以, P(x) = 1.1752+ 1.1036 x+ 0.3578(|X2 -|) +0.0705(|x3 -|x) =0.9963+0.9978 x + 0.5367 x2 + 0.1762 x3 7.在著名的高次插值的龙格反例中，在区间［_5, 5］上的10次拉格朗日插 5 + x 值出现振荡现象。为了使插值余项极小化，可以利用切比雪夫多项式的极性。试推导11次 切比雪夫多项式零点所对应［-5, 5］的上的插值结点。 解：山11次切比雪夫多项式零点，得 x” =5c0S(2 + l?)(k = 0, 1, 2, , 10) k 11 2 二、例题 1.