機(jī)械社區(qū)

 找回密碼
 注冊(cè)會(huì)員

QQ登錄

只需一步,快速開(kāi)始

搜索
樓主: 子子61961
打印 上一主題 下一主題

在日本做機(jī)械設(shè)計(jì)的感受23:一道計(jì)算題的回帖

  [復(fù)制鏈接]
11#
發(fā)表于 2016-2-15 12:48:43 | 只看該作者
關(guān)注
12#
發(fā)表于 2016-2-15 14:51:40 | 只看該作者
非常感謝你的分享,上來(lái)就學(xué)了一手單變量求解。!0.329784444989256
13#
發(fā)表于 2016-2-15 15:05:59 | 只看該作者
謝謝分享
14#
發(fā)表于 2016-2-15 15:16:40 | 只看該作者
大俠說(shuō)到的回帖方式就是一種網(wǎng)絡(luò)交流禮儀啊。如果大家都能注意,大家都會(huì)受益。
" |/ D' @( m! g6 h! V也謝謝大俠的單變量求解。
15#
發(fā)表于 2016-2-15 15:27:18 | 只看該作者
大俠,學(xué)習(xí)了。
" ]4 a- |# j1 ?- n2 l+ |關(guān)于中國(guó)文化,您提到的事情,實(shí)實(shí)在在存在著,也一直困擾著我。起初我一直懷疑是管理的問(wèn)題,現(xiàn)在明白了一點(diǎn),太深了。謝謝。
' I% c& B) n1 f& A9 |$ J$ D6 c溫習(xí)并熟練了單變量求解。. C! K& R9 d$ u# K
附圖為我的計(jì)算結(jié)果:精度好像不太高,正在想方法解決。

單變量求解(方法).JPG (51.36 KB, 下載次數(shù): 77)

單變量求解(方法).JPG

評(píng)分

參與人數(shù) 1威望 +1 收起 理由
子子61961 + 1 謝謝您的認(rèn)真回復(fù)。

查看全部評(píng)分

16#
 樓主| 發(fā)表于 2016-2-15 15:47:01 | 只看該作者
點(diǎn)滴積累 發(fā)表于 2016-2-15 16:27 $ d! |( }3 v$ Z' D+ w
大俠,學(xué)習(xí)了。
. Q% b# W; |! ?2 T: \  B關(guān)于中國(guó)文化,您提到的事情,實(shí)實(shí)在在存在著,也一直困擾著我。起初我一直懷疑是管理的問(wèn) ...

