矩阵论课外报告.docVIP

—、报告 在己知曲线大致模型的情况下，运用曲线拟合最小二乘法，使得观测数据与 曲线模型数据之间的误差平方和最小。进而求得曲线的模型参数，并由所求的曲 线模型进行分析预测。 二、题目内容 一颗导弹从敌国发射，通过雷达我们观测到了它的飞行轨迹，具体有如下数 据： 水平距离/m 0 250 500 750 1000 高度/m 0 8 15~~19~~20 我国军情处分析得出该导弹沿抛物线轨道飞行。 问题：预测该导弹在什么水平距离着地。 三、基本术语 1.内积 设V是实数域R上的线性空间，如果V屮任意两个向量4,0都按某一个确 定的法则对应于惟一确定的实数，记作(％0),并且(％0)满足 对任意的 有(G,0) = (0,Q) 对任意的a.p./g V ,有(G+0，7)= (a,/) + (0,y) 对任意的 k = R,a,/3eV 有(ka,0) = k(a,0) 对任意的ae V ,有(a,a)0 o当且仅当q = 0时，(a,a) = 0 则称(％0)为向量的内积。如无特殊说明的，我们认为对任意向量 心=（坷卫2，???,色）,0 =（勺,优,???,?），其内积（G,0）为 （% 0） = a}b} + a2b2 + ? ?? + anbn 2.范数 如果V是数域K上的线性空间，且对于V的任以向量力，对应于一个实数函 数惻，它满足如下三个条件。 i.非负性当力H0时||/|| 0 ；当力=0时，||/|| = 0 ； ii? 齐次性 ||?z|| = |^|||^||,/eV ； iii?三角不等式|力+ J ||力|| + ||了 ||,力,^eV； 则称||力||为V上力的范数。 可以证明对于向量“(—6)的长度 kll=』纤+|纤+???+￡『 是一种范数，我们称为2■范数，记为||Z||2 o 3.线性方程组 设有n个未知数m个方程的线性方程组 q內+再2兀2+??? + %￡二勺 a2lx^a22x2+-^a2nxn=b2 0”+佥2兀2+???+ 4””￡=4 可以写成以向量尢为未知元的向量方程 Ax = h 则A为该方程的系数矩阵，B = (A,b)为增广矩阵。该线性方程有解的条件如下 i. 当A的秩A)和B的秩R(B)满足R(A) R(B)时，该方程无解 ii. 当R(A) = R(B)= 〃时，该方程有唯一解。 iii. 当R（A） = R（B） b时，该方程有无穷解。 四、基本原理 对于一组给定的实验数据（兀，y, ） （i = 0丄…，m ）,要求出自变量兀与因 变量y的函数关系y = S（x）,由于观测数据总有误差，所以不要求y = S（x）通过 已知点（壬，x ）（心0,1,…，加），而只要求在给定点兀?上的误差 Q = Sa）-x（i = O,l,…，加）的平方和工疔最小。 1=0 若已知 S （兀）=绳％ （兀）+ 4% （兀）+…+ 叱（兀） 这里（p°g（P\g ???,（/）是线性无关的函数族，假定有一组数据 ｛（兀,必），，=0丄?「加｝，要求y = S（x）使/ =（绳卫］,…最小，其中 I =(6f0,6Zp---,^,) = ^[S(X/-X)]2/=! 这就是最小二乘逼近，得到的拟合曲线为y = SCr）,这种方法称为曲线拟合的最 小二乘法。 / = （°0,4,???,%）实际上是关于兔卫］,…卫”的多元函数，求/ = 如…，陽） 的最小值就是求多元函数/ = /（d（）,q,???“）的极值，由极值的必要条件，可得 /Uo) /Uo)_ 0_ 我们令％ = ■ ■ ■ ，尸 1 ■ ■ _(pM_ ? _x,_ （Z = 0,1,…,ni）,即%是将实验数据Oo,…兀“）带入函 31 =2 工 )% (兀)+ …+ an(pn (兀)一 ;炳(%,) = 0 /=0 (k = 0,1,???,/!) 数0心）所得的列向量，y是实验数据（y°,y,…，儿）的列向量。则上式可改写为 (%,0?do+(%，輕M+…+(%，件)Q” =(y，0? 伙=0,1,…，72) 这是关于参数…4的线性方程组，用矩阵表示为 （%,%） （%,%）…（％,%）、 Z 、 0（）、 ?%）、 （俗如（％,%）…（％?） ? ? ? ? ? ? ? a\ ■ ■ — （y，%） ■ ■ ? ? @,%） （%,%）…（％,%）丿 ■ 5 ■ 力 该线性方程称为法方程。对于该方程的系数矩阵A和增广矩阵B,当R (A) =R (B)二n时，该方程有唯一解。记方程解为％ =g； (k二0,1,…,n),从而得到最小 二乘拟合曲线 y = S* (兀)=a；% (x) + % (x) + …+ an(pn (x) 可以证明对任意(伽?,…,“，G R卄 ,有/(G；,G：y,d：) /(q