作者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