购房贷款的利率.pptVIP

* end 一、引例 ：购房贷款的利率 不难算出，你向银行总共借了25.2万，30年内共要还51.696万，约为当初借款的两倍。这个案例中贷款年利率是多少呢？ 下面是《新民晚报》2000年3月30日第七版上的一则房产广告： 解 设xk为第k个月的欠款数， a为月还款数, r为月利率。 xk+1 = (1+r) xk- a 那么 xk = (1+r) xk-1- a = (1+r)2 xk-2 – (1+r)a – a =…… = (1+r)k x0 – a[1+(1+r)+……+(1+r)k-1] = (1+r)k x0 – a[(1+r)k-1]/r 根据 a=0.1436, x0=25.2, x360=0 得到 25.2(1+r)360 – 0.1436[(1+r)360-1]/r=0 二、数学理论复习：非线性方程（组） 一元非线性方程的一般形式为 f (x) = 0 若对于数?有f (?) = 0, 则称?为方程的解或根，也称为函数f (x)的零点 若对于数?有f (?) = 0, f ’(?)?0则?称为单根若有k 1, f (?) = f ’(?) = …= f (k-1)(?) = 0 但f (k)(?)?0 , 称为k重根 非线性方程（组）求解通常用数值方法求近似解，常见的有二分法、牛顿法等 1 图解法 适用于求一元或二元方程(组)低精度解或找迭代初值 例1 解方程 sin(x) = 0.1x 解 -1? sinx ?1, 得| x |10, 作出sin(x), 0.1x 在[-10,10]范围内的图, 可看出根的大致位置 ? fplot(sin,[-10,10]);hold on; ? fplot(0.1*x,[-10,10]) ；grid; zoom 三、数值解法：图解法和迭代法 2 迭代法 迭代法是从解的初始近似值x0(简称初值)开始，利用某种迭代格式x k+1 = g (x k )，求得一近似值序列x1, x2, …, xk, xk+1, …逐步逼近于所求的解?(称为不动点)。最常用的迭代法是牛顿迭代法，其迭代格式为 例2 求方程 x 2 - 3 x + e x = 2 的正根 (要求精度? = 10 -6) 解 令f (x) = x 2 - 3 x + e x - 2, f(0)=-1,当x 2, f (x) 0, f ’(x) 0即f (x)单调上升，所以根在[0,2]内。 先用图解法找初值， ?fplot(x^2-3*x+exp(x)-2,[0,2]) ,grid 取x0 = 1, 迭代格式 四、使用MATLAB 1 多项式的根 roots(p)可求得多项式p的所有复根 MATLAB中一个多项式用系数降幂排列向量来表示。 例如 求多项式x3 + 2 x2 - 5的根 用MATLAB命令 ? roots([1 2 0 -5]) 2 一元函数零点 fzero (F,X,tol),返回函数F的一个零点,其中F为字符串表示的函数或M函数名。X为标量时,为迭代初值； X为向量[a,b]时, 返回函数F在[a,b]中的一个零点,要求F在a,b异号。tol为可变精度(缺省值1e-4)。 例1可求解如下 ?fzero(sin(x)-0.1*x,6) ?fzero(sin(x)-0.1*x,[2,6]) 3 非线性方程组求解 fsolve可求非线性方程或多元非线性方程组的实根， 用法与fzero类似。 例3 解方程组 解法一：本题可通过作二元函数图求解 解法二：先写一个M函数eg2_2fun.m， 然后用 ?x=fsolve(‘eg2_2fun’,[0,0]) 求解 解法三：直接求解 ?x=fsolve(‘[4*x(1)-x(2)+exp(x(1))/10-1,… -x(1) +4*x(2)+x(1).^2/8]’,[0,0]) 4 解析求解solve solve是符号数学工具箱中的解方程组命令。可求各种类型方程组的解析解。当找不到解析解时，solve会自动寻求一个近似解，并且精度很高 ?solve(‘a*x^2+b*x+c’, ‘x’) 可得二次多项式求根公式 ?solve(‘x^2-3*x+exp(x)-2’,‘x’) 可解例2 ?[x,y]=solve(‘4*x-y+exp(x)/10=1’,… ‘-x+4*y+x^2/8=0’, ‘x,y’) 可解例3 五、实验例题 例4 (引例) 常识上，r应比当时活期存款月利率略高一些。我们用活期存款月利率0.0198/12 作为迭代初值，用fzero求解 ?r=fzero(25.2*(1+x)^360-((1+x)^360-1)/... x*0.1436,0.