在日本做機械設計的感受23：一道計算題的回帖

laotounihao

關注

eastking1111

非常感謝你的分享，上來就學了一手單變量求解！！！0.329784444989256

bbk1

謝謝分享

鋼鐵俠2014

大俠說到的回帖方式就是一種網絡交流禮儀啊。如果大家都能注意，大家都會受益。
也謝謝大俠的單變量求解。

點滴積累

大俠，學習了。
關于中國文化，您提到的事情，實實在在存在著，也一直困擾著我。起初我一直懷疑是管理的問題，現在明白了一點，太深了。謝謝。
' N7 E, a) e7 u* p0 x  ]( v$ H! S溫習并熟練了單變量求解。
附圖為我的計算結果：精度好像不太高，正在想方法解決。

子子61961

點滴積累 發表于 2016-2-15 16:27
大俠，學習了。
3 D' B3 k7 D* [/ r' @9 v/ l關于中國文化，您提到的事情，實實在在存在著，也一直困擾著我。起初我一直懷疑是管理的問 ...

感謝您的認真回帖。
0 k- k, D( t8 O; T/ @, v關于管理方面的問題，其實可以考慮為兩個方面。
一方面是從上向下的監視和修正，
4 [1 t0 x" v0 i5 x) \/ }6 [) P( H: k一方面是從下層開始的自省和改善。
中國提到的管理，往往更偏向于前者，
2 l; t; ]& G  A6 u3 H4 a而日本的管理，其實是后者做得更好一些，
% K# H" ~  k6 l4 \4 y所以在日本，從上層向下的管理會輕松一些。
以前在第二話的文章里面寫過一點，以后考慮也繼續寫一些。
3 U) k% N; _! t% R
關于單變量求解，它有計算的設定。
在Excel的菜單里面，工具--〉選項--〉計算方法--〉反復計算
& G% _9 o, S* S% q9 |
$ H% C1 \3 w6 J+ f! ?- K
5 z/ v, m. c, V( u* g9 o（抱歉手邊現在只有日文電腦，可以看里面的漢字參考）
3 C( ~2 Q/ b+ ]0 z1 j0 W( p
里面有反復計算的設定值。最大反復次數100次，變化的最大值 0.001。
將這個次數改大一些，或者變化值改的更小一些，則會得到更精確一些的結果。

4 V2 D, [" O# K# `* C1 |3 R5 s另外，如何選取可變單元格也是一個因素。
如果選擇使用弧度，就是在0.32左右每次變化0.001進行求解。
如果選擇使用角度，就是在18.8左右每次變化0.001進行求解。
其收斂的速度是不同的。

點滴積累

子子61961 發表于 2016-2-15 15:47
  A  y  u! {" `感謝您的認真回帖。
關于管理方面的問題，其實可以考慮為兩個方面。
一方面是從上向下的監視和修正，

大俠，中文版，呵呵

andyany

: V! X0 e9 r( a* U' ~. ]- h9 n
! O4 L+ ^' c7 E8 K- Z3 f8 ?據我所碰到的，國內的交流能力都很低。
/ E2 d) r0 \0 b: u( Y6 ~3 P) U7 t5 Z1
當你問問題時，其實對方也不會，但他會說”這你都不會？！“ 讓人備受羞辱。
但有段時間和老外一起工作，老外不會時，一般也不會說”我不會“導致你干等或自己一通找，而是說”你找xxx，他會幫助你“，給你指路。讓人倍感溫暖。
( A+ _, G' r- N0 o# ]  C7 U2
& l/ h: \( X# G  @+ n: ?/ h( d  @4 U你覺得國內領導應該經歷過，于是去求助。領導說“這個你應該自己搞定”，你也沒話說。而你若是求助外方的領導，即便很大的官，也會給出你答案甚至理由和引申，讓你佩服得不要不要的。

408465954

我想樓主想說的是 不要想著什么工具 軟件 用筆頭算吧！！ 這是一個做設計人員對基本定理規律的理解 這是根本 不要太依賴工具軟件 否則就失去了根本的理論基礎 這樣談何開發

ngsxngtd

AC-arctan(AC/80)*80=1是數學問題。
# S- |5 w" w7 @子子大俠既然讀過VB，那么讀VBA代碼就應該沒什么困難，核心是數值計算。
' X' X* M1 D- ?2 K. M' F# S---------------------------------------------------------------------------------------------------------------------------------------
定義待求解函數：
Public Function QesFun(ByVal Var_AC As Double) As Double
2 b  m- c* v9 ~
               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1

End Function
9 x1 m# G$ e. _1 Y; x* W--------------------------------------------------------------------------------------------------
1.        二分法
0 r! ]5 S. a3 c* J5 G1.1  由      Arctan(AC/80)=(AC-1)/80
) _: T& S( _1 [2 `      知      -PI()/2<(AC-1)/80< PI()/2
' {# f6 |, p) C2 F1 _, n/ T- T' @     即     1 -80*PI()/2<AC<1+ 80*PI()/2
& T9 }7 y( T+ l# R0 a9 [! n6 D! w++++++++++++++++++++++++++++++++++++++
1.2定義求解函數：
4 d4 M0 C. K4 P# y5 m' s# w( [Public Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
0 q0 T+ S! o3 s. V( w# E
Dim Res#, VarAdj#
1 R5 g9 C- M1 z! R+ Z$ D
+ C7 u8 N; k7 a, d, iVarAdj = 10 ^ -6

If QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then
$ |1 Y. d1 \+ a. D8 Z
7 J' H: f8 T4 j+ k/ J      Do While (1)

- J. t2 h, G# R7 o  H' J              Res = (MaxLim + MinLim) / 2

% A. x% E. w3 c  m' |, n$ b4 ~               If Abs(QesFun(Res)) <= 10 ^ -12 Then

                       SolFunDic = Res: Exit Do

              ElseIf (QesFun(Res) < 0) Then

2 e$ _) d" {$ f                      MinLim = Res

$ Q9 e6 c& I+ T1 y" M1 r             Else

, z. T2 K! i/ y                     MaxLim = Res
1 u& j! v1 Q" v
             End If
            
      Loop
! E+ x9 J7 G  ~) V+ [1 L" Y
& @( f& s9 X$ C; t/ p3 m4 SElse

. w( y' t4 X1 x% t# h% i; `, F5 }     Do While (1)
     
/ n/ [; D! S+ S: S. F               Res = (MaxLim + MinLim) / 2
) u" t% F0 T/ i
             If Abs(QesFun(Res)) <= 10 ^ -12 Then
9 j& P  P% |$ I3 j
; i6 }$ q- f0 u5 Y. x. l7 i1 r                        SolFunDic = Res: Exit Do
& \& `3 X9 R+ g( d. d
              ElseIf (QesFun(Res) > 0) Then

                        MinLim = Res
( T8 O3 e% u9 {1 o0 l  H: ?- e
              Else
* ^; k9 B$ E9 y8 B) `: X
                         MaxLim = Res
7 I" _7 O( M- N' D! B
% S- p, _; Q* Q4 @7 I' @& Z$ A3 F            End If
            
      Loop
' U8 A$ z  I& u# v- N5 o! u      
: r. V9 P/ l4 R- }/ o( lEnd If
End Function
$ U0 m2 w' _* l+ Z( i--------------------------------------------------------------
, E4 h: N: `; x% E& I" i; F2.        牛頓法
2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
/ s8 [# n: Z7 Z        求導     f(AC)’=1/(1+(AC/80)^2)-1
% r* M1 E9 e( F: V        即      AC_1=AC_0- f(AC)/ f(AC)’
. |7 G+ P- x+ H0 \9 l" O--------------------------------------------------------------
9 ?9 p$ }: x! ]* N0 n7 Q# z2.2定義迭代函數：
" F5 F$ l. S* B$ KPublic Function QesFunNew(ByVal Var_AC As Double) As Double
1 X& _$ @% Q. ?6 q, ?
              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
4 H5 P5 k8 y# E" R1 @& A$ n              
' {, [) Q5 _2 B# ?End Function
$ Q0 m  P2 I' b/ l$ d$ W---------------------------------------------------------------
/ |: H% n" O) `6 _, ?1 i" C2.3定義求解函數
" w2 D: E2 j/ W7 ?- c& LPublic Function SolFunNew(ByVal IniAC As Double) As Double

Dim Res#

/ k- O2 U) [; C: p5 Z% YDo While (1)

4 x# k9 E7 @, w7 t                  Res = QesFunNew(IniAC)

                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
  
5 _7 l6 f+ D5 h' q8 H) N9 G! a" h                                   SolFunNew = Res: Exit Do
  z; ]! j6 h: }8 `% k$ O  
                   Else
                  
                                     IniAC = Res
+ @0 O; N1 U1 q$ t+ y* A- N+ I; F$ T                  
                  End If
+ f- I- j( t1 g9 ~
7 T8 [( ^: O0 E. C) jLoop
----------------------------------------------------------------------------------------------------------
4 k. C! a, n6 Y: A. j$ ~" k
6 {4 I% z* n7 O這樣做可能有點麻煩，但涉及到循環，迭代時，可自由調用自定義VBA函數和工作表函數（矩陣計算連桿機構），還可控制輸出表格，便于插圖。計算冷卻塔時，積分得用辛普遜；解汽水比，得解非線性方程。對這些問題，EXCEL 工作表自身好像很吃力。

. d9 [! v2 n" G3 [/ n! U9 w