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

逍遙處士

% t* B4 a* L; S$ o1 U

9 U5 z( k" j" A1 y8 b+ Y+ Y2 G
2 Q/ j( _" ?4 H& W. u$ ~這是用MathCAD做的計(jì)算,，程序跟一樓一樣,，可見第5個(gè)也是29，就這來看樓主的程序可能有缺陷。
$ _3 Y+ z! Q  w假如齒數(shù)為57,，則360*3600/57不能整除,，問題是否出在這兒？
( {5 ]1 R2 [/ G- Z* c' @# Z* T
假設(shè)所給齒數(shù)能整除的話,，不妨只考慮將第一個(gè)角度依次翻倍,，代碼如下:

int a=0, b=0, c=0;
while ( a<360 )
{
  c=(c+s)%60; if ( (c+s)>=60 )  b=(c+s)/60;
1 f2 Q! ]" G8 L9 |  b=(b+f)%60; if ( (b+f)>=60 )  a=(b+f)/60;
+ o7 e$ o4 O) _0 y0 o3 O  a=a+j;
% y) S+ C* `" F- _  printf("%d°%d′%d″\n",a,b,c);
5 G& t$ w4 e* g1 Q}
7 [7 Z3 j, u; \/ N
" ?+ m) E/ ^+ `4 j# v. G) c/ ]不妨試一試，看下結(jié)果,。

fmdd

隱含的升位運(yùn)算和強(qiáng)制降位運(yùn)算都會(huì)遭遇精度損失，如果取值范圍夠用的話,，參與運(yùn)算的數(shù)據(jù)最好都統(tǒng)一為一種類型,，比如全部為float型。

樓主程序中的

* e" Y" X5 L4 N, a; U. p4 Z: B    int  m=22734;

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

昨日明月

fmdd 發(fā)表于 2013-9-15 08:20
: c* X, n- _7 g& D' E5 J; V3 ^! B隱含的升位運(yùn)算和強(qiáng)制降位運(yùn)算都會(huì)遭遇精度損失，如果取值范圍夠用的話,，參與運(yùn)算的數(shù)據(jù)最好都統(tǒng)一為一種類 ...

取的值,，可以變化的。是6°18′54″轉(zhuǎn)化為秒數(shù)
$ X) A) b- o. u( p8 }+ ]' s

fmdd

$ O% O$ V9 d6 v$ e7 _

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

這個(gè)取值本身就損失了精度

6°18′54″是指一個(gè)齒的轉(zhuǎn)角是么？這個(gè)值的精度不夠,，計(jì)算出來是6.31578947……°

這里的精度損失造成了最后累加數(shù)值偏差,。
4 I9 c% n' v# G. A/ e+ L  P0 J+ L0 a

昨日明月

fmdd 發(fā)表于 2013-9-15 08:31
這個(gè)取值本身就損失了精度
9 S8 o) D" y/ j* `
# e* z7 f2 c0 C+ W- J' k6°18′54″是指一個(gè)齒的轉(zhuǎn)角是么？這個(gè)值的精度不夠,，計(jì)算出來是6.3157894 ...

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

昨日明月

逍遙處士 發(fā)表于 2013-9-14 23:11
這是用MathCAD做的計(jì)算，程序跟一樓一樣,，可見第5個(gè)也是29,，就這來看樓主的程序可能有缺陷。
假如齒數(shù) ...

* q% ^0 U" R' C按照逍遙前輩的代碼運(yùn)行還是有問題呢？

zmq

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

逍遙處士

把那兩行改成這樣試一試：
& U7 u- S6 l: Sif ( (c+s)>=60 )  b=(c+s)/60; c=(c+s)%60;
if ( (b+f)>=60 )  a=(b+f)/60;  b=(b+f)%60;

昨日明月

逍遙處士 發(fā)表于 2013-9-15 10:41
把那兩行改成這樣試一試：
$ r* ~- V! g+ w' a) x. o* Pif ( (c+s)>=60 )  b=(c+s)/60; c=(c+s)%60;
if ( (b+f)>=60 )  a=(b+f)/60;  ...

改后是這個(gè)結(jié)果，也不行,。

逍遙處士

昨日明月 發(fā)表于 2013-9-15 11:42
$ \( D. D2 f2 m0 ?; h$ ~9 r5 p改后是這個(gè)結(jié)果,，也不行,。

4 A: T( A6 i3 J' u$ Y; u+ m
漏掉了一個(gè) ＋ 號(hào)，哈哈,！
5 b9 m+ C; Z4 F& G/ r吾寶刀未老也,！
- S" h, e9 m: j