3 |! {) f  m- y% ^6 F: j: ^0 {感謝您的認(rèn)真回帖。8 f# E! _, z8 z' m
關(guān)于管理方面的問(wèn)題,其實(shí)可以考慮為兩個(gè)方面。9 k* y9 v. u1 B1 p  ^5 F( t
一方面是從上向下的監(jiān)視和修正,
8 w) k2 c! f# k  C# G一方面是從下層開(kāi)始的自省和改善。9 @0 b2 I% `2 y
中國(guó)提到的管理,往往更偏向于前者,
2 q9 ^: [1 H- U1 Z8 ]9 y9 q而日本的管理,其實(shí)是后者做得更好一些,2 X) e( w" ?  h2 }( e# n1 x
所以在日本,從上層向下的管理會(huì)輕松一些。
$ V$ o- l. v3 v" K' s9 v  \以前在第二話的文章里面寫(xiě)過(guò)一點(diǎn),以后考慮也繼續(xù)寫(xiě)一些。
( K- @" Y$ H, u
$ T. k4 z0 I: G: @關(guān)于單變量求解,它有計(jì)算的設(shè)定。
: o4 ^! k  u" Z* _在Excel的菜單里面,工具--〉選項(xiàng)--〉計(jì)算方法--〉反復(fù)計(jì)算
8 p% K1 _; _5 `1 D- w1 W
5 s$ d6 c8 [" B  B! I
. f6 E# g0 w) U! l  [(抱歉手邊現(xiàn)在只有日文電腦,可以看里面的漢字參考)- o2 \" h3 U; c# {" x2 x
) {/ s" O, e3 T* i# O$ X
里面有反復(fù)計(jì)算的設(shè)定值。最大反復(fù)次數(shù)100次,變化的最大值 0.001。0 c  s2 M3 P$ ]' j0 x6 x
將這個(gè)次數(shù)改大一些,或者變化值改的更小一些,則會(huì)得到更精確一些的結(jié)果。# D/ V9 f5 _# y; z' ]+ ^

6 T2 M; Q! u; c+ K" q6 v另外,如何選取可變單元格也是一個(gè)因素。
* `; {  ^8 |1 \& l3 O7 `; x如果選擇使用弧度,就是在0.32左右每次變化0.001進(jìn)行求解。
# `% `7 y0 h8 @, w: R! n, u1 r如果選擇使用角度,就是在18.8左右每次變化0.001進(jìn)行求解。. o( E6 L. _, W4 e; K0 |4 \3 z4 w
其收斂的速度是不同的。
) _3 a& k; z! E- u. r

點(diǎn)評(píng)

500,500,500,300,200  發(fā)表于 2016-2-18 13:19
17#
發(fā)表于 2016-2-15 16:06:37 | 只看該作者
子子61961 發(fā)表于 2016-2-15 15:47
" \3 P0 t0 m1 T4 i3 o+ A感謝您的認(rèn)真回帖。- w, u+ q0 g8 \1 D! i
關(guān)于管理方面的問(wèn)題,其實(shí)可以考慮為兩個(gè)方面。9 y  Y9 \' a+ w# a3 ]7 h1 y
一方面是從上向下的監(jiān)視和修正,

7 n, ~  G# X) c) k$ y/ s' U1 s6 W大俠,中文版,呵呵

單變量求解(方法)精度問(wèn)題.JPG (111.1 KB, 下載次數(shù): 83)

單變量求解(方法)精度問(wèn)題.JPG

點(diǎn)評(píng)

THANKS~  發(fā)表于 2016-2-15 16:23
18#
發(fā)表于 2016-2-15 17:06:51 | 只看該作者
本帖最后由 andyany 于 2016-2-15 17:08 編輯
/ C0 c; t5 V; X% _
7 r% B) k' x, w. i; S- ]據(jù)我所碰到的,國(guó)內(nèi)的交流能力都很低。
: n# H& V& w2 ~* q4 a0 q1, w& m0 ^8 ^# ^4 {) H! i
當(dāng)你問(wèn)問(wèn)題時(shí),其實(shí)對(duì)方也不會(huì),但他會(huì)說(shuō)”這你都不會(huì)?!“ 讓人備受羞辱。; m0 k( s. L. E# h. T/ ]
但有段時(shí)間和老外一起工作,老外不會(huì)時(shí),一般也不會(huì)說(shuō)”我不會(huì)“導(dǎo)致你干等或自己一通找,而是說(shuō)”你找xxx,他會(huì)幫助你“,給你指路。讓人倍感溫暖。
# J4 E2 ]5 S3 s9 C2
* M; f& O, q2 ^你覺(jué)得國(guó)內(nèi)領(lǐng)導(dǎo)應(yīng)該經(jīng)歷過(guò),于是去求助。領(lǐng)導(dǎo)說(shuō)“這個(gè)你應(yīng)該自己搞定”,你也沒(méi)話說(shuō)。而你若是求助外方的領(lǐng)導(dǎo),即便很大的官,也會(huì)給出你答案甚至理由和引申,讓你佩服得不要不要的。

評(píng)分

參與人數(shù) 1威望 +1 收起 理由
子子61961 + 1 謝謝分享!

查看全部評(píng)分

19#
發(fā)表于 2016-2-15 19:15:28 | 只看該作者
我想樓主想說(shuō)的是 不要想著什么工具 軟件 用筆頭算吧! 這是一個(gè)做設(shè)計(jì)人員對(duì)基本定理規(guī)律的理解 這是根本 不要太依賴(lài)工具軟件 否則就失去了根本的理論基礎(chǔ) 這樣談何開(kāi)發(fā)
20#
發(fā)表于 2016-2-15 19:51:33 | 只看該作者
AC-arctan(AC/80)*80=1是數(shù)學(xué)問(wèn)題。
2 B; \# A& G% K% N子子大俠既然讀過(guò)VB,那么讀VBA代碼就應(yīng)該沒(méi)什么困難,核心是數(shù)值計(jì)算。3 e! r3 \: s0 J3 S. r
---------------------------------------------------------------------------------------------------------------------------------------
/ p% H) o( E9 {$ V3 A- T7 q" n定義待求解函數(shù):' G. e/ |) B" }
Public Function QesFun(ByVal Var_AC As Double) As Double" r3 u3 ]3 H' z' M7 ~
# \. a$ M- O9 T" N4 `: P8 u- z
               QesFun = Var_AC - Atn(Var_AC / 80) * 80 - 15 |% Y% J6 L& Y+ Y0 g) n
- j3 t- \5 b# ?
End Function
4 f* H( Y- T/ ^# I: i--------------------------------------------------------------------------------------------------
4 P: [: b4 z' j) p% J0 X1.        二分法
! D. A3 q7 f! S. @( C) Q. P1.1  由      Arctan(AC/80)=(AC-1)/80
5 h: R+ g& u" n      知      -PI()/2<(AC-1)/80< PI()/2. [, m/ g6 W; v* q/ ]% O
     即     1 -80*PI()/2<AC<1+ 80*PI()/2+ c5 I6 k' N: J+ p
++++++++++++++++++++++++++++++++++++++
/ Q" _5 f8 c9 P1.2定義求解函數(shù):
9 ]8 i: r. e- }- lPublic Function SolFunDic(ByVal MaxLim As Double, ByVal MinLim As Double) As Double
$ }. H' x1 l6 k' k; S. A7 ]2 H/ _& Q9 A9 X) A+ ]
Dim Res#, VarAdj#
( d/ t+ D3 g  f  g6 R1 `! [# K! \& f) H, V! m1 I
VarAdj = 10 ^ -6/ |# P4 s" E+ M  N
$ \. Z  i# ?+ Z
If QesFun(MinLim + VarAdj) < QesFun(MaxLim - VarAdj) Then
! j" \# [# r0 e9 _  U6 I9 `/ j2 q$ _9 R& F; ]
      Do While (1)
