司老师第九讲QQ讨论部分.doc

海航司守奎 19:08:24 clc,?clearfr=@(r)50000*r*(1+r)^36-1637*(1+r)^36+1637;?%定义代数方程左边项的匿名函数r=fsolve(fr,0.5)??%求月利率的第一种方法r?=fzero(fr,[0.0001,0.1])??%求月利率的第二方法 海航司守奎 19:09:01 海航司守奎 19:11:51 a=zeros(1,181);a(1)=1;?a(61)=-12.3;?a(181)=11.3;b=roots(a) 201307029(785144557) 19:12:52 fx=@(x)x^180-12.3*x^120+11.3 海航司守奎 19:13:13 fx=@(x)x^180-12.3*x^120+11.3;x=fzero(fx,[1.0001,1.1]) 海航司守奎 19:15:36 clc,?clearformat?long?g?fr=@(r)80000*r*(1+r)^180-660.88*(1+r)^180+660.88;?%定义代数方程左边项的匿名函数f=fzero(fr,[0.0001,0.1])??%求月利率format?%恢复到短小数的显示 海航司守奎 19:20:09 clc,?clearformat?long?g?fr=@(r)80000*r*(1+r)^180-660.88*(1+r)^180+660.88;?%定义代数方程左边项的匿名函数r=fzero(fr,[0.0001,0.1])??%求月利率N=181;?A0=80000;b=N*A0*r/2??%计算还款总利息delta=A0*r/N??%计算每月的还款差额xk=A0/N+(N+1)*A0*r/N/2??%计算每月还款额format?%恢复到短小数的显示 海航司守奎 19:33:11 clc,?clearp=[0.08?0.9?0.02;0.08?0.9?0.02;0.24?0.7?0.06];a=[p-eye(3);?ones(1,3)];?%构造非线性方程组的系数阵b=[zeros(3,1);1];?%构造非线性方程组的常数项列plimit=a\b 海航司守奎 19:37:02 clc,?clearp=[0.08?0.9?0.02;0.08?0.9?0.02;0.24?0.7?0.06];a=[p-eye(3);?ones(1,3)];?%构造非线性方程组的系数阵b=[zeros(3,1);1];?%构造非线性方程组的常数项列plimit=a\b[vec,val]=eigs(p,1);vec=vec/sum(vec);[plimit,vec] 海航司守奎 19:38:44 [1]许哲，李号雷，基于灰色模型和Bootstrap理论的大规模定制质量控制方法研究，数学的实践与认识，第42卷第21期，121~127，2012年11月。 海航司守奎 19:47:27 cumsum(x) 海航司守奎 19:47:33 就是累加运算 海航司守奎 20:00:08 clc,?clearx0=[2.5320,2.6470,2.6290,2.5840,2.6090,2.6010,2.5280,2.5630,2.6540,2.6190];n=length(x0);me=quantile(x0,0.5)??%计算中位数[h,p,stat]=runstest(x0,me)??%进行游程检验x1=cumsum(x0);??%求累加序列zk=(x1(1:end-1)+x1(2:end))/2??%求累加序列的均值序列B=[-zk,?ones(size(zk))];?yn=x0(2:end);ab=B\yn??%拟合参数a,bx=dsolve(Dx+a*x=b,x(0)=x00)?%求微分方程的符号解x=subs(x,{a,b,x00},{ab(1),ab(2),x0(1)});xx=vpa(x,6)yuce=subs(x,t,[0:n+5]);??%求累加序列的预测值yuce0=[x0(1),diff(yuce)]??%求原始数据的预测值c=std(yuce0(1:n))/std(x0)??%求后验差比值cnyuce=yuce0(n+1:end)??%提取6个新的预测值nyb=[x0,?nyuce];??%构造新的样本数据nnyb=reshape(nyb,[4,4])mu