线性代数几何背景及应用幻灯片.pptVIP

Matlab程序ea110为 A=[1,-1,0,0;0,1,-1,0;0,0,1,-1;-1,0,0,1] b=[160;-40;210;-330] U0=rref([A,b]) 可以得出其最简行阶梯形矩阵 由于U0的最后一行为全零，也就是说，四个方程中实际上只有三个独立方程，所以该方程组为欠定方程，存在无穷多组解。 若以 为自由变量，方程组的解可以表示为： 如果有一些车围绕十字路的矩形区反时针绕行，流量 。都会增加，但并不影 响出入十字路的流量。这就是方程组有无穷多解的原因。 人口迁徙问题 例5 假设一个城市的总人口数是固定不变，但人口的分布情况变化如下：每年都有5％的市区居民搬到郊区；而有15％的郊区居民搬到市区。若开始有700000人口居住在市区，300000人口居住在郊区。请利用分析： （1）10年后市区和郊区的人口各是多少？ （2）30年后、50年后市区和郊区的人口各是多少？ （3）分析（2）中数据相似的原因。 解：这个问题可以用矩阵乘法来描述。令人口变量 其中 为市区人口所占比例， 为郊区人口所占比例。在n+1年的人口分布状态为： 用矩阵乘法可写成： 开始市区和郊区的人口数为 可以得到n年后市区和郊区的人口分布： 因此10年后的人口可用程序计算如下： A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; X10=A^10*X0 程序运行的结果为： 市区和郊区人口数约为：744630和255370。 无限增加时间n，市区和郊区人口之比将趋向一组常数0.25/0.75。为了弄清为什么它趋向于一个稳态值，可以将A对角化。令 ，其中Λ为对角矩阵，则有 对角矩阵的幂次可以化为元素的幂次 所以，它就很容易计算。 程序la24 % 分析n年后城市人口分布 A=[0.95,0.15;0.05,0.85]; X0=[700000;300000]; [P,lamda]=eig(A); syms n % 定义符号变量n Xn=P*lamda.^n*inv(P)*X0 % .^n对矩阵lamda中所有元素进行幂运算 计算结果为： 随n增大后一项(4/5)^n趋近于零。 多项式插值与拟合 例6 下表给出了平面坐标系中五个点的坐标。 （1）请过这五个点作一个四次多项式函数， 并求x=5时的函数值。用MATLAB绘制多项式函数曲线、通过已知点及插值点。 （2）请根据这五个点，拟合一个二次多项式函数， 并用MATLAB绘制多项式函数曲线及已知的五个点。 x 0 1 2 3 4 y -27 0 21 0 -75 解：（1）根据已知条件，把五个点的坐标值分别代入四次多项式函数，可以得到如下线性方程组： 其中矩阵： 系数矩阵A的行列式为范德蒙(Vandermonde)行列式，且五个坐标点的横坐标各不相同，则该行列式不等于零，所以方程组有唯一解。 写出程序： x=[0;1;2;3;4]; % 输入已知点坐标 y=[-27;0;21;0;-75]; % 构造范德蒙矩阵，也可用内置的vander函数 A=[x.^0,x.^1,x.^2,x.^3,x.^4]; a=A\y; % 得到适定方程组的唯一解a， 运行程序，得到 a(1)=-27,a(2)=12,a(3)=26,a(4)=-12,a(5)= 1 多项式拟合要解一个超定方程 把五个点的坐标值分别代入二次多项式函数，可以得到如下线性方程组： 其中， 该方程组有三个未知数，但有五个方程，进一步分析可以得到该方程组无解，即不存在一个二次多项式曲线刚好能过已知的五个点。MATLAB软件提供了一个利用最小二乘法解决超定方程组解的方法。求系数的公式也是a=A\y，以找到一条二次曲线来近似地描述已知5点的变化情况。 对比插值和拟合的曲线如下图 刚体的平面运动 例7 用平面坐标系中的一个闭合图形来描述刚体，用一个矩阵X来表示它。X的一列表示刚体一个顶点的坐标。为了使图形闭合，X的最后一列和第一列相同；为了实现刚体的平移运算，给矩阵X添加元素值都为1的一行，使矩阵X的形状为3×n。 若有矩阵： 可以证明，矩阵是刚体X沿x轴正方向平移，沿y轴正方向平移后的结果；矩阵是刚体X以坐标原点为中心逆时针转动t弧度的结果。 计算实例 x 0 4 6 10 8 5 3.5 6.1 6.5 3.2 2 0 y 0 14 14 0