下面是源代码(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

用matlab预测双色球的方法相关推荐

  1. caffe matlab 重启,Matlab caffe 具体使用方法

    Matlab caffe 具体使用方法 caffe 的matlab接口网络上的资料比较少,这里把几个基本的操作简单介绍,主要包括: classification_demo的简单介绍 卷积核的显示 各层 ...

  2. Matplotlib可视化预测双色球

    步骤 1.获得数据 2.清洗保存数据 3.读取数据并可视化 4.得出结论 首先,简单说下国民彩票双色球的规则:选6位数字作为红球和1位数字作为篮球组成一注彩票,红球的区间是1~33,蓝球的区间是1 ~ ...

  3. 写一个预测双色球的代码

    预测双色球的代码需要考虑以下几个方面: 定义一个函数来生成双色球号码.双色球号码由 6 个红球和 1 个蓝球组成,红球号码从 1 到 33 中选取,蓝球号码从 1 到 16 中选取. 在生成号码的同时 ...

  4. 【 MATLAB 】DFT的性质讨论(二)序列的循环移位及其 MATLAB 实现(频域方法)

    上篇博文:[ MATLAB ]DFT的性质讨论(二)序列的循环移位及其 MATLAB 实现(时域方法) 提到了对序列x(n)做循环移位后的DFT形式为: 上篇博文已经讨论过了第一种实现循环移位的方法, ...

  5. matlab中 三种方法计算 Ax b,在MATLAB中,方程Ax=B的解可以用哪个命令求得? matlab 求助 解方程组...

    matlab中解方程组还是很方便的,例如,对于代数方程组Ax=b(A为系数矩阵,非奇异)的求解,MATLAB中有两种方法: (1)x=inv(A)*b - 采用求逆运算解方程组: (2)x=A\B - ...

  6. Paper:《A Unified Approach to Interpreting Model Predictions—解释模型预测的统一方法》论文解读与翻译

    Paper:<A Unified Approach to Interpreting Model  Predictions-解释模型预测的统一方法>论文解读与翻译 导读:2017年11月25 ...

  7. ML之回归预测之BE:利用BE算法解决回归(实数值评分预测)问题—线性方法解决非线性问题

    ML之回归预测之BE:利用BE算法解决回归(实数值评分预测)问题-线性方法解决非线性问题 目录 输出结果 设计思路 代码实现 输出结果 设计思路 代码实现 for row in xList:newRo ...

  8. Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动

    Algorithm之MC:基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 目录 输出结果 设计代码 输出结果 基于Matlab实现通过蒙特卡洛方法模拟二维布朗运动 设计代码 %基于Matlab实 ...

  9. matlab中怎么画五维图,使用MatLab绘制散点图的方法

    使用MatLab绘制散点图的方法: 1.首先从给出的excel表格中采用数据筛选方法,将所有某类型地区的相关数据全部筛选 出来.例如第1类区域的数据如下: 4043 1895 14 1 9.17 24 ...

最新文章

  1. 提供openssl -aes-256-cbc兼容加密/解密的简单python函数
  2. 主成分分析计算机如何实现,主成分分析(PCA)原理及R语言实现及分析实例
  3. opencv利用矩形框选中某一区域_【从零学习OpenCV】4Ubuntu系统中安装OpenCV 4
  4. Programming Protocol-independent Packet Processors (P4)
  5. 91Android万能驱动最新版,NVIDIA万能驱动电脑版下载|NVIDIA万能驱动官方安装版下载 -爱安卓...
  6. 定时任务getScheduler
  7. Android实现友盟分享小程序
  8. 微信公众服务号如何快速申请注册并认证开通支付功能
  9. 最好用的共享文件服务器,文件共享有哪些方式,哪种比较好?
  10. 愿天下有情人都是失散多年的兄妹 (25 分)
  11. 默认关闭idea2020的注释doc的rendered view模式
  12. 浏览器主页劫持的解决办法
  13. kakfka2.2源码分析之handleFetchRequest
  14. html怎么查看cad文件,怎么将电脑CAD文件在手机上查看?原来怎么简单!
  15. 计算机夏令营英语自我介绍,夏令营英语自我介绍
  16. 计算机基础 电子工业出版,电子工业出版社21世纪计算机基础教育系列教材2吴功宜吴英编着.ppt...
  17. 万物云开启路演:拟募资8亿美元 淡马锡与高瓴是基石投资者
  18. 数字调制系列:IQ调制基本理论
  19. onvif 视频28181 1400 区别
  20. matlab工具箱拟合函数

热门文章

  1. 计算机速记速录的发展,速录发展趋势及技巧
  2. Oracle密码策略设置指引
  3. 解决导入obj模型时出现模型镂空的问题
  4. 区块链解决方案-最新全套文件
  5. 2023.3.5-TensorFlow学习
  6. 快速给视频添加背景和水印的步骤
  7. 《薄冰实用英语语法详解》独家连载之十六:名词性从句
  8. AWS S3 V4签名实现(nodejs)
  9. Mac电脑使用:通过Mac电脑导出iPhone手机里面的照片(不是iTunes导出照片的方法)
  10. A Thorough Examination of the CNN/Daily Mail Reading Comprehension Task