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

状態空間表現,バネマスダンパー系,状態ベクトルの変換

状態空間表現

{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]

状態ベクトルの変換

x=T1x=12ξ[γξ1γξ1][xv]T1=12ξ[γξ1γξ1]ξ=γ2ω02=12ξ[(γξ)xv(γξ)x+v]=[γxξxv2ξγxξx+v2ξ]=[γx+ξx+v2ξγx+ξxv2ξ]=[x2+v+γx2ξx2v+γx2ξ] x(0)=[x02+v0+γx02ξx02v0+γx02ξ]C1,C2
Tで戻して確認 Tx=[11γ+ξγξ][x2+v+γx2ξx2v+γx2ξ]=[1(x2+v+γx2ξ)+1(x2v+γx2ξ)(γ+ξ)(x2+v+γx2ξ)+(γξ)(x2v+γx2ξ)]=[x2+v+γx2ξ+x2v+γx2ξ(γ+ξ)x2ξξ+(γ+ξ)(v+γx)2ξ+(γξ)x2ξξ(γξ)(v+γx)2ξ]=[x2+x2γξx+ξ2xγv+ξvγ2x+γξx2ξ+γξxξ2x+γv+ξv+γ2x+γξx2ξ]=[xξv2ξ+ξv2ξ]=[xv]=x
Tx=TT1x=[11γ+ξγξ]12ξ[γξ1γξ1][xv]=12ξ[11γ+ξγξ][γξ1γξ1][xv]=[1001][xv]=[xv]=x

0 件のコメント:

コメントを投稿