matlab高等工程数学作业实践报告.doc

目 录 非线性方程求根的数值解法 房贷年利率………………………………………1 2.线性方程组的数值解法 配置指定成分合金………………………………3 3.估计与检验 铝合金板的批次检验……………………………8 4插值与拟合算法 晶粒尺寸与退火时间的关系…………………10 回归分析 镁合金析氢与时间关系………………………14 数值积分法 人口统计模型…………………………………17  1.非线性方程求根的数值解法 ……………房贷年利率 问题描述 某人欲购置一套总价为600000元的房产，首付30%，余下的70%按揭并共20年还清。已知每月还款额为3000元，求银行的年利率。 建立模型 设xk为第k年的欠款数，a为每年还款数，r为年利率。可得一下关系式： xk+1=(1+r)xk-a 由该递推关系可得： xk=(1+r)kx0-a[(1+r)k-1] 其中a=36000,x0=420000,x20=0 算法的选择 这是一个典型的非线性方程的求解问题，在此使用Newton迭代法求近似解。即运用迭代公式： 计算过程及MATLAB程序 %HouseLoanNewton x0=0.2 x1=x0-(420*x0*(1+x0)^20-60*((1+x0)^20-1))/(420*(1+x0)^20+8400*x0*(1+x0)^19-1200*(1+x0)^19); n=1; while (abs(x1-x0)=1.0e-6)(n=100) x0=x1; x1=x0-(420*x0*(1+x0)^20-60*((1+x0)^20-1))/(420*(1+x0)^20+8400*x0*(1+x0)^19-1200*(1+x0)^19); n=n+1; end x1 N MATLAB程序截图 计算结果与分析 由计算结果可以看出，Newton迭代收敛速度很快，仅通过7次计算便得到了所需进度的近似解。  2．线性方程组的数值解法 ……………配置指定成分合金 问题描述 先要求配置某种合金，其成分为Al-Hg-Ge-Pb，要求利用已有的Al-Hg-Ge-Pb合金、Al-Hg-Ge合金、Al-Pb合金和纯Al配置。其原料成分和目标成分要求如下表： 元素 中间合金（wt%） 目标合金（wt%) Al-Hg-Ge-Pb Al-Hg-Ge Al-Pb Al 40 80 50 67 Hg 30 10 \ 15 Ge 15 10 \ 8 Pb 15 \ 50 10 若现在欲配置100kg目标成分的Al-Hg-Ge-Pb合金，需要准备多少上述中间合金。 (2)建立模型 为了搭配合理，我们假设3种不同的中间合金的质量分别为x1、x2、x3，纯Al的质量为x4，单位为kg。得到如下方程组： 0.4x1+0.8x2+0.5x3+x4=67 0.3x1+0.1x2=15 0.15x1+0.1x2=8 0.15x1+0.5x3=10 将上面的方程组写成矩阵的形式为：aX=b 其中， ， ， 模型假设：以上数据真实有效 由此该问题的模型可以建立如下： 可见解上述问题就需要对上面这个线性方程组进行求解。 算法选择 在此选择的算法为Jacobi迭代算法，同时与矩阵直接除法即精确解做对比。Jacobi迭代公式为： 计算流程与MATLAB程序 %Jacobi迭代法 x0=[0,0,0,0]; x1=[1,1,1,1]; i=1; while norm(x1-x0)=0.001i100 x0=x1; x1(1)=1/0.3*(-0.1*x0(2)+15); x1(2)=1/0.1*(-0.15*x0(1)+8); x1(3)=1/0.5*(-0.15*x0(1)+10); x1(4)=1/1*(-0.4*x0(1)-0.8*x0(2)-0.5*x0(3)+67); i=i+1; end x1 i %直接除法 a=[0.4 0.8 0.5 1 0.3 0.1 0 0 0.15 0.1 0 0 0.15 0 0.5 0]; b=[67 15 8 10]’; x=a\b MATLAB程序截图： 计算结果与分析 由结果可见，经过37次迭代，Jacobi迭代法产生了足够精确的结果（误差小于0.001）。与精确解对比，已经足以完成原先的合金配置要求。同时为了确保Jacobi迭代矩阵的收敛，在编程的时候将方程组中的各方程顺序做了适当的调整。3.估计与检