Java 通过sin(x)和cos(x)幂级数 编程计算其结果

最新推荐文章于 2024-04-14 17:32:12 发布
最新推荐文章于 2024-04-14 17:32:12 发布
阅读量688 收藏
点赞数
分类专栏： 笔记 Algorithm Java
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xiuye2015/article/details/110441106
版权
Java 同时被 3 个专栏收录
297 篇文章 1 订阅
订阅专栏
笔记
124 篇文章 0 订阅
订阅专栏
Algorithm
111 篇文章 0 订阅
订阅专栏 
package com.xiuye.util.algorithm.basics;

import com.xiuye.util.X;

import java.math.BigDecimal;
import java.math.BigInteger;

public class Sin {
    /**
     * y = 0.987862x - 0.155271x^3 + 0.00564312x^5
     *
     * @param c
     * @return
     */
    public static double sin1(double cn) {
        X.lg("sin1:", cn);
        double a = 0.987862;
        double b = 0.155271;
        double c = 0.00564312;


        return a * cn - b * cn * cn * cn + c * cn * cn * cn * cn * cn;
    }


    //sin(x) 幂级数
    public static double sin2(double x) {
        //把大于 -2PI 和 2PI 的值放入区间内

        int rSign = 1;
        if(x<0){
            rSign = -rSign;
            x = -x;
        }
        // >2PI
        x %= 2*Math.PI;
        // >

        X.lg("sin2:", x);
        int cnt = 0;
        double derr = 1e-15;

        long exp = 1;
        int sign = 1;
        double r = x;

        for (long n = 1; ; n++) {

            sign = -sign;//-1 or +1
            double xm = x;
            long m = 2 * n + 1;
            exp *= m * (m - 1);
            xm = Math.pow(x, m);
//            X.lg("m:",m);
//            X.lg("exp:",exp);
//            X.lg("xm:",xm);
            if (exp < 0 || xm / exp < derr) {
                break;
            }
            r += sign * xm / exp;
            cnt++;


        }
        X.lg("sin2.cnt:", cnt);
        return rSign*r;

    }

    //cos(x) 幂级数
    public static double cos2(double x) {
        //把大于 -2PI 和 2PI 的值放入区间内

        int rSign = 1;
        if(x<0){
//            rSign = -rSign;
            x = -x;
        }
        // >2PI
        x %= 2*Math.PI;
        // >

        X.lg("cos2:", x);
        int cnt = 0;
        double derr = 1e-15;

        long exp = 1;
        int sign = 1;
        double r = 1;

        for (long n = 1; ; n++) {

            sign = -sign;//-1 or +1
            double xm = x;
            long m = 2 * n;
            exp *= m * (m - 1);
            xm = Math.pow(x, m);
//            X.lg("m:",m);
//            X.lg("exp:",exp);
//            X.lg("xm:",xm);
            if (exp < 0 || xm / exp < derr) {
                break;
            }
            r += sign * xm / exp;
            cnt++;


        }
        X.lg("cos2.cnt:", cnt);
        return rSign*r;

    }




