第四章MATLAB线性代数.pptVIP

*/16 小行星的轨道方程 特征值问题及应用 离散数据的多项式拟合 第四章 线性代数 ? ? ? 例4.2 小行星轨道方程 以太阳为坐标原点,测得小行星坐标 x 4.5596 5.0816 5.5546 5.9636 6.2756 y 0.8145 1.3685 1.9895 2.6925 3.5265 a1x12 + 2a2x1y1 + a3 y12 +2a4 x1 + 2a5 y1 = –1 a1x22 + 2a2x2y2 + a3 y22 +2a4 x2 + 2a5 y2 = –1 a1x32 + 2a2x3y3 + a3 y32 +2a4 x3 + 2a5 y3 = –1 a1x42 + 2a2x4y4 + a3 y42 +2a4 x4 + 2a5 y4 = –1 a1x52 + 2a2x5y5 + a3 y52 +2a4 x5 + 2a5 y5 = –1 椭圆二次曲线方程 a1x2 + 2a2xy + a3 y2 +2a4 x + 2a5 y + 1 = 0 ? Az = b ? MATLAB 求解方程组方法：A\b 创建方程组系数矩阵方法： A=[X.^2, 2*X.*Y, Y.^2, X, Y] X=[4.5596;5.0816;5.5546;5.9636;6.2756]; Y=[0.8145;1.3685;1.9895;2.6925;3.5265]; A=[X.*X,2*X.*Y,Y.*Y,2*X,2*Y]; E=[-1;-1;-1;-1;-1]; z=A\E; a1=z(1);a2=z(2);a3=z(3);a4=z(4);a5=z(5); syms x y F=a1*x^2+2*a2*x*y+a3*y^2+2*a4*x+2*a5*y+1; ezplot(F,[-1,6.5,-1.5,6]) hold on,plot(X,Y,ro) 程序文件 mlab42.m a1x2 + 2a2xy + a3 y2 +2a4 x + 2a5 y + 1 = 0 或 XTAX + 2XTb + 1 = 0 平移变换,X = X0 + Z, 其中X0(椭圆中心)待定 (X0 + Z)TA(X0 + Z) + 2(X0 + Z)Tb + 1 = 0 ZTAZ + 2ZT(AX0 + b) +( X0TAX0 + 2X0Tb + 1) = 0 F = X0TAX0 + 2X0Tb + 1 0 消去方程中一次项,令 二次方程化简: ZTAZ + F = 0 AX0 + b = 0 ? 设?1，?2 是 A 的特征值,对应特征向量为 令 ? 令 ? ? 其中 变量 t 的离散数据 ? 轨道数据 原坐标下的轨道离散数据 ? ? MATLAB解算特征值问题方法 lamda=eig(A) —— 计算A的特征值,这里lamda是A的全部特征值构成的列向量。 [P,D]=eig(A) ——计算出A的全部特征值和对应的特征向量. 其中, D是对角矩阵,保存矩阵A的全部特征值; P是满阵, P的列向量构成对应于D的特征向量组。 矩阵特征值问题 A是n阶方阵,求非零向量 和数 使得 称 为特征向量,称 为特征值. x=[4.5596 5.0816 5.5546 5.9636 6.2756]; y=[0.8145 1.3685 1.9895 2.6925 3.5265]; D=[x.^2 2*x.*y y.^2 2*x 2*y]; E=[1;1;1;1;1];e=-D\E; A=[e(1) e(2);e(2) e(3)];b=[e(4);e(5)]; X0=-A\b; F=X0*A*X0+2*X0*b+1; [U d]=eig(A);a=sqrt(-F/d(1,1));b=sqrt(-F/d(2,2)); t=linspace(0,2*pi,2000); u=a*cos(t);v=b*sin(t); X=U*[u;v]+X0*ones(1,2000);xt=X(1,:);yt=X(2,:); plot(xt,yt,x,y,b*,0,0,ro),hold on comet(xt,yt) 例4.5 出租汽车问题。 出租汽车公司在仅有A城和B城的海岛上,设了A,B两营业部。如果周一A城有120辆可出租汽车，而B城有150辆。统计数据表明，平均每天A城营业部汽车的10%被顾客租用开到B城?,B城营业部汽车的12%被开到了A城。假设所有汽车正常，试计算一周后两城的汽车数量。寻找方案使每天汽车正常流动而A城和B城的汽车数量不增不减。 设第n天A城营业部汽车数为x1(n)，B城营业部汽车数为x2(n)。 则有 营业部汽车