《C语言及程序设计》实践参考——sin泰勒展式中的错误

返回:贺老师课程教学链接  项目要求

【项目1-sin泰勒展式中的错误】

下面是sin函数的泰勒展式:
(注:x取弧度值,而非角度值)
编写了double mysin(double x)用于求sin值,却“死”在了123°上(122°度的结果已经出来了)。剧透一下,循环没有问题(当然问题会表现在循环中)。试着用调试工具找出问题出现在哪里,然后给出解决问题的方案。

#include<stdio.h>
#define pi 3.1415926
double mysin(double x);
double myabs(double x);
int main( )
{
    double angle;
    for(angle=0; angle<=180; angle++)
        printf("sin(%.0f°) = %.3f\n", angle, mysin((angle/180)*pi));
    return 0;
}

//下面定义mysin函数,求sin值
double mysin(double x)
{
    double sum=x,x_pow=x,item;
    int n=1,fact=1, sign=1;     //定义变量时赋初值,已经将第一项考虑到累加和sum中
    do
    {
        fact=fact*(n+1)*(n+2);  //fact用于表示阶乘,在公式中作分母
        x_pow*=x*x;             //x_pow是分子中用于表示阶乘,在公式中作分母
        sign=-sign;             //确定即将要累加的这一项的符号
        item =x_pow/fact*sign; //计算出要累加的项
        sum+=item;              //将该项累加上去
        n+=2;
    }while(myabs(item)>1e-5);
    return sum;
}

//下面定义myabs函数
double myabs(double x)
{
    return ((x>=0)?x:-x);
}

提示:请进入到mysin中后,注意各变量的变化,看通项是否会收敛,从而使循环能够结束。
[参考解答]

(若需要参考,下面的锦囊逐个找开。你要是将所有锦囊全看了再干,……老贺会伤心的:每个锦囊里都有心血,一个一个做出来不容易。)

(若你在实施中,还有意外阻碍了你,请在评论中说明,帮助老贺细化锦囊。)

锦囊1:跟踪要进到mysin函数中,注意用step into。

锦囊2:跟踪mysin函数的执行,离不了进循环,你要是一直用next line(用step into不遇到自定义函数时效果也一样),点鼠标很单调,还容易分散注意力。请在循环中某语句上设置断点,用Debug/Continue按钮“跨越式”跟踪。

锦囊3:因为问题出在123°,你要是从0°开始跟踪mysin中的循环,我相信你看到这个锦囊时,大概angle不超过10(这已经说明你有足够的耐心了)。我们需要直接进入到对当角度是123°时对mysin的调用。有两种方法:

第一种:改一下main函数,例如(还可以有很多方式,只要能直接调用mysin(123°)即可):

int main( )
{
    double angle;
    printf(" %.3f\n", mysin((123/180)*pi));
    //for(angle=0; angle<=180; angle++)
    //    printf("sin(%.0f°) = %.3f\n", angle, mysin((angle/180)*pi));
    return 0;
}

第二种:在观察窗口(Watches)中,还可以在跟踪中改变变量的值,以便看到对应的执行结果。所以如下图,可以在进入到mysin函数之前,在观察窗口中加入angle(需要在系统自动出现的变量的下面,自己再输入变量名),然后在后面直接将想要的值输入。

其实还可以step into到mysin中后(或者通过断点直接进入到mysin中后),同样的办法修改x的值。

锦囊4:很可能你听了老贺的指点,用锦囊3中的第一种方法,却发现进到mysin后,x的值是0,不是123°对应的弧度值。
对不起,我故意挖个坑,你也就进来了。进来了,坐会儿再走。
求123°的sin值,不是调用mysin((123/180)*pi),而是mysin((123.0/180)*pi)!mysin((123/180)*pi)的确就是mysin(0),注意数据类型。
进了这个坑,就要知道,程序的修改,可能会引入新的错误。这是软件工程中的一个规律。

锦囊5:现在进入关键时刻(建议还是设好断点跟踪),你会发现问题在于,item的绝对值本来是逐渐递减趋于0的(这体现了泰勒公式的收敛性,实际上,当角度值没有达到123°时,这种收敛是能保证的),但是,某个时候,item的绝对值却又大了起来,退出循环的希望,逐渐渺茫,以致于失去了希望……你可以观察fact的值,这里也发生了一些似乎不可思议的事情。

锦囊到此,现象都有了,下面就需要你的诊断了。请试着解决这个问题。解决了,或者实在想不出来了,再看下面的“真相”。

真相(倒着看,一来你得活动活动身体了,二来,实在不想让你很容易地放弃自己给出解决方案的探索历程):

时间: 2024-11-10 10:52:18

《C语言及程序设计》实践参考——sin泰勒展式中的错误的相关文章

《C++语言基础》实践参考—— 链表类

