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

laotounihao

关注

eastking1111

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

bbk1

谢谢分享

钢铁侠2014

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

点滴积累

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

子子61961

点滴积累 发表于 2016-2-15 16:27
% v+ D: b# T' {% T: M! \  E- o0 ]大侠，学习了。
关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问 ...

" v5 v; C# f9 E% h1 y感谢您的认真回帖。
关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，
( Y, N2 k# \$ x, T, f! ~+ q+ ^& p一方面是从下层开始的自省和改善。
1 V" X' g  \/ E5 i% w, I& I# x0 E4 g& G中国提到的管理，往往更偏向于前者，
) N& i) L7 a2 l/ |  W6 W3 C而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
以前在第二话的文章里面写过一点，以后考虑也继续写一些。

关于单变量求解，它有计算的设定。
& \& R! J( N& }0 _2 j7 }" x在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算
% c! E$ x# \: ~" s# @

（抱歉手边现在只有日文电脑，可以看里面的汉字参考）
  G) ^4 L% N9 ?3 H
里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
: C/ G, d( T& ~6 g" m; f( w5 G将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。
7 g  P; ]' r8 W' q) _
1 S- e& C9 t& b' c* q# u% c另外，如何选取可变单元格也是一个因素。
6 z9 o* ^7 ?' z+ d# W* A2 B6 _如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
$ C) Z( [  v2 f' g# k# S0 c1 p- g其收敛的速度是不同的。
' E: A6 J1 i5 |5 H( Y

点滴积累

子子61961 发表于 2016-2-15 15:47
感谢您的认真回帖。
4 p. w: U4 J9 O4 m2 S7 c) X关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

据我所碰到的，国内的交流能力都很低。
! [2 f2 H1 [- ~1
0 _  B, L/ y" ?+ ^) [" O当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
4 E( F7 p( l- n$ @- f3 A: e3 M但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
2
  i( `& I9 [+ G( K" i  U! O- {你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

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

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
( \+ q4 b& @0 i% \6 B4 o$ H7 J/ n" ^1 Z---------------------------------------------------------------------------------------------------------------------------------------
: o, O5 r, h' V# ^% @1 {$ k定义待求解函数：
# c# {1 B% }/ S; H) S5 D5 iPublic Function QesFun(ByVal Var_AC As Double) As Double

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

- V$ u4 p! c9 e+ B% e8 XEnd Function
--------------------------------------------------------------------------------------------------
* l0 i* T. K  ^7 z7 T8 n1.        二分法
1.1  由      Arctan(AC/80)=(AC-1)/80
      知      -PI()/2<(AC-1)/80< PI()/2
# c" ^' u6 e" \; c& F4 E+ F     即     1 -80*PI()/2<AC<1+ 80*PI()/2
++++++++++++++++++++++++++++++++++++++
1.2定义求解函数：
( y: A7 a# C* _# APublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double

Dim Res#, VarAdj#

  X4 t# T8 j; X+ Q0 s. b# K: mVarAdj = 10 ^ -6
% B& W, z9 [" N
3 s: u& d( i1 l6 }4 X! LIf QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then

! p2 @2 p- z) G, [0 w# I% ^      Do While (1)
; q- m- [( G2 f3 t8 ^
8 i7 d# {- w! L              Res = (MaxLim + MinLim) / 2

9 Y+ D2 c) n; X( u3 d" ~1 _0 M! R               If Abs(QesFun(Res)) <= 10 ^ -12 Then

, d0 `& s6 `% H                       SolFunDic = Res: Exit Do

% S+ m5 t6 X2 H  T' Y2 K              ElseIf (QesFun(Res) < 0) Then

                      MinLim = Res

+ q2 \' G+ J# t, K/ f             Else
9 ]& G7 z- f* b" X8 P
6 h. `* g# y; P0 a* M6 J% j6 y                     MaxLim = Res
. e4 y9 h" [2 f: U
             End If
+ G7 J: m/ X& c! p" ^. F# K            
      Loop
3 Z8 L# x$ @2 H4 ]
Else

. H) ^6 f  T5 x: x+ k! ~# Q+ z" d     Do While (1)
     
3 V" J- Q1 e- @) \, _$ c2 C               Res = (MaxLim + MinLim) / 2

' \/ l: u+ p/ d  T" |             If Abs(QesFun(Res)) <= 10 ^ -12 Then
8 S" t8 O0 p9 H
- g; _* w* f* `) X! J. v* O3 w+ }                        SolFunDic = Res: Exit Do

              ElseIf (QesFun(Res) > 0) Then
0 @) l- Q/ q4 ?8 u- ]) d6 C% x- Y
* L3 M/ `7 o# s* g0 e& x                        MinLim = Res

+ y* p9 N' R! w8 |( _              Else
# \. [4 \( s1 m3 D$ K: H
" I' h7 j/ i* E8 d& b& Q3 b7 F& R  x                         MaxLim = Res
: {; k; B6 g" z8 V* n9 \8 H& P
$ e- H; E4 h0 C' M            End If
+ X# F; P& W0 ?% n9 r            
0 o" r7 e' b! i( M6 U      Loop
! m' \9 s/ B5 r, H! S      
) b5 q; a& `6 |. ^/ j% cEnd If
$ O5 Q9 T- _4 f) g" b' nEnd Function
--------------------------------------------------------------
+ E) P. j- e/ M2.        牛顿法
2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
: Z/ c( C' a  ?3 a4 K6 w6 K% o        求导     f(AC)’=1/(1+(AC/80)^2)-1
$ {/ W; m0 @- E/ h, E% P; X/ |% Q        即      AC_1=AC_0- f(AC)/ f(AC)’
--------------------------------------------------------------
! h6 K8 |' A5 t1 C2.2定义迭代函数：
- T! s3 X" u: \6 R- TPublic Function QesFunNew(ByVal Var_AC As Double) As Double

              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
) Z' _: S1 b" }1 L. F              
* g/ I5 I( F* b: i$ nEnd Function
2 G' \0 Q7 D4 G  W1 ~---------------------------------------------------------------
; d, E$ d3 o" M+ d$ ^  G2.3定义求解函数
Public Function SolFunNew(ByVal IniAC As Double) As Double
% v6 S# S6 m1 T8 f& K9 z! m
* ~" o3 z. m; `; r! v& }' j6 GDim Res#
& G( V1 F8 N4 \
) p" L+ d& U% ~6 J- W# \Do While (1)
0 a  w  B4 l9 t& |! t
) ]# C4 D0 Y) A5 r9 q                  Res = QesFunNew(IniAC)
6 p1 A2 n- Q6 j* i, q( B
" G0 A8 _1 U) c                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
5 m  o. q$ u0 O+ e) o% w5 j  
4 Q# X& G# P& f* i6 Y8 O/ k                                   SolFunNew = Res: Exit Do
* y5 [. @! I$ u$ q! a$ o4 t1 s9 F. w  
                   Else
                  
                                     IniAC = Res
                  
# h/ i# K7 u3 g9 Y3 K/ J$ t& N% d2 D! }8 ~                  End If

Loop
----------------------------------------------------------------------------------------------------------
8 d, A+ J. i; G7 A* w# o- `6 z$ ^- o
+ `; k1 W0 A3 I& {* L/ J这样做可能有点麻烦，但涉及到循环，迭代时，可自由调用自定义VBA函数和工作表函数（矩阵计算连杆机构），还可控制输出表格，便于插图。计算冷却塔时，积分得用辛普逊；解汽水比，得解非线性方程。对这些问题，EXCEL 工作表自身好像很吃力。
4 d8 `9 m2 A$ G% Y4 C% S. P3 x, i
; k% n( g6 J4 a" |: I# g! x" c  o( u