实验报告样例2

PAGEPAGE6数值实验报告Ⅱ实验名称三对角方程求解问题实验时间2012姓名吴海洋班级数学10-1学号10072116成绩一、实验目的，内容1．理解三对角方程的概念2．掌握三对角方程的解法，主要是追赶法二、算法描述追赶法的步骤及思想：三对角矩阵的追赶法分解实质上也是基于对矩阵的LU分解，思想上大致是一致的，但三对角也有其很好的性质，其具体分解如下：=根据公式即可得到三对角矩阵LU分解的计算公式如下：方程组求解的计算公式为：以上便是三对角矩阵的追赶法分解以及求解过程。三．程序代码三对角矩阵的追赶法：n=100;%n=500;%方程组的真解x0=ones(1,n);a1=-1*ones(1,n-1);a3=a1;a2=2*ones(1,n);b=zeros(1,n);b(1)=1;b(n)=1;%a1为对角下向量%a2为对角向量%a3为对角上向量%b为方程常数项r1=size(a1);m=r1(2);r2=size(a2);n=r2(2);ifsize(a1)~=size(a2)|m~=n-1|size(a2)~=size(b)error('变量不匹配，检查变量输入情况！');endp=ones(1,m);Y=ones(1,n);x=Y;p(1)=a1(1)/a2(1);Y(1)=b(1)/a2(1);t=0;fori=2:mt=a2(i)-a1(i-1)*p(i-1);p(i)=a3(i)/t;Y(i)=(b(i)-a1(i-1)*Y(i-1))/t;endY(n)=(b(n)-a1(n-1)*Y(n-1))/(a2(n)-a1(n-1)*p(n-1));x(n)=Y(n);fori=n-1:-1:1x(i)=Y(i)-p(i)*x(i+1);end%方程组的解为disp('方程组的解为')x%数值解与真值的误差为disp('数值解与真值的误差为:')c=norm(x-x0,2)四.数值结果方程组的解为：（1）n=100时x=Columns1through71.00001.00001.00001.00001.00001.00001.0000Columns8through141.00001.00001.00001.00001.00001.00001.0000Columns15through211.00001.00001.00001.00001.00001.00001.0000Columns22through281.00001.00001.00001.00001.00001.00001.0000Columns29through351.00001.00001.00001.00001.00001.00001.0000Columns36through421.00001.00001.00001.00001.00001.00001.0000Columns43through491.00001.00001.00001.00001.00001.00001.0000Columns50through561.00001.00001.00001.00001.00001.00001.0000Columns57through631.00001.00001.00001.00001.00001.00001.0000Columns64through701.00001.00001.00001.00001.00001.00001.0000Columns71through771.00001.00001.00001.00001.00001.00001.0000Columns78through841.00001.00001.00001.00001.00001.00001.0000Columns85through911.00001.00001.00001.00001.00001.00001.0000Columns92through981.00001.00001.00001.00001.00001.00001.0000Columns99through1001.00001.0000数值解与真值的误差为:c=7.3995e-014（2）n=500时方程组的解为x=Columns1through71.00001.00001.00001.00001.00001.00001.0000Columns8through141.00001.00001.00001.00001.00001.00001.0000Columns15through211.00001.00001.00001.00001.00001.00001.0000Columns22through281.00001.00001.00001.00001.00001.00001.0000Columns29through351.00001.00001.00001.00001.00001.00001.0000Columns36through421.00001.00001.00001.00001.00001.00001.0000Columns43through491.00001.00001.00001.00001.00001.00001.0000Columns50through561.00001.00001.00001.00001.00001.00001.0000Columns57through631.00001.00001.00001.00001.00001.00001.0000Columns64through701.00001.00001.00001.00001.00001.00001.0000Columns71through771.00001.00001.00001.00001.00001.00001.0000Columns78through841.00001.00001.00001.00001.00001.00001.0000Columns85through911.00001.00001.00001.00001.00001.00001.0000Columns92through981.00001.00001.00001.00001.00001.00001.0000Columns99through1051.00001.00001.00001.00001.00001.00001.0000Columns106through1121.00001.00001.00001.00001.00001.00001.0000Columns113through1191.00001.00001.00001.00001.00001.00001.0000Columns120through1261.00001.00001.00001.00001.00001.00001.0000Columns127through1331.00001.00001.00001.00001.00001.00001.0000Columns134through1401.00001.00001.00001.00001.00001.00001.0000Columns141through1471.00001.00001.00001.00001.00001.00001.0000Columns148through1541.00001.00001.00001.00001.00001.00001.0000Columns155through1611.00001.00001.00001.00001.00001.00001.0000Columns162through1681.00001.00001.00001.00001.00001.00001.0000Columns169through1751.00001.00001.00001.00001.00001.00001.0000Columns176through1821.00001.00001.00001.00001.00001.00001.0000Columns183through1891.00001.00001.00001.00001.00001.00001.0000Columns190through1961.00001.00001.00001.00001.00001.00001.0000Columns197through2031.00001.00001.00001.00001.00001.00001.0000Columns204through2101.00001.00001.00001.00001.00001.00001.0000Columns211through2171.00001.00001.00001.00001.00001.00001.0000Columns218through2241.00001.00001.00001.00001.00001.00001.0000Columns225through2311.00001.00001.00001.00001.00001.00001.0000Columns232through2381.00001.00001.00001.00001.00001.00001.0000Columns239through2451.00001.00001.00001.00001.00001.00001.0000Columns246through2521.00001.00001.00001.00001.00001.00001.0000Columns253through2591.00001.00001.00001.00001.00001.00001.0000Columns260through2661.00001.00001.00001.00001.00001.00001.0000Columns267through2731.00001.00001.00001.00001.00001.00001.0000Columns274through2801.00001.00001.00001.00001.00001.00001.0000Columns281through2871.00001.00001.00001.00001.00001.00001.0000Columns288through2941.00001.00001.00001.00001.00001.00001.0000Columns295through3011.00001.00001.00001.00001.00001.00001.0000Columns302through3081.00001.00001.00001.00001.00001.00001.0000Columns309through3151.00001.00001.00001.00001.00001.00001.0000Columns316through3221.00001.00001.00001.00001.00001.00001.0000Columns323through3291.00001.00001.00001.00001.00001.00001.0000Columns330through3361.00001.00001.00001.00001.00001.00001.0000Columns337through3431.00001.00001.00001.00001.00001.00001.0000Columns344through3501.00001.00001.00001.00001.00001.00001.0000Columns351through3571.00001.00001.00001.00001.00001.00001.0000Columns358through3641.00001.00001.00001.00001.00001.00001.0000Columns365through3711.00001.00001.00001.00001.00001.00001.0000Columns372through3781.00001.00001.00001.00001.00001.00001.0000Columns379through3851.00001.00001.00001.00001.00001.00001.0000Columns386through3921.00001.00001.00001.00001.00001.00001.0000Columns393through3991.00001.00001.00001.00001.00001.00001.0000Columns400through4061.00001.00001.00001.00001.00001.00001.0000Columns407through4131.00001.00001.00001.00001.00001.00001.0000Columns414through4201.00001.00001.00001.00001.00001.00001.0000Columns421through4271.00001.00001.00001.00001.00001.00001.0000Columns428through4341.00001.00001.00001.00001.00001.00001.0000Columns435through4411.00001.00001.00001.00001.00001.00001.0000Columns442through4481.00001.00001.00001.00001.00001.00001.0000Columns449through4551.00001.00001.00001.00001.00001.00001.0000Columns456through4621.00001.00001.00001.00001.00001.00001.0000Columns4