示例

原创 2013年12月03日 18:03:09


1 聚类分析

  1. 样题

    8 个样品, 每个样品有2 个指标(m=2, 2 维变量), 它们的量纲相同, 试用绝对距离对以上样本进行聚类。

    编号 1 2 3 4 5 6 7 8
    x1 2 2 4 4 -4 -2 -3 -1
    x2 5 3 4 3 3 2 2 -3
    x=[2,2,4,4,-4,-2,-3,-1;
    5,3,4,3,3,2,2,-3]';
    y=pdist(x,'cityblock');
    obslabel=num2str(1:8);
    Z=squareform(y);
    Z=linkage(y,'single')
    
    Z =
    
         6     7     1
         3     4     1
         5     9     2
         1     2     2
        10    12     2
        11    13     5
         8    14     6
    
    H = dendrogram(Z);
    

    matlabclust.png

2 Markov链稳定分布

  1. 设Markov 链的转移概率矩阵为
    0.7 0.1 0.2
    0.1 0.8 0.1
    0.05 0.05 0.9

    求它的平稳分布.

  2. 马尔科夫稳定性求解
    P=matrix(c(0.7,0.1,0.2,0.1,0.8,0.1,0.05,0.05,0.9),ncol=3,byrow=TRUE)
    A=rbind((t(P)-diag(rep(1,3)))[1:2,],rep(1,3));
    a<-solve(A,c(0,0,1))
    
    c(3/17,4/17,10/17)
    
    [1] 0.176 0.235 0.588
    

3 相似系数

  1. 样题

    已知三个三维样本 X1=(5,2,8) , X2=(11,7,10) , X3=(7,6,11), 求它们的相似系数矩阵

    -5 2 8
    -11 7 10
    7 6 11
    ρX,Y=i=1n(xix¯)(yiy¯)i=1n(Xix¯)2i=1n(yiy¯)2
  2. 求解过程

    分别计算第一个变量的标准化

    x1=c(-5,2,8);
    x2=c(-11,7,10);
    x3=c(7,6,11);
    cor(cbind(x1,x2,x3))
    
          x1    x2    x3
    x1 1.000 0.940 0.726
    x2 0.940 1.000 0.449
    x3 0.726 0.449 1.000
    
    x1=[-5,2,8];
    x2=[-11,7,10];
    cor(x1,x2);
    

4 回归分析

  1. 例题

    为考察碳含量(x)和合金强度(y)的关系,收集了 12 种合金,分别测得其碳含量和强度如下

    序号 1 2 3 4 5 6 7 8 9 10.0 11 12
    x 0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.2 0.21 0.23
    y 42.0 43.50 45 45.50 45 47.50 49 53 50 55.0 55 60

    试根据上述数据进行一元回归分析,找出碳含量和合金强度之间的经验回归模型

  2. 求解

    计算可得

    x¯=0.1583,y¯=49.21i=112x2i=0.3194,i=112xiyi=95.925Sxx=i=112x2in(x¯)2=0.3194120.15832=0.0186,Sxy=i=112xiyinx¯y¯=2.4292β^1=SxySxx=130.83,β^0=y¯x¯β^1=28.49
  3. R 求解
    thl<-c(0.1,0.11,0.12,0.13,0.14,0.15,0.16,0.17,0.18, 0.2,0.21,0.23)
    qd<-c(42.0,43.50,45.00,45.50,45.00,47.50,49.00,53.00,50.00,55.0,55.00,60.00)
    lm(qd~thl)
    
    Call:
    lm(formula = qd ~ thl)
    
    Coefficients:
    (Intercept)          thl
          28.49       130.83
    


相关文章推荐

c++ primer plus6th 示例源代码

  • 2017年11月03日 20:26
  • 363KB
  • 下载

matplotlib基本画图示例

在这里我介绍一下plot方法的函数的使用。 Series.plot方法的函数: 参数 说明 label 用于图例的标签 ax 要在其上进行绘制的matplotlib...

echarts动态关系图示例

  • 2017年11月02日 17:38
  • 12KB
  • 下载

android studio 使用NDK和swig编译c++示例

  • 2017年10月30日 10:16
  • 10.41MB
  • 下载

Exec Maven Plugin全面解析和使用示例

本文原文连接: http://blog.csdn.net/bluishglc/article/details/7622286 ,转载请注明出处!1.为什么使用exec?现在的工程往往依赖众多的jar包...

DELPHI HL7示例

  • 2017年10月27日 11:11
  • 2.83MB
  • 下载

Swiper 3 的特色功能、示例

无需加载公用库:Swiper无需加载任何公共库(如jquery)即可运行,这保证了Swiper的轻量和运行速度。Swiper也可以在加载了公共库的环境下安全的运行,如jQuery, Zepto, jQ...
  • henrypt
  • henrypt
  • 2016年12月06日 09:10
  • 707

spring+springmvc+mybatis 示例

  • 2017年10月25日 23:58
  • 33.93MB
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:示例
举报原因:
原因补充:

(最多只允许输入30个字)