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

laotounihao

关注

eastking1111

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

bbk1

谢谢分享

钢铁侠2014

大侠说到的回帖方式就是一种网络交流礼仪啊。如果大家都能注意，大家都会受益。
) r2 S# `/ g! S9 S9 z6 Q' d也谢谢大侠的单变量求解。

点滴积累

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

子子61961

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

# Q$ d3 ?8 `% p感谢您的认真回帖。
3 x. B3 y* T: C7 p关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，
) v: X* K. M7 O; K一方面是从下层开始的自省和改善。
8 e* O" v' K6 h& w" J% U中国提到的管理，往往更偏向于前者，
而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
' s3 C- j3 J9 D0 D# @! D7 j以前在第二话的文章里面写过一点，以后考虑也继续写一些。

关于单变量求解，它有计算的设定。
. ^8 n% P1 W1 O$ F8 G/ q2 g1 b在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算
5 V! n$ i" \$ v) [* g
( O5 f" Q# n4 f/ Q& G: }
（抱歉手边现在只有日文电脑，可以看里面的汉字参考）
/ l8 A2 J3 [. X/ m( C
里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
3 P5 e) H8 M( R9 L+ l% `; _将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。
1 ^/ |  q# O. u4 w. ?
# w5 f) S8 h% \( D2 o) X1 U+ z另外，如何选取可变单元格也是一个因素。
+ w/ I) N0 |2 V0 k& k% U% r如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
* {, I1 ]& ^6 X- b* P; c; [如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
$ u. Y. W8 @1 J* z; y' Q其收敛的速度是不同的。
1 S) B- y9 z+ p9 S

点滴积累

子子61961 发表于 2016-2-15 15:47
* a6 R+ f- ^2 J' M5 [* ?+ m& @感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

" W, x2 y4 {6 w
据我所碰到的，国内的交流能力都很低。
3 j6 B( y, L$ w3 ^2 C8 Z+ D8 y1
3 b2 y- G* ]* S& i当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
# g/ s( S2 L6 z+ U5 L% _! J2
你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

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

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
---------------------------------------------------------------------------------------------------------------------------------------
1 V3 Z2 q" N# q  `/ J# H: P定义待求解函数：
Public Function QesFun(ByVal Var_AC As Double) As Double

1 ?8 `8 z+ e: [) k) z; e               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1

( v3 c3 z- a/ S7 O, P) mEnd Function
, r% ?5 n" M' g4 d6 M, `; z--------------------------------------------------------------------------------------------------
: r" r; X$ X3 f$ W" y1.        二分法
8 J4 a6 y! P7 D& ?, m/ W1 b1.1  由      Arctan(AC/80)=(AC-1)/80
) H' z+ h2 U, @- q- s      知      -PI()/2<(AC-1)/80< PI()/2
     即     1 -80*PI()/2<AC<1+ 80*PI()/2
; a! J0 i, C' i4 F& M2 }2 G$ N++++++++++++++++++++++++++++++++++++++
1.2定义求解函数：
/ G2 ?# R/ B6 K% z, t/ t; }, ~Public Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
2 K2 u1 s4 P: f6 g2 O% Z; S
  j; G* b; O9 i9 i0 ~! S: D* oDim Res#, VarAdj#
1 [3 k  z! h/ {: P4 w& F
& ^) G+ H& q; ^! MVarAdj = 10 ^ -6
7 G, m. k# D& X, a' c4 s+ K1 \
2 ^8 }4 g, X- M% @If QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then
+ p5 |: P! d& }
      Do While (1)

8 i4 |. l. L0 k" ^              Res = (MaxLim + MinLim) / 2
6 A" q4 Q8 L2 R( }6 `1 y
               If Abs(QesFun(Res)) <= 10 ^ -12 Then
( Q1 C* f7 p+ D7 ]  S' ^, F8 z& \$ Y
                       SolFunDic = Res: Exit Do

              ElseIf (QesFun(Res) < 0) Then
% K+ F# d: g5 r* G) h
8 {- x1 i$ U" O- E& e; r                      MinLim = Res

             Else

                     MaxLim = Res
# x5 N- n0 p8 |+ n: j- y$ I, V: l
  L. I/ J" K8 P             End If
2 l) I% p  \" ?            
      Loop

Else

  a7 T1 e3 m) T2 H  T" [5 M! j     Do While (1)
     
               Res = (MaxLim + MinLim) / 2

+ g" M- M, G% Q# @" W             If Abs(QesFun(Res)) <= 10 ^ -12 Then
& r* }# z* W& c* v
                        SolFunDic = Res: Exit Do

              ElseIf (QesFun(Res) > 0) Then

                        MinLim = Res
( S3 x: J9 Y8 H% ]9 W
/ c! e3 {, Z2 T1 r              Else
7 f4 _, g0 _1 P; p- k
1 }% ~! d, w& |                         MaxLim = Res

            End If
            
      Loop
      
End If
End Function
. T- @! N# M3 P. Y9 l5 T/ W--------------------------------------------------------------
4 N0 I, L1 `2 c8 e9 {! j. V# _& H2.        牛顿法
4 [6 p+ p4 i# Z3 T3 @2 @2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
        求导     f(AC)’=1/(1+(AC/80)^2)-1
( @5 H. j( f( ~$ e  _) I5 Q        即      AC_1=AC_0- f(AC)/ f(AC)’
5 V3 i+ e( [! G+ n. S7 {--------------------------------------------------------------
6 M. Y( F8 U  j" h- _2 B2.2定义迭代函数：
Public Function QesFunNew(ByVal Var_AC As Double) As Double
0 j9 D: v% Y/ u8 w6 ^; _
              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
. _$ C) u% H8 d$ O, |* r4 u              
End Function
---------------------------------------------------------------
2.3定义求解函数
7 F5 E' z& j7 p! T/ q; aPublic Function SolFunNew(ByVal IniAC As Double) As Double
; ]; b. [7 @/ @/ H1 _: `" D9 M# z
Dim Res#
4 I1 I/ h# [- A7 w5 `# T
Do While (1)
" X7 ?+ }+ e8 q* I& M5 q) [
" ^! j# T* x) \$ Z$ z  n                  Res = QesFunNew(IniAC)
" O+ J2 n5 K' G9 c/ I$ f" r# \
                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
  
                                   SolFunNew = Res: Exit Do
  
$ g' R6 l+ T8 }7 \, u                   Else
                  
                                     IniAC = Res
                  
                  End If
" s: ~. Z+ [( r1 X. J' R, B+ x6 Q
Loop
9 G* A, i5 _5 \# x$ @----------------------------------------------------------------------------------------------------------

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

- n. z% @) m; W! {  |