4 B. s( G5 t, \% V: D6 I( E* [
; S% |7 \3 _: S4 W              Res = (MaxLim + MinLim) / 27 J9 B4 L3 u& |* h
6 O1 \' y7 I2 `( v( c; }
               If Abs(QesFun(Res)) <= 10 ^ -12 Then
' v/ n% J. x3 x, k! c" D  Z( B
' ~. s) r, K$ o6 e# `' x                       SolFunDic = Res: Exit Do+ e9 e8 k& w, x. d1 K3 f* L5 ~
: L6 ~. x2 Z5 L( b( y
              ElseIf (QesFun(Res) < 0) Then
7 Y7 p0 H* E  m2 l+ o7 ]8 E" E, v2 C/ R9 B! j9 u' G
                      MinLim = Res1 V+ j. Y6 M/ y4 d) o

* r. @8 r( y* b1 }4 d7 K             Else0 P+ ], f3 M; [8 K

" Y* l7 s' [- W4 `0 N9 `                     MaxLim = Res, N2 K# U; B  q- i+ K

/ j# y8 {3 N+ v; o, J! r             End If8 p  k* e1 P: [+ G9 g- h
            
- ^2 D% [7 \; G7 x6 N      Loop
% C% d3 A* J. Z( q5 S4 V2 G2 {
- G, l/ p& e9 ~* ?7 hElse' i% D) k% k% d$ ~/ T% Z; {) h
+ G# ^9 c: Q* l6 F0 x- O7 _
     Do While (1)" m2 H! f, I' W5 p2 R
     % x5 F2 T, }" H2 M. j! N
               Res = (MaxLim + MinLim) / 2
8 h2 k! P8 A0 W3 S
7 u0 P& I/ C" G4 e             If Abs(QesFun(Res)) <= 10 ^ -12 Then
3 {' Y/ r' J7 E7 P/ ]1 Q) l8 K" _4 Z! Q
                        SolFunDic = Res: Exit Do+ X* f2 H: l0 P% U% S% Y

$ V# i. R! b: [& N8 U2 l% ^              ElseIf (QesFun(Res) > 0) Then
& B/ ]$ j3 f' @# `% x6 B0 ~
6 s' D  A( e- V2 c! |+ \                        MinLim = Res
4 o2 A  k8 `$ p' ~4 o# ^5 d) O) |" Z1 f! y& \) T2 h
              Else; ?5 i( ~: p  Q' W7 N) R7 R
2 ~/ w; c5 k$ [3 E
                         MaxLim = Res
- r& d' o. ?' m9 G
/ W4 p7 p; j. [+ w            End If1 ~& @% M7 m/ B' h& x% y+ ~
            
, U3 F+ j$ _, X% J$ ~      Loop2 O) t4 _1 @+ t* ]& F1 i
      * Q% I* B+ r0 g8 U0 H1 O+ Q
End If/ `# ?4 z+ r( F* ^% T" n9 F, S
End Function* Q4 k6 M" `( E. o! n
--------------------------------------------------------------  r% y( g4 [. Q
2.        牛頓法
. M; v( x% i8 t- f% U% B2.1     由      f(AC)=arctan(AC/80)*80+1-AC  u. y- F8 V+ L0 z: |8 Z
        求導(dǎo)     f(AC)’=1/(1+(AC/80)^2)-1$ o$ [3 D4 N% u
        即      AC_1=AC_0- f(AC)/ f(AC)’
3 ~; e& R4 K9 ?8 r( w--------------------------------------------------------------
" Z. m" Y. ]2 i; K0 t; ?2.2定義迭代函數(shù):; G. i1 h" o/ @8 u8 Z  G3 t
Public Function QesFunNew(ByVal Var_AC As Double) As Double: {+ y; D$ k4 F  R2 [9 g# q2 V

, n5 w4 h1 b% M0 K4 L7 C4 n              QesFunNew = Var_AC - (Atn(Var_AC / 80) * 80 + 1 - Var_AC) / (1 / (1 + (Var_AC / 80) ^ 2) - 1)
6 h+ Z) {  w) P* d+ G7 w- ]              
3 [' z8 h# Q6 h( V# M# iEnd Function5 V! n6 G! \) T* r  R
---------------------------------------------------------------
  t: D$ X0 P' {* k. e4 k$ w! j2.3定義求解函數(shù): r& ^( W; q! z5 B  N9 @! n
Public Function SolFunNew(ByVal IniAC As Double) As Double
. H% H5 f4 {, Z8 |/ X. t+ G/ R; V7 u0 l) A! v' i7 T) }1 [
Dim Res#
' u* `5 q1 x0 Y6 m+ m8 ^* C# [2 G. v2 U3 T" M. D, |
Do While (1), g1 i* Z( n" r7 ?1 V/ p
4 r+ d' v: H$ S6 J0 S
                  Res = QesFunNew(IniAC)7 J- p. ]/ }2 d5 J
9 S. I1 Y; v! u' P
                   If Abs(QesFun(Res)) <= 10 ^ -12 Then4 v6 E7 f: u0 }
  
5 z3 q- ]- r( F$ U                                   SolFunNew = Res: Exit Do
" k" p( a% v& Z8 e0 O- e; R  : Y4 P6 {4 z+ E' @9 A, y' R3 U
                   Else) Y' a1 @( I, h# N) A7 c# w4 A" b
                  ) n' A- H5 d% M6 L, U
                                     IniAC = Res- z5 W: e2 s; x& n; r" D, J
                  
" i# W; c7 k+ p2 I" P: x                  End If, Y% y. P1 N! B0 c
: M: U6 l6 E6 H6 p
Loop8 l) K+ L& d# V
----------------------------------------------------------------------------------------------------------
+ C) Z5 T* J  s
) c% D0 G+ w$ o/ a' S6 k: d  t( ?" L7 `6 v這樣做可能有點(diǎn)麻煩,但涉及到循環(huán),迭代時(shí),可自由調(diào)用自定義VBA函數(shù)和工作表函數(shù)(矩陣計(jì)算連桿機(jī)構(gòu)),還可控制輸出表格,便于插圖。計(jì)算冷卻塔時(shí),積分得用辛普遜;解汽水比,得解非線性方程。對(duì)這些問(wèn)題,EXCEL 工作表自身好像很吃力。( r/ |& _" N0 J! o
: r2 L1 v, ]/ F0 j! T2 \

點(diǎn)評(píng)

我的天啦,我的圈子里沒(méi)有一個(gè)會(huì)寫(xiě)這樣完整計(jì)算程序的人。  發(fā)表于 2021-5-10 17:28
佩服佩服!  發(fā)表于 2016-2-22 10:14

評(píng)分

參與人數(shù) 3威望 +3 收起 理由
mbampa1234 + 1 思想深刻,見(jiàn)多識(shí)廣!
一朵嬌粉嫩的花 + 1 不覺(jué)明歷,俺當(dāng)時(shí)學(xué)Fortran 90就是0分。。
子子61961 + 1 專(zhuān)業(yè)精湛,謝謝分享!

查看全部評(píng)分

本版積分規(guī)則

小黑屋|手機(jī)版|Archiver|機(jī)械社區(qū) ( 京ICP備10217105號(hào)-1,京ICP證050210號(hào),浙公網(wǎng)安備33038202004372號(hào) )

GMT+8, 2024-9-23 03:32 , Processed in 0.060778 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4 Licensed

© 2001-2017 Comsenz Inc.

快速回復(fù) 返回頂部 返回列表