第五讲非线性方程模型实验购房贷款的利率教学幻灯片.pptVIP

第五讲 非线性方程模型实验  购房贷款的利率 实验目的 实验内容 2、学会用Matlab代数方程的数值解. 1、学会用Matlab求代数方程的解析解. 1、求代数方程的解析解. 4、实验作业. 2、求代数方程的数值解. 问题：如下是一则房产广告。不难算出，你向银行共借了25.2万，30年内共要还51.696万，约为当初借款的两倍，这个案例中贷款年利率是多少？ 建筑面积 总价 30%首付 70%按揭 月还款 86.98m2 36万 10.8万 30年 1436元 分析 有人可能会这样算 年利率=(51.696-25.2)/30/25.2=3.5% 错的，因为你并不是等到30年后一次性还款。 根据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命令 roots 求多项式的根 fsolve 方程（组）数值解 fzero 求一元函数实根 solve 符号方程（组）求解 1. 多项式的根 roots(p) 多项式p的所有复根。例 x3+2x2-5的根 roots([1 2 0 -5]) ans = -1.6209 + 1.1826i -1.6209 - 1.1826i 1.2419 2. 一元函数零点 fzero(F,X,tol) F为字符串表示的函数或M函数名； x为标量时，作为迭代初值；X为向量[a,b]时，返回F在[a,b]中的一个零点，这时要求F在a,b两点异号；tol为精度（缺损值1e-4). 例: y=sin(x)-0.1x fzero(sin(x)-0.1*x,6) ans = 7.0682 fzero(sin(x)-0.1*x,[2,6]) ans = 2.8523 注：fzero 只能求零点附近变号的根，试用fzero求解(x-1)2=0, 看看发生了什么？ 3. 非线性方程组求解 fsolve 用法与fzero类似,例：解方程组 写M函数eg2_1fun.m function y=fun(x) y(1)=4*x(1)-x(2)+exp(x(1))/10-1; y(2)=-x(1)+4*x(2)+x(1)^2/8; 然后用 [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))/1