[从头学数学] 第175节 统计

最新推荐文章于 2023-01-10 21:39:00 发布
最新推荐文章于 2023-01-10 21:39:00 发布
阅读量3.1k 收藏 1
点赞数 1
分类专栏： 从头学数学 从头学数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mwsister/article/details/51134240
版权
在星历2016年的某一天，银河系厄尔斯星球上的工程师阿伟与机器小伟共同探讨统计学。小伟设计了一个随机样本工具，并运用最小二乘法进行数据拟合。他们对数据的深入分析得到了[人叫板老师]的认可，揭示了统计学在解决实际问题中的应用。
摘要由CSDN通过智能技术生成
剧情提要：
[机器小伟]在[工程师阿伟]的陪同下进入了结丹中期的修炼，
这次要修炼的目标是[统计]。

正剧开始：

星历2016年04月12日 16:00:33, 银河系厄尔斯星球中华帝国江南行省。
[工程师阿伟]正在和[机器小伟]一起研究[统计]。






小伟想了一下，做了这样一个取随机样本的工具：

<span style="font-size:18px;">>>> ================================ RESTART ================================
>>> 
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
[10, 1, 9]
>>> ================================ RESTART ================================
>>> 
[1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
[5, 4, 8]


#从N个样本中抽取n个样本
def sample(N, n):
    a = list(range(1, N+1));
    print(a);
    random.shuffle(a);
    return a[0:n];

if __name__ == '__main__':
    print(sample(10, 3));</span>

还是能玩一玩的。 









<span style="font-size:18px;">	if (1) {
		var config = new PlotConfiguration();        
		config.init();    
		  
		config.setPreference();    
		var r = 20;  
		//config.setSector(1,1,1,1);    
		//config.graphPaper2D(0, 0, r);  
		//config.axis2D(0, 0,180);    
	  
		var stat = new Statistic();          
		var data = [4,8,15,22,25,14,6,4,2];  
		var text = [];  
		  
		//初始化各组频数  
		//0到60岁，5岁一组  
		for (var i = 0; i < 5; i+=0.5) {  
			text.push(i.toFixed(1));  
		}     

		  
		stat.init(data, '月用水量(t)', '人数', 1);   
		stat.histogram(text, 0, 0);   
	
	}
	</span>


 


<span style="font-size:18px;">	if (1) {
		var config = new PlotConfiguration();        
		config.init();    
		  
		config.setPreference();    
		var r = 20;  
		//config.setSector(1,1,1,1);    
		//config.graphPaper2D(0, 0, r);  
		//config.axis2D(0, 0,180);    
	  
		var stat = new Statistic();          
		var data = [4,8,15,22,25,14,6,4,2];  
		for (var i = 0; i < data.length; i++) {
			data[i]/=100;
		}
		
		var text = [];  
		  
		//初始化各组
		for (var i = 0; i < 5; i+=0.5) {  
			text.push(i.toFixed(1));  
		}     

		  
		stat.init(data, '月用水量(t)', '频率', 1);   
		stat.rateHistogram(text, 0, 0);   
	
	}
	
	//频率分布直方图
	this.rateHistogram = function(lableArray, xOffset, yOffset, yMin, yMax) {
		lableArray = lableArray ? lableArray : [];
		var lables = lableArray.length;
		xOffset = xOffset ? xOffset : 0;
		yOffset = yOffset ? yOffset : 0;
		
		var colorArray = ['red', 'orange', '#0088FF', 'green', 'cyan', 'blue', '#FF00FF',
			'#FF8844', 'purple'];
		var colors = colorArray.length;
			
		var height = 380, width = 580;
		
		plot.save()
			.translate(xOffset+60, yOffset+50);

		plot.setLineWidth(2)
			.setTextAlign('right');
			

		var max = yMax ? yMax : this.max();  
		var min = yMin ? yMin : 0; 		

		//最大值的末位为0的近似数，比如最大值25，最合适的近似数为30
		var dy = Math.round((max-min)/10*100)/100;

		
		var size = this.size();  
		var perH = Math.round((height-100) / 10);  
		var perW = Math.round((width-100) / (size+2));
		
		//宽和高度边界
		var wBound = (size+2)*perW, hBound = 10*perH;
		plot.setLineWidth(5)
			.strokeRect(0, 0, wBound, hBound);
			
		this.axis2D(0, hBound, wBound+20, hBound+20, this.xLabel, this.yLabel);
			
		plot.setLineWidth(2);
		var count = 0;
		for (var i = hBound; i >-1; i -= perH) {
			plot.fillText((min+dy*count).toFixed(2), -10, i+10, 30);
			count++;
			
			if (i > 0) {
				plot.beginPath()
					.moveTo(0, i)
					.lineTo(wBound, i)
					.closePath()
					.stroke();
			}
		}
		
		for (var i = 0; i < size+2; i++) {
			plot.beginPath()
				.moveTo(i*perW, 0)
				.lineTo(i*perW, hBound)
				.closePath()
				.stroke();
		}
		

		
		var xpos, xpos2;
		
		for (var i = 0; i < size; i++) { 
			var h = -(this.statisticalSample[i]-min)/dy*perH;

			xpos = perW*(1+i);
			xpos2 = xpos + 0.5*perW;
			
			//对于值为0的组，就不多加操作了
			if (h < 0) {
				plot.setFillStyle(colorArray[i%colors]);  
				plot.fillRect(perW*(1+i), hBound, perW, h); 
				
				
				if (i == 0) {
					xpos += 0.75*perW;
				}
				else {
					xpos += 0.5*perW;
				}
			
				plot.setFillStyle('blue');
				plot.fillText(this.statisticalSample[i].toString(), xpos, h+hBound-5, 100);
			}
			
			plot.setFillStyle('blue');
			plot.setTextAlign('center');
			if (i < lables) {
				plot.fillText(lableArray[i], xpos2-0.5*perW, 
					hBound+30, 100);  

			}
			
			//增加x轴的最后一个标注
			if (i == size-1) {
				var n1 = parseFloat(lableArray[i]);
				var n2 = parseFloat(lableArray[i-1]);
				var n3 = (2*n1-n2);
				var s;
				if (Math.abs(n3-Math.round(n3))<0.001) {
					s = n3.toFixed(0);
				}
				else {
					s = n3.toFixed(1);
				}
				plot.fillText(s, xpos2+0.5*perW, hBound+30, 100); 
			}
			//plot.fillText(this.statisticalSample[i].toFixed(0), xpos2, hBound+40, 100);  
		}  

		
		plot.restore();
	
	}</span>


 













<span style="font-size:18px;">	//一组点    
    if (1) {  
        var r = 20;        
        config.setSector(1,1,1,1);          
        config.graphPaper2D(0, 0, r);        
        config.axis2D(0, 0,180);          
            
        //坐标轴设定    
        var scaleX = 2*r, scaleY = 2*r;      
        var spaceX = 10, spaceY = 10;       
        var xS = -10, xE = 80;      
        var yS = -10, yE = 50;      
        config.axisSpacing(xS, xE, spaceX, scaleX, 'X');        
        config.axisSpacing(yS, yE, spaceY, scaleY, 'Y');     
          
          
        var X=[ 23,27,39,41,45,49,50,53,54,56,57,58,60,61],   
            Y=[ 9.5,17.8,21.2,25.9,27.5,26.3,28.2,29.6,30.2,31.4,30.8,33.5,35.2,34.6];  
              
        var array = [];  
        var size = X.length;  
          
        for (var i = 0; i < size; i++) {  
            array.push([X[i], Y[i]]);  
        }  
              
        var transform = new Transform();  
        var tmp = [];  
          
        array = transform.scale(transform.translate(array, 0, 0), scaleX/spaceX, scaleY/spaceY);      
        tmp = [].concat(array);        
        shape.pointDraw(tmp, 'green');        
        //tmp = [].concat(array);        
        //shape.multiLineDraw(tmp, '#CCFF22');      
              
      
    }
	</span>

进行拟合，这个就是用了最小二乘法的公式了：

<span style="font-size:18px;">>>> 
X= [23, 27, 39, 41, 45, 49, 50, 53, 54, 56, 57, 58, 60, 61]
Y= [9.5, 17.8, 21.2, 25.9, 27.5, 26.3, 28.2, 29.6, 30.2, 31.4, 30.8, 33.5, 35.2, 34.6]
拟合结果: y =    0.57648 x +   -0.44780 , r=   0.97074

def linefit(x , y):  
    N = float(len(x))  
    sx,sy,sxx,syy,sxy=0,0,0,0,0  
    for i in range(0,int(N)):  
        sx  += x[i]  
        sy  += y[i]  
        sxx += x[i]*x[i]  
        syy += y[i]*y[i]  
        sxy += x[i]*y[i]  
    a = (sy*sx/N -sxy)/( sx*sx/N -sxx)  
    b = (sy - a*sx)/N  
    r = abs(sy*sx/N-sxy)/math.sqrt((sxx-sx*sx/N)*(syy-sy*sy/N))  
    return a,b,r  
  
  
def tmp():  
    X=[ 23,27,39,41,45,49,50,53,54,56,57,58,60,61];
    Y=[ 9.5,17.8,21.2,25.9,27.5,26.3,28.2,29.6,30.2,31.4,30.8,33.5,35.2,34.6];  
    a,b,r=linefit(X,Y)  
    print("X=",X)  
    print("Y=",Y)  
    print("拟合结果: y = %10.5f x + %10.5f , r=%10.5f" % (a,b,r) )  </span>

结果： 


<span style="font-size:18px;">    if (1) {  
        var r = 20;        
        config.setSector(1,1,1,1);          
        config.graphPaper2D(0, 0, r);        
        config.axis2D(0, 0,180);          
            
        //坐标轴设定    
        var scaleX = 2*r, scaleY = 2*r;      
        var spaceX = 10, spaceY = 10;       
        var xS = -10, xE = 80;      
        var yS = -10, yE = 50;      
        config.axisSpacing(xS, xE, spaceX, scaleX, 'X');        
        config.axisSpacing(yS, yE, spaceY, scaleY, 'Y');     
          
          
        var X=[ 23,27,39,41,45,49,50,53,54,56,57,58,60,61],   
            Y=[ 9.5,17.8,21.2,25.9,27.5,26.3,28.2,29.6,30.2,31.4,30.8,33.5,35.2,34.6];  
              
        var array = [];  
        var size = X.length;  
          
        for (var i = 0; i < size; i++) {  
            array.push([X[i], Y[i]]);  
        }  
              
        var transform = new Transform();  
        var tmp = [];  
          
        array = transform.scale(transform.translate(array, 0, 0), scaleX/spaceX, scaleY/spaceY);      
        tmp = [].concat(array);        
        shape.pointDraw(tmp, 'green');       

		array = [];
        for (var i = 0; i < size; i++) {  
            array.push([X[i], taskFun(X[i])]);  
        }  		
		array = transform.scale(transform.translate(array, 0, 0), scaleX/spaceX, scaleY/spaceY);  
        tmp = [].concat(array);        
        shape.multiLineDraw(tmp, 'red');     

		plot.setFillStyle('blue');
		plot.fillText('拟合结果: y = 0.57648 x - 0.44780', -270, -170, 200);
              
      
    }
	
	
}

