1、第一题：1.有 MATLAB进行灰色模型求解。clear all;clc;a = 2750 2450 2050 1800 1300;n = length(a);a = a;a1 = cumsum(a); %一次累加数列B=-0.5*(a1(1:end-1)+a1(2:end),ones(n-1,1);Y=a(2:end);r=BY;y=dsolve(Dy+a*y=b,y(0)=y0);y=subs(y,a,b,y0,r(1),r(2),a1(1);yuce1=subs(y,t,0:2*n-1);digits(6),y=vpa(y)%微分方程yuce=diff(yuce1)yuce=a(1),y
2、uce %预测值2.有 MATLAB对灰色模型进行检验。ex=ec./x;r=0;rou=0.5;for k=1:nr=r+rou* s(ec(k)+rou*max(a (ec);endr=r/n;%原始序列的标准差s1=std(x);%计算残差的标准差s2=std(ec) %计算 CC=s2/s1;%计算后验概率deta=ec-mean(ec);index=fineta)0.6745*s1);P=length(index)/n;3.画出预测的长江江豚数量的变化的图。x=1:5:71y=2700 2422 1991 1636 1342 1099 899 734 597 485 391 315
3、251 199 155;plot(x,y,+b)hold on;plot(x,y,r)x=xlabel(时间)y=ylabel(江豚的数量)hh=title(预测长江江豚数量的变化)4.利用 MATLAB拟合指数函数：在 M文档中输入下面的命令：function f=cur(x,data)f=x(1)+x(2)*exp(x(3)*data) 再在命令窗口输入如下的命令： x=1:5:71y=2700 2422 1991 1636 1342 1099 899 734 597 485 391 315 251 199 155;x0=1,-0.0001,0.01;x=lsqcurvefit (curvefun1,x0,tdata,cdata)5.用 SPSS对 MATLAB拟合指数函数的检验：* 曲线估计. TSET NEWVAR=NONE. CURVEFIT /VARIABLES=长江江豚数量 WITH 时间/CONSTANT /MODEL=EXPONENTIAL /PLOT FIT.可得到如下图表：卓金武，魏永生，秦健，李必文 MATLAB在数学建模中的应用，北京：航空航天大学出版社，2010 年。肖 文，张先锋，长江江及其支流长江江豚种群数量及分布/中国科学院水生生物研究所，2001 年 04 期，
