Scilab-Branch-5.3-GIT
- Scilab help
- CACSD
- abcd
- abinv
- arhnk
- arl2
- arma
- arma2p
- armac
- armax
- armax1
- arsimul
- augment
- balreal
- bilin
- black
- bode
- bstap
- cainv
- calfrq
- canon
- ccontrg
- chart
- cls2dls
- colinout
- colregul
- cont_frm
- cont_mat
- contr
- contrss
- copfac
- csim
- ctr_gram
- dbphi
- dcf
- ddp
- des2ss
- des2tf
- dhinf
- dhnorm
- dscr
- dsimul
- dt_ility
- dtsi
- equil
- equil1
- evans
- feedback
- findABCD
- findAC
- findBD
- findBDK
- findR
- findx0BD
- flts
- fourplan
- frep2tf
- freq
- freson
- fspecg
- fstabst
- g_margin
- gainplot
- gamitg
- gcare
- gfare
- gfrancis
- gtild
- h2norm
- h_cl
- h_inf
- h_inf_st
- h_norm
- hallchart
- hankelsv
- hinf
- imrep2ss
- inistate
- invsyslin
- kpure
- krac2
- lcf
- leqr
- lft
- lin
- linf
- linfn
- linmeq
- lqe
- lqg
- lqg2stan
- lqg_ltr
- lqr
- ltitr
- m_circle
- macglov
- markp2ss
- minreal
- minss
- mucomp
- narsimul
- nehari
- nicholschart
- noisegen
- nyquist
- obs_gram
- obscont
- observer
- obsv_mat
- obsvss
- p_margin
- parrot
- pfss
- phasemag
- ppol
- prbs_a
- projsl
- reglin
- repfreq
- ric_desc
- ricc
- riccati
- routh_t
- rowinout
- rowregul
- rtitr
- sensi
- sgrid
- show_margins
- sident
- sm2des
- sm2ss
- sorder
- specfact
- ss2des
- ss2ss
- ss2tf
- st_ility
- stabil
- svplot
- sysfact
- syssize
- tf2des
- tf2ss
- time_id
- trzeros
- ui_observer
- unobs
- zeropen
- zgrid
- nyquistfrequencybounds
Please note that the recommended version of Scilab is 2025.0.0. This page might be outdated.
See the recommended documentation of this function
evans
エバンス根軌跡
呼び出し手順
evans(H [,kmax])
パラメータ
- H
リスト (線形システム
syslin
)- kmax
実数 (プロットに用いる最大ゲイン )
説明
状態空間または伝達形式 H(s)の線形システム
(syslin
リスト)のエバンス根軌跡を指定します.
これは,複素平面における
1+k*H(s)=1+k*N(s)/D(s)
の根の軌跡です.
選択された点のゲインk <= kmax
に関して,
D(s)+k*N(s)
の根の虚部が実部に対してプロットされます.
根軌跡の指定した点におけるゲインを得るには,
以下の命令を実行します:
k=-1/real(horner(h,[1,%i]*locate(1)))
を実行し,
根軌跡の好きな点をクリックします.
選択された点の座標が,
実数の2 x 1 ベクトルP=locate(1)
の場合,
このk
は
方程式 k*N(w) + D(w) =0
を
w=P(1)+%i*P(2)=[1,%i]*P
の条件の元で
解いたものです.
例
H=syslin('c',352*poly(-5,'s')/poly([0,0,2000,200,25,1],'s','c')); evans(H,100) P=3.0548543 - 8.8491842*%i; //P=selected point k=-1/real(horner(H,P)); Ns=H('num');Ds=H('den'); roots(Ds+k*Ns) //contains P as particular root // Another one clf();s=poly(0,'s');n=1+s; d=real(poly([-1 -2 -%i %i],'s')); evans(n,d,100); // clf();n=real(poly([0.1-%i 0.1+%i,-10],'s')); evans(n,d,80);
<< equil1 | CACSD | feedback >> |