    public static void main(String[] args) {
        X.lg("PI/2", sin1(Math.PI / 2));
        X.lg("PI/6", sin1(Math.PI / 6));
        X.lg("PI/4", sin1(Math.PI / 4));
        X.lg(1 / Math.sqrt(2));
        X.lg("PI/2", sin2(Math.PI / 2));
        X.lg("PI/2+2PI", sin2(Math.PI / 2 + Math.PI * 2));//有精度损失
        X.lg("PI/6", sin2(Math.PI / 6));
        X.lg("PI/6+2PI", sin2(Math.PI / 6 + Math.PI * 2));//有精度损失
        X.lg("PI/4", sin2(Math.PI / 4));
        X.lg("PI/4+2PI", sin2(Math.PI / 4 + Math.PI * 2));//有精度损失
        X.lg("-PI/6", sin2(-Math.PI / 6 + Math.PI * 2));//有精度损失
        X.lg("-PI/6+2PI", sin2(-Math.PI / 6 + Math.PI * 2));//有精度损失
        X.lg("2PI", sin2(Math.PI * 2));//有精度损失
        X.lg("-2PI", sin2(Math.PI * 2));//有精度损失

        X.lg(11.0 % 88);
        X.lg(11.0 / 88);
        X.lg(11.0 / 88 * 88);
        X.lg(Math.PI / Math.PI);
        X.lg(2 * Math.PI / (2 * Math.PI));
        X.lg(-Math.PI / Math.PI);
        X.lg(-2 * Math.PI / (2 * Math.PI));

        X.lg("PI/2",cos2(Math.PI/2));
        X.lg("PI",cos2(Math.PI));
        X.lg("-PI",cos2(-Math.PI));
        X.lg("PI/6",cos2(Math.PI/6));
        X.lg("PI/3",cos2(Math.PI/3));
        X.lg("PI/4",cos2(Math.PI/4));
        X.lg("2PI",cos2(2*Math.PI));
        X.lg("0",cos2(0));
//        BigInteger bigInteger = BigInteger.valueOf(2432902008176640000L);
//        X.lg(bigInteger.multiply(BigInteger.valueOf(2432902008176640000L)));
//        X.lg(bigInteger.multiply(BigInteger.valueOf(2432902008176640000L)).doubleValue());

    }
}


sin1: 1.5707963267948966
PI/2 1.0038988590395785
sin1: 0.5235987755982988
PI/6 0.4951766393046648
sin1: 0.7853981633974483
PI/4 0.7023268134010443
0.7071067811865475
sin2: 1.5707963267948966
sin2.cnt: 9
PI/2 1.0
sin2: 1.5707963267948966
sin2.cnt: 9
PI/2+2PI 1.0
sin2: 0.5235987755982988
sin2.cnt: 6
PI/6 0.5
sin2: 0.5235987755982991
sin2.cnt: 6
PI/6+2PI 0.5000000000000003
sin2: 0.7853981633974483
sin2.cnt: 7
PI/4 0.7071067811865475
sin2: 0.7853981633974483
sin2.cnt: 7
PI/4+2PI 0.7071067811865475
sin2: 5.759586531581287
sin2.cnt: 9
-PI/6 -0.5001706202630702
sin2: 5.759586531581287
sin2.cnt: 9
-PI/6+2PI -0.5001706202630702
sin2: 0.0
sin2.cnt: 0
2PI 0.0
sin2: 0.0
sin2.cnt: 0
-2PI 0.0
11.0
0.125
11.0
1.0
1.0
-1.0
-1.0
cos2: 1.5707963267948966
cos2.cnt: 10
PI/2 6.08176268469988E-17
cos2: 3.141592653589793
cos2.cnt: 10
PI -0.9999999999243493
cos2: 3.141592653589793
cos2.cnt: 10
-PI -0.9999999999243493
cos2: 0.5235987755982988
cos2.cnt: 7
PI/6 0.8660254037844386
cos2: 1.0471975511965976
cos2.cnt: 8
PI/3 0.5000000000000004
cos2: 0.7853981633974483
cos2.cnt: 8
PI/4 0.7071067811865475
cos2: 0.0
cos2.cnt: 0
2PI 1.0
cos2: 0.0
cos2.cnt: 0
0 1.0



精度有问题 用BigDecimal解决

 

__BMGT()
关注
专栏目录
Maple拥有优秀的符号计算和数值计算能力
weixin_34226706的博客
03-26 7852
    https://www.maplesoft.com/products/maple/     Maple高级应用和经典实例:　　https://wenku.baidu.com/view/f24696210722192e4536f65d.html Maple高级应用和经典实例: http://vdisk.weibo.com/s/dbLrQxb6KthZA     Maple是目前世界上最为通用...
