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

laotounihao

关注

eastking1111

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

bbk1

谢谢分享

钢铁侠2014

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

点滴积累

大侠，学习了。
3 y; Q" q  S$ T0 _0 O0 q: }" ]关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问题，现在明白了一点，太深了。谢谢。
温习并熟练了单变量求解。
- b5 y+ I7 a# Y1 m( o7 u" h附图为我的计算结果：精度好像不太高，正在想方法解决。

子子61961

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

感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
, k1 @6 @: u* E' s( v一方面是从上向下的监视和修正，
4 x% K8 N' l# `. s8 u% v! y一方面是从下层开始的自省和改善。
中国提到的管理，往往更偏向于前者，
& d! e1 ]2 r" I7 Q# S: `7 O而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
以前在第二话的文章里面写过一点，以后考虑也继续写一些。
) O2 W7 ^' _* V3 U) P& t2 ?( Q3 Y
关于单变量求解，它有计算的设定。
. _9 d5 O) }, g2 @- x在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算


（抱歉手边现在只有日文电脑，可以看里面的汉字参考）
) v7 C7 z: }1 Y% i6 ?
$ ^$ U0 M+ _- K6 D里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。
. E; v2 y: f- h; {& X% ]* Q8 N
另外，如何选取可变单元格也是一个因素。
如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
) u/ }3 {2 I) }1 B& x如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
8 z' v% G3 e7 h( Y. }- {9 d5 X其收敛的速度是不同的。
; B# r- J4 X& ^: ?6 E. U  z! |

点滴积累

子子61961 发表于 2016-2-15 15:47
感谢您的认真回帖。
3 g" }3 W) Y3 n4 C9 K* L$ D关于管理方面的问题，其实可以考虑为两个方面。
6 ~5 A+ K% m. F% k* j) I一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

' ~, [2 x  M& Q% ~& B# O
据我所碰到的，国内的交流能力都很低。
. p/ s+ P8 y  I1
当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
2
你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

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

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
---------------------------------------------------------------------------------------------------------------------------------------
定义待求解函数：
Public Function QesFun(ByVal Var_AC As Double) As Double
; y) E: f5 e  x- D7 M. @, M- S
               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1

End Function
3 J1 w. |7 N( O) [( \9 n--------------------------------------------------------------------------------------------------
1.        二分法
2 l) K6 ?7 [. P* A1 o1.1  由      Arctan(AC/80)=(AC-1)/80
      知      -PI()/2<(AC-1)/80< PI()/2
     即     1 -80*PI()/2<AC<1+ 80*PI()/2
+ X6 \! y7 k* D; y  z( c++++++++++++++++++++++++++++++++++++++
* }* g) `( i: ?. K& N1.2定义求解函数：
5 ]9 {& z9 {7 v# {/ p/ ~0 gPublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
) s3 [4 e- x" Z$ \+ W% B5 }
Dim Res#, VarAdj#

VarAdj = 10 ^ -6
) ]& X; j* P: O$ i5 ]+ _
If QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then

* ^" U# G) F/ v( Q2 v: @/ n9 I      Do While (1)

              Res = (MaxLim + MinLim) / 2

; V- v1 U5 G2 V  S/ V% C               If Abs(QesFun(Res)) <= 10 ^ -12 Then

( Q( f1 g5 f6 |& D4 c6 G# X$ ~                       SolFunDic = Res: Exit Do

              ElseIf (QesFun(Res) < 0) Then
9 L+ s! B, z0 f: i
" [4 g5 O- @5 l2 P8 U                      MinLim = Res
* o; B5 ~4 q4 S4 F% v6 w
( I% p3 i* T( ?  h             Else

                     MaxLim = Res
  a/ w: C& o" }+ h/ C
0 o& B8 q# Z8 A, j; i' Z             End If
            
! {4 T7 o1 M, e6 n; l9 H& l      Loop

Else
$ \$ z+ j8 F- q  T' ?; a3 [
     Do While (1)
     
4 q- R% s2 k; k7 U+ o8 V               Res = (MaxLim + MinLim) / 2

             If Abs(QesFun(Res)) <= 10 ^ -12 Then

" w# T: {9 l- H8 z% T+ e. ?4 q( d                        SolFunDic = Res: Exit Do

5 g. N, F5 A; ]/ P2 j7 X( G$ v, Y              ElseIf (QesFun(Res) > 0) Then
7 g8 R" S% c. L2 @5 |) s
                        MinLim = Res

              Else

                         MaxLim = Res
; d( i* M$ A5 {) F5 y- L" r* U
            End If
            
( `  F  ]2 A9 z  X1 t! p, B      Loop
      
* t+ c  m) u% QEnd If
9 q5 v( V8 A4 j9 |End Function
--------------------------------------------------------------
2.        牛顿法
2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
9 w' W& C1 w/ a( `# X8 ~- S# E; S        求导     f(AC)’=1/(1+(AC/80)^2)-1
  Q+ L, R( _$ s' X        即      AC_1=AC_0- f(AC)/ f(AC)’
: d; O; c+ C/ l. a4 ~6 O# R--------------------------------------------------------------
2.2定义迭代函数：
9 C" o# f$ l3 X7 A" D7 f( fPublic Function QesFunNew(ByVal Var_AC As Double) As Double
4 I6 Y8 x! d2 Y
              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
$ j1 m8 {5 q1 f              
End Function
---------------------------------------------------------------
2.3定义求解函数
Public Function SolFunNew(ByVal IniAC As Double) As Double

Dim Res#
) M3 D% w, @; D! M3 Y- u+ |: \
  v! ^( [3 g& Q4 K0 x1 G/ ADo While (1)

                  Res = QesFunNew(IniAC)

( P* a( e& j. e- @                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
  
                                   SolFunNew = Res: Exit Do
' b* E# `" i9 s5 R$ i" F- b1 u  
5 a8 k" b+ G6 L. N$ Z0 j3 B) V/ H                   Else
                  
                                     IniAC = Res
                  
4 W- T2 }! z1 r+ P: @0 B. B) z                  End If

Loop
8 ~) ?# F# p7 w0 E9 ?----------------------------------------------------------------------------------------------------------

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

! W+ @( u' K( m& d( _