西安电子科技大学高等代数机算与应用作业题答案.doc

高等代数机算与应用作业题学号： 姓名： 成绩： 一、机算题1利用函数rand和函数round构造一个55的随机正整数矩阵A和B。 a=round(rand(5)a = 0 0 1 1 1 1 1 0 1 0 1 0 1 0 1 0 1 1 0 0 0 0 1 0 1 b=round(rand(5)b = 0 0 0 0 0 0 1 1 0 1 0 1 1 1 0 1 1 0 1 0 0 1 1 1 0(1)计算AB，AB和6A a+bans = 0 0 1 1 1 1 2 1 1 1 1 1 2 1 1 1 2 1 1 0 0 1 2 1 1 a-bans = 0 0 1 1 1 1 0 -1 1 -1 1 -1 0 -1 1 -1 0 1 -1 0 0 -1 0 -1 1 6*aans = 0 0 6 6 6 6 6 0 6 0 6 0 6 0 6 0 6 6 0 0 0 0 6 0 6(2)计算，和 (a*b)ans = 1 1 0 0 0 3 2 2 2 2 2 1 2 2 2 3 1 2 1 2 0 1 0 1 0 b*aans = 1 1 0 0 0 3 2 2 2 2 2 1 2 2 2 3 1 2 1 2 0 1 0 1 0 (a*b)100ans = 1.0e+078 * 1.4732 7.6495 6.1764 5.5225 2.1271 1.0117 5.2535 4.2418 3.7927 1.4608 0.9229 4.7921 3.8692 3.4596 1.3325 0.9229 4.7921 3.8692 3.4596 1.3325 0.9229 4.7921 3.8692 3.4596 1.3325(3)计算行列式，和 det(a)ans = 1 det(b)ans = 0 det(a*b)ans = 0(4)若矩阵A和B可逆，计算和 inv(a)ans = 0 0 1.0000 0 -1.0000 -1.0000 1.0000 -1.0000 0.0000 2.0000 1.0000 -1.0000 1.0000 1.0000 -2.0000 1.0000 0.0000 -0.0000 -0.0000 -1.0000 -1.0000 1.0000 -1.0000 -1.0000 3.0000b不存在逆矩阵(5)计算矩阵A和矩阵B的秩。 rank(a)ans = 5 rank(b)ans = 3(6)生成一个6行5列秩为3的矩阵，并求其最简阶梯形。a=1,1,1,1,1;0,1,1,1,1;0,0,1,1,1;0,0,0,0,0;0,0,0,0,0;0,0,0,0,0a = 1 1 1 1 1 0 1 1 1 1 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 rref(a)ans = 1 0 0 0 0 0 1 0 0 0 0 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 02求解下列方程组(1)求非齐次线性方程组的唯一解。 a=2,1,2,4;-14,17,-12,7;7,7,6,6;-2,-9,21,-7a = 2 1 2 4 -14 17 -12 7 7 7 6 6 -2 -9 21 -7 b=5;8;5;10b = 5 8 5 10 inv(a)*bans = -0.8341 -0.2525 0.7417 1.3593(2)求非齐次线性方程组的通解。 a=5,9,7,2,8;4,22,8,25,23;1,8,1,8,8;2,6,6,9,7a = 5 9 7 2 8 4 22 8 25 23 1 8 1 8 8 2 6 6 9 7 b=4;9;1;7b = 4 9 1 7 c=a,bc = 5 9 7 2 8 4 4 22 8 25 23 9 1 8 1 8 8 1 2 6 6 9 7 7 d=rref(c)d = 1.0000 0 0 -4.1827 -0.8558 -1.6635 0 1.0000 0 1.3269 1.0577 0.1346 0 0 1.0000 1.5673 0.3942 1.5865 0 0 0 0 0 0对应齐次方程组的基础解系为a1= 4.1827 a2= 0.8558 -1.3269 -1.0577 -1.5673 -0.3942 1.0000 0 0 1.0000非齐次线性方程组的特解为b= -1.6635 0.1346 1.5865 0 0原方程组的通解为x=k1a1+k2a2+b3已知向量组，求出它的最大无关组，并用该最大无关组来线性表示其它向量。 a=3,1,2,9,0;4,1,3,3,8;0,0,0,2,-2;8,2,6,1,21;3,2,1,2,10a = 3 1 2 9 0 4 1 3 3 8 0 0 0 2 -2 8 2 6 1 21 3 2 1 2 10 rref(a)ans = 1 0 1 0 2 0 1 -1 0 3 0 0 0 1 -1 0 0 0 0 0 0 0 0 0 0最大无关组为(a1,a2,a4).a3=a1-a2; a5=2a1+3a2-a3.4求向量空间中向量在基下的坐标 a=1,2,3;0,1,2;0,0,1a = 1 2 3 0 1 2 0 0 1 b=3;2;5b = 3 2 5 inv(a)*bans = 4 -8 5a的坐标为(4,-8,5).5求下列矩阵的特征值和特征向量，并判断其正定性。(1)；(2) a=1,2,3;2,5,6;3,6,25a = 1 2 3 2 5 6 3 6 25 v,d=eig(a)v = 0.9357 0.3279 0.1303 -0.3518 0.8961 0.2706 -0.0280 -0.2990 0.9538d = 0.1582 0 0 0 3.7297 0 0 0 27.1121a为正定矩阵. b=-20,3,1;3,-10,-6;1,-6,-22b = -20 3 1 3 -10 -6 1 -6 -22 v,d=eig(b)v = -0.3810 0.9059 0.1850 0.4005 -0.0186 0.9161 0.8334 0.4231 -0.3557d = -25.3404 0 0 0 -19.5947 0 0 0 -7.0649b为负定矩阵.6用正交变换法将下列二次型化为标准形。其中“”为自己学号的后三位。k1=0;k2=0;k3=1;二次型的矩阵形式为 1 0 0 0 2 0.5 0 0.5 3 a=1,0,0;0,2,0.5;0,0.5,3a = 1.0000 0 0 0 2.0000 0.5000 0 0.5000 3.0000 v,d=eig(a)v = 1.0000 0 0 0 -0.9239 0.3827 0 0.3827 0.9239d = 1.0000 0 0 0 1.7929 0 0 0 3.2071 c=orth(v)c = 1.0000 0 0 0 -0.9239 0.3827 0 0.3827 0.9239作正交变换 x=cy,则得f的标准型:f=y12+1.7929y22+3.2071y32二、应用题1在钢板热传导的研究中，常常用节点温度来描述钢板温度的分布。假设下图中钢板已经达到稳态温度分布，上下、左右四个边界的温度值如图所示，而表示钢板内部四个节点的温度。若忽略垂直于该截面方向的热交换，那么内部某节点的温度值可以近似地等于与它相邻四个节点温度的算术平均值，如。请计算该钢板的温度分布。4T1-T2-T3=70-T1+4T2-T4=50-T1+4T3-T4=50-T2-T3+4T4=30 A=4,-1,-1,0;-1,4,0,-1;-1,0,4,-1;0,-1,-1,4A = 4 -1 -1 0 -1 4 0 -1 -1 0 4 -1 0 -1 -1 4 B=70;50;50;30B = 70 50 50 30 inv(A)*Bans = 30.0000 25.0000 25.0000 20.00001处温度为30,2处温度为25,3处温度为25，4处温度为20.2 下表给出了平面坐标系中六个点的坐标。x012345y260262941302请过这六个点作一个五次多项式函数，并求当时的函数值。 a=polyfit(0,1,2,3,4,5,2,6,0,26,294,1302,5)a = 1.0000 -3.0000 -0.0000 1.0000 5.0000 2.0000 y=polyval(a,6)y = 3.9560e+0033. 李博士培养了一罐细菌，在这个罐子里存放着A、B、C三类不同种类的细菌，最开始A、B、C三种细菌分别有、2、3个。但这些细菌每天都要发生类型转化，转化情况如下：A类细菌一天后有5的变为B类细菌、15的变为C细菌；B类细菌一天后有30的变为A类细菌、10的变为C类细菌；C类细菌一天后有30的变为A类细菌、20的变为B类细菌。请利用MATLAB软件分析：(1)一周后李博士的A、B、C类细菌各有多少个？(2)两周后和三周后李博士的A、B、C类细菌各有多少个？(3)分析在若干周后，李博士的各种细菌的个数几乎不发生变化的原因。第0天s0=a b c第一天s1=0.85a+0.30b+0.30c 0.60b+0.05a+0.20c 0.50c+0.15a+0.10b则其变换矩阵为c=0.85 0.05 0.150.30 0.60 0.100.30 0.20 0.50 c=0.85,0.05,0.15;0.30,0.60,0.10;0.30,0.20,0.50c = 0.8500 0.0500 0.1500 0.3000 0.6000 0.1000 0.3000 0.2000 0.5000 s0=1,2,3s0 = 1 2 3 s0*c7ans =4.4048 1.2320 1.4726 s0*c14ans = 5.4928 1.4867 1.8309 s0*c21ans = 6.8098 1.8426 2.2699单位:108个4. 一个混凝土生产企业可以生产出三种不同型号的混凝土，它们的具体配方比例如表所示。表11.3混凝土的配方型号1混凝土型号2混凝土型号3混凝土水101010水泥222618砂323129石子536450灰058(1)分析这三种混凝土是否可以用其中的两种来配出第三种？(2)现在有甲、乙两个用户要求混凝土中含水、水泥、砂、石子及灰的比例分别为：24，52，73，133，12和36，75，100，185，20。那么，能否用这三种型号混凝土配出满足甲和乙要求的混凝土？如果需要这两种混凝土各500吨，问三种混凝土各需要多少？ . a=10,22,32,53,0;10,26,31,64,5;10,18,29,50,8a = 10 22 32 53 0 10 26 31 64 5 10 18 29 50 8 rank(a)ans = 3不能用其中两种配出第三种. a=aa = 10 10 10 22 26 18 32 31 29 53 64 50 0 5 8 b=24;52;73;133;12b = 24 52 73 133 12 c=a,bc = 10 10 10 24 22 26 18 52 32 31 29 73 53 64 50 133 0 5 8 12 rank(a)ans = 3 rank(c)ans = 3能满足甲的要求. rref(c)ans = 1.0000 0 0 0.6000 0 1.0000 0 0.8000 0 0 1.0000 1.0000 0 0 0 0 0 0 0 00.6吨1号混凝土,0.8吨2号混凝土,1吨3号混凝土可配成2.4吨甲型混凝土.需1号混凝土500/2.4*0.6=125吨,2号混凝土500/2.4*0.8=167吨,3号混凝土500/2.4*1=208吨. b=36;75;100;185;20b = 36 75 100 185 20 c=a,bc = 10 10 10 36 22 26 18 75 32 31 29 100 53 64 50 185 0 5 8 20 rank(a)ans = 3 rank(c)ans = 4不能配成乙型混凝土.5. 在某网格图中，每个中间节点的值与其相邻的上、下、左、右四个节点的值有如下关系：。其中系数，网格图的上下、左右四个节点的值如图所示。如，请计算该网格节点1，2，3，4的值。2134a aC Cbbdd令A=a;B=b;C=c;D=d;a=;b=;c=;d=T1=*A+*T3+*B+*T2T2=*A+*T4+*T1+*DT3=*T1+*C+*B+*T4T4=*T2+*C+*T3+*D即T1-*T3-*T2=*A+*BT2-*T4-*T1=*A+*DT3-*T1-*T4=*C+*BT4-*T2-*T3=*C+*D x=sym(1,-d,-b,0;-c,1,0,-b;-a,0,1,-d;0,-a,-c,1)x = 1, -d, -b, 0 -c, 1, 0, -b -a, 0, 1, -d 0, -a, -c, 1 y=sym(a*A+c*B;a*A+d*D;b*C+c*B;b*C+d*D)y = a*A+c*B a*A+d*D b*C+c*B b*C+d*D inv(x)*yans =-(-1+d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+c*B)+d*(1-d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+d*D)-b*(-d*c-1+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+c*B)+2*b*d/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+d*D) c*(1-d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+c*B)-(-1+d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+d*D)+2*b*c/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+c*B)-b*(-d*c-1+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+d*D) -a*(-d*c-1+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+c*B)+2*d*a/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+d*D)-(-1+d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+c*B)+d*(1-d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+d*D) 2*c*a/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+c*B)-a*(-d*c-1+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(a*A+d*D)+c*(1-d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+c*B)-(-1+d*c+a*b)/(1-2*d*c-2*a*b+d2*c2-2*c*a*b*d+a2*b2)*(b*C+d*D)6. 假如一个数据库包含以下10种图书：B1：高等代数，B2:线性代数，B3：工程线性代数，B4：初等线性代数，B5：线性代数及其应用，B6：MATLAB在数值线性代数中应用，B7：矩阵代数及其应用，B8：矩阵理论，B9：线性代数及MATLAB入门，B10：基于MATLAB的线性代数及其应用。而检索的6个关键词按拼音字母次序排列为：“代数，工程，矩阵，MATLAB，数值，应用”。读者1的检索策略为：“代数，MATLAB”；读者2的检索策略是：“代数，应用”。请用矩阵运算来为这两位读者检索图书。以关键词为行,图书类别为列,建立数据库矩阵:1 1 1 1 1 1 1 0 1 10 0 1 0 0 0 0 0 0 00 0 0 0 0 0 1 1 0 00 0 0 0 0 1 0 0 1 10 0 0 0 0 1 0 0 0 00 0 0 0 1 1 1 0 0 1读者1的检索矩阵:1 0 0 1 0 0a=1,1,1,1,1,1,1,0,1,1;0,0,1,0,0,0,0,0,0,0;0,0,0,0,0,0,1,1,0,0;0,0,0,0,0,1,0,0,1,1;0,0,0,0,0,1,0,0,0,0;0,0,0,0,1,1,1,0,0,1a = 1 1 1 1 1 1 1 0 1 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 1 0 0 1 1 0 0 0 0 0 1 0 0 0 0 0 0 0 0 1 1 1 0 0 1 b=1,0,0,1,0,0b = 1 0 0 1 0 0 b*aans = 1 1 1 1 1 2 1 0 2 2找到图书B6,B9,B10.读者2的检索矩阵:1 0 0 0 0 1b = 1 0 0 0 0 1 b*aans = 1 1 1 1 2 2 2 0 1 2找到图书B5,B6,B7,B10.7 某城市有如图所示的9节点交通图，每一条道路都是单行道，图中数字表示某一个时段该路段的车流量。若针对每一个十字路口，进入和离开的车辆数相等。请计算每两个相邻十字路口间路段上的交通流量。若已知AB段和FO段在修路，即：，又已知，。求此时各个路段的交流流量。单行道9节点交通流图x1+x8=200x1-x2+x10=300x2-x3=166-x7+x8+x9=350-x9-x10+x11+x12=0x3-x4-x12= -500x6+x7=150x5-x6-x11= -400x4-x5= 134x8=0x12=0x10=300x11=660 a1= 1 0 0 0 0 0 0 1 0 0 0 0 200a1 = 1 0 0 0 0 0 0 1 0 0 0 0 200 a2=1 -1 0 0 0 0 0 0 0 1 0 0 300a2 = 1 -1 0 0 0 0 0 0 0 1 0 0 300 a3=0 1 -1 0 0 0 0 0 0 0 0 0 166a3 = 0 1 -1 0 0 0 0 0 0 0 0 0 166 a4=0 0 0 0 0 0 -1 1 1 0 0 0 350a4 = 0 0 0 0 0 0 -1 1 1 0 0 0 350 a5=0 0 0 0 0 0 0 0 -1 -1 1 1 0a5 = 0 0 0 0 0 0 0 0 -1 -1 1 1 0 a6=0 0 1 -1 0 0 0 0 0 0 0 -1 -500a6 = 0 0 1 -1 0 0 0 0 0 0 0 -1 -500 a7=0 0 0 0 0 1 1 0 0 0 0 0 150a7 = 0 0 0 0 0 1 1 0 0 0 0 0 150 a8=0 0 0 0 1 -1 0 0 0 0 -1 0 -400a8 = 0 0 0 0 1 -1 0 0 0 0 -1 0 -400 a9=0 0 0 1 -1 0 0 0 0 0 0 0 134a9 = 0 0 0 1 -1 0 0 0 0 0 0 0 134 a10=0 0 0 0 0 0 0 1 0 0 0 0 0a10 = 0 0 0 0 0 0 0 1 0 0 0 0 0 a11=0 0 0 0 0 0 0 0 0 0 0 1 0a11 = 0 0 0 0 0 0 0 0 0 0 0 1 0 a12=0 0 0 0 0 0 0 0 0 1 0 0 300a12 = 0 0 0 0 0 0 0 0 0 1 0 0 300 a13=0 0 0 0 0 0 0 0 0 0 1 0 660a13 = 0 0 0 0 0 0 0 0 0 0 1 0 660 a=a1;a2;a3;a4;a5;a6;a7;a8;a9;a10;a11;a12;a13a = 1 0 0 0 0 0 0 1 0 0 0 0 200 1 -1 0 0 0 0 0 0 0 1 0 0 300 0 1 -1 0 0 0 0 0 0 0 0 0 166 0 0 0 0 0 0 -1 1 1 0 0 0 350 0 0 0 0 0 0 0 0 -1 -1 1 1 0 0 0 1 -1 0 0 0 0 0 0 0 -1 -500 0 0 0 0 0 1 1 0 0 0 0 0 150 0 0 0 0 1 -1 0 0 0 0 -1 0 -400 0 0 0 1 -1 0 0 0 0 0 0 0 134 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 300 0 0 0 0 0 0 0 0 0 0 1 0 660 rank(a)ans = 12 rref(a)ans = 1 0 0 0 0 0 0 0 0 0 0 0 200 0 1 0 0 0 0 0 0 0 0 0 0 200 0 0 1 0 0 0 0 0 0 0 0 0 34 0 0 0 1 0 0 0 0 0 0 0 0 534 0 0 0 0 1 0 0 0 0 0 0 0 400 0 0 0 0 0 1 0 0 0 0 0 0 140 0 0 0 0 0 0 1 0 0 0 0 0 10 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 360 0 0 0 0 0 0 0 0 0 1 0 0 300 0 0 0 0 0 0 0 0 0 0 1 0 660 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0此时各路段流量为200,200,34,534,400,140,10,0,360,300,660,0.8用MATLAB软件实现以下操作：(1)构造一个直角三角形刚体矩阵X；(2)先对刚体逆时针转动，然后再向下移动20，向右移动20； (3)先对刚体向下移动20，向右移动20，然后对刚体逆时针转动。逆时针转动45的变换矩阵为cos(45) -sin(45) 0sin(45) cos(45) 00 0 1向下移动20,向右移动20的变换矩阵为1 0 00 1 020 -20 1构造一个直角三角形刚体矩阵x: a=0 0 1;0 4 1;3 0 1a = 0 0 1 0 4 1 3 0 1 b=0.7071 -0.7071 0;0.7071 0.7071 0;0 0 1b = 0.7071 -0.7071 0 0.7071 0.7071 0 0 0 1.0000 c=1 0 0;0 1 0;20 -20 1c = 1 0 0 0 1 0 20 -20 1先对刚体逆时针转动，然后再向下移动20，向右移动20: a*b*cans = 20.0000 -20.0000 1.0000 22.8284 -17.1716 1.0000 22.1213 -22.1213 1.0000先对刚体向下移动20，向右移动20，然后对刚体逆时针转动: a*c*bans = 0 -28.2840 1.0000 2.8284 -25.4556 1.0000 2.1213 -30.4053 1.00009 (1)甲方收到与之有秘密通信往来的乙方的一个密文信息，密文内容如下： WOWUYSBACPGZSAVCOVKPEWCPADKPPABUJCQLYXQEZAACPP已知密钥为能否知道这段密文的意思？ a=1 2;0 3a = 1 2 0 3 A=det(a)A = 3 alpha=mod(inv(a).*27,26)alpha = 1 8 0 9 m=23 15;23 21;25 19;2 1;