北京科技大学数学实验报告---第4次.doc

文档介绍：《数学实验》报告实验名称MATLAB基础知识学院专业班级姓名学号2013年5月2一、【实验目的】1、掌握MATLAB线性代数相关运算;2、掌握矩阵的基本运算、特征值、特征向量和线性方程组的求解;能熟练运用数值方法求定积分3、运用MATLAB求积分。二、【实验任务】完成P114第12题、第14题、第21(1)题;P167第17(2)题、第18题三、【实验程序】:a=magic(6)a’det(a)rank(a)rref(a):A=[2,1,1;1,2,1;1,1,2][V,D]=eig(A)p=poly(A)poly2str(p,'X')(1)题:A=[1121;-1130;2-34-1];rref(A)(2)题x=0:pi/50:pi;y=(x.*(sin(x)))./(1+(cos(x)).^2);s1=sum(y(1:50))*(pi/50)s2=sum(y(2:51))*(pi/50)s3=trapz(x,y)s4=quad('(x.*(sin(x)))./(1+(cos(x)).^2)',0,pi):x=0:pi/50:pi/4;3y=1./(1-sin(x));s1=sum(y(1:))*(pi/50)s2=sum(y(2:))*(pi/50)s3=trapz(x,y)s4=quad('1./(1-sin(x))',0,pi/4)u1=s1-sqrt(2)u2=s2-sqrt(2)u3=s3-sqrt(2)u4=s4-sqrt(2)四、【实验结果】=1626192433272123253192222720828331710153053412141643629131811转置a'=353318304132928536672333429262122171213192327101418242520151611秩:5最简形:10000-201000-;;b=X^3-6X^2+9X-4特征向量:--