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

1、司老师第九讲QQ讨论部分海航司守奎(94339146) 19:08:24clc,clearfr=(r)50000*r*(1+r)36-1637*(1+r)36+1637;%定义代数方程左边项的匿名函数r=fsolve(fr,0.5)%求月利率的第一种方法r=fzero(fr,0.0001,0.1)%求月利率的第二方法海航司守奎(94339146) 19:09:01海航司守奎(94339146) 19:11:51a=zeros(1,181);a(1)=1;a(61)=-12.3;a(181)=11.3;b=roots(a)201307029(785144557) 19:12:52fx=(x)x1

2、80-12.3*x120+11.3海航司守奎(94339146) 19:13:13fx=(x)x180-12.3*x120+11.3;x=fzero(fx,1.0001,1.1)海航司守奎(94339146) 19:15:36clc,clearformatlonggfr=(r)80000*r*(1+r)180-660.88*(1+r)180+660.88;%定义代数方程左边项的匿名函数f=fzero(fr,0.0001,0.1)%求月利率format%恢复到短小数的显示海航司守奎(94339146) 19:20:09clc,clearformatlonggfr=(r)80000*r*(1+r)

3、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%恢复到短小数的显示海航司守奎(94339146) 19:33:11clc,clearp=0.080.90.02;0.080.90.02;0.240.70.06;a=p-eye(3);ones(1,3);%构造非线性方程组的系数阵b=zeros(3,1);1;%构造非线性方程组的

4、常数项列plimit=ab海航司守奎(94339146) 19:37:02clc,clearp=0.080.90.02;0.080.90.02;0.240.70.06;a=p-eye(3);ones(1,3);%构造非线性方程组的系数阵b=zeros(3,1);1;%构造非线性方程组的常数项列plimit=abvec,val=eigs(p,1);vec=vec/sum(vec);plimit,vec海航司守奎(94339146) 19:38:441许哲，李号雷，基于灰色模型和Bootstrap理论的大规模定制质量控制方法研究，数学的实践与认识，第42卷第21期，121127，2012年11月。

5、海航司守奎(94339146) 19:47:27cumsum(x)海航司守奎(94339146) 19:47:33就是累加运算海航司守奎(94339146) 20:00:08clc,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=-z