c语言实现 三角函数,关于数学：快速实现C ++三角函数
weixin_28683689的博客
05-23 2410
简短版本：我想知道是否有标准三角函数的实现比math.h中包含的函数更快。较长的版本：我有一个程序，它非常注重数值(这是物理模拟)，并且需要调用三角函数，主要是sin和cos。 当前，我只是在使用math.h中包含的实现。 分析显示，对这些函数的调用花费比我期望的要多(希望)。虽然在代码的其他部分中肯定有很多优化的余地，但是更快的sin和cos可能会给我带来一些额外的收益。.那么，你们有什么建议吗...
求sinx和cosx的近似值代码
12-07
利用泰勒展式，求sinx和cosx的近似值
java语言判定sinx的近似值
mengxin2123的博客
03-06 1882
import java.util.Scanner; public class Main { //main函数,程序入口 public static void main(String[] args) { //定义变量sum用于保存结果 double sum=0.0,x; //定义变量n用于循环,flag置正负标志位 int n,flag=1; //输出提示 System.out.prin...
java三角函数
牢骚太盛防肠断，风物长宜放眼量
07-30 2149
package com.xiaoz.classTest; public class TestAngle { public static void main(String[] args) { /* sin()、cos()、tan() 传入的是角度(角度是弧度制),输出的是正弦值、余弦值、正切值 弧度制 = (PI/180) * 度数 */ System.out.println(Math.sin(Math.
Java输入小数和整数求COSx_Java SE考试--编程题总结（优秀范文五篇）
weixin_34558565的博客
02-27 393
线程和进程的区别：(1)进程是运行中的程序，拥有自己独立的内存空间和资源；(2)一个进程可以有一个或多个线程组成，且至少有一个线程称为主线程； (3)线程是最小的处理单位，多个线程共享一块内存和资源；(4)当一个线程改变了所属进程的变量时，其它线程下次访问该变量时得到这种改变。 简述什么是多线程。线程的最主要功能是多任务处理，即多线程。多线程也就是在主线程中有多个线程在运行，多个线程的执行是并发的...
java中如何表达sin_在Java中使用“sincos”
weixin_33539366的博客
02-27 2196
在很多情况下,我不仅需要正弦,还要求相同参数的余弦.对于C,常见的unix m数学库中存在sincos函数.实际上,至少在i386上,这应该是一个单一的汇编指令fsincos.sincos, sincosf, sincosl – calculate sin and cos simultaneously我猜这些好处是存在的,因为计算正弦和余弦有一个明显的重叠：sin(x)^ 2 cos(x)^ 2 ...
程序设计语言（LOGO/Scratch/VB/C）
最新发布
weixin_44525733的博客
04-14 1281
（3）某城市的出租车计费标准为起步价10元(3公里内)，此外，在3公里到10公里之间每公里 2.1元，超过 10公里部分每公里3元，输入行车距离x，输出车费y。用100元钱买100只鸡，公鸡每只5元，母鸡每只3元，小鸡3只1元，要求每种都买，问能买公鸡、母鸡、小鸡各几只?Int(Rnd(0) * 10)： [0,1)→[0,10)，题目可能会问正整数的范围就要注意了！Rnd(x)：生成一个[0, 1)之间的单精度随机数，其中(x)可选。（2）输人学生成绩(百分制)，判断该成绩的等级(优良、及格、不及格)
maplesim matlab,maple
weixin_28197941的博客
04-08 1226
maple(工程计算软件)编辑锁定讨论Maple是目前世界上最为通用的数学和工程计算软件之一，在数学和科学领域享有盛誉，有“数学家的软件”之称。Maple 在全球拥有数百万用户，被广泛地应用于科学、工程和教育等领域，用户渗透超过96%的世界主要高校和研究所，超过81%的世界财富五百强企业。Maple系统内置高级技术解决建模和仿真中的数学问题，包括世界上最强大的符号计算、无限精度数值计算、创新的互联...
big-math:使用任意精度的高级Java BigDecimal数学函数（pow，sqrt，log，sin等）
04-18
使用任意精度的高级Java BigDecimal数学函数（ pow ， sqrt ， log ， sin ，...）。 另请参阅官方的。 大十进制数学 BigDecimalMath类为以下各项提供了有效且准确的实现： log(BigDecimal, MathContext) exp(BigDecimal, MathContext) pow(BigDecimal, BigDecimal, MathContext)计算x ^ y sqrt(BigDecimal, MathContext) root(BigDecimal, BigDecimal, MathContext)计算x的第n个根 sin(BigDecimal, MathContext) cos(BigDecimal, MathContext) tan(BigDecimal, MathContext) asin(BigD
BigDecimal运算
10-19
BigDecimal运算Jar包。支持sin、cos、tan、asin、acos、atan、pow、sqrt、cbrt、root、log10、log、ln、exp、sinh、cosh、tanh、asinh、acosh、atanh、deg、rad
java中如何表达sin,使用“sincos”和“sincos”。在Java中
weixin_33321992的博客
02-27 733
In a lot of situations I not only need the sine, but also the cosine of the same parameter.For C, there is the sincos function in the common unix m math library. And actually, at least on i386, this s...
用sin 和cos 的泰勒展开式编程序，求出sin( π/2)、sin(56°)、cos(87°)、cos( π/3)
u012332009的专栏
11-13 2828
/* *Corpyright (c)2013,烟台大学计算机学院 *All right reseved. *作者：高古尊 *完成日期：2013年11月9日 *版本号：v1.0 *输入描述： *问题描述：用sin 和cos 的泰勒展开式编程序， 求出sin( π/2)、sin(56°)、cos(87°)、cos( π/3) *程序输出： *问题分析： *算法设计： *
深入学习java源码之Math.sin()与 Math.sqrt()
悦分享
01-05 2268
深入学习java源码之Math.sin()与 Math.sqrt() native关键字 凡是一种语言，都希望是纯。比如解决某一个方案都喜欢就单单这个语言来写即可。Java平台有个用户和本地C代码进行互操作的API，称为JNI native关键字告诉编译器（其实是JVM）调用的是该方法在外部定义，这里指的是C。 Modifier and Type Method and Des...
Java 使用BigDecimal，正弦函数sin(x)级数计算其值
修也
12-02 1255
private static BigDecimal bigDecimal(String val){ return new BigDecimal(val); } public static BigDecimal sin3(BigDecimal x) { //把大于 -2PI 和 2PI 的值放入区间内 BigDecimal rSign = bigDecimal("1"); // int rSign = 1; i.
java 级数_Java练习： 数列/级数求和
weixin_28705087的博客
02-15 557
级数通常指无穷级数，而有穷级数通常称为数列。数列/级数求和或求近似值是常见的循环编程练习。下面将给出一些题目。1.调和数调和级数Harmonic numbers，H(n)= 1/1 + 1/2 + 1/3 + 1/4 + ... + 1/n。求调和级数前n项的和/*** 调和级数Harmonic numbers,* H(n)= 1/1 + 1/2 + 1/3 + 1/4 + ... + 1/n*/...
求导程序编写(oo-java编程)
diding0909的博客
03-23 710
本单元的任务为求导。 即将一个含自变量x的多项式F求导成为另外一个含自变量x的多项式f。使得 dF/dx = f 为降低我们的难度，这个任务被分解成了三个阶段： （1）对幂函数进行求导（不允许嵌套） （2）对幂函数和三角函数进行求导（不允许嵌套，三角函数中只能有x） （3）对幂函数和三角函数进行求导（允许嵌套，三角函数中只能有因子） 一、聊聊思路 1、字符串...
Java 大数BigDecimal常用处理
富士山下
04-28 2231
import java.math.BigDecimal; import java.util.Scanner; public class Main { public static void main(String[] args) { Scanner cin=new Scanner(System.in); while(cin.hasNext())
C/C++编程实践：绘制cos(x)和sin(x)图形
这两个实例都是C语言实践教学的重要组成部分，它们不仅涉及基本的控制结构（如for循环）、数学函数的使用（如acos()），还展示了如何通过控制字符输出创建图形，这对于初学者掌握C语言的控制流和数值计算非常有帮助...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值