在日本做機(jī)械設(shè)計(jì)的感受23：一道計(jì)算題的回帖

laotounihao

關(guān)注

eastking1111

非常感謝你的分享，上來(lái)就學(xué)了一手單變量求解�。�！0.329784444989256

bbk1

謝謝分享

鋼鐵俠2014

大俠說(shuō)到的回帖方式就是一種網(wǎng)絡(luò)交流禮儀啊。如果大家都能注意，大家都會(huì)受益。
3 b" G& J- _% q  p+ n也謝謝大俠的單變量求解。

點(diǎn)滴積累

大俠，學(xué)習(xí)了。
6 B: A6 Z% h, l& N/ i關(guān)于中國(guó)文化，您提到的事情，實(shí)實(shí)在在存在著，也一直困擾著我。起初我一直懷疑是管理的問(wèn)題，現(xiàn)在明白了一點(diǎn)，太深了。謝謝。
5 t+ L. R( a3 b. u9 a溫習(xí)并熟練了單變量求解。
# J8 W5 Q% M: d0 E' B" r9 E* v附圖為我的計(jì)算結(jié)果：精度好像不太高，正在想方法解決。

子子61961

點(diǎn)滴積累 發(fā)表于 2016-2-15 16:27
大俠，學(xué)習(xí)了。
關(guān)于中國(guó)文化，您提到的事情，實(shí)實(shí)在在存在著，也一直困擾著我。起初我一直懷疑是管理的問(wèn) ...

3 J4 Q# e# z9 w2 \, }感謝您的認(rèn)真回帖。
關(guān)于管理方面的問(wèn)題，其實(shí)可以考慮為兩個(gè)方面。
7 ^' |& O6 I5 n/ f8 z% U2 F5 F一方面是從上向下的監(jiān)視和修正，
; g* ~# _6 j: k  `4 z* t" L一方面是從下層開(kāi)始的自省和改善。
中國(guó)提到的管理，往往更偏向于前者，
' R! d6 z5 n3 ]4 _2 l5 T而日本的管理，其實(shí)是后者做得更好一些，
4 K% J  E+ B+ k3 Q) K所以在日本，從上層向下的管理會(huì)輕松一些。
; ^: j& N0 L. f. R! s/ }% s以前在第二話的文章里面寫過(guò)一點(diǎn)，以后考慮也繼續(xù)寫一些。
8 h  m) S4 o) `# X  g" t5 @! X+ f& Q
) A% i- I- G  u: k$ J& P" `關(guān)于單變量求解，它有計(jì)算的設(shè)定。
在Excel的菜單里面，工具--〉選項(xiàng)--〉計(jì)算方法--〉反復(fù)計(jì)算
  U; @/ }7 f, S' X, F, i, H. u
5 J8 _3 V  s, r4 w+ `, ]
- [2 |  h) R( w8 z& \（抱歉手邊現(xiàn)在只有日文電腦，可以看里面的漢字參考）
% _8 X" o1 G3 |/ }
6 C+ [6 g. t: b4 m$ m7 a$ I里面有反復(fù)計(jì)算的設(shè)定值。最大反復(fù)次數(shù)100次，變化的最大值 0.001。
將這個(gè)次數(shù)改大一些，或者變化值改的更小一些，則會(huì)得到更精確一些的結(jié)果。

: T7 E# _( t+ p. |5 K% Z另外，如何選取可變單元格也是一個(gè)因素。
如果選擇使用弧度，就是在0.32左右每次變化0.001進(jìn)行求解。
如果選擇使用角度，就是在18.8左右每次變化0.001進(jìn)行求解。
其收斂的速度是不同的。

點(diǎn)滴積累

子子61961 發(fā)表于 2016-2-15 15:47
感謝您的認(rèn)真回帖。
  ^: G- S5 i+ S) _0 D/ w: w關(guān)于管理方面的問(wèn)題，其實(shí)可以考慮為兩個(gè)方面。
一方面是從上向下的監(jiān)視和修正，

大俠，中文版，呵呵

andyany

9 e: {7 }+ a1 u$ H) b據(jù)我所碰到的，國(guó)內(nèi)的交流能力都很低。
' k1 o* R5 U4 }! X1
當(dāng)你問(wèn)問(wèn)題時(shí)，其實(shí)對(duì)方也不會(huì)，但他會(huì)說(shuō)”這你都不會(huì)？！“ 讓人備受羞辱。
但有段時(shí)間和老外一起工作，老外不會(huì)時(shí)，一般也不會(huì)說(shuō)”我不會(huì)“導(dǎo)致你干等或自己一通找，而是說(shuō)”你找xxx，他會(huì)幫助你“，給你指路。讓人倍感溫暖。
' y: k+ H9 |5 S/ B2
你覺(jué)得國(guó)內(nèi)領(lǐng)導(dǎo)應(yīng)該經(jīng)歷過(guò)，于是去求助。領(lǐng)導(dǎo)說(shuō)“這個(gè)你應(yīng)該自己搞定”，你也沒(méi)話說(shuō)。而你若是求助外方的領(lǐng)導(dǎo)，即便很大的官，也會(huì)給出你答案甚至理由和引申，讓你佩服得不要不要的。

408465954

我想樓主想說(shuō)的是 不要想著什么工具 軟件 用筆頭算吧！！ 這是一個(gè)做設(shè)計(jì)人員對(duì)基本定理規(guī)律的理解 這是根本 不要太依賴工具軟件 否則就失去了根本的理論基礎(chǔ) 這樣談何開(kāi)發(fā)

ngsxngtd

AC-arctan(AC/80)*80=1是數(shù)學(xué)問(wèn)題。
子子大俠既然讀過(guò)VB，那么讀VBA代碼就應(yīng)該沒(méi)什么困難，核心是數(shù)值計(jì)算。
! J. x$ k6 u1 R7 c" ?& F, H3 m5 W---------------------------------------------------------------------------------------------------------------------------------------
( X' Q% z  y  z9 N" J定義待求解函數(shù)：
Public Function QesFun(ByVal Var_AC As Double) As Double
! ~3 J$ V8 U/ ~5 \  k& [
( F& x+ j& ~9 ~0 E               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1
3 L9 v  {* O6 f7 [
End Function
) ~/ c- V+ W5 T* R' e7 h$ ?--------------------------------------------------------------------------------------------------
1.        二分法
1.1  由      Arctan(AC/80)=(AC-1)/80
      知      -PI()/2<(AC-1)/80< PI()/2
, Q  W3 i$ s; a8 ^     即     1 -80*PI()/2<AC<1+ 80*PI()/2
++++++++++++++++++++++++++++++++++++++
% r* f" j& a: M6 a0 h1.2定義求解函數(shù)：
/ _7 a% a" U9 Q( m8 W! lPublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double

% ^  f+ y: f9 K' v9 Q+ \, e* I( ]Dim Res#, VarAdj#
$ Z2 h# x% ~$ E* ^7 @. L0 Q
3 w: O+ V4 v# p9 n$ V! j3 sVarAdj = 10 ^ -6
' D: ]; S1 K6 d' ]* ~; v
4 Z9 {0 @; d/ G' ZIf QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then

      Do While (1)
; Q) u4 @: t! N$ j3 a2 ?. H
              Res = (MaxLim + MinLim) / 2
# l8 a! B$ a) ^% ?/ R6 L# W6 ?
  r1 I$ C( f& l/ E  A8 U! ]4 W               If Abs(QesFun(Res)) <= 10 ^ -12 Then
: j) \+ [# E) q9 p; @
- ]5 |- {# ]4 p( z; Q. f0 W, \% s* V                       SolFunDic = Res: Exit Do

. m" g6 @5 k& z; @4 l              ElseIf (QesFun(Res) < 0) Then
2 m3 T3 b6 G4 u
, z# s; b7 A5 P; q                      MinLim = Res
; O. T% a' t+ m: `, C8 B
             Else

                     MaxLim = Res
