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

laotounihao

关注

eastking1111

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

bbk1

谢谢分享

钢铁侠2014

大侠说到的回帖方式就是一种网络交流礼仪啊。如果大家都能注意，大家都会受益。
也谢谢大侠的单变量求解。

点滴积累

大侠，学习了。
关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问题，现在明白了一点，太深了。谢谢。
温习并熟练了单变量求解。
附图为我的计算结果：精度好像不太高，正在想方法解决。

子子61961

点滴积累 发表于 2016-2-15 16:27
大侠，学习了。
关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问 ...

感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
( n- H, ^( ^( Z" g* m1 l一方面是从上向下的监视和修正，
1 `) }+ i' t! w- ~一方面是从下层开始的自省和改善。
中国提到的管理，往往更偏向于前者，
而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
以前在第二话的文章里面写过一点，以后考虑也继续写一些。
) h$ g. o5 j5 I
" T* v/ f+ ^, |. P" t, G0 @关于单变量求解，它有计算的设定。
在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算
/ m  M% ?9 f4 b. Y' t

（抱歉手边现在只有日文电脑，可以看里面的汉字参考）
6 h6 b4 J$ e: |. t* A7 [
里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。
9 d/ z& [- `. q4 w8 O
另外，如何选取可变单元格也是一个因素。
' q, s( M+ p4 ^/ F如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
: n9 M2 J; a8 a6 @5 O4 b( L) e如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
7 s) k0 n/ K0 X# B) S0 M其收敛的速度是不同的。
/ z% w  |; _0 y% A) i

点滴积累

子子61961 发表于 2016-2-15 15:47
- A0 D+ O* d: p, f感谢您的认真回帖。
; [, W- {) B! s. n5 ^7 F& H关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

& v, O. k: J: l0 u2 b3 i/ d3 ~
据我所碰到的，国内的交流能力都很低。
$ S9 R  k  j; O, a+ z1 ^3 n1
& `: r, X% i+ Z9 [当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
2
你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

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

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
! X5 p1 m, h; B6 c: ]$ q9 K8 V子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
---------------------------------------------------------------------------------------------------------------------------------------
; r' v& v# n! @3 p7 h* e3 f0 I定义待求解函数：
Public Function QesFun(ByVal Var_AC As Double) As Double

               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1

6 ^$ l4 z+ y7 @" ~4 y' GEnd Function
--------------------------------------------------------------------------------------------------
" z% [7 P: B' G  E7 N1.        二分法
1.1  由      Arctan(AC/80)=(AC-1)/80
      知      -PI()/2<(AC-1)/80< PI()/2
     即     1 -80*PI()/2<AC<1+ 80*PI()/2
0 {2 U* i* d& @. M- y++++++++++++++++++++++++++++++++++++++
: t% [+ q8 M# s1.2定义求解函数：
Public Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
4 z% j0 Q" i- t1 M4 o0 ], i
% B7 ^- s2 x) Q/ f1 ]+ |* GDim Res#, VarAdj#

: g4 m* q  l1 RVarAdj = 10 ^ -6
, l! i1 i% B  W9 K1 y& L% |
3 ^+ p8 P% C. sIf QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then
. B/ @# \' i) D5 @, R' t, Z, r0 J  j
6 p7 T9 [- F4 i) z: b2 S      Do While (1)

              Res = (MaxLim + MinLim) / 2

               If Abs(QesFun(Res)) <= 10 ^ -12 Then
$ J6 c  f  `2 p7 `0 F% o9 j# L
                       SolFunDic = Res: Exit Do

              ElseIf (QesFun(Res) < 0) Then
$ G3 j/ D1 j+ E) s* f. j2 F  \
                      MinLim = Res
1 f' L3 p: N) s  A
2 k) D5 z5 L  j% y8 q             Else

( g' P! b% H  O+ I+ N/ |5 _/ R                     MaxLim = Res
7 M& K$ s: D9 C$ ]/ [% B: N" `4 c. a
             End If
9 T$ N* A5 M- r- R; F0 a            
      Loop
' O4 R$ H$ S, `4 }% ~( ~
7 ?' u/ Y5 U) @( b9 [Else

     Do While (1)
/ y% I1 \$ Q, H9 G     
               Res = (MaxLim + MinLim) / 2
4 u. _' {: H$ E  ^+ L# s
             If Abs(QesFun(Res)) <= 10 ^ -12 Then

9 N% S; ^% E) Y" G                        SolFunDic = Res: Exit Do
: z2 t: G- }9 D+ P1 \/ a$ B- v, j
' U5 m/ P' H. x! T              ElseIf (QesFun(Res) > 0) Then

/ _, ~, l% }* d* s  z5 o* Z# U                        MinLim = Res

: _, J8 f& T4 j              Else

6 Z! V2 J; P7 ^- C' E# e% Y                         MaxLim = Res
7 b  ?8 k5 \( @6 l
            End If
            
( W7 C, V' p! T0 J      Loop
      
End If
End Function
--------------------------------------------------------------
9 E/ a9 D! |) p- q( |+ N2.        牛顿法
: @7 m6 F0 O8 i+ S1 V6 y2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
! m& R4 D+ s9 j. T        求导     f(AC)’=1/(1+(AC/80)^2)-1
        即      AC_1=AC_0- f(AC)/ f(AC)’
8 T* z9 E) Z4 V  v--------------------------------------------------------------
2.2定义迭代函数：
Public Function QesFunNew(ByVal Var_AC As Double) As Double

3 V3 }' d0 i- r, ?% J! P              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
              
End Function
+ t/ }/ M& m  F1 i+ t. L---------------------------------------------------------------
+ {$ _' [1 T0 C& }2.3定义求解函数
$ M6 c5 g2 E+ I2 S* DPublic Function SolFunNew(ByVal IniAC As Double) As Double

Dim Res#

Do While (1)

                  Res = QesFunNew(IniAC)

                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
  
: i* b* ]& ~' D( e                                   SolFunNew = Res: Exit Do
  
                   Else
                  
1 ^* Y. M1 z" I( y. v3 f4 b1 _6 G                                     IniAC = Res
                  
                  End If
; S* ~* j3 F* M6 V) \9 D  i
' v) u5 k/ _  w! c( I& _- ^Loop
+ e; n: ~. g2 k----------------------------------------------------------------------------------------------------------

. s) m" B' h! e/ M* ~+ K& z- r) N这样做可能有点麻烦，但涉及到循环，迭代时，可自由调用自定义VBA函数和工作表函数（矩阵计算连杆机构），还可控制输出表格，便于插图。计算冷却塔时，积分得用辛普逊；解汽水比，得解非线性方程。对这些问题，EXCEL 工作表自身好像很吃力。