灰色预测水果销量

clearsyms a b;c=[a b]';H1 = [97209375  25245141    18823858    14280912    6860331     8328256     2211484     359329 98147156    28009225    20595451    13286870    6962323     9023196     2471999     431855 111824181   32678720    23644462    14726781    7994102     9652477     2922246     366416 110304098   31279460    23330802    14471040    8342154     9716920     2470808     632805 114710579   31009584    24121896    15341226    8494085     10863106    2880975     733774 122267474   33867597    25546753    16743196    9566318     11122704    2956350     898121 126608205   33773914    27165180    17098139    11211999    11728772    2622915     945448 ];H = H1';title_name = {'水果','苹果','柑橘','梨','葡萄','桃','杏','猕猴桃'};row_name = {'A1:A14','B1:B14','C1:C14','D1:D14','E1:E14','F1:F14','G1:G14','H1:H14'};for m = 1:8A = H(m,:);row_str = row_name(:,m);B=cumsum(A);  % 原始数据累加n=length(A);for i=1:(n-1)    C(i)=(B(i)+B(i+1))/2;  % 生成累加矩阵end% 计算待定参数的值D=A;D(1)=[];D=D';E=[-C;ones(1,n-1)];c=inv(E*E')*E*D;c=c';a=c(1);b=c(2);% 预测后续数据F=[];F(1)=A(1);for i=2:(n+7)    F(i)=(A(1)-b/a)/exp(a*(i-1))+b/a ;endG=[];G(1)=A(1);for i=2:(n+7)    G(i)=F(i)-F(i-1); %得到预测出来的数据end t1=2007:2013;t2=2007:2020;Gxlswrite('E:\陕科大金专\数学建模\真题练习\2014年真题\水果产量.xls',G','Sheet2',row_str{1,1});plot(t1,A,'ko', 'LineWidth',2)hold onplot(t2,G,'k', 'LineWidth',2)xlabel('年份', 'fontsize',12)ylabel('产量（万吨）','fontsize',12)set(gca,'LineWidth',2);title(title_name(:,m))figureend