五章节 非线性方程模型实验 购房贷款利率.pptVIP

第五讲 非线性方程模型实验  购房贷款的利率;实验目的;问题：如下是一则房产广告。不难算出，你向银行共借了25.2万，30年内共要还51.696万，约为当初借款的两倍，这个案例中贷款年利率是多少？ ;分析 ; 设xk—第k个月的欠款数；a—月还款数；r—为月利率，我们得到迭代关系式 xk+1=(1+r)xk-a (2.1) 那么 xk=(1+r)xk-1-a=(1+r)2xk-2-(1+r)a-a=… =… =(1+r)kx0-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 (2.2) 关于月利率r的高次代数方程。 年利率R=12r.; 非线性方程（组）简介 若方程是未知量x的多项式，称为高次代数方程；若方程包含x的超越函数，称为超越方程。 一元非线性方程的一般形式为 f(x)=0 (2.3) 若对于数a有f(a)=0,则称 a 为方程（2.3）的解或根，也称为函数f(x)的零点。方程的根可能是实数也可能是复数。相应地称为实根和复根。如果对于数a有f(a)=0,f? (a)?0,则a称为单根，如果有k1,f(a)=f ?(a)=…=f(k-1)(a)=0 但f(k)(a)≠0,称为k重根，对于高次代数方程，其根的个数与其次数相同（包括重数），至于超越方程，其界可能是一个或几个甚至无穷多，也可能无解。 ;常见的求解问题有如下两重要求：一种是要求定出在给定范围内的某个解，而解的粗略位置事先从问题的物理背景或应用（作图等）其他方法得知；另一种是定出方程的全部解，或者给定区域内的所有解，而解的个数未知。除少数特殊的方程可以利用公式直接求解（如4次以下代数方程），一般都没有解析求解方法，只能靠数值方法求得近似解。常见的数值方法有二分法等。 n元非线性方程组的一般形式为 fi(x1,x2,…,xn)=0, i=1,…,m (2.4) 非线性方程组的解极少能用解析法求得。常用的数值方法是Newton法、拟Newton法和最优化方法等。 ;解方程和方程组的MATLAB命令;1. 多项式的根;2. 一元函数零点; fzero(sin(x)-0.1*x,6) ans = 7.0682 fzero(sin(x)-0.1*x,[2,6]) ans = 2.8523;3. 非线性方程组求解;然后用 [x,y,f]=fsolve(eg2_2fun,[0,0]) x = 0.2326 0.0565 y = 1.0e-006 * 0.0908 0.1798 f = 1 注：X返回解向量,y返回误差向量，f0则解收敛。;或直接用 [x,y,f]=fsolve([4*x(1)-x(2)+exp(x(1))/10-1,-x(1)+4*x(2)+x(1).^2/8],[0,0]) x = 0.2326 0.0565 y = 1.0e-006 * 0.0908 0.1798 f = 1 注意：fsolve采用最小二乘优化法，稳定性比fzero好，但fsolve 可能陷入局部极小。试用fsolve解x2+x+1=0,看会发生什么？不要完全相信计算机。;4.解析求解solve; [x,y]=solve(4*x-y+exp(x)/10=1,-x+4*y+y^2/8=0,x,y) x = .23297580773115396971569236570313 y = .58138324907069742242891748561961e-1 注意所得的解与fsolve的不同。 注意：虽然solve可用于求数值解，但速度很慢，且有很大的局限性，不提倡使用。;数值解法：图解法和迭代法;可知±8.5,±7,±3,0 附近各有一解。 ?;2、迭代法（牛顿法，切线法）;例 求如下方程的正根（要求精度?=10-6）x2-3x+ex=2;唯一正根在1附近，取x0=1,迭代格式;贷款利率问题求解;练习;补充：混沌 线性迭代要么收敛于它的不动点，要么趋于无穷大；而不收敛的非线性迭代可能会趋于无穷大，也可能趋于一个周期解，但也可能在一个有限区域内杂乱无章地动弹，由确定性运动导致的貌似