返回:贺老师课程教学链接  项目要求 [项目 - 链表类]动态链表也是程序设计中的一种非常有用的数据结构.可以说,是否能够理解有关操作的原理,决定了你是否有资格称为"科班"出身.在后续的专业基础课中,相关的内容还会从不同的角度,反复地认识,反复地实践.不过,在现阶段多些体验,也是很有必要的了.(1)阅读下面的程序,回顾一下动态链表,阅读程序过程中,请用笔画一画形成链表的过程中指针值的变化. #include <iostream> using namespace std; s

《C++语言基础》实践参考——旱冰场造价

返回:贺老师课程教学链接 [项目1 - 旱冰场造价] 有一个圆形的旱冰场地,场地内抹水泥,造价为每平方米M元,围栏用木条围成,每米造价N元,输入旱冰场半径,用程序计算出其造价. 一级提示: 先定义类,确定其数据成员和成员函数.在这里要定义的旱冰场类,实际我们只关心其面积和周长,根据面积计算场地内抹水泥的造价,根据周长计算围栏的造价.有了类的定义,在main函数中定义该类对象,调用成员函数即可以完成求解.二级提示: 实际上,对于旱冰场,我们只关心它就是一个圆.这就是设计中抓住了本质.于是,设计一个

《C++语言基础》实践参考——指向学生类的指针

返回:贺老师课程教学链接 [项目4-指向学生类的指针] 设计一个学生类Student,数据成员包括学号(num)和成绩(score),成员函数根据需要自行设计(建议配备需要的set.get函数,以及必要的输入或输出,给出的代码中也可以找到需要成员函数的线索).在main函数中,要做到: 建立一个对象数组,通过初始化,设置5个学生的数据,要求: 用指针指向数组首元素,输出第1.3.5个学生的信息: 设计一个函数int max(Student *arr);,用指向对象的指针作函数参数,在max函数中

《C++语言基础》实践参考——复数模板类

返回:贺老师课程教学链接 [项目6-复数模板类]    阅读教材例10.1.该例实现了一个复数类,但是美中不足的是,复数类的实部和虚部都固定只能是double型的.可以通过模板类的技术手段,设计Complex,使实部和虚部的类型为定义对象时指定的实际类型.    (1)要求类成员函数在类外定义.    (2)在此基础上,再实现减法.乘法和除法    你可以使用的main()函数如下. int main( ) { Complex<int> c1(3,4),c2(5,-10),c3; //实部和虚

《C++语言基础》实践参考——Josephus(约瑟夫环)问题

返回:贺老师课程教学链接  项目要求 [项目-Josephus(约瑟夫环)问题]n个小孩子围成一圈,从第一个小孩子开始顺时针方向数数字,到第m个小孩子离开,这样反反复复,最终只剩下一个小孩子,求第几个小孩子留下?    提示:约瑟夫环即是一个首尾相连的链表,在建立好这个环以后,从头结点开始,每次间隔m孩子删除一个结点,直至只余下一个结点(删除了n-1个).     参考下面的代码,也可以自行设计类. //链表结点kid,其中number为这个人的编号 struct kid { int numbe

《C++语言基础》实践参考——友元类

返回:贺老师课程教学链接 [项目5-友元类]定义下面两个类的成员函数(为体验友元类,实际上本例并不一定是一个好的设计,将两个类的合并为一个DateTime,日期.时间都处理更好) class Date; //对Date类的提前引用声明 class Time { public: Time(int,int,int); void add_a_second(Date &); //增加1秒,1秒后可能会到了下一天,乃到下一月.下一年 void display(Date &); //显示时间,格式:月

《C++语言基础》实践参考——分数类中的运算符重载

返回:贺老师课程教学链接 [项目3-分数类中的运算符重载] (1)实现分数类中的运算符重载,在分数类中可以完成分数的加减乘除(运算后再化简).比较(6种关系)的运算.可以在第4周分数类代码的基础上开始工作. class CFraction { private: int nume; // 分子 int deno; // 分母 public: //构造函数及运算符重载的函数声明 }; //重载函数的实现及用于测试的main()函数 [参考解答] #include <iostream> #inclu

《C++语言基础》实践参考——点、圆的关系

返回:贺老师课程教学链接  项目要求 [项目4 - 点.圆的关系](1)先建立一个Point(点)类,包含数据成员x,y(坐标点):(2)以Point为基类,派生出一个Circle(圆)类,增加数据成员(半径),基类的成员表示圆心:(3)编写上述两类中的构造.析构函数及必要运算符重载函数(本项目主要是输入输出):(4)定义友元函数int locate,判断点p与圆的位置关系(返回值<0圆内,==0圆上,>0 圆外): int main( ) { Circle c1(3,2,4),c2(4,5,

《C++语言基础》实践参考——三角形类1

返回:贺老师课程教学链接 [项目1 - 三角形类1]下面设计一个三角形类,请给出各成员函数的定义  #include<iostream> #include<Cmath> using namespace std; class Triangle { public: void setABC(double x, double y, double z);//置三边的值,注意要能成三角形 void getABC(double *x, double *y, double *z);//取三边的值