作者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