5 _. N* z. q0 z6 [' D
             End If
            
! c1 L9 |) ^  o+ C) M4 ]* n      Loop
. f( m+ t; m/ D( o3 z$ E
4 a4 |7 C; i( z' ^7 \* `Else
, J0 b5 {% X! d; y
     Do While (1)
& F. L7 e! U  W8 ~     
               Res = (MaxLim + MinLim) / 2

             If Abs(QesFun(Res)) <= 10 ^ -12 Then
9 p" g& Z! \" g$ M/ u
$ L/ o# @9 u, l6 X6 p                        SolFunDic = Res: Exit Do
! r# e7 y! o* I0 a$ e7 J
5 n: M- d) j7 ]% }/ D              ElseIf (QesFun(Res) > 0) Then
* P! u; f! e; h2 _
7 ]7 M3 [7 f7 e, }! U! l! l/ \                        MinLim = Res

              Else

! k# T+ L' x; q1 ?0 q                         MaxLim = Res
# I. j  O. z8 B/ @2 |0 X+ J5 H
8 H& V8 F5 B& U9 m8 s  o% i+ h            End If
            
      Loop
      
4 L( \* k9 K( T4 N5 \End If
End Function
--------------------------------------------------------------
/ e1 ^* m' O6 h2.        牛頓法
! g. \  S: m, x: I2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
        求導(dǎo)     f(AC)’=1/(1+(AC/80)^2)-1
        即      AC_1=AC_0- f(AC)/ f(AC)’
--------------------------------------------------------------
2.2定義迭代函數(shù)：
Public Function QesFunNew(ByVal Var_AC As Double) As Double
' z% k( w2 T0 _, @8 W9 _$ \$ g
/ G& T" D9 f7 ^" `9 E# M; s( @              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
) v8 }8 z1 B' q! f5 d5 {              
; d; j; ~3 s8 d( B+ ?; _, GEnd Function
7 `  |1 b: g% j# g# W---------------------------------------------------------------
! Q8 l' Q4 ?6 w6 B9 {9 U  x: l7 c2.3定義求解函數(shù)
9 K% h, d! t: [Public Function SolFunNew(ByVal IniAC As Double) As Double

; ~) d3 i2 g5 V$ k! u" k+ m: eDim Res#
$ x9 W! c- C6 g; ]/ |
$ X! e9 k2 b7 W7 j) W4 L( bDo While (1)
2 y* [8 t; Y) n5 j4 k; j
                  Res = QesFunNew(IniAC)

2 J+ a9 r; W/ N7 D: X                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
  
                                   SolFunNew = Res: Exit Do
9 ?* o8 a8 e" s& s+ p  
1 w; w  g4 z0 {1 p0 y  g! I                   Else
1 S; y4 u+ l% g9 G                  
                                     IniAC = Res
                  
0 b6 M8 K; k6 {                  End If
9 J  F. p- I' J; X! B3 N( H8 ~
7 A* K2 ?, V* z& k7 ?* LLoop
4 T& i; o- k6 F/ _----------------------------------------------------------------------------------------------------------
; l! i: v' j) {8 {
這樣做可能有點(diǎn)麻煩，但涉及到循環(huán)，迭代時(shí)，可自由調(diào)用自定義VBA函數(shù)和工作表函數(shù)（矩陣計(jì)算連桿機(jī)構(gòu)），還可控制輸出表格，便于插圖。計(jì)算冷卻塔時(shí)，積分得用辛普遜；解汽水比，得解非線性方程。對(duì)這些問(wèn)題，EXCEL 工作表自身好像很吃力。