作者fjf1980 (Yes We Can!)
看板C_and_CPP
标题[课业] 请问各位一个国考程式题目
时间Mon Oct 26 20:01:58 2009
※ [本文转录自 Examination 看板]
作者: fjf1980 (Yes We Can!) 看板: Examination
标题: [课业] 96地特/程式语言/泰勒展开
时间: Mon Oct 26 19:33:45
4.想问的内容：
题目：f(x)=1- x/2! + x/4! -x/6!...
每项算法为si = (-1)^i * x^2i / 2i!
题目要求设计把引数x先模数到 0<= x <2*pi, 然後用级数计算f(x)相加
准确到项数|Sn|<10^-8
5.想法：下面是我在DEV-C++实际跑程式，可是不管我x换多少，
我的答案都是1.000000 = = 买尬！ 为什麽？
#include <stdio.h>
#include <stdlib.h>
#define pi 3.14159
double sn(int x)
{
double result=0.0;
double temp=0.0;
int i=0;
while(temp<0.00000001 && temp >-0.00000001)
{
temp = pow(-1,i)*pow(x*pi,2*i)/fact(2*i);
i++;
result += temp;
}
return result;
}
int fact(int x)
{
if(x==0 || x==1)
return 1;
else
return x*fact(x-1);
}
int main(int argc, char *argv[])
{
double x=1.1;
double input_x=0.0;
x = x-((int)(x/2.0)*2.0);
printf("%f\n",sn(x));
system("PAUSE");
return 0;
}
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.230.197.214
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.230.197.214
1F：→ fjf1980:请教高手们 我这样写有问题吗？ 10/26 20:02
※ 编辑: fjf1980 来自: 61.230.197.214 (10/26 20:03)
2F：推 POSIX:sn的argument 是int? main里面给double? 10/26 20:09
3F：→ fjf1980:我改double sn(double x) 结果还是一样耶 10/26 20:12
4F：推 VictorTom:实作演算法有没有错呢??i从0开始代的话, -1^0是1, x*pi 10/26 20:18
5F：→ VictorTom:只要不是0, 代入2*i也相当於是^0, fact传入0回传1, 10/26 20:18
6F：→ VictorTom:(上面的^暂当作次方描述, 不是C/C++的xor op^^||) 10/26 20:19
7F：→ VictorTom:temp在第一轮回圈会算出1, 超过你的while执行范围, 就算 10/26 20:19
8F：→ VictorTom:出1跳回去了. 10/26 20:19
9F：→ VictorTom:既然是误差小於0.00000001, 所以应该是它和某个值相减以 10/26 20:20
10F：→ VictorTom:後的绝对值小於它, 而您temp的值看起来只是一项(第一项) 10/26 20:20
11F：→ VictorTom:所以您可能要重新检查您的实作与演算法:) 10/26 20:21
12F：→ VictorTom:试着step by step/debug几个测试资料就会发现这问题了:) 10/26 20:21
13F：→ VictorTom:另外, 像2F大大说的型别问题最好也注意一下, 还有printf 10/26 20:22
14F：→ VictorTom:要印double时要用%lf印:) 10/26 20:22
15F：推 PsMonkey:这样推文成章甘好？ 10/26 20:24
16F：推 VictorTom:小弟的推文常常是想到一点就推一点, 上面其实是分开想到 10/26 20:27
17F：→ VictorTom:的三件事Orz 要是一开始就想完了, 就会回文骗P币了XD 10/26 20:28
18F：→ fjf1980:甘温唷！ 10/26 20:32
19F：推 cismjmgoshr:fact(int x)也会有问题..x超过13的话会造成整数溢位 10/26 20:42
20F：→ cismjmgoshr:不过我觉得直接用math.h里面的cos(x)比较快 :p 10/26 21:01
21F：推 VictorTom:(大惊)原来这是算cos用的展开示啊....Orz 10/26 21:06
22F：→ fjf1980:越来越复杂了 数学果然烦人 10/26 21:31