間違いしかありません.コメントにてご指摘いただければ幸いです(気が付いた点を特に断りなく頻繁に書き直していますのでご注意ください).

状態空間表現,バネマスダンパー系,状態行列の固有ベクトル

状態空間表現

  {ddtx=Ax+BuA:,B:,x:,u:y=Cx+DuC:,D:,y:

バネマスダンパー系(入力なし)

{ddt[xv]=[01kmcm][xv]k:,m:,c:,x:,v=dxdt:=[01ω022γ][xv]γ=c2m,ω0=km[y]=[10][xv]

固有値λ1,2

λ1,2=γ±γ2ω02γ=c2m,ω0=km=γ±ω0ζ21ζ=γω0=γ±ξξ=ω0ζ21=γ2ω02,

固有値λ1に対応する固有ベクトルv1

(λ1IA)v1=0([λ100λ1][01ω022γ])v1=0[λ11ω02λ1+2γ]v1=0[(γ+ξ)1ω02(γ+ξ)+2γ]v1=0λ1=γ+ξ[γ+ξ1ω02γ+ξ]v1=0 {(γ+ξ)v11v12=0ω02v11+(γ+ξ)v12=0 v11=t(γ+ξ)tv12=0v12=(γ+ξ)tv12=(γ+ξ)t v1=[v11v12]=t[1γ+ξ]
2番目の式で確認 ω02v11+(γ+ξ)v12=0ω02t+(γ+ξ)v12=0v11=t(γ+ξ)v12=ω02tv12=ω02γ+ξt=(γ2ξ2)γ+ξtξ=γ2ω02,ξ2=γ2ω02,ω02=γ2ξ2=(γ+ξ)(γξ)γ+ξt=(γξ)t=(γ+ξ)t


固有値λ2に対応する固有ベクトルv2

(λ2IA)v2=0([λ200λ2][01ω022γ])v2=0[λ21ω02λ2+2γ]v2=0[γξ1ω02(γξ)+2γ]v2=0λ2=γξ[γ+ξ1ω02γξ]v2=0 {(γξ)v21v22=0ω02v21+(γξ)v22=0 v21=t(γξ)tv22=0v22=(γξ)tv22=(γξ)t v2=[v21v22]=t[1γξ]
2番目の式で確認 ω02v21+(γξ)v22=0ω02t+(γξ)v22=0v21=t(γξ)v22=ω02tv22=ω02γξt=(γ2ξ2)γξtξ=γ2ω02,ξ2=γ2ω02,ω02=γ2ξ2=(γ+ξ)(γξ)γξt=(γ+ξ)t=(γξ)t

0 件のコメント:

コメントを投稿