北京科技大學(xué)數(shù)學(xué)實(shí)驗(yàn)報告-第4次

PAGEPAGE2《數(shù)學(xué)實(shí)驗(yàn)》報告實(shí)驗(yàn)名稱

MATLAB基礎(chǔ)知識

學(xué)院專業(yè)班級

姓

名

學(xué)

號

2013年5月【實(shí)驗(yàn)?zāi)康摹空莆誐ATLAB線性代數(shù)相關(guān)運(yùn)算；掌握矩陣的基本運(yùn)算、特征值、特征向量和線性方程組的求解；能熟練運(yùn)用數(shù)值方法求定積分3、運(yùn)用MATLAB求積分?！緦?shí)驗(yàn)任務(wù)】完成P114第12題、第14題、第21（1）題；P167第17（2）題、第18題【實(shí)驗(yàn)程序】P114第12題：a=magic(6)a’det(a)rank(a)rref(a)P114第14題：A=[2,1,1;1,2,1;1,1,2][V,D]=eig(A)p=poly(A)poly2str(p,'X')第21（1）題：A=[1121;-1130;2-34-1];rref(A)P167第17（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)P167第18題：x=0:pi/50:pi/4;y=1./(1-sin(x));s1=sum(y(1:12.5))*(pi/50)s2=sum(y(2:13.5))*(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)【實(shí)驗(yàn)結(jié)果】1.a=1626192433272123253192222720828331710153053412141643629131811轉(zhuǎn)置a'=353318304132928536672333429262122171213192327101418242520151611秩：5最簡形:10000-201000-20010010001020000120000002.特征值;1.00001.0004.0000特征多項(xiàng)式;b=X^3-6X^2+9X-4特征向量：0.40820.70710.57740.4082-0.70710.5774-0.816500.57743.ans=1000.480100.6001-0.04則原方程等價為：方程的個數(shù)<未知量個數(shù)有無窮多個解取=1，則=-0.48，=-0.6，=0.04所以，方程的通解為：其中k為任意常數(shù)。4.s1=2.4669s2=2.4669s3=2.4669s4=2.46745.s1=1.2448s2=1.3811s3=1.3129s4=1.4142u1=-0.16945u2=-0.033098u3=-0.10127u4=3.