南京工业大学研究生工程数学应用B第一次作业..docxVIP

工程应用数学作业（一）专业：安全科学与工程 姓名：** 学号：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