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

laotounihao

关注

eastking1111

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

bbk1

谢谢分享

钢铁侠2014

大侠说到的回帖方式就是一种网络交流礼仪啊。如果大家都能注意，大家都会受益。
% _1 }/ ]& w: u2 R* n( t: N也谢谢大侠的单变量求解。

点滴积累

大侠，学习了。
# d1 i$ z2 Z8 h( Z( k9 h4 `1 {& ^关于中国文化，您提到的事情，实实在在存在着，也一直困扰着我。起初我一直怀疑是管理的问题，现在明白了一点，太深了。谢谢。
% g/ M+ W! r6 ~8 Y" c5 H温习并熟练了单变量求解。
附图为我的计算结果：精度好像不太高，正在想方法解决。

子子61961

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

9 h6 ^, A2 @. H' T$ Z% N感谢您的认真回帖。
5 o5 ?% G* i. I& Z, k3 D关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，
" I5 y+ B/ u. `+ k一方面是从下层开始的自省和改善。
* J1 G) e8 A# i  q8 \中国提到的管理，往往更偏向于前者，
6 m3 W8 W: R) \. G4 t$ _8 ^+ m而日本的管理，其实是后者做得更好一些，
所以在日本，从上层向下的管理会轻松一些。
! p7 H$ k/ \; m5 C1 P  G* }以前在第二话的文章里面写过一点，以后考虑也继续写一些。

关于单变量求解，它有计算的设定。
% f9 U5 V- c$ S" z- e* |在Excel的菜单里面，工具--〉选项--〉计算方法--〉反复计算
' ?' d/ W% C/ t

1 S, {* N+ P5 B& e4 ^. L! C（抱歉手边现在只有日文电脑，可以看里面的汉字参考）
- }" p$ P! N* w  Y
里面有反复计算的设定值。最大反复次数100次，变化的最大值 0.001。
2 ^6 W6 {6 F; \$ W将这个次数改大一些，或者变化值改的更小一些，则会得到更精确一些的结果。
) `  T8 A) {7 g( V4 N2 z
/ _% z/ w5 p5 P: s% m另外，如何选取可变单元格也是一个因素。
如果选择使用弧度，就是在0.32左右每次变化0.001进行求解。
/ H# d$ v3 Q1 L! w如果选择使用角度，就是在18.8左右每次变化0.001进行求解。
5 o% L$ H, }4 f; J! n1 h+ S其收敛的速度是不同的。
2 a  g+ e4 n* l+ e2 g

点滴积累

子子61961 发表于 2016-2-15 15:47
1 G: D! Q$ N6 B% n5 E感谢您的认真回帖。
1 s: e+ {% n9 l! |关于管理方面的问题，其实可以考虑为两个方面。
一方面是从上向下的监视和修正，

大侠，中文版，呵呵

andyany

据我所碰到的，国内的交流能力都很低。
1
0 E& r4 e, [1 m5 o; c当你问问题时，其实对方也不会，但他会说”这你都不会？！“ 让人备受羞辱。
但有段时间和老外一起工作，老外不会时，一般也不会说”我不会“导致你干等或自己一通找，而是说”你找xxx，他会帮助你“，给你指路。让人倍感温暖。
2
! S( ~5 Y. P; A( |) A( X" U你觉得国内领导应该经历过，于是去求助。领导说“这个你应该自己搞定”，你也没话说。而你若是求助外方的领导，即便很大的官，也会给出你答案甚至理由和引申，让你佩服得不要不要的。

408465954

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

ngsxngtd

AC-arctan(AC/80)*80=1是数学问题。
1 W7 @" G# h$ R子子大侠既然读过VB，那么读VBA代码就应该没什么困难，核心是数值计算。
---------------------------------------------------------------------------------------------------------------------------------------
7 W/ U/ b) E! P8 S3 \, b' J. P定义待求解函数：
0 H* M5 @) ~0 H+ T. HPublic Function QesFun(ByVal Var_AC As Double) As Double
# W( t% S) ~9 L* T
2 E5 g3 t' z& \$ T3 N' \) a               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 1

9 R) _: y" f3 u! Y; REnd Function
--------------------------------------------------------------------------------------------------
4 [, F* c6 n; J. ?* q" J- A/ j: D1.        二分法
1.1  由      Arctan(AC/80)=(AC-1)/80
8 Z1 h6 A5 e$ _; n6 A" p7 R      知      -PI()/2<(AC-1)/80< PI()/2
. ?- W: y3 i* u- F% u" @     即     1 -80*PI()/2<AC<1+ 80*PI()/2
++++++++++++++++++++++++++++++++++++++
1.2定义求解函数：
; {9 j* H$ x, r1 J5 }+ [7 bPublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double

Dim Res#, VarAdj#

7 R1 l3 j: X! N+ V- qVarAdj = 10 ^ -6
2 P0 H2 @$ V/ m
! l9 n2 P! u, }( _; y1 i* Y: r* ZIf QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then

0 T' ^2 F) t1 f! a' L4 N      Do While (1)
0 Z* `! N/ e+ |( h
              Res = (MaxLim + MinLim) / 2
+ l4 q! T% D1 A$ @1 P  W
, C4 ^: }$ Y6 |8 L" V               If Abs(QesFun(Res)) <= 10 ^ -12 Then
: @& l3 \# I5 N7 ~
                       SolFunDic = Res: Exit Do
