JavaScript实现计算圆周率到小数点后100位的方法示例
作者:ncs2002 发布时间:2024-02-27 02:38:58
标签:JavaScript,计算,圆周率
本文实例讲述了JavaScript实现计算圆周率到小数点后100位的方法。分享给大家供大家参考,具体如下:
浮点数的有效数位是16位,我自己做了一个大数类,能存储100位有效数位,并实现了大数类的基本运算。我用它来计算圆周率(割圆法,即多边形逼近),得到了小数点后一百位有效数字,比对了Machin 公式的计算结果,没有误差。用时约2秒。
完整示例如下:
<!DOCTYPE html>
<html>
<head>
<meta charset="utf-8">
<title>www.aspxhome.com js计算圆周率</title>
</head>
<body>
<script>
<!--
function BigNum(str, n, b)
{
/*
BigNum -- 大数类
私有成员:
data -- 119 位数字,放在长度为 17 的数组里,每个数组元素存放 7 位数字。
decimal_place -- 小数点的位置,从最左位开始算。
positive -- 是否是正数。
recalc() -- 为了尽可能存放最多的有效数位,去除前缀的 0,并重新计算小数点位置。
init() -- 部分初始化工作。
公有成员:
BigNum( String str, INT n, BOOL b) --
构造函数。参数:str -- 字符串,各个有效数位;n -- 整数,小数点位置,从最左位开始算,比如 BigNum("123", 2) = 12.3; BigNum("123", 0) = 0.123; BigNum("123", -2) = 0.00123;b -- 布尔值,是否是正数。
Add( BigNum num ) -- 加法。
Subtract( BigNum num ) -- 减法。
Multiply( BigNum num ) -- 乘法。
Divide( BigNum num ) -- 除法。
SquareRoot() -- 平方根。
toString() -- 转换为字符串(包括小数点),以便以文本形式输出计算结果。
Clone() -- 复制。
*/
this.recalc = function() /* 去除前缀的 0,并重新计算小数点位置 */
{
for(var i = 0; i < 17; i ++)
{
if(this.data[0] != 0) break;
this.data.shift();
this.data.push(0);
this.decimal_place --;
}
}
this.init = function() /* 部分初始化工作 */
{
this.decimal_place = Math.ceil( n / 7 ); //小数点位置
this.data = new Array(17); //保存有效数位的数组
if(n % 7 > 0)
{
var arr = new Array( 8 - n % 7 );
}
else
{
var arr = new Array( 1 - n % 7 );
}
str = arr.join("0") + str;
if(str.length > 119)
{
str = str.substr(0, 119);
}
else if(str.length < 119)
{
var arr = new Array(120 - str.length);
str += arr.join("0");
}
for( var i = 0; i < 17; i ++ )
{
this.data[i] = parseInt( str.substr(i * 7, 7) , 10 );
}
}
/* 初始化开始 */
this.positive = b;
if( ! /^0*$/.test(str) )
{
this.init();
this.recalc();
}
else
{
this.data = new Array(17);
for( var i = 0; i < 17; i ++ )
{
this.data[i] = 0;
}
this.decimal_place = 0;
}
/* 初始化结束 */
this.Add = function(num) /* 加法 */
{
if(this.positive && !num.positive)
{
num.positive = true;
var result = this.Subtract(num);
num.positive = false;
return result;
}
else if(num.positive && !this.positive)
{
this.positive = true;
var result = num.Subtract(this);
this.positive = false;
return result;
}
var result = new BigNum("", 0, this.positive);
var num1,num2;
if(this.decimal_place >= num.decimal_place)
{
num1 = this;
num2 = num;
}
else
{
num1 = num;
num2 = this;
}
result.decimal_place = num1.decimal_place;
if(num1.decimal_place - num2.decimal_place >= 17)
{
for(var i = 0; i < 17; i ++)
{
result.data[i] = num1.data[i];
}
return result;
}
var nOffDec = num1.decimal_place - num2.decimal_place;
var nTmp = 0;
for( var i = 16; i >= 0; i -- )
{
var nTmp1 = i - nOffDec;
var nTmp2 = 0;
if(nTmp1 >= 0)
{
nTmp2 = num1.data[i] + num2.data[nTmp1];
}
else
{
nTmp2 = num1.data[i];
}
nTmp2 += nTmp;
nTmp = Math.floor(nTmp2 / 10000000);
result.data[i] = nTmp2 % 10000000;
}
if(nTmp > 0)
{
result.data.unshift(nTmp);
result.decimal_place ++;
}
return result;
}
this.Subtract = function(num) /* 减法 */
{
if(this.positive && !num.positive)
{
num.positive = true;
var result = this.Add(num);
num.positive = false;
return result;
}
else if(!this.positive && num.positive)
{
this.positive = true;
var result = this.Add(num);
result.positive = false;
this.positive = false;
return result;
}
else
{
var num1 = num2 = null;
var bPositive;
if(this.decimal_place > num.decimal_place)
{
num1 = this;
num2 = num;
bPositive = this.positive;
}
else if(this.decimal_place < num.decimal_place)
{
num1 = num;
num2 = this;
bPositive = !this.positive;
}
else
{
for( var i = 0; i < 17; i ++ )
{
if(this.data[i] > num.data[i])
{
num1 = this;
num2 = num;
bPositive = this.positive;
break;
}
else if(this.data[i] < num.data[i])
{
num1 = num;
num2 = this;
bPositive = !this.positive;
break;
}
}
}
if( num1 == null)
{
return new BigNum("", 0, true);
}
else
{
if(num1.decimal_place - num2.decimal_place >= 17)
{
var result = new BigNum("", 0, bPositive);
for(var i = 0; i < 17; i ++)
{
result.data[i] = num1.data[i];
}
result.decimal_place = num1.decimal_place;
return result;
}
var result = new BigNum("", 0, bPositive);
result.decimal_place = num1.decimal_place;
var nOffDec = num1.decimal_place - num2.decimal_place;
var nTmp = 0;
for( var i = 16; i >= 0; i -- )
{
var nTmp1 = i - nOffDec;
var nTmp2 = 0;
if(nTmp1 >= 0)
{
nTmp2 = 10000000 + nTmp + num1.data[i] - num2.data[nTmp1];
}
else
{
nTmp2 = 10000000 + nTmp + num1.data[i];
}
if(nTmp2 >= 10000000)
{
result.data[i] = nTmp2 - 10000000;
nTmp = 0;
}
else
{
result.data[i] = nTmp2;
nTmp = -1;
}
}
result.recalc();
return result;
}
}
}
this.Multiply = function(num) /* 乘法 */
{
var bPositive;
var nDecimalPlace = this.decimal_place + num.decimal_place - 1;
if(this.positive == num.positive)
{
bPositive = true;
}
else
{
bPositive = false;
}
var result = new BigNum("", 0, bPositive);
var nTmpData = 0;
for( var i = 16; i >= 0; i -- )
{
for( var j = 16; j >= 0; j -- )
{
if(isNaN(result.data[j + i]))
result.data[j + i] = 0;
result.data[j + i] += this.data[j] * num.data[i];
if(result.data[j + i] >= 10000000)
{
if( j + i -1 >= 0 )
{
result.data[j + i -1] += Math.floor(result.data[j + i] / 10000000);
}
else
{
nTmpData += Math.floor(result.data[j + i] / 10000000);
}
result.data[j + i] = result.data[j + i] % 10000000;
}
}
}
if(nTmpData > 0)
{
result.data.unshift(nTmpData);
result.data.pop();
nDecimalPlace ++;
}
result.decimal_place += nDecimalPlace;
return result;
}
this.Divide = function(num) /* 除法 */
{
var bPositive;
var nDecimalPlace = this.decimal_place - num.decimal_place + 1;
if(this.positive == num.positive)
{
bPositive = true;
}
else
{
bPositive = false;
}
var result = new BigNum("", 0, bPositive);
var arrTemp = new Array(17);
for( var i = 0; i < 17; i ++ )
{
arrTemp[i] = this.data[i];
}
var bTest = true;
var nTest = 0;
for( var i = 0; i < 17; i ++ )
{
if(bTest)
{
nTest = Math.floor( ( arrTemp[0] * 10000000 + arrTemp[1] ) / ( num.data[0] * 10000000 + num.data[1] ) );
}
else
{
bTest = true;
}
if(nTest == 0)
{
result.data[i] = 0;
arrTemp[1] += arrTemp[0] * 10000000;
arrTemp.shift();
arrTemp.push(0);
continue;
}
var arrTemp1 = new Array(17);
for( var j = 0; j < 17; j ++ )
{
arrTemp1[j] = 0;
}
for( var j = 16; j >= 0; j -- )
{
arrTemp1[j] += nTest * num.data[j];
if(arrTemp1[j] >= 10000000)
{
if(j != 0)
{
arrTemp1[j - 1] += Math.floor( arrTemp1[j] / 10000000);
arrTemp1[j] = arrTemp1[j] % 10000000;
}
}
}
for( var j = 0; j < 17; j ++ )
{
if(arrTemp[j] < arrTemp1[j])
{
bTest = false;
break;
}
else if(arrTemp[j] > arrTemp1[j])
{
break;
}
}
if(bTest)
{
result.data[i] = nTest;
for( var j = 16; j >= 0; j -- )
{
if(arrTemp[j] >= arrTemp1[j])
{
arrTemp[j] -= arrTemp1[j];
}
else
{
arrTemp[j] = 10000000 + arrTemp[j] - arrTemp1[j];
arrTemp[j - 1] --;
}
}
}
else
{
nTest --;
i --;
continue;
}
arrTemp[1] += arrTemp[0] * 10000000;
arrTemp.shift();
arrTemp.push(0);
}
result.decimal_place = nDecimalPlace;
result.recalc();
return result;
}
this.SquareRoot = function() /* 平方根 */
{
var result = new BigNum("", 0, true);
nDecimalPlace = Math.ceil(this.decimal_place / 2);
var arrTemp = new Array(17);
for(var i = 0; i < 17; i ++)
{
arrTemp[i] = this.data[i];
}
var bTest = true;
for(var i = 0; i < 17; i ++)
{
if( i == 0 )
{
if(this.decimal_place % 2 == 0)
{
var nTemp = arrTemp[0] * 10000000 + arrTemp[1];
var nTemp1 = Math.floor( Math.sqrt( nTemp ) );
var nTemp2 = nTemp - nTemp1 * nTemp1;
arrTemp[0] = 0;
arrTemp[1] = nTemp2;
arrTemp.shift();
arrTemp.push(0);
}
else
{
var nTemp1 = Math.floor( Math.sqrt( arrTemp[0] ) );
var nTemp2 = arrTemp[0] - nTemp1 * nTemp1;
arrTemp[0] = nTemp2;
}
}
else
{
if(bTest)
{
if( i == 1 )
{
nTemp1 = Math.sqrt( (arrTemp[0] * 10000000 + arrTemp[1]) + 100000000000000 * Math.pow(result.data[0], 2) ) - 10000000 * result.data[0];
nTemp1 = Math.floor(nTemp1);
}
else
{
nTemp = result.data[0] * 10000000 + result.data[1];
nTemp1 = Math.floor( ( arrTemp[0] * 10000000 + arrTemp[1] ) / ( 2 * nTemp ) );
}
}
else
{
bTest = true;
}
var arrTemp1 = new Array(17);
var nTemp3 = 0
for( var j = i - 1; j >= 0; j -- )
{
arrTemp1[j] = result.data[j] * 2 + nTemp3;
if( arrTemp1[j] >= 10000000 && j > 0 )
{
nTemp3 = 1;
arrTemp1[j] = arrTemp1[j] % 10000000;
}
else
{
nTemp3 = 0;
}
}
arrTemp1[i] = nTemp1;
nTemp3 = 0;
for( var j = i; j >= 0; j -- )
{
arrTemp1[j] = arrTemp1[j] * nTemp1 + nTemp3;
if( arrTemp1[j] >= 10000000 && j > 0 )
{
nTemp3 = Math.floor( arrTemp1[j] / 10000000 );
arrTemp1[j] = arrTemp1[j] % 10000000;
}
else
{
nTemp3 = 0;
}
}
var arrTemp2 = new Array(17);
for( var j = 0; j < 17; j ++ )
{
arrTemp2[j] = arrTemp[j];
}
for( var j = i; j >= 0; j -- )
{
if( arrTemp2[j] >= arrTemp1[j] )
{
arrTemp2[j] -= arrTemp1[j];
}
else
{
if(j > 0)
{
arrTemp2[j] = arrTemp2[j] + 10000000 - arrTemp1[j];
arrTemp2[j - 1] --;
}
else
{
bTest = false;
break;
}
}
}
if(bTest)
{
arrTemp = arrTemp2;
}
else
{
nTemp1 --;
i --;
continue;
}
}
result.data[i] = nTemp1;
arrTemp[1] += arrTemp[0] * 10000000;
arrTemp.shift();
arrTemp.push(0);
}
result.decimal_place = nDecimalPlace;
result.recalc();
return result;
}
this.toString = function() /* 转换为字符串(包括小数点),以便以文本形式输出计算结果 */
{
var szData = "";
var szOutPut = "";
this.recalc();
for( var i = 0; i < 17; i ++ )
{
var szTmpData = this.data[i].toString()
var arr = new Array(8 - szTmpData.length);
szData += arr.join("0") + szTmpData;
}
if( /^0*$/.test(szData) )
{
return "0";
}
var n = this.decimal_place * 7;
for(var i = 0; i < 7; i++)
{
if(szData.substr(i, 1) != 0) break;
n --;
}
szData = szData.replace(/^0+/g,"");
szData = szData.substr(0, 101);
szData = szData.replace(/0+$/g,"");
if( n < 1 )
{
szOutPut = szData.substr(0, 1) +
( ( szData.length > 1 ) ? "." : "" ) +
szData.substr(1) + "e" + ( n - 1 ).toString();
}
else if(n == szData.length)
{
szOutPut = szData;
}
else if(n > szData.length)
{
szOutPut = szData.substr(0, 1) + "." + szData.substr(1) + "e+" + (n - 1).toString();
}
else
{
szOutPut = szData.substr(0, n) + "." + szData.substr(n);
}
return ( this.positive ? "" : "-" ) + szOutPut;
}
this.Clone = function() /* 复制 */
{
var result = new BigNum("", 0, true);
for( var i = 0; i < 17; i ++)
{
result.data[i] = this.data[i];
}
result.decimal_place = this.decimal_place;
result.positive = this.positive;
return result;
}
}
var a = new BigNum("1", 1, true)
var count = 168;
var nTwo = new BigNum("2", 1, true);
function loop(intTmpvar,intCount)
{
if(intCount == 0) return intTmpvar;
var v1 = intTmpvar.Divide( nTwo );
var v11 = v1.Clone();
var nTemp = v1.Multiply( v11 );
var a1 = a.Clone();
a1 = a.Multiply(a1);
var nTemp1 = a1.Subtract( nTemp )
v2 = nTemp1.SquareRoot();
v3 = a.Subtract( v2 );
var v31 = v3.Clone();
var nTemp2 = v3.Multiply( v31 );
var nTemp3 = nTemp2.Add(nTemp);
v4 = nTemp3.SquareRoot();
return loop( v4 , -- intCount )
}
var a1 = a.Clone();
var nTemp = a.Multiply(a1);
var nTemp1 = nTemp.Clone();
nTemp = nTemp.Add(nTemp1);
nTemp = loop(nTemp.SquareRoot(), count);
var nFour = new BigNum("4", 1, true);
nTemp = nTemp.Multiply( nFour );
nTemp1 = new BigNum("2", 1, true);
var nTemp2 = new BigNum("2", 1, true);
for(var i = 0; i < count - 1; i ++)
{
nTemp1 = nTemp1.Multiply( nTemp2 );
}
nTemp = nTemp.Multiply( nTemp1 );
nTemp = nTemp.Divide( nTwo );
nTemp = nTemp.Divide( a );
document.write( nTemp )
//-->
</script>
</body>
</html>
运行结果:
3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679
希望本文所述对大家JavaScript程序设计有所帮助。
来源:https://blog.csdn.net/ncs2002/article/details/429821
0
投稿
猜你喜欢
- 一、存储过程存储过程(Stored Procedure)是在大型数据库系统中,一组为了完成特定功能的SQL语句集,经编译后存储在数据库中,用
- 前言使用np.savetxt()方法可以将数据保存为txt文件或者是csv文件。1 np.savetxt()存储txt文件1-1 基础参数n
- pytorch在数据中采用Dataset的数据保存方式,需要继承data.Dataset类,如果需要自己处理数据的话,需要实现两个基本方法。
- #常用的几个代码--查询临时表空间select name from v$tempfile;--查询表空间select name from v
- 前言最近学完Python,写了几个爬虫练练手,网上的教程有很多,但是有的已经不能爬了,主要是网站经常改,可是爬虫还是有通用的思路的,即下载数
- 在数据库开发方面,通过单表所表现的实现,有时候需要组合查询来找到我们需要的记录集,这时候我们就会用到连接查询。连接查询主要包括以下几个方面:
- 英文文档:classmethod(function)Return a class method for function.A class m
- 通过Python脚本批量生成插入数据的SQL语句原始SQL语句:INSERT INTO system_user (id, login_nam
- 先来看查看效果:在代码连接数据库后,并且执行三条sql后,将mysql直接重启掉,故我们的连接池连接均是不ok的,所以,它会全部删除再抓新的
- 本例子程序展示了长白山火山气体地球化学2002年观测数据中CO2和He两种气体元素深度的时间序列。程序中用到了常用的时间序列python数据
- go fmt命令主要是用来帮你格式化所写好的代码文件【很多第三方集成软件都是使用了go fmt命令】一、使用:go fmt <文件名&
- 首先,在models.py中创建UserModels类from django.db import modelsfrom django.con
- # 比较两个字符串,如果不同返回第一个不相同的位置# 如果相同返回0def cmpstr(str1, str2): &
- 1. python三维图表绘制方法简介python三维图表的绘制算是二维图表的一个进阶版本,本质上和二维图表的绘制并无差别,唯一的区别在于使
- 介绍本篇文章主要介绍如何爬取麦子学院的课程信息(本爬虫仍是单线程爬虫),在开始介绍之前,先来看看结果示意图怎么样,是不是已经跃跃欲试了?首先
- xml即可扩展标记语言,它可以用来标记数据、定义数据类型,是一种允许用户对自己的标记语言进行定义的源语言。从结构上,很像HTML超文本标记语
- 引言日常开发中,我们经常会使用到group by。亲爱的小伙伴,你是否知道group by的工作原理呢?group by和having有什么
- 从98年某月某日我的第一个个人猪页诞生.到2008年的今天.宣告了我从事设计行业整整十年.十年,很多变化,从摆弄个人猪页到现在以做设计为生,
- 一、图像直方图1.1 定义图像直方图是图像的基本属性之一,也是反映图像像素数据分布的统计学特征,其横坐标代表了图像像素点在[0,255]范围
- 存储函数也是过程式对象之一,与存储过程相似。他们都是由SQL和过程式语句组成的代码片段,并且可以从应用程序和SQL中调用。然而,他们也有一些