[課業] 請問各位一個國考程式題目

看板C_and_CPP作者 (Yes We Can!)時間14年前 (2009/10/26 20:02), 編輯推噓6(6016)
留言22則, 5人參與, 最新討論串1/1
※ [本文轉錄自 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

10/26 20:02, , 1F
請教高手們 我這樣寫有問題嗎?
10/26 20:02, 1F
※ 編輯: fjf1980 來自: 61.230.197.214 (10/26 20:03)

10/26 20:09, , 2F
sn的argument 是int? main裡面給double?
10/26 20:09, 2F

10/26 20:12, , 3F
我改double sn(double x) 結果還是一樣耶
10/26 20:12, 3F

10/26 20:18, , 4F
實作演算法有沒有錯呢??i從0開始代的話, -1^0是1, x*pi
10/26 20:18, 4F

10/26 20:18, , 5F
只要不是0, 代入2*i也相當於是^0, fact傳入0回傳1,
10/26 20:18, 5F

10/26 20:19, , 6F
(上面的^暫當作次方描述, 不是C/C++的xor op^^||)
10/26 20:19, 6F

10/26 20:19, , 7F
temp在第一輪迴圈會算出1, 超過你的while執行範圍, 就算
10/26 20:19, 7F

10/26 20:19, , 8F
出1跳回去了.
10/26 20:19, 8F

10/26 20:20, , 9F
既然是誤差小於0.00000001, 所以應該是它和某個值相減以
10/26 20:20, 9F

10/26 20:20, , 10F
後的絕對值小於它, 而您temp的值看起來只是一項(第一項)
10/26 20:20, 10F

10/26 20:21, , 11F
所以您可能要重新檢查您的實作與演算法:)
10/26 20:21, 11F

10/26 20:21, , 12F
試著step by step/debug幾個測試資料就會發現這問題了:)
10/26 20:21, 12F

10/26 20:22, , 13F
另外, 像2F大大說的型別問題最好也注意一下, 還有printf
10/26 20:22, 13F

10/26 20:22, , 14F
要印double時要用%lf印:)
10/26 20:22, 14F

10/26 20:24, , 15F
這樣推文成章甘好?
10/26 20:24, 15F

10/26 20:27, , 16F
小弟的推文常常是想到一點就推一點, 上面其實是分開想到
10/26 20:27, 16F

10/26 20:28, , 17F
的三件事Orz 要是一開始就想完了, 就會回文騙P幣了XD
10/26 20:28, 17F

10/26 20:32, , 18F
甘溫唷!
10/26 20:32, 18F

10/26 20:42, , 19F
fact(int x)也會有問題..x超過13的話會造成整數溢位
10/26 20:42, 19F

10/26 21:01, , 20F
不過我覺得直接用math.h裡面的cos(x)比較快 :p
10/26 21:01, 20F

10/26 21:06, , 21F
(大驚)原來這是算cos用的展開示啊....Orz
10/26 21:06, 21F

10/26 21:31, , 22F
越來越複雜了 數學果然煩人
10/26 21:31, 22F
文章代碼(AID): #1AvO-uor (C_and_CPP)