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

状態空間表現,バネマスダンパー系,状態行列の対角化

状態空間表現

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

状態行列の対角化(基底の変換)とそれに伴う出力行列の変換

A¯=T1AT=12ξ[γξ1γξ1][01ω022γ][11γ+ξγξ]T=[11γ+ξγξ]T1=12ξ[γξ1γξ1]ξ=γ2ω02,ξ2=γ2ω02,ω02=γ2ξ2=12ξ[(γ+ξ)1γξ1][01(γ2ξ2)2γ][11γ+ξ(γ+ξ)]=12ξ[(γ+ξ)1γξ1][01+1(γ+ξ)01+1(γ+ξ)(γ2ξ2)1+(2γ)(γ+ξ)(γ2ξ2)1+(2γ)(γ+ξ)]=12ξ[(γ+ξ)1γξ1][γ+ξ(γ+ξ)γ2+ξ2+2γ22γξγ2+ξ2+2γ2+2γξ]=12ξ[(γ+ξ)1γξ1][γ+ξ(γ+ξ)ξ2+γ22γξξ2+γ2+2γξ]=12ξ[(γ+ξ)1γξ1][γ+ξ(γ+ξ)(γξ)2(γ+ξ)2]=12ξ[(γ+ξ)(γ+ξ)+(1)(γξ)2(γ+ξ){(γ+ξ)}+(1)(γ+ξ)2(γξ)(γ+ξ)+(1)(γξ)2(γξ){(γ+ξ)}+(1)(γ+ξ)2]=12ξ[(ξ2γ2)(γξ)2(γ+ξ)2(γ+ξ)2(γξ)2+(γξ)2(γ2ξ2)+(γ+ξ)2]=12ξ[ξ2+γ2(γ22γξ+ξ2)00γ2+ξ2+(γ2+2γξ+ξ2)]=12ξ[ξ2+γ2γ2+2γξξ200γ2+ξ2+γ2+2γξ+ξ2]=12ξ[2γξ2ξ2002γξ+2ξ2]=[2γξ2ξ22ξ002γξ+2ξ22ξ]=[γ+ξ00γξ]=[λ100λ2]C¯=CT=[10][11γ+ξγξ]T=[11γ+ξγξ]=[11+0(γ+ξ)11+0(γξ)]=[11]

0 件のコメント:

コメントを投稿