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

目录

1.非线性方程求根的数值解法

房贷年利率………1

2.线性方程组的数值解法

配置指定成分合金………………3

3.估计与检验

铝合金板的批次检验……………8

4插值与拟合算法

晶粒尺寸与退火时间的关系…10

5.回归分析

镁合金析氢与时间关系………14

6.数值积分法

人口统计模型…………………17

1.非线性方程求根的数值解法

……………房贷年利率

（1）问题描述

某人欲购置一套总价为600000元的房产，首付30%，余

下的70%按揭并共20年还清。已知每月还款额为3000元，

求银行的年利率。

（2）建立模型

k

设x为第k年的欠款数，a为每年还款数，r为年利率。

可得一下关系式：

k+1k

x=(1+r)x-a

由该递推关系可得：

kk

k0

x=(1+r)x-a[(1+r)-1]

020

其中a=36000,x=420000,x=0

（3）算法的选择

这是一个典型的非线性方程的求解问题，在此使用

Newton迭代法求近似解。即运用迭代公式：

f(x)

k

xx

k1kf(x)

k

（4）计算过程及MATLAB程序

%HouseLoanNewton

1

x0=0.2

x1=x0-(420*x0*(1+x0)^20-60*((1+x0)^20-1))/(420*(1+x

0)^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+x

0)^20+8400*x0*(1+x0)^19-1200*(1+x0)^19);

n=n+1;

end

x1

N

MATLAB程序截图

（5）计算结果与分析

由计算结果可以看出，Newton迭代收敛速度很快，仅通

过7次计算便得到了所需进度的近似解。

2

2．线性方程组的数值解法

……………配置指定成分合金

（1）问题描述

先要求配置某种合金，其成分为Al-Hg-Ge-Pb，要求利用已有的

Al-Hg-Ge-Pb合金、Al-Hg-Ge合金、Al-Pb合金和纯Al配置。其原

料成分和目标成分要求如下表：

元素中间合金（wt%）目标合金（wt%)

Al-Hg-Ge-PbAl-Hg-GeAl-Pb

Al40805067

Hg3010\15

Ge1510\8

Pb15