(10)--实验六 方程及方程组的解.pptVIP

*编写函数文件f.mfunctiony=f(r)y=80*r*(1+r)^36-3*(1+r)^36+3;确定根的隔离区间fplot(f,[0.01,0.5]),gridon↙由图知根的隔离区间［0.01,0.5］.*调用二分法的程序,在命令窗口输入:erfen(f,0.01,0.51.0e-4)↙于是得月利率年利率为所以分期付款买车要付的利率高达20.64%，与3年贷款的利率比较，最好选择贷款买车.n=10ans=0.0172*例11某人向银行贷款30000元购房,这笔贷款在15年内分180次按月等额偿还，每月还款450元.房地产开发商向他提出更“诱人”的的承诺：在每月还款额基本不变的情况下,提前两年还清贷款.他们的条件是:(1)首期先预付半年的还款;(2)改每月还款450元为每半月还款225元.在这两种情况下,分别求这笔贷款应付的年利率.解（1）按月等额偿还.设借款的月利率为r，由公式*整理得于是有记编写函数文件f.mfunctiony=f(r)y=1000*r*(1+r)^180-15*(1+r)^180+15;*计算函数f(r)的值:feval(f,0.05)↙feval(f,0.01)↙有根.ans=2.2812e+05ans=-14.9790*编写导函数文件df.mfunctiony=df(r)y=1000*(1+r)^180+180000*r*(1+r)^179-15*180*(1+r)^179;在牛顿迭代法程序中,取初值在命令窗口输入：newton（f,df,0.01）↙r1=0.0137n=9*(2)首期先预付半年的还款,其余按半月等额偿还.设借款的半月利率为r，由公式(1)整理得*根据经验,方程(2)的最小正根在0附近(但不等于0),分别取*编写导函数文件f.mfunctiony=f(r)y=10000*r*(1+12*r)*(1+r)^300-900*r*(1+r)^300-75*(1+12*r)*(1+r)^300+75*(1+12*r);计算f(r)在r0=0.01和0.001的值:feval(f,0.01)↙feval(f,0.001)↙ans=459.9890ans=-14.0951*内有根.编写导函数文件df.mfunctiony=df(r)y=10000*(1+12*r)*(1+r)^300+120000*r*(1+r)^300+3000000*r*(1+12*r)*(1+r)^299-1800*(1+r)^300-270000*r*(1+r)^299-22500*(1+12*r)*(1+r)^299+900;*由此可见，按房地产开发商提供的方式买房比到银行贷款买房年利率要高出0.84%.在命令窗口输入：newton（f,df,0.01）↙r1=0.0072n=4*离区间,重复上述做法，*（2）MATLAB实现用二分法求方程根的近似值，一般分为两步：第一步根的隔离用fplot函数作出y=f(x)的图形，它与x轴的交点就是方程f(x)=0的根.确定根的存在区间[a,b](即根的隔离区间)使方程第二步根的逐次逼近若函数f(x)在根的隔离区间[a,b]上单调连续,且则在(a,b)内方程有唯一根.我们可以反复用将根的隔离区间一分为二，其中含有根的区间为新的隔离区间的方法，最终搜索到达到要求精度的根的近似值.*functiony=erfen(fun,a,b,esp)iffeval(fun,a).*feval(fun,b)0n=1;c=(a+b)/2;whileabs(b-a)espiffeval(fun,a).*feval(fun,c)0b=c;c=(a+b)/2;elseiffeval(fun,c).*feval(fun,b)0a=c;c=(a+b)/2;elsey=c;esp=10000;end打开M-edit窗口,编写二分法程序：erfen·m*n=n+1;endy=c;elseiffeval(fun,a)==0y=a;elseiffeval(fun,b)==0y=b;elsedisp(these,maynotbearoo