function taskFun(x) {
	return 0.57648*x -0.44780;
}
</span>







这位勤奋的小店主，小伟决定要好好算一下他的这些数据。







<span style="font-size:18px;">>>> 
X= [-5, 0, 4, 7, 12, 15, 19, 23, 27, 31, 36]
Y= [156, 150, 132, 128, 130, 116, 104, 89, 93, 76, 54]
拟合结果: y =   -2.35170 x +  147.76710 , r=   0.98113


	//一组点    
    if (1) {  
        var r = 20;        
        config.setSector(1,1,1,1);          
        config.graphPaper2D(0, 0, r);        
        config.axis2D(0, 0,180);          
            
        //坐标轴设定    
        var scaleX = 2*r, scaleY = 2*r;      
        var spaceX = 10, spaceY = 50;       
        var xS = -10, xE = 50;      
        var yS = -10, yE = 200;      
        config.axisSpacing(xS, xE, spaceX, scaleX, 'X');        
        config.axisSpacing(yS, yE, spaceY, scaleY, 'Y');     
          
          
        var X=[-5, 0, 4, 7, 12, 15, 19, 23, 27, 31, 36],   
            Y=[ 156, 150, 132, 128, 130, 116, 104, 89, 93, 76, 54];  
              
        var array = [];  
        var size = X.length;  
          
        for (var i = 0; i < size; i++) {  
            array.push([X[i], Y[i]]);  
        }  
              
        var transform = new Transform();  
        var tmp = [];  
          
        array = transform.scale(transform.translate(array, 0, 0), scaleX/spaceX, scaleY/spaceY);      
        tmp = [].concat(array);        
        shape.pointDraw(tmp, 'green');       

		
		array = [];
        for (var i = 0; i < size; i++) {  
            array.push([X[i], taskFun(X[i])]);  
        }  		
		array = transform.scale(transform.translate(array, 0, 0), scaleX/spaceX, scaleY/spaceY);  
        tmp = [].concat(array);        
        shape.multiLineDraw(tmp, 'red');     

		plot.setFillStyle('blue');
		plot.fillText('拟合结果: y = -2.35170 x +  147.76710', -270, -170, 300);
              
      
    }
	
	
	if (0) {
		var config = new PlotConfiguration();        
		config.init();    
		  
		config.setPreference();    
		var r = 20;  
		//config.setSector(1,1,1,1);    
		//config.graphPaper2D(0, 0, r);  
		//config.axis2D(0, 0,180);    
		
		var X=[-5, 0, 4, 7, 12, 15, 19, 23, 27, 31, 36],   
            Y=[ 156, 150, 132, 128, 130, 116, 104, 89, 93, 76, 54];  
	  
		var stat = new Statistic();          
		var data = Y; 

		var text = [];  
		  
		//横坐标
		for (var i = 0; i < X.length; i++) {  
			text.push(X[i].toFixed(0));  
		}     

		  
		stat.init(data, '气温', '卖出热饮（杯）', 1);   
		stat.barChart(text, 0, 0);   
	
	}
	
	
}

