matlab文件 南邮数学实验3答案.docVIP

第三次练习题 书上习题：(实验九) 2，3，4，9，10，12，14，16 练习2 对，，求出的通项. 输入： A=[4,2;1,3]; A=sym([4,2;1,3]); [P,D]=eig(A);inv(P); Xn=P*D.^x*inv(P)*[1;2] 输出： Xn = -2^x+2*5^x 2^x+5^x 练习3 对于练习1中的，，求出的通项. 输入： A=[0.4,0.2;0.1,0.3]; A=sym([4,2;1,3]); [P,D]=eig(A);inv(P); Xn=P*D.^x*inv(P)*[1;2] 输出： Xn = -2^x+2*5^x 2^x+5^x 练习4 对随机给出的，观察数列.该数列有极限吗? 输入：temp.m clear all; clc; x1 = rand(1); x2 = rand(1); syms a m n x; a = sym([4,2;1,3]); [m, n] = eig(a); for i = 1:10 b = eval(limit([x1,x2]*m*n.^x*inv(m),10*i)); b(1)/b(2) end 输出： ans = 1.0000 ans = 1.0000 ans = 1.0000 ans = 1 ans = 1 ans = 1 ans = 1 ans = 1 ans = 1 ans = 1 极限为1 练习9 对上面的例子（A=[2.1,3.4,-1.2,2.3;0.8,-0.3,4.1,2.8;2.3,7.9,-1.5,1.4;3.5,7.2,1.7,-9.0]，x0=[1;2;3;4]），继续计算.观察及的极限是否存在. 输入：temp.m clear all; clc; % x1 = rand(1); % x2 = rand(1); syms a p d x; x0 = [1,2,3,4]; a = sym([2.1,3.4,-1.2,2.3;0.8,-0.3,4.1,2.8;2.3,7.9,-1.5,1.4;3.5,7.2,1.7,-9.0]); [p, d] = eig(a); for i = 1:10 b = eval(limit(x0*p*d.^x*inv(p),10*i)) end 输出： b = 1.0e+010 * -2.1642 -5.8793 0.6730 9.6298 b = 1.0e+021 * -0.7141 -1.9127 0.1887 3.0332 b = 1.0e+031 * -2.2616 -6.0565 0.5963 9.6004 b = 1.0e+042 * -0.7159 -1.9171 0.1888 3.0388 b = 1.0e+052 * -2.2660 -6.0681 0.5975 9.6189 b = 1.0e+063 * -0.7173 -1.9207 0.1891 3.0447 b = 1.0e+073 * -2.2703 -6.0798 0.5986 9.6374 b = 1.0e+084 * -0.7186 -1.9244 0.1895 3.0505 b = 1.0e+094 * -2.2747 -6.0915 0.5998 9.6559 b = 1.0e+105 * -0.7200 -1.9281 0.1898 3.0564 极限不存在 练习10 求出的所有特征值与特征向量，并与练习9的结论作对比. 输入： [P,D]=eig(A) 输出： P = -0.3779 -0.8848 -0.0832 -0.3908 -0.5367 0.3575 -0.2786 0.4777 -0.6473 0.2988 0.1092 -0.7442 -0.3874 -0.0015 0.9505 0.2555 D = 7.2300 0 0 0 0 1.1352 0 0 0 0 -11.2213 0 0