在日本做机械设计的感受23：一道计算题的回帖

laotounihao

关注

eastking1111

非常感谢你的分享，上来就学了一手单变量求解！！！0.329784444989256

bbk1

谢谢分享

钢铁侠2014

大侠说到的回帖方式就是一种网络交流礼仪啊。如果大家都能注意，大家都会受益。
1 q' F, V. ]2 v6 Z4 T也谢谢大侠的单变量求解。

点滴积累

大侠，学习了。
/ R7 N- |3 }* k3 b* w关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问题，现在明白了一点，太深了。谢谢。
# m- Q$ c8 W2 j- B6 s0 {% V8 Y温习并熟练了单变量求解。
* _$ }# v6 h  a/ ]附图为我的计算结果：精度好像不太高，正在想方法解决。

子子61961

点滴积累 发表于 2016-2-15 16:27
) s% ^- i9 j. C( M, S大侠，学习了。
$ g# R1 ]2 k- x关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问 ...

感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，
' a* O$ N1 }, m$ ^3 |2 g& N一方面是从下层开始的自省和改善。
& f. P. e' s4 }1 E" W中国提到的管理，往往更偏向于前者，
而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
以前在第二话的文章里面写过一点，以后考虑也继续写一些。

关于单变量求解，它有计算的设定。
在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算
: E& R; G  `9 t+ S

（抱歉手边现在只有日文电脑，可以看里面的汉字参考）

里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。
+ m: D6 h: ]) [: o: J5 C+ N; w
另外，如何选取可变单元格也是一个因素。
( m0 e9 T! X8 Q/ Y- m2 o如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
其收敛的速度是不同的。

点滴积累

子子61961 发表于 2016-2-15 15:47
, x9 N# ]# v/ x/ B' k) f+ |% d感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

8 ]+ }' b5 X* G$ U& i
# M0 P* Y2 S( c7 B1 s据我所碰到的，国内的交流能力都很低。
1
当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
5 n) `% a, ^! ]9 a( a但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
  r, }# t; [% m- b: i. `2
3 I7 t  O  D1 n4 S; Z* G, [你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

我想楼主想说的是 不要想着什么工具 软件 用笔头算吧！！ 这是一个做设计人员对基本定理规律的理解 这是根本 不要太依赖工具软件 否则就失去了根本的理论基础 这样谈何开发

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
/ Y0 ?5 B5 H6 H) p6 U0 D, P; Z---------------------------------------------------------------------------------------------------------------------------------------
定义待求解函数：
0 v1 i* n( O4 M: |: W! qPublic Function QesFun(ByVal Var_AC As Double) As Double
8 G3 G& F8 V0 i) e. u) U/ @
5 v' B/ U3 x; G" j& l  w               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1

8 ]) G  `; g, I: l) CEnd Function
--------------------------------------------------------------------------------------------------
1.        二分法
1.1  由      Arctan(AC/80)=(AC-1)/80
      知      -PI()/2<(AC-1)/80< PI()/2
     即     1 -80*PI()/2<AC<1+ 80*PI()/2
2 g/ ?, d! h4 C% a  ?: G++++++++++++++++++++++++++++++++++++++
$ B' O- C! r' g) c+ e1.2定义求解函数：
' k7 N9 V9 |; {: j4 z( _3 [6 VPublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
7 m" w% I5 |7 n6 |
% T: V. O6 j5 x4 FDim Res#, VarAdj#

VarAdj = 10 ^ -6

# M  T' E$ S3 s6 ?" Y( x7 d3 ?3 lIf QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then
/ }2 j9 l  ~( o) Q# Y1 U
      Do While (1)

  g$ c& f. G1 O0 @0 u5 W              Res = (MaxLim + MinLim) / 2
0 F. B4 O% {  p' q0 U
               If Abs(QesFun(Res)) <= 10 ^ -12 Then
: I* Q9 ^2 g/ |& c
' M& l! S+ w& \' R* q/ \* Q                       SolFunDic = Res: Exit Do

# ]$ h' y: U) F              ElseIf (QesFun(Res) < 0) Then
0 j' W% W7 n$ T' l# }% d2 @  @. f
                      MinLim = Res

% B6 F- d0 n  T8 V& `1 L7 J             Else
8 b: e! u. V) X1 y
9 p2 \2 L- D0 `2 s7 z                     MaxLim = Res
0 \) M$ a+ w+ g6 ?
# `8 F2 s8 N* g3 K1 f3 p9 ^             End If
            
      Loop

8 T. c* n7 y. FElse
# g& J- V4 q/ F+ j* @6 c% ^
/ n  l" L. U+ m2 U4 Q8 M4 j  D     Do While (1)
     
               Res = (MaxLim + MinLim) / 2
% D2 F' |- g7 K( {3 `) _, A- E
             If Abs(QesFun(Res)) <= 10 ^ -12 Then
" l; p- S1 K: F! b8 c7 ?1 e
                        SolFunDic = Res: Exit Do
3 `4 {, D7 ]9 _+ H
# M& N4 T  C5 W0 n              ElseIf (QesFun(Res) > 0) Then

3 |: F( ^+ v) ~  _" k                        MinLim = Res

2 Y/ d4 b- q/ g6 m              Else

. P( V& U) C: o1 K5 j                         MaxLim = Res

            End If
            
      Loop
6 M$ o$ H" C. e' [! o' C* v* ]' d      
2 J+ n7 G4 P4 ]5 n' W& VEnd If
* y: ?7 N, ?+ ]4 Z* TEnd Function
( e8 E7 z9 ~- V/ y--------------------------------------------------------------
2.        牛顿法
' O3 C2 R/ M: j& B2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
        求导     f(AC)’=1/(1+(AC/80)^2)-1
        即      AC_1=AC_0- f(AC)/ f(AC)’
: @2 A9 n+ Y2 _( ?& A9 \--------------------------------------------------------------
" Y9 y' @. G' }  P: \2.2定义迭代函数：
Public Function QesFunNew(ByVal Var_AC As Double) As Double

8 T7 q3 J/ ?' p1 w              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
              
# I8 G* U& b$ Y: f$ i3 Z% tEnd Function
4 p8 \/ x, S+ i8 \7 t0 E4 _$ t---------------------------------------------------------------
( B) l( Z+ ?) }0 o- a3 @: K# k2.3定义求解函数
Public Function SolFunNew(ByVal IniAC As Double) As Double

Dim Res#

Do While (1)
0 ?/ C5 N8 _3 f; S( O* ]
                  Res = QesFunNew(IniAC)
3 K( t# D% h/ t  ?! r* N, |
! j3 U0 K. q9 B2 ?: E: o                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
  
$ O; d( L2 }! ?. v. U; s  i                                   SolFunNew = Res: Exit Do
  
                   Else
! ~7 x" ~8 i5 P! o6 Z( O                  
* M& f9 N) r2 W: l) f/ y# ^                                     IniAC = Res
                  
4 \$ U) K! s5 P! @* ^                  End If
) b6 M# q& S0 S7 z! z8 {
Loop
----------------------------------------------------------------------------------------------------------
" @% o# z. A' [6 c
( G. J! c" N" S# P4 A3 w这样做可能有点麻烦，但涉及到循环，迭代时，可自由调用自定义VBA函数和工作表函数（矩阵计算连杆机构），还可控制输出表格，便于插图。计算冷却塔时，积分得用辛普逊；解汽水比，得解非线性方程。对这些问题，EXCEL 工作表自身好像很吃力。