function taskFun(x) {
	return -2.35170*x +  147.76710;
}</span>

看来[人叫板老师]也认为是这个答案。





本节到此结束，欲知后事如何，请看下回分解。

mwsister
关注
专栏目录
输入 n , 找 range(1, n+1) 的最小公倍数 - 华为技术面手撕代码（Python）
学习,你不是一个人在战斗 ～ 961302305（QQ,微信）
03-14 122
输入 n , 找 range(1, n+1) 的最小公倍数
【ChatGPT核心原理实战】手动求解 Transformer：分步数学示例 | Solving Transformer by Hand: A Step-by-Step Math Example
程序员光剑
12-22 3011
手动求解 Transformer：分步数学示例Understanding Transformers: A Step-by-Step Math Example — Part 1了解 Transformer：分步数学示例 — 第 1 部分I understand that the transformer architecture may seem scary, and you might have encountered various explanations on…我知道变压器架构可能看起来很可怕，并且
用python中的for循环语句，写1+2+3+4+....+N的代码
weixin_42588877的博客
01-10 4949
可以使用Python的range()函数来生成1到N之间的整数，然后使用for循环语句来遍历这些整数。代码如下： N = 10 for i in range(1, N+1): print(i) 这个代码会输出1到10之间的整数。如果要打印其它范围内的整数，只需要更改range(1, N+1)中的1和N+1的值即可。 range(1, N+1)中的1是起始值，N+1是结束值,左闭右开区间。...
（数据结构python版）第1章python入门
weixin_42194337的博客
08-23 388
python入门 python概述 python解释器 python是一种解释语言，命令在python解释器的软件中执行 python解释器可以通过再命令行中输入“python”启动，可以使用一个额外的-i标志进入交互模式（python -i demo.py） python程序预览 python语法依赖，缩进 一条命令写在多行，利用反斜杠（\） 注释：# python对象 标识符，对象和赋值语句 ...
Python学习笔记
XHQ1999的博客
03-29 2251
第一章 Python的基础语法 1.1print()函数及变量引用 print()的用途： （1）用来打印数字 （2）用来打印字符串，需要加上引号’’ （3）用来计算 注意：打印多个字符时，需要使用逗号隔开。多个字符串相连可以使用逗号或者+号相连。 常用的转义字符： （1）\n，表示新起一行显示后面的内容 （2）\t,空格 （3）\r，对\r前的所有内容进行覆盖 （4）\b，对\b前按删除键一次，退格按一次 （5）\，保留后面的内容输出 （6）r，输出原字符，让我们转义字符失效 标识符不能是保留字，不能以数
python-用range()函数快速生成列表
qq_37143903的博客
06-01 1万+
任务1描述 本关任务： 编程实现：输入一个整数n，用range()函数产生一个从n开始的10个整数组成的列表listx;将列表listx中的每个元素的值乘以2，形成一个新的列表listy。 编程要求 输入一个整数n，用range()函数产生一个从n开始的10个整数组成的列表listx; 将列表listx中的每个元素的值乘以2，形成一个新的列表listy。 根据提示，在右侧编辑器补充代码 。 测试说明 平台会对你编写的代码进行测试： 测试输入：1； 预期输出： [1, 2, 3, 4,
"控制工程数学学习支持及评估在英国谢菲尔德大学的应用
第九届国际会计师联合会控制教育进展国际自动控制联合会，俄罗斯下诺夫哥罗德，2012年帮助控制工程师学习数学J.A. 罗西特公司帕特尔英国谢菲尔德大学自动控制与系统工程系，地址：Mappin Street，Sheffield，S。约克...
二元混合物的量子化学统计热力学计算软件Replemaker 2的更新版本及应用
CA 92697，USAb波恩大学，Mulliken理论化学中心，53115 Bonn，Germanyar t i cl e i nf o文章历史记录：2018年9月3日收到收到修订版，2018年10月30日接受，2018年关键词：分子团簇二元混合物液体量子化学统计热力学a...
"乳腺癌组织学图像分类的迁移学习比较研究
⃝⃝可在www.sciencedirect.com在线ScienceDirectICT Express 4（2018）247www.elsevier.com/locate/icte乳腺癌组织学图像分类：从头开始训练还是迁移学习？拉杰什？梅赫拉？沙卢？印度昌迪加尔国家技术教师培训研究...
"乳腺癌组织学图像分类：迁移学习与传统训练网络比较
⃝⃝可在www.sciencedirect.com上在线ScienceDirectICT Express 4（2018）247www.elsevier.com/locate/icte乳腺癌组织学图像分类：从头开始训练还是迁移学习？拉杰什？梅赫拉？沙卢？印度昌迪加尔国家技术教师培训...
python生成一个1到n的队列
热门推荐
welcome to daijiguo's blog
08-05 5万+
n = 5 nlist = range(1,n) print nlist 结果： [1, 2, 3, 4]
通过range函数快速产生一个list
jason的笔记
06-24 4207
Python中range()函数用法
qq_38737191的博客
12-07 4198
Python中range()函数的用法 函数原型：range(start, stop, step) 这里的start和step都是可以缺省的，start缺省状态下默认取值为0，step缺省状态下默认取值为1 几种常见的用法：print(list(range(0,5)) 输出为 [0,1,2,3,4] print(list(range(5)) 输出为[0,1,2,3,4] 这里的range()函数所返回的对象在python3中并不是list类型的，也不继承list类型。 python3中的range
python的range函数
狂想曲
07-22 3205
python range() 可以用于for循环 ，和列表推导式。 使用方式: range(start, stop[, step]) start: 计数从 start 开始。默认是从 0 开始。range(4) 相当于 range(0,4) stop: 计数到 stop 结束，不包括 stop。range(0,4) 表达 [0, 1, 2, 3] step：步长，默认为1。range(0,3) 相当于 range(0, 3, 1) 得到从2开始，小于100 的递增12的列表： t=list(.
Python range() 函数用法及易错点
08-19 1万+
python range() 函数可创建一个整数列表，一般用在 for 循环中。 函数语法 range(start, stop[, step]) 参数说明： start: 计数从 start 开始。默认是从 0 开始。例如range（5）等价于range（0， 5）; stop: 计数到 stop 结束，但不包括 stop。例如：range（0， 5） 是[0, 1, 2, 3, 4]...
Python的range（n）的用法
weixin_30954607的博客
11-08 8329
Python的range（n） 方法就是： API定义：　　If you do need to iterate（迭代） over a sequence（一系列） of numbers, the built-in function range() comes in handy（方便的）.It generates arithmetic progressions 　　如果确实需要迭代一组数字，那么...
初学python,range()函数的用法
liufang12的专栏
09-08 1186
range()函数的用法： http://www.cnblogs.com/buro79xxd/archive/2011/05/23/2054493.html 大神用python写的打印三角形 n = 5 for i in range(n, 0, -1):     for s in range(i, 0, -1):         print ((n-s)*' '+(2*s -1)*
【Python】一元线性回归的分析
计算机小白努力学习
06-18 1702
在科学实验中，人们研究各物理量之间的关系，通常是通过实验测得一批实验数据，然后求得自变量 x 与因变量y之间的近似函数关系 y=f(x) , 这就是按给定的 n 个点的数据拟合的问题。 线性回归（最小二乘方法） 1、今有实验数据：(共有N= 6组数据 ) | 序号 | 0 |1 | 2 |3 |4 |5| |Xi |1 |2 |3 4 5 6 Yi 2.5 3.51 4.45 5.52 6.47 7.51 序号 0 1 2 3 4 5 Xi 1 2 3 4 5 6 Yi 2.5 3.
python将多重列表中的成绩按大小输出_零基础入门Python：基本命令、函数、数据结构...
weixin_39542936的博客
10-24 288
导读：本文对Python的基本使用做一个简单的介绍。限于篇幅，本文不可能详细讲解Python的使用，只是针对本书涉及的数据挖掘案例所用到的代码进行基本讲解。如果读者是初步接触Python，并且使用Python的目的就是数据挖掘，那么相信本文的介绍对你来说是比较充足的了。作者：张良均 谭立云 刘名军 江建明来源：大数据DT(ID：hzdashuju)01 运行方式本文示例代码使用的Pyth...
如何从头学习python
最新发布
12-28
从头学习Python可以通过以下步骤进行： 1. 学习基本概念和语法：首先，你需要学习Python的基本概念和语法。这包括变量、数据类型、运算符、条件语句、循环语句等。你可以通过阅读Python的官方文档、参考书籍或在线...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值