求解一個C語言程序功能是每隔6°18′54″就增加6°18′54″并輸出結(jié)果的程序問題。

逍遙處士

1 z, }9 Y$ v  {, N4 F. ?) m5 r+ I! S
6 S( t: P- P  i, y

這是用MathCAD做的計算，程序跟一樓一樣，可見第5個也是29，就這來看樓主的程序可能有缺陷。
假如齒數(shù)為57，則360*3600/57不能整除，問題是否出在這兒？

! }1 ^: ]9 w' o6 ^假設(shè)所給齒數(shù)能整除的話，不妨只考慮將第一個角度依次翻倍，代碼如下:

int a=0, b=0, c=0;
3 U! s* J" ?! [& w) A3 ewhile ( a<360 )
{
7 l2 T: a- P: o. @- B% W- M# P  c=(c+s)%60; if ( (c+s)>=60 )  b=(c+s)/60;
3 `/ B4 b7 Y. T  [' H7 [6 F' E  b=(b+f)%60; if ( (b+f)>=60 )  a=(b+f)/60;
  a=a+j;
  printf("%d°%d′%d″\n",a,b,c);
) q" h. b- `  {: Q5 T& \}

不妨試一試，看下結(jié)果。

fmdd

隱含的升位運(yùn)算和強(qiáng)制降位運(yùn)算都會遭遇精度損失，如果取值范圍夠用的話，參與運(yùn)算的數(shù)據(jù)最好都統(tǒng)一為一種類型，比如全部為float型。
( S8 ?" K5 r8 L! f
$ [9 H& d1 `3 H. p樓主程序中的
% a, D! ^2 O9 o
8 P. S& o/ M1 |% w. O    int  m=22734;

該值是做什么用的？該值如何來的？

昨日明月

fmdd 發(fā)表于 2013-9-15 08:20
隱含的升位運(yùn)算和強(qiáng)制降位運(yùn)算都會遭遇精度損失，如果取值范圍夠用的話，參與運(yùn)算的數(shù)據(jù)最好都統(tǒng)一為一種類 ...

( a5 V# I1 O* I  S取的值，可以變化的。是6°18′54″轉(zhuǎn)化為秒數(shù)
8 p1 P- D2 x% ?. Z1 Y

fmdd

昨日明月 發(fā)表于 2013-9-15 08:22
取的值，可以變化的。是6°18′54″轉(zhuǎn)化為秒數(shù)

& l* s" ~6 v+ Y3 T) n這個取值本身就損失了精度

4 P* G9 t0 `* C9 e. ^6°18′54″是指一個齒的轉(zhuǎn)角是么？這個值的精度不夠，計算出來是6.31578947……°

這里的精度損失造成了最后累加數(shù)值偏差。

昨日明月

fmdd 發(fā)表于 2013-9-15 08:31
這個取值本身就損失了精度

  k" u5 l2 C9 Q5 f5 Z  a' [# }: I6°18′54″是指一個齒的轉(zhuǎn)角是么？這個值的精度不夠，計算出來是6.3157894 ...

老式銑床精度要求沒那么高。。。 達(dá)到秒數(shù)就夠了。。。

昨日明月

逍遙處士 發(fā)表于 2013-9-14 23:11
! i# E2 S6 U( }* u+ ~5 b這是用MathCAD做的計算，程序跟一樓一樣，可見第5個也是29，就這來看樓主的程序可能有缺陷。
  f: @7 D) h( U  ?假如齒數(shù) ...

1 U6 \% i! O- w( r1 O. A( D按照逍遙前輩的代碼運(yùn)行還是有問題呢？

zmq

這是一個很有意思的問題。我也試了以下，結(jié)果當(dāng)i = 5的時候，c 確實等于29。為此，當(dāng)i = 5的時候，檢查 c 取整以后的余數(shù) =1 ，也是一個整數(shù)。它就是 c 所差的數(shù)，由此，我就懷疑程序本身的內(nèi)部函數(shù)在處理“余數(shù)”的問題？！為此，我把6°18′54″化成22734″的數(shù)值改成了22734.0001″，結(jié)果一切正常。這說明你所編的程序并沒有錯。請試試看。

逍遙處士

把那兩行改成這樣試一試：
" Q/ w' k- r: t2 W' t9 i8 dif ( (c+s)>=60 )  b=(c+s)/60; c=(c+s)%60;
' k3 M9 @  ^' w9 {- {) Xif ( (b+f)>=60 )  a=(b+f)/60;  b=(b+f)%60;

昨日明月

逍遙處士 發(fā)表于 2013-9-15 10:41
( f" V, o( I  g0 j/ h把那兩行改成這樣試一試：
0 y3 u7 o1 b5 w0 Rif ( (c+s)>=60 )  b=(c+s)/60; c=(c+s)%60;
0 l# s7 u3 l; Tif ( (b+f)>=60 )  a=(b+f)/60;  ...

# m, e0 M2 n/ N1 u' J- l) q# w改后是這個結(jié)果，也不行。

逍遙處士

昨日明月 發(fā)表于 2013-9-15 11:42
: f4 E- Z/ I& B0 a# ~改后是這個結(jié)果，也不行。

4 D% R% u6 d$ r2 j9 G  Y
漏掉了一個 ＋ 號，哈哈！
1 i' u7 u1 ^: k5 D吾寶刀未老也！
6 V! {1 Y3 ]* @5 w- Y
1 {' ?' u  J) G# Q3 M" `" g9 A