数学建模课程实验报告.doc

《数学建模课程实验报告.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、《数学建模》课程实验报告学院班级姓名学号上课时间12/12一、输入下面的矩阵A、B并完成相应的运算；1、分别求出矩阵A和矩阵B的逆矩阵、矩阵的秩、矩阵所对应的行列式的值及A^2;B^2解：>>inv(A)ErrorusinginvMatrixmustbesquare.>>det(A)ErrorusingdetMatrixmustbesquare.>>rank(A)ans=3>>A.^2ans=191699251619940916411111>>inv(B)ErrorusinginvMatrixmustbesquare.>>det(B

2、)ErrorusingdetMatrixmustbesquare.>>rank(B)ans412/12>>B.^2ans=9191404164649094916911112.分别求出矩阵A和矩阵B的伴随矩阵、矩阵的特征值及特征向量、并将矩阵化为最简的阶梯型矩阵；>>compan(A)Errorusingcompan(line13)Inputargumentmustbeavector.>>eig(A)ErrorusingeigMatrixmustbesquare.>>[v,d]=eig(A)ErrorusingeigMatrixmus

3、tbesquare.>>rref(A)ans=1.0000001.400001.000000.4000001.0000-0.800000000000>>compan(B)12/12Errorusingcompan(line13)Inputargumentmustbeavector.>>eig(B)ErrorusingeigMatrixmustbesquare.>>[v,d]=eig(B)ErrorusingeigMatrixmustbesquare.>>rref(B)ans=100001000010000100001、完成下列矩阵的

4、运算A*B、A/B、AB、A.*B、A./B>>A*BErrorusing*Innermatrixdimensionsmustagree.>>A/Bans=2.1622-1.594600.5135-2.75682.0000-1.000001.0000-4.00000.9595-0.351400.6216-2.31080.8784-0.054100.8649-2.93240.0000-0.00000-0.00001.000012/12>>ABWarning:Rankdeficient,rank=3,tol=3.845925e-15

5、.ans=0000-0.30681.53401.05850.5082-0.10770.53860.99770.81971.4145-1.0726-1.0562-0.3279>>A.*Bans=-3312-3-608-4424-60928-8-31111>>A./Bans=-0.33333.00001.3333-3.0000-1.5000Inf2.0000-0.25001.00000.3750-0.6667NaN1.00000.5714-0.5000-0.33331.00001.00001.00001.0000；二、解下面的线性方程组

6、；12/12（1）（2）解：>>C=[21-51;1-30-6;02-12;14-76]C=21-511-30-602-1214-76>>d=[89-50]'d=89-50>>E=[123;225;351]E=12/12123225351>>f=[123]'f=123>>x=Cdx=3.0000-4.0000-1.00001.0000>>x=Efx=1.00000.00000.0000三、解决下列高等数学中的问题（一）求出下列极限的值；12/12解：>>symsabx;f=((a*x+4)/(a*x+1))^b*x;limit(

7、f,x,inf)ans=Inf-(imag(b)*(real(a)*((2*real(a)^2)/(imag(a)^2+real(a)^2)^2-1/(imag(a)^2+real(a)^2))-(4*real(a))/(imag(a)^2+real(a)^2)+(2*imag(a)^2*real(a))/(imag(a)^2+real(a)^2)^2)*i)/(imag(a)^2/(imag(a)^2+real(a)^2)+real(a)^2/(imag(a)^2+real(a)^2))-(real(b)*((real(a)*((

8、2*real(a)^2)/(imag(a)^2+real(a)^2)^2-1/(imag(a)^2+real(a)^2))+imag(a)*((2*real(a)^2)/(imag(a)^2+real(a)^2)^2-1/(imag(