6、k,ones(size(zk);yn=x0(2:end);ab=Byn%拟合参数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=mea

7、n(nnyb)%分别求4个子样本的均值jc=range(nnyb)%分别求4个子样本的极差xlswrite(hb.xls,nnyb;mu;jc)%把数据写到Excel文件中，便于做表使用b=rand(16,250);%产生16行250列的随机数矩阵h=floor(n*b)+1;%把每一列的随机数映射为相应的编号mx=nyb;bb=mx(h);bb=reshape(bb,4,1000);mmu=mean(bb);mjc=range(bb);%计算1000个子样本的均值和极差smu=sort(mmu);sjc=sort(mjc);%把均值和极差按照从小到大的次序排列alpha=0.01;k1=fl

8、oor(1000*alpha/2),k2=floor(1000*(1-alpha/2)mqj=smu(k1),smu(k2)%显示均值的置信区间jqj=sjc(k1),sjc(k2)%显示极差的置信区间subplot(1,2,1),plot(mu,*-),holdonplot(1,4,mqj(1),mqj(1),plot(1,4,mqj(2),mqj(2)ylabel(样本均值)subplot(1,2,2),plot(jc,*-),holdonplot(1,4,jqj(1),jqj(1),plot(1,4,jqj(2),jqj(2)ylabel(极差)海航司守奎(94339146) 20:12

9、:51clc,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=Byn%拟合参数a,bx=dsolve(Dx+a*x=b,x(0)=x00)%求微分方程的符号解x=su

10、bs(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)%求原始数据的预测值yuce0=double(yuce0);c=std(yuce0(1:n)/std(x0)%求后验差比值c海航司守奎(94339146) 20:14:0382.9 92 17.1 94 8.488.0 93 21.3 96 9.699.9 96 25.1 97 10.4105.3 94 29.0 97 11.4117.7 100 34.0 100 12.2131.0 101 40.0 10

11、1 14.2148.2 105 44.0 104 15.8161.8 112 49.0 109 17.9174.2 112 51.0 111 19.6184.7 112 53.0 111 20.8185.3 108 48.6 103 18.9175.9 99 52 99 20.5166.8 106 48.8 102 19.2156.8 105 49.7 98.8 18.6148.9 93 50.6 96.5 17.9136.5 96 47.9 97.9 17.6海航司守奎(94339146) 20:20:49clc,cleard=load(huigui2.txt);xs=ones(16,1),

12、d(:,1:4);cs=xsd(:,5)cs2,cs2int,r,rint,stats=regress(d(:,5),xs)cs,cs2海航司守奎(94339146) 20:36:35clc,cleard=load(huigui2.txt);xs=ones(16,1),d(:,1,3);cs1,cs1int,r,rint,stats=regress(d(:,5),xs)海航司守奎(94339146) 20:40:04clc,cleard=load(huigui2.txt);a=d(:,1:4);b=d(:,5);rstool(a,b)海航司守奎(94339146) 20:45:30clc,cl

13、eara=load(huigui.txt);AA=a(:,1:end-1);B=a(:,end);rstool(AA,B)n=size(a,1);xs=ones(n,1),AA,AA(:,1).*AA(:,2),AA(:,1).*AA(:,3),AA(:,1).*AA(:,4),AA(:,2).*AA(:,3),AA(:,2).*AA(:,4),AA(:,3).*AA(:,4),AA.2cs=xsB%为了验证完全二次式的回归答案海航司守奎(94339146) 19:37:02clc,clearp=0.080.90.02;0.080.90.02;0.240.70.06;a=p-eye(3);on

14、es(1,3);%构造非线性方程组的系数阵b=zeros(3,1);1;%构造非线性方程组的常数项列plimit=abvec,val=eigs(p,1);vec=vec/sum(vec);plimit,vec海航司守奎(94339146) 19:38:441许哲，李号雷，基于灰色模型和Bootstrap理论的大规模定制质量控制方法研究，数学的实践与认识，第42卷第21期，121127，2012年11月。海航司守奎(94339146) 19:47:27cumsum(x)海航司守奎(94339146) 20:57:13clc,cleara=load(huigui3.txt);x0=a(:,1:en

15、d-1);y0=a(:,end)stepwise(x0,y0,1:8)%初始时使用全部8个自变量海航司守奎(94339146) 21:02:28海航司守奎(94339146) 21:02:51clc,cleara=load(huigui3.txt);x0=a(:,1:end-1);y0=a(:,end)stepwise(x0,y0)%初始时不使用任何变量海航司守奎(94339146) 21:04:24下面开始答疑201304001(903528428) 21:05:46马尔科预测的适用范围是什么？微分霸天虎(1355097263) 21:07:44“匿名函数”格式是什么？201304001(9

16、03528428) 21:07:44预测极限状态海航司守奎(94339146) 21:08:34函数格式：function返回值=函数名（输入参数）海航司守奎(94339146) 21:09:10匿名函数格式：返回值=(输入参数)函数表达式201307031(820014430) 21:10:32逐步回归向前选择，向后删除，选择哪个方法更好？微分霸天虎(1355097263) 21:11:56何时用“匿名函数”？201304001(903528428) 21:12:34岭回归的含义是什么？海航司守奎(94339146) 21:14:28微分霸天虎(1355097263) 21:15:01b=z

17、eros(3,1);1具体是怎样的内容？海航司守奎(94339146) 21:15:02海航司守奎(94339146) 21:16:25201304001(903528428) 21:18:32如果有多个自变量,但自变量间可能相关，这用主成分回归或者偏最小二乘回归比较好，对吗？201304001(903528428) 21:22:33下次课讲？微分霸天虎(1355097263) 21:22:47建议老师尽量多分析程序来源海航司守奎(94339146) 21:22:48最后一次课讲随机模拟海航司守奎(94339146) 21:23:00什么意思？201304001(903528428) 21:2

18、3:18讲讲预测方法各自的特点和适用范围微分霸天虎(1355097263) 21:24:19就是matlab程序的构建201304001(903528428) 21:24:21谢谢老师！海航司守奎(94339146) 21:27:40海航司守奎(94339146) 21:28:12DescriptionY=wblpdf(X,A,B)computestheWeibullpdfateachofthevaluesinXusingthecorrespondingscaleparameter,Aandshapeparameter,B.X,A,andBcanbevectors,matrices,ormul

19、tidimensionalarraysthatallhavethesamesize.Ascalarinputisexpandedtoaconstantarrayofthesamesizeastheotherinputs.TheparametersinAandBmustbepositive.TheWeibullpdfisSomereferencesrefertotheWeibulldistributionwithasingleparameter.ThiscorrespondstowblpdfwithA=1.Examples海航司守奎(94339146) 21:28:37海航司守奎(9433914

20、6) 21:29:25画图的Matlab程序如下%数学与Matlab参数：b=tau;a=beta(-1/b);subplot(121)fx=(x,tau,beta)tau*beta*x.(tau-1).*exp(-beta*x.tau);fplot(x)fx(x,2,0.5),0,5,.)text(1.6,0.5,$beta=0.5,tau=2$,Interpreter,Latex)xlabel($x$,Interpreter,Latex)ylabel($f(x)$,Interpreter,Latex)title(数学上的参数表示)subplot(122)fx2=(x,a,b)wblpdf(

21、x,a,b);fplot(x)fx2(x,0.5(-1/2),2),0,5,.)str=a=,num2str(0.5(-1/2),b=2;text(1.6,0.5,str)xlabel($x$,Interpreter,Latex)ylabel($f(x)$,Interpreter,Latex)title(Matlab的参数表示)海航司守奎(94339146) 21:30:21201307014-程小乐(646092267) 09:16:21这是是不是不能在matlab7.0上运行啊201307019(347667369) 09:19:13为啥201307014-程小乐(646092267) 0

22、9:20:02201307019(347667369) 09:24:08呃，那倒真有可能201307019(347667369) 09:24:17工具箱里没有。201307004(82985573) 09:49:07用2013a吧201307014-程小乐(646092267) 09:49:33我们机房安装的是7.0的内民大-刘京(1374413450) 16:25:10net=newff(minmax(P),8,5,logsig,logsig,traingdx)这个BP网络里面的minmax（P）是什么意思？201307019(347667369) 16:32:45minmax不是最小值和最

23、大值的意思么。201307019(347667369) 16:35:21xmminmax（x）求x的最小值和最大值内民大-刘京(1374413450) 16:36:04内民大-刘京(1374413450) 16:36:15这步变化是怎么得到的？201307019(347667369) 16:37:10p是个啥？内民大-刘京(1374413450) 16:37:21P是我定义的矩阵201307019(347667369) 16:37:25矩阵？内民大-刘京(1374413450) 16:37:32对201307019(347667369) 16:38:59呃，对矩阵做minmax我就不清楚了，要

24、不help一下看看201307029(785144557) 16:40:54是不是第一行最小值0，最大值1？第二行也是？内民大-刘京(1374413450) 16:41:10我DOC了一下。内民大-刘京(1374413450) 16:41:17这个例子。完全看不明白。201307019(347667369) 16:41:26按行返回结果？201307029(785144557) 16:41:34doc哪个内民大-刘京(1374413450) 16:41:37按列返回结果。内民大-刘京(1374413450) 16:41:42DOCminmax201307029(785144557) 16:42

25、:04那就不知道了201307019(347667369) 16:42:12哦201307029(785144557) 16:42:28bp是什么201307005(1445891025) 16:42:39minmax(A)即是求解A矩阵中每个行向量中最小值和最大值201307005(1445891025) 16:42:58bp里面用的不是mapminmax吗内民大-刘京(1374413450) 17:24:16哦哦哦哦哦。内民大-刘京(1374413450) 17:24:19我去。内民大-刘京(1374413450) 17:24:23这么简单内民大-刘京(1374413450) 17:24:40内民大-刘京(1374413450) 17:24:46就是没想明白。内民大-刘京(1374413450) 17:24:49谢谢内民大-刘京(1374413450) 17:25:18那mapminmax是啥意思？201307017(122192130) 17:25:47标准化？201307005(1445891025) 18:22:08mapminmax是归一化处理，防止bp网络里面如果有些参量过大可能影响bp网络的结果。.海航司守奎(94339146) 21:31:5388