. j8 H0 [5 A" |8 B, ^& g7 B( t
              ElseIf (QesFun(Res) < 0) Then
! [$ F5 k2 V/ |! M% Z5 e. d
                      MinLim = Res
/ w1 s# o$ m8 i- x; ]
             Else

                     MaxLim = Res
: k) R2 F; Z/ i* ^8 ]
& _7 Q0 N# f. b* ^* M( w$ c             End If
            
) O/ g; ~: Z5 |; ]$ |% C) o      Loop
7 i/ B% F2 I1 [+ c( D
9 X+ y- R* M( r- ?) A) L2 YElse

& v+ f: l, B+ d$ F9 J     Do While (1)
     
' Z& h6 t! z7 f               Res = (MaxLim + MinLim) / 2
, |, ~) X5 f; p3 H; u. J6 P' m
             If Abs(QesFun(Res)) <= 10 ^ -12 Then
; V$ I, j% k$ S, a
                        SolFunDic = Res: Exit Do

- x" S/ q1 e- F3 C  O0 E. o              ElseIf (QesFun(Res) > 0) Then
# c' |$ m/ e; F$ g8 S# V5 c" k
8 u# @1 U' \" W+ [                        MinLim = Res
0 w' [  f3 e' h" r" L" ~
5 S: R$ M. ~2 `( q              Else
3 `* f. y+ M8 T. p$ u
0 J- Y  U% R# A7 N6 l. S: x                         MaxLim = Res

            End If
            
      Loop
2 x$ n- t8 J& V8 K6 ]; Z      
8 u. |! o1 n# j1 G3 \! M5 KEnd If
  `. Y* D- O7 UEnd Function
--------------------------------------------------------------
2.        牛顿法
2.1     由      f(AC)＝arctan(AC/80)*80+1-AC
        求导     f(AC)’=1/(1+(AC/80)^2)-1
1 n3 a, S5 p% q' _' K- k! l        即      AC_1=AC_0- f(AC)/ f(AC)’
--------------------------------------------------------------
2.2定义迭代函数：
Public 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)
! J2 ~) d8 V) Y' V* ?& C              
End Function
---------------------------------------------------------------
, T: x* ]8 t- S/ ]  K- p2.3定义求解函数
Public Function SolFunNew(ByVal IniAC As Double) As Double
. N- w$ k0 u2 p; {; O- z
Dim Res#

Do While (1)

                  Res = QesFunNew(IniAC)
: k6 P. b! S/ S4 K: s+ E
$ _7 q6 {1 z8 z3 v6 D8 z  L( p                   If Abs(QesFun(Res)) <= 10 ^ -12 Then
3 R. ]. I1 U9 k" I1 N: M' X- i5 r  
6 h+ B/ ]0 ]4 H7 N: r& \$ W; Z- p+ L                                   SolFunNew = Res: Exit Do
  
                   Else
                  
7 S( g7 B- }: O* v1 [5 h                                     IniAC = Res
                  
                  End If
1 y+ R5 i9 @4 H  q
' c5 Z5 Y* J) O- P# KLoop
----------------------------------------------------------------------------------------------------------

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