资源描述：
《桂林理工大学数值分析上机实验作业》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。
1、〔值分析上机实验报告学校：桂林理工大学学院:环境科学与工程学院专业：市政工程学号：102015184姓名：邢梦龙P45-46实验题2(1)(3)，5(A为12阶改为5阶),8,10(n=300改为n=10)；第三章迭代法P71-72实验题2,3;习题2,1234第四章数据建模-P106实验题1,2,4;习题43第五章数值微积分P135实验题1(1)-(4)4,5,8;第六章数值分析及其MATLAB实验•53P166实验题1(1),5,6第二章数值代数实验2求下列矩阵的行列式、逆、特征值、特征向量、各种范数、和条件数:(1)"41-1"32-61-53(1)»A=[41-1;32-6;1-53
2、];a=det(A),B=inv(A),[V，D]=eig(A)，t=eig(A)%矩阵的行列式a、逆B、特征值t、特征向量V[norm(A),norm(A,1),norm(A,inf)]%分别为矩阵A的2,1，范数[cond(A),cond(A,1),cond(A,inf)]%分别为矩陈A的2,1，条什数-940.2553-0.02130.04260.1596-0.1383-0.22340.1809-0.2234-0.05320.0185-0.9009-0.3066-0.7693-0.1240-0.7248-0.6386-0.41580.6170D二-3.05273.67608.3766a
3、ns=8.608910.000011.0000ans=3.74945.95745.7340(2)"5765710876810957910»A=[5765;71087;68109;57910];a=det(A),B=inv(A),[V,D]=eig(A),t=eig(A)%矩阵的行列式a、逆B、特征值t、特征向量V[norm(A),norm(A,1),norm(A,inf)]%分别为矩陈A的2,1，范数[cond(A),cond(A,1),cond(A,inf)]%分别为矩陈A的2,1，条件数1.0000B二68.0000-41.0000-17.000010.0000-41.000025.00
4、0010.0000-6.0000-17.000010.00005.0000-3.000010.0000-6.0000-3.00002.00000.83040.09330.39630.3803-0.5016-0.30170.61490.5286-0.20860.7603-0.27160.55200.1237-0.5676-0.62540.52090.010200000.843100003.8581000030.28870.01020.84313.858130.288730.288733.000033.0000ans=1.0e+03
5、2.98414.48804.4880实验5Hilbert矩阵是
6、著名的病态矩阵,n阶Hilbert矩阵定义为A=(ad),其中aij=l/(i+j-l)。设A为5阶Hilbert矩阵，计算cond(A),A-1,norm(A^A-E)及
7、A
8、
9、A"1卜1，并分析结果的精度。再比较MATLAB求解Hilbert矩阵及其逆函数hilb(5),invhilb(5)。»A=hilb(5)1.00000.50000.33330.25000.20000.50000.33330.25000.20000.16670.33330.25000.20000.16670.14290.25000.20000.16670.14290.12500.20000.16670.14290.
10、12500.1111cond(A)ans=4.7661e+005»inv(A)ans=1.0e+005*0.0002-0.00300.0105-0.01400.0063-0.00300.0480•0.18900.2688-0.12600.0105-0.18900.7938-1.17600.5670-0.01400.2688-1.17601.7920-0.88200.0063-0.12600.5670-0.88200.44104.8038e-012»det(A)*det(inv(A))-lans=2.4158e-013»det(hilb(5))ans=3.7493e-012»det(inv(h
11、ilb(5)))ans=2.6672e+011实验8分别用顺序Gauss消去法（程序2.1）和选列主元Gauss消去法（程序2.2）求解下列方程组，验证计算结果，并分析误差产生的原因：0.3xl0"1559.143r"59.17_5.291-6.13-12x246.7811.29521_1211_%4__2_顺序Gauss消去法程序functionx=nagauss(a,b,flag)ifnargin<3,f
