下面是源代码(fc.m)：
function y=fc(isPloy,factorial);
%function y=fc(isPloy,factorial);
%factorial表示预测函数的阶乘个数,默认是5
%isPloy判断参数的个数，来决定是否，每个阶都进行画图,默认是0(不进行画图)
%close %关闭图形界面
%clear %清理系统变量
if nargin==1 %判断参数的个数，来决定是否，每个阶都进行画图
isPloy=1;
factorial=5;
elseif nargin>=2
else
isPloy=0;
factorial=5;
end
%n=factorial;
fd=fopen('fcout.dat','w+'); %打开结果的保存文件
load fcin.dat; %加载历史数据
x1=fcin(:,1); %加载第一列为x1
y1=fcin(:,2); %加载第二列为y1，红球1
y2=fcin(:,3); %加载第三列为y2，红球2
y3=fcin(:,4); %红球3
y4=fcin(:,5); %红球4
y5=fcin(:,6); %红球5
y6=fcin(:,7); %红球6
y7=fcin(:,8); %篮球
nextPoint=x1(length(x1))+1; %从第一列的最后一个值中,获取下一个预期x轴值
k=nextPoint;
for n=1:1:factorial %n表示预测函数的阶数
%%计算出预测函数的多项式参数
p1=polyfit(x1, y1, n);
p2=polyfit(x1, y2, n);
p3=polyfit(x1, y3, n);
p4=polyfit(x1, y4, n);
p5=polyfit(x1, y5, n);
p6=polyfit(x1, y6, n);
p7=polyfit(x1, y7, n);
xi=linspace(1,k , k); %产生预测的x轴向量
%%计算出根据这些函数，预测出来的与x1对应的y轴值
z1=polyval(p1, xi);
z2=polyval(p2, xi);
z3=polyval(p3, xi);
z4=polyval(p4, xi);
z5=polyval(p5, xi);
z6=polyval(p6, xi);
z7=polyval(p7, xi);
%打印预测向量的最后一个值，为预测值
%fprintf('阶数为:%3d, 红球%3d,%3d,%3d,%3d,%3d,%3d,篮球为:%3d/n',n,round(z1(k)),round(z2(k)),round(z3(k)),round(z4(k)),round(z5(k)),round(z6(k)),round(z7(k)));
fprintf(fd,'阶数为:%3d, 红球%3d,%3d,%3d,%3d,%3d,%3d,篮球为:%3d/n',n,round(z1(k)),round(z2(k)),round(z3(k)),round(z4(k)),round(z5(k)),round(z6(k)),round(z7(k)));
if isPloy>0
val=input('stop,enter');%起到暂停作用
plot(x1, y1,xi,z1, '+', x1, y2,xi,z2,'+', x1, y3,xi,z3,'+', x1, y4,xi,z4,'+', x1, y5,xi,z5 ,'+', x1, y6,xi,z6,'+', x1, y7,xi,z7 ,'o' );
end
end
fclose(fd);
close
type fcout.dat
%画图
plot(x1, y1,xi,z1, '+', x1, y2,xi,z2,'+', x1, y3,xi,z3,'+', x1, y4,xi,z4,'+', x1, y5,xi,z5 ,'+', x1, y6,xi,z6,'+', x1, y7,xi,z7 ,'o' );
%plot(x1, y1,xi,z1, x1, y2,xi,z2, x1, y3,xi,z3, x1, y3,xi,z3, x1, y4,xi,z4 , x1, y5,xi,z5, x1, y6,xi,z6 ,x1, y7,xi,z7 );
下面是 fcint.dat的数据：
01 03 05 18 20 27 33 01 2007-125
02 09 10 19 23 26 31 09 2007-126
03 06 09 13 16 24 28 11 2007-127
04 09 10 19 21 27 31 05 2007-128
05 05 07 20 21 22 30 08 2007-129
06 03 05 09 11 27 31 04 2007-130
07 03 05 07 16 22 27 05 2007-131
08 01 09 16 21 22 23 05 2007-132
09 03 06 07 11 13 33 10 2007-133
10 01 04 10 13 18 25 15 2007-134
11 01 11 16 26 31 33 16 2007-135
12 01 02 18 21 25 29 14 2007-136
13 03 07 08 18 20 22 03 2007-137
14 02 03 15 17 19 25 16 2007-138
15 06 10 12 14 16 22 06 2007-139
16 01 05 16 21 22 26 11 2007-140
17 02 03 04 06 17 31 08 2007-141
18 11 20 25 26 27 30 08 2007-142
19 01 06 22 23 24 26 04 2007-143
20 08 14 23 25 28 32 16 2007-145
21 08 09 11 12 25 31 11 2007-145
22 04 18 22 24 26 30 09 2007-146
23 03 07 18 24 26 27 04 2007-147
24 03 09 16 17 23 28 07 2007-148
25 01 17 19 22 28 30 03 2007-149
26 03 05 11 13 19 24 05 2007-150
27 01 06 10 11 23 25 02 2007-151
28 11 17 21 29 30 33 08 2007-152
29 01 04 19 20 25 31 15 2007-153
30 02 04 07 09 14 29 03 2008-1
31 03 04 18 22 25 29 09 2008-2
|
|
发表于 2018-11-20 17:38:11
楼主
