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

laotounihao

关注

eastking1111

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

bbk1

谢谢分享

钢铁侠2014

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

点滴积累

大侠，学习了。
! x% q5 X  c4 `7 L5 R7 g关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问题，现在明白了一点，太深了。谢谢。
8 ]# x" L  E$ M温习并熟练了单变量求解。
' V2 k3 r5 B$ D4 B附图为我的计算结果：精度好像不太高，正在想方法解决。

子子61961

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

感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
, ^: K. n4 \0 f( w% @& t一方面是从上向下的监视和修正，
一方面是从下层开始的自省和改善。
- y9 O0 G) C- n+ e/ z* Q: L中国提到的管理，往往更偏向于前者，
而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
4 R0 S) z  b( J) V# Y5 I$ X" k以前在第二话的文章里面写过一点，以后考虑也继续写一些。
; K% J! H8 v" k
4 e4 B. l7 ^  B& L# I8 r关于单变量求解，它有计算的设定。
在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算
. O( T5 B. O0 s

# @* S- i) w* M  b/ Z, @' Q（抱歉手边现在只有日文电脑，可以看里面的汉字参考）

里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
2 S; T8 c4 ?/ j" R将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。

( b6 T* M4 P: B" u8 K" O$ z) h另外，如何选取可变单元格也是一个因素。
" I. l* S+ Y+ U$ S: @4 w如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
& V  Z! ~: L5 i" B# k8 G5 c0 O5 z如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
其收敛的速度是不同的。
+ Z* |( D8 m% @/ F& ~

点滴积累

子子61961 发表于 2016-2-15 15:47
0 m" ]* m; }  c3 T5 h, D感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

) j- i' f' g. p) p
据我所碰到的，国内的交流能力都很低。
3 A. c; B0 g% Y4 W1
当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
2
& X5 l' T; v8 I7 f你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

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

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
. I) s5 L' }( e+ ]/ \- s2 d---------------------------------------------------------------------------------------------------------------------------------------
' M- \% z! h; f/ x# I& U% w定义待求解函数：
Public Function QesFun(ByVal Var_AC As Double) As Double
5 F4 h/ F9 p, D
               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1
" D0 \( `9 t6 Q
9 i- l$ h; p* NEnd Function
--------------------------------------------------------------------------------------------------
1.        二分法
1.1  由      Arctan(AC/80)=(AC-1)/80
5 \- F1 o" }; `8 }* |      知      -PI()/2<(AC-1)/80< PI()/2
     即     1 -80*PI()/2<AC<1+ 80*PI()/2
: A$ Y4 L/ W% P9 ^7 L6 [++++++++++++++++++++++++++++++++++++++
1.2定义求解函数：
! \0 ~: A' z$ V9 a* ~4 UPublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
: j2 |9 g, d8 _5 D/ \, V
1 u* W8 ]8 N2 L4 P1 ^9 m! j0 zDim Res#, VarAdj#
" U7 S+ X7 Z/ E" k
+ G! N" r4 |3 E; B# o% m  a9 yVarAdj = 10 ^ -6

# S% m- H$ c2 O/ s7 XIf QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then

      Do While (1)

) ^! T* M6 I" G8 Z              Res = (MaxLim + MinLim) / 2
. _. B9 |8 b8 A9 h: |# F- H8 w. I
/ T' N& |( L( _7 Z# s6 b" M               If Abs(QesFun(Res)) <= 10 ^ -12 Then

                       SolFunDic = Res: Exit Do

1 {+ @# ~4 @, W# O; i              ElseIf (QesFun(Res) < 0) Then
: ~- v! J" s2 ?: S' E
; @4 ?. v( G3 i                      MinLim = Res
8 d0 Y+ H( v7 t# P7 Q
             Else

+ D2 g9 B- j* |& l                     MaxLim = Res

5 e5 f- A0 Y( V( y$ ]             End If
; {) ^2 j$ H7 q4 f( h/ E3 v- _            
: V) n% i9 Y$ O+ I) U1 b      Loop
. L+ F+ ?5 F3 |% p) T1 S
Else
, R; }% z9 v1 e  C5 Y
     Do While (1)
. B5 Z1 I9 l# S% P     
! W2 P$ B0 c, n0 _$ B2 U1 T               Res = (MaxLim + MinLim) / 2

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

& D7 s! X4 p' ]1 r' G7 q" F0 j1 ^                        SolFunDic = Res: Exit Do
$ S9 |/ M; ?- y; c# A# t
! a- r7 w! J2 a              ElseIf (QesFun(Res) > 0) Then

                        MinLim = Res
1 W- q& E/ B. N; V5 R
3 f: I3 ^" V" c5 `              Else
) b8 U5 d3 O( e/ y! z/ G- I; ^
6 k4 d" b* \# f7 X# @                         MaxLim = Res
# Y7 h7 ?, z1 g( X, o# A
8 E0 `6 |, R  T, q% H% t5 |            End If
$ h6 s0 s% O4 \0 Z) i9 g. \* c& z  y            
/ b6 A- n, j; }: i8 v% O3 [      Loop
* x/ u7 a) m( m      
End If
End Function
--------------------------------------------------------------
3 w! ^0 W$ R% r1 d* Z/ w7 N  h2.        牛顿法
2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
) c0 v8 u6 w/ N9 [6 i* b        求导     f(AC)’=1/(1+(AC/80)^2)-1
        即      AC_1=AC_0- f(AC)/ f(AC)’
--------------------------------------------------------------
2.2定义迭代函数：
Public Function QesFunNew(ByVal Var_AC As Double) As Double
% A. O$ @3 Y- ]1 N. f& F
7 r- _. w* {# |+ g              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
. b" U1 B1 x5 I; r/ t2 r; a0 e2 c              
End Function
---------------------------------------------------------------
3 E. a- z! u0 M' e0 o2.3定义求解函数
* X! F; h. M8 V8 r$ KPublic Function SolFunNew(ByVal IniAC As Double) As Double

Dim Res#
7 k) D% |' w3 r6 `' C
3 L  Y9 i2 A, G$ j/ XDo While (1)

                  Res = QesFunNew(IniAC)
4 T  ~) C6 K! w) Z# t! M
                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
) }: j3 H$ v2 r  
9 P# {/ k/ `% {! x) Q                                   SolFunNew = Res: Exit Do
6 z$ G( a$ o! a! |( B$ U2 ^, F1 D  
                   Else
                  
                                     IniAC = Res
                  
* N2 {3 j2 ?1 u* @% V: b1 Q! }/ O% Y2 _                  End If

Loop
1 b" K2 c& H: V; W! P) N) k5 y----------------------------------------------------------------------------------------------------------

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