作者kerent (MacGyver)
看板graduate
標題[教學]casio fx-991ms解一元四次以上的方程式
時間Sat Feb 26 18:07:56 2005
範例..張碩那一本自動控制5-14頁..71年台大電機
畫出以下的根軌跡
G(S)=K(S+6)/S(S+4)(S^2+4S+8)
Sol:
1.極點(pole)根零點(zero)
pole = 0,-2+i2,-2-i2,-4.....n=4
zero = -6....n=1
2.所以實軸上的根軌跡在..負無限大到負六與負四到零
3.漸進線與實軸的交點
sigmaA=-2/3
漸進線角度
thetaA=+60度,-60度,-180度
4.極點-2+2i的離開角
這也可以用casio fx-991ms計算
算法如下
(1) 按下mode鍵..轉到complex(複數模式)
mode->2
(2) 將-2+2i存為一個變數A
-2+2i->shift->sto->A
(3) 計算角度..叫出A的方式為..ALPHA->A
如果連上面這都不會的話..那你工程計算機是白買的..連說明書也不會看
shift->arg((A+6)/S(S+4)->=
恩..算出來不是答案是吧..還要扣掉90度..為什麼要扣90度..自己想吧
就可以得到phi1=-153.4349488-90=-243.4349488
所以離開角phiP=-63.4349488度
5.分離點..這就是重點了
如果是三次方或二次方的話casio fx-991ms可由內建程式計算
但是這題是四次方..所以你就掛了..慢慢try and error嗎?
是可以..如果你時間多的話..可以慢慢算..不過如果你買了casio fx-991ms
還慢慢算的話..那真的是笨蛋了..有工具可以用還不會用(當然不只這台可以計算)
以上是廢話
-----------------------------------
分離點要對G(S)微分分子為零可以得到
微分可得下式
3S^4+40S^3+168S^2+288S+192=0
四次方以上要用牛頓法求近似值
原理不解釋..以下是直接
這題由上面的的推論可以知道分離點在負無限大到負六與負四到零之間
所以我先假設這點為-10與-1
先將上式一次積分算出來..可得
12S^3+120S^2+336S+288=0
先算負無限大到負六之間的這一點(mode->1..先轉回一般形式)
-10->shift->sto->A
A-((3A^4+40A^3+168A^2+288A+192)/(12A^3+
120A^2+336A+288))
這邊的括號要注意..不小心會括錯..多一個或少一個
重點來了..上面那長長的式子在計算機上打完後->shift->sto->A
把它再存回A..就變成了一個迴圈
這時候螢幕上可以得到-8.66145833
這是牛頓法逼近第一次的結果..然後你只要再按下=這個按鍵
就會在逼近第二次..通常按個四五次就會得到一個不會再改變的答案..就是結果
以這題為例..會後回得到-7.297457982
與書上的-7.3相當接近..甚至更準
-----------------------------------
同理可得得到另一點..將初始A的值改為-1同樣的步驟便可以得到結果
注意..如果初始值沒設好..沒設好的意思是你故意去找沒解的地方
以這題來說你如果設A=-5就會有問題..設-100就沒問題
只是要多按一兩次=
-----------------------------------
6.剩下的不解了..如果還不會解接下來的地方的話..那你今年應該沒希望了....
-----------------------------------
結論:
1.此方法可以在四次方以上的方程式使用..不過步驟要熟練一點
上面的步驟習慣後應該可以在一兩分鐘內打完得到答案
2.這台計算機的功能不只這些
其實在算自動控制的時候有很多節省時間的小技巧..那本說明書真的很不錯
有空的話請多翻翻..技巧有待你們自己發覺..我在這邊懶的說了
3.以上方法是我自己發現使用的..覺得有用的話就推一下文吧
--
http://macgyver.demian.idv.tw/main.asp
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.120.25.225
1F:推 afei0912:狂推阿!!!!好文章~~~謝啦~~ 210.85.18.50 02/26
2F:推 afei0912:推暴~~!!希望這台考試可以帶進去 控制之寶 210.85.18.50 02/26
3F:推 cool1985:如果還不會解接下來的地方的話那你今年應該沒望 218.170.61.88 02/26
4F:推 lcjk:真感謝你...將那麼好的方法告訴大家,我算B.C Kuo 61.230.115.93 02/26
5F:→ lcjk:全卡在求分離點....Orz 61.230.115.93 02/26
※ 編輯: kerent 來自: 140.120.25.225 (02/26 19:17)
6F:推 kimuralin:這一台應該是不能帶進考場吧...140.115.205.173 02/26
7F:推 lcjk:可以吧..之前我麼適用5500,但他的功能太強大了,考 61.230.115.93 02/26
8F:→ lcjk:試過的學長都推這一台 61.230.115.93 02/26
9F:推 taitel:推991 謝謝原作者大公無私 219.68.110.131 02/26
10F:推 afei0912:再推原po 讓我靠著計算機上準台大 Orz 210.85.18.50 02/26
11F:推 newlly:好心動...不過真的可以用嗎? 222.157.53.52 02/26
12F:推 aaaaa99999:好讚的文 感謝分享~ 219.84.28.122 02/26
13F:推 aaaaa99999:這個方法350也可用真是太讚辣 350大都可帶巴 219.84.28.122 02/26
14F:推 llat:感謝.這個方法好像也可以求三個極點的相位教越頻率 61.231.80.56 02/26
15F:推 devilstick:純推不下,大感激 211.74.216.139 02/26
16F:推 PerWang:屋屋屋...Orz...我是笨蛋.... 61.228.89.154 02/27
17F:推 chaosteam:嗚嗚 太感動了啦 感謝經驗分享 211.73.187.1 02/27
18F:推 guteres:如果有人再用這台我就直接向監考老師舉發y 59.104.43.26 04/18
※ lpissonice:轉錄至某隱形看板 12/13 18:01
19F:推 jason77626:受益良多 前輩分離點這招真的太強大了 02/13 01:06
20F:推 tkufc: 受益推 42.70.39.229 02/02 02:39