南京工业大学生工程数学应用B第一次作业教程范本.docx

工程应用数学作业（一）专业：安全科学与工程 姓名：** 学号：652083700***一、第三章练习题说明：解题所需函数均提前保存为m文件，方便直接调用，以下解题过程不包括保存m文件过程。练习题1：设两点边值问题的精确解为现以h为步长划分区间为100等份，用差分近似代替微分，将微分方程离散化为线性方程组，代入初始条件后，得到如下的方程组问题其中，，。(1) 分别用J迭代法，G-S迭代法和SOR迭代法求解，并与精确解进行比较；(2) 如果，，再求解该问题解：（1）1）输入：t=1;a=1/2;h=1/100;A=(-h-2*t)*eye(99)+diag((t+h)*ones(98,1),1)+diag((t)*ones(98,1),-1);C(99,1)=zeros;C(99,1)=t+h;b=(a*h^2)*ones(99,1)-C;x0=zeros(99,1);e=2e-10;y=jacobi(A,b,x0,e)进行jacobi迭代。2）输入：t=1;a=1/2;h=1/100;A=(-h-2*t)*eye(99)+diag((t+h)*ones(98,1),1)+diag((t)*ones(98,1),-1);C(99,1)=zeros;C(99,1)=t+h;b=(a*h^2)*ones(99,1)-C;x0=zeros(99,1);e=2e-10;for i=1:15 w=0.95+i/20;N=8; [y,r,n]=sor(a,b,x0,w,1e-15,N); h(i)=r;end[s,j]=min(h);w=0.95+j/20得松弛因子w=1，再输入：[x,r,n]=sor(A,b,x0,w,e,6900)完成sor算法。3）输入：t=1;a=1/2;h=0.001;A=(-h-2*t)*eye(99)+diag((t+h)*ones(98,1),1)+diag((t)*ones(98,1),-1);C(99,1)=zeros;C(99,1)=t+h;b=(a*h^2)*ones(99,1)-C;x0=zeros(99,1);e=2e-10; [y,r,n]=seidel(A,b,x0,e)完成sediel迭代。输出结果如下表所示：Jacobi迭代Sediel迭代Sor迭代0.01290.01050.01280.0256 0.02090.02560.0383 0.03140.03820.0510 0.04190.05080.0635 0.05230.06340.0760 0.06270.07580.0884 0.07310.08820.1007 0.08350.10050.1130 0.09390.11270.1251 0.10430.12490.1372 0.11470.13700.1493 0.12510.14900.1613 0.13540.16090.1732 0.14580.17280.1850 0.15610.18460.1968 0.16650.19630.2085 0.17680.20800.2201 0.18710.21960.2317 0.19740.23120.2432 0.20770.24270.2546 0.21800.25410.2660 0.22820.26540.2773 0.23850.27670.2885 0.24870.28800.2997 0.25900.29910.3108 0.26920.31020.3219 0.27940.32130.3329 0.28960.33230.3438 0.29980.34320.3547 0.31000.35410.3656 0.32020.36490.3763 0.33040.37570.3870 0.34060.38640.3977 0.35070.39700.4083 0.36090.40760.4188 0.37100.41820.4293 0.38110.42860.4398 0.39120.43910.4501 0.40130.44940.4605 0.41140.45980.4707 0.42150.47000.4810 0.43160.48030.4911 0.44170.49040.5013 0.45170.50060.5113 0.46180.51060.5214 0.47180.52060.5313 0.48180.53060.54