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

状態空間表現,バネマスダンパー系,対角正準形式

状態空間表現

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

対角正準形式

{ddtx=A¯xy=C¯x {ddt[x2+v+γx2ξx2v+γx2ξ]=[γ+ξ00γξ][x2+v+γx2ξx2v+γx2ξ]y=[11][x2+v+γx2ξx2v+γx2ξ]

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

状態空間表現

{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

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

状態空間表現

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

状態空間表現,バネマスダンパー系,状態行列の基底の変換行列

状態空間表現

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

基底の変換行列[v1v2]

T=[v1v2]=[[v11v12][v21v22]]=[v11v21v12v22]=[11γ+ξγξ]v1=[v11v12]=t[1γ+ξ]v2=[v21v22]=t[1γξ]ξ=γ2ω02,

基底の変換行列の逆行列

T1=1|T|[v22v21v12v11]A=[abcd],A1=1|A|[dbca]=11(γξ)1(γ+ξ)[γξ1(γ+ξ)1]=1γξ+γξ[γξ1γξ1]=12ξ[γξ1γξ1]

状態空間表現,対角正準形式

状態空間表現

{ddtx=Ax+Buy=Cx+Du

対角正準形式

状態ベクトルxが正則行列Txを用いて表せるとする. x=Tx これを代入すると状態空間表現は以下のようになる. {ddtTx=ATx+Buy=CTx+Du この時,一つ目の式に対してTの逆行列であるT1を左から掛けると以下のようになる. {ddtT1Tx=T1ATx+T1Buy=CTx+Du T1T=Iなので以下のようになる. {ddtx=T1ATx+T1Buy=CTx+Du 各行列部をA¯,B¯,C¯で置き直すことで以下のようになる. {ddtx=A¯x+B¯uA¯=T1AT,B¯=T1By=C¯x+DuC¯=CT この時にTとしてAの固有ベクトル(基底)を並べた行列(変換行列)を用いればT1AT=A¯は対角化され,対角行列となり,このような状態空間表現を対角正準形と呼ぶ. T=[v1v2vn] 対角正準形では状態ベクトルxは状態ベクトルxへと変換されている. x=TxT1x=T1Txx=T1x

状態空間表現,バネマスダンパー系(入力なし)

状態空間表現

{ddtx=Ax+BuA:,B:,x:,u:y=Cx+DuC:,D:,y: {ddt[x1x2]=[a11a12a21a22][x1x2]+[b11b12b21b22][u1u2]2[y1y2]=[c11c12c21c22][x1x2]+[d11d12d21d22][u1u2]

バネマスダンパー系

mx¨+cx˙+kx=0k:,m:,c:,x:,v=dxdt:d2xd2t+cmdxdt+kmx=0 d2xd2t=cmdxdtkmxdvdt=cmvkmxv=dxdtd2xd2t=ddtddtx=ddtv {ddt[xv]=[01kmcm][xv]k:,m:,c:[y]=[10][xv] 入力行列Bと直達行列Dはともに0とする(入力なし).

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

状態空間表現

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

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

{ddt[xv]=[01kmcm][xv]k:,m:,c:,x:,v=dxdt:[y]=[10][xv]

状態行列の固有値

|λIA|=0|λ01kmλ(cm)|=λ{λ(cm)}(km)=λ(λ+cm)+km=λ2+cmλ+kmλ2+cmλ+km=0λ1,2=cm±(cm)24km2=c2m±c24m2km=γ±γ2ω02γ=c2m,ω0=km=γ±ω0ω0γ2ω02=γ±ω01ω02(γ2ω02)=γ±ω0γ2ω02ω02ω02=γ±ω0γ2ω021=γ±ω0(γω0)21=γ±ω0ζ21ζ=γω0=γ±ξξ=ω0ζ21

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

状態空間表現

  {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

e^xとマクローリン展開と双曲線凾数

exとマクローリン展開と双曲線凾数

ex=1+x+12!x2+13!x3+14!x4+15!x5+16!x6+17!x7+18!x8+=(1+12!x2+14!x4+16!x6+18!x8+)+(x+13!x3+15!x5+17!x7+)=cosh(x)+sinh(x)x=cosh(x)=1+12!x2+14!x4+16!x6+18!x8+sinh(x)=x+13!x3+15!x5+17!x7+=ex+ex2+exex2=ex+ex+exex2=2ex2=ex

バネマスダンパー系,過減衰の式からの他の式の導出

バネマスダンパー系

運動方程式

md2xd2t+cdxdt+kx=0m:,c:,d2xd2t+cmdxdt+kmx=0d2xd2t+2γdxdt+ω02x=0γ=c2m,ω0=km

単振動,振動減衰,臨界減衰,過減衰

γ=c2m,ω0=km,x0=x(0),v0=dx(t)dt|t=0
γ ζ=γω0 ξ=ω0ζ21 eγt[x0cosh(ξt)+v0+γx0ξsinh(ξt)] x(t)
単振動 0 0 ω0i e0t[x0cosh(ω0it)+v0+0x0ω0isinh(ω0it)] x0cos(ω0t)+v0ω0sin(ω0t)
振動減衰 0<γ<ω0 0<ζ<1 ωi
0<ω<1
eγt[x0cosh(ωit)+v0+γx0ωisinh(ωit)] eγt[x0cos(ωit)+v0+γx0ωsin(ωt)]
臨界減衰 ω0 1 0 eγt[x0(1+12!(ξt)2+)+v0+γx0ξ(ξt+13!(ξt)3+)] eγt[x0+(v0+γx0)t]
過減衰 ω0<γ 1<ζ ω
1<ω
eγt[x0cosh(ωt)+v0+γx0ωsinh(ωt)] eγt[x0cosh(ωt)+v0+γx0ωsinh(ωt)]

過減衰の式

x(t)=eγt[x0cosh(ξt)+v0+γx0ξsinh(ξt)]()

ω0<γ

x(t)=eγt[x0cosh(ξt)+v0+γx0ξsinh(ξt)]|1<γ,ξ=ωω0<γ,1<ζ,ξ=ω0ζ21=ωω=ω0|ζ21|=eγt[x0cosh(ωt)+v0+γx0ωsinh(ωt)]

0<γ<ω0

x(t)=eγt[x0cosh(ξt)+v0+γx0ξsinh(ξt)]|0<γ<1,ξ=ωi0<γ<ω0,0<ζ<1,ξ=ω0ζ21=ωiω=ω0|ζ21|=eγt[x0cosh(ωit)+v0+γx0ωisinh(ωit)]cosh(ix)=cos(x),sinh(ix)=isin(x)=eγt[x0cos(ωt)+v0+γx0ωiisin(ωt)]=eγt[x0cos(ωt)+v0+γx0ωsin(ωt)]

γ=0

x(t)=eγt[x0cosh(ξt)+v0+γx0ξsinh(ξt)]|γ=0,ξ=ω0iγ=0,ζ=0,ξ=ω001=ω0iω=ω0|ζ21|=e0t[x0cosh(ω0it)+v0+0x0ω0isinh(ω0it)]=x0cos(ω0t)+v0ω0iisin(ω0t)cosh(ix)=cos(x),sinh(ix)=isin(x)e0=1=x0cos(ω0t)+v0ω0sin(ω0t)

γω0

x(t)=limγω0,ξ0eγt[x0cosh(ξt)+v0+γx0ξsinh(ξt)]γω0,ζ1,ξ0=limγω0,ξ0eγt[x0(1+12!(ξt)2+)+v0+γx0ξ(ξt+13!(ξt)3+)]cosh(ξt)=10!(ξt)0+12!(ξt)2+,sinh(ξt)=11!(ξt)1+13!(ξt)3+=limγω0,ξ0eγt[x0(1+12!(ξt)2+)+(v0+γx0)(t+13!ξ2t3+)]=eω0t[x0(1)+(v0+ω0x0)(t)]=eω0t[x0+(v0+ω0x0)t]=eγt[x0+(v0+γx0)t]

sinh(x)のマクローリン展開

a点まわりのテイラー展開

f(x)=k=0f(k)(a)x!(xa)ka=10!f(0)(a)(xa)0+11!f(1)(a)(xa)1+12!f(2)(a)(xa)2+f(n)(x):f(x)n

マクローリン展開(0点まわりのテイラー展開)

f(x)=k=0f(k)(a)x!(xa)k|a=0=k=0f(k)(0)x!(x0)k=k=0f(k)(0)x!(x)k=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+

f(x)=sinh(x)のマクローリン展開(0点まわりのテイラー展開)

sinh(x)=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+=10!{sinh(0)}x0+11!{cosh(0)}x+12!{sinh(0)}x2+13!{cosh(0)}x3+14!{sinh(0)}x4+15!{cosh(0)}x5+16!{sinh(0)}x6+17!{cosh(0)}x7+18!{sinh(0)}x8+ddxcosh(x)=sinh(x),ddxsinh(x)=cosh(x)=10!0x0+11!1x+12!0x2+13!1x3+14!0x4+15!1x5+16!0x6+17!1x7+18!0x8+cosh(0)=e0+e02=1+12=1,sinh(0)=e0e02=112=0=11!x1+13!x3+15!x5+17!x7+=x+13!x3+15!x5+17!x7+

cosh(x)のマクローリン展開

a点まわりのテイラー展開

f(x)=k=0f(k)(a)x!(xa)ka=10!f(0)(a)(xa)0+11!f(1)(a)(xa)1+12!f(2)(a)(xa)2+f(n)(x):f(x)n

マクローリン展開(0点まわりのテイラー展開)

f(x)=k=0f(k)(a)x!(xa)k|a=0=k=0f(k)(0)x!(x0)k=k=0f(k)(0)x!(x)k=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+

f(x)=cosh(x)のマクローリン展開(0点まわりのテイラー展開)

cosh(x)=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+=10!{cosh(0)}x0+11!{sinh(0)}x+12!{cosh(0)}x2+13!{sinh(0)}x3+14!{cosh(0)}x4+15!{sinh(0)}x5+16!{cosh(0)}x6+17!{sinh(0)}x7+18!{cosh(0)}x8+ddxcosh(x)=sinh(x),ddxsinh(x)=cosh(x)=10!1x0+11!0x+12!1x2+13!0x3+14!1x4+15!0x5+16!1x6+17!0x7+18!1x8+cosh(0)=e0+e02=1+12=1,sinh(0)=e0e02=112=0=10!x0+12!x2+14!x4+16!x6+18!x8+=1+12!x2+14!x4+16!x6+18!x8+

cosh, sinhの微分

cosh(x)の微分

ddxcosh(x)=ddxex+ex2=12(ddxex+ddxex)ddxCf(x)=Cddxf(x)ddx(f(x)+g(x))=ddxf(x)+ddxg(x)=12(ddxex+ddxex)=12{ex+(1)ex}=12(exex)=sinh(x)

sinh(x)の微分

ddxsinh(x)=ddxexex2=12(ddxexddxex)ddxCf(x)=Cddxf(x)ddx(f(x)+g(x))=ddxf(x)+ddxg(x)=12(ddxexddxex)=12{ex(1)ex}=12(ex+ex)=cosh(x)

cosh(i x), sinh(i x) (純虚数に対するcosh, sinh)

cosh(ix)

cosh(x)=ex+ex2cosh(ix)=eix+eix2={cos(x)+isin(x)}+{cos(x)+isin(x)}2eix=cos(x)+isin(x)={cos(x)+isin(x)}+{cos(x)isin(x)}2cos(x)=cos(x),sin(x)=sin(x)=cos(x)+isin(x)+cos(x)isin(x)2=cos(x)+cos(x)2=2cos(x)2=cos(x)

sinh(ix)

sinh(x)=exex2sinh(ix)=eixeix2={cos(x)+isin(x)}{cos(x)+isin(x)}2eix=cos(x)+isin(x)={cos(x)+isin(x)}{cos(x)isin(x)}2cos(x)=cos(x),sin(x)=sin(x)=cos(x)+isin(x)cos(x)+isin(x)2=isin(x)+isin(x)2=2isin(x)2=isin(x)

cos(ix),sin(ix) (純虚数に対するcos,sin)

cos(ix),sin(ix) (純虚数に対するcos,sin)

e^(ix)のマクローリン展開

eixのマクローリン展開(0点まわりのテイラー展開)

eix=1+(ix)+12!(ix)2+13!(ix)3+14!(ix)4+15!(ix)5+16!(ix)6+17!(ix)7+18!(ix)8+ex=1+x+12!x2+13!x3+14!x4+15!x5+16!x6+17!x7+18!x8+=1+ix+12!i2x2+13!i3x3+14!i4x4+15!i5x5+16!i6x6+17!i7x7+18!i8x8+=1+ix+12!(1)x2+13!(i)x3+14!x4+15!ix5+16!(1)x6+17!(i)x7+18!x8+=1+ix12!x2i13!x3+14!x4+i15!x516!x6i17!x7+18!x8+=112!x2+14!x416!x6+18!x8++i(x13!x3+15!x517!x7+)=cos(x)+isin(x)cos(x)=112!x2+14!x416!x6+18!x8+sin(x)=x13!x3+15!x517!x7+

eを底とした指数凾数のマクローリン展開

a点まわりのテイラー展開

f(x)=k=0f(k)(a)x!(xa)ka=10!f(0)(a)(xa)0+11!f(1)(a)(xa)1+12!f(2)(a)(xa)2+f(n)(x):f(x)n

マクローリン展開(0点まわりのテイラー展開)

f(x)=k=0f(k)(a)x!(xa)k|a=0=k=0f(k)(0)x!(x0)k=k=0f(k)(0)x!(x)k=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+

f(x)=exのマクローリン展開(0点まわりのテイラー展開)

ex=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+=10!(e0)x0+11!(e0)x+12!(e0)x2+13!(e0)x3+14!(e0)x4+15!(e0)x5+16!(e0)x6+17!(e0)x7+18!(e0)x8+=10!1x0+11!1x+12!1x2+13!1x3+14!1x4+15!1x5+16!1x6+17!1x7+18!1x8+=10!x0+11!x1+12!x2+13!x3+14!x4+15!x5+16!x6+17!x7+18!x8+=1+x+12!x2+13!x3+14!x4+15!x5+16!x6+17!x7+18!x8+

sin(x)のマクローリン展開

a点まわりのテイラー展開

f(x)=k=0f(k)(a)x!(xa)ka=10!f(0)(a)(xa)0+11!f(1)(a)(xa)1+12!f(2)(a)(xa)2+f(n)(x):f(x)n

マクローリン展開(0点まわりのテイラー展開)

f(x)=k=0f(k)(a)x!(xa)k|a=0=k=0f(k)(0)x!(x0)k=k=0f(k)(0)x!(x)k=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+

f(x)=sin(x)のマクローリン展開(0点まわりのテイラー展開)

sin(x)=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+=10!{sin(0)}x0+11!{cos(0)}x+12!{sin(0)}x2+13!{cos(0)}x3+14!{sin(0)}x4+15!{cos(0)}x5+16!{sin(0)}x6+17!{cos(0)}x7+18!{sin(0)}x8+=10!0x0+11!1x+12!0x2+13!1x3+14!0x4+15!1x5+16!0x6+17!1x7+18!0x8+=11!x113!x3+15!x517!x7+=x13!x3+15!x517!x7+

cos(x)のマクローリン展開

a点まわりのテイラー展開

f(x)=k=0f(k)(a)x!(xa)ka=10!f(0)(a)(xa)0+11!f(1)(a)(xa)1+12!f(2)(a)(xa)2+f(n)(x):f(x)n

マクローリン展開(0点まわりのテイラー展開)

f(x)=k=0f(k)(a)x!(xa)k|a=0=k=0f(k)(0)x!(x0)k=k=0f(k)(0)x!(x)k=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+

f(x)=cos(x)のマクローリン展開(0点まわりのテイラー展開)

cos(x)=10!f(0)(0)x0+11!f(1)(0)x+12!f(2)(0)x2+=10!{cos(0)}x0+11!{sin(0)}x+12!{cos(0)}x2+13!{sin(0)}x3+14!{cos(0)}x4+15!{sin(0)}x5+16!{cos(0)}x6+17!{sin(0)}x7+18!{cos(0)}x8+=10!1x0+11!0x+12!1x2+13!0x3+14!1x4+15!0x5+16!1x6+17!0x7+18!1x8+=10!x012!x2+14!x416!x6+18!x8+=112!x2+14!x416!x6+18!x8+

二階微分のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

二階微分のラプラス変換

f(t)=d2g(t)d2tL[f(t)]=0f(t)estdt=0d2g(t)d2testdt=[dg(t)dtest]00dg(t)dt(s)estdtabf(t)g(t)dt=[f(t)g(t)]ababf(t)g(t)dt(f:f,g:g)=[dg()dtes{dg(t)dt|t=0}es0](s)0dg(t)dtestdt=[0g(0)]+s0dg(t)dtestdtdg(t)dt|t=0=g(0)=s0dg(t)dtestdtg(0)=s{sL[g(t)]g(0)}g(0)0dg(t)dtestdt=L[dg(t)dt]=sL[g(t)]g(0)=s2L[g(t)]sg(0)g(0)

バネマスダンパー系,γ,ω0による比較

バネマスダンパー系

γ,ω0による比較

x¨(t)+cmx˙(t)+kmx(t)=0d2d2tx(t)+cmddtx(t)+kmx(t)=0 γ=c2mω02=kmω2=|γ2ω02|=|ω02ω02(γ2ω02)|=ω02|(γω0)21|=ω02|ζ21|ζ=γω0ω=ω02|ζ21|=ω0|ζ21| x(t)=x0cos(ω0t)+v0ω0sin(ω0t)γ=0ζ=0=eγt[x0cos(ωt)+v0+γx0ωsin(ωt)]0<γ<ω00<ζ<1=eγt[x0+(v0+γx0)t]γ=ω0ζ=1=eγt[x0cosh(ωt)+v0+γx0ωsinh(ωt)]ω0<γ1<ζ

バネマスダンパー系,臨界減衰(γ=ω0)の場合

バネマスダンパー系

γ=ω0

λ1,2=γ±γ2ω02|γ=ω0λ1,2=γ±γ2ω02=ω0±ω02ω02=ω0 x(t)=C(t)eω0tx˙=C(t)eω0tω0C(t)eω0t(fg)=fg+fgx¨={C(t)eω0tω0C(t)eω0t}+{ω0C(t)eω0t+ω02C(t)eω0t}(fg)=fg+fg=C(t)eω0t2ω0C(t)eω0t+ω02C(t)eω0t x¨+2γx˙+ω02x|γ=ω0=0x¨+2ω0x˙+ω02x=0[C(t)eω0t2ω0C(t)eω0t+ω02C(t)eω0t]+2ω0[C(t)eω0tω0C(t)eω0t]+ω02C(t)eω0t=0C(t)eω0t2ω0C(t)eω0t+ω02C(t)eω0t+2ω0C(t)eω0t2ω02C(t)eω0t+ω02C(t)eω0t=0C(t)eω0t=0C(t)=0eω0t>0 C(t)=0C(t)=C1C(t)=C1t+C2x(t)=(C1t+C2)eω0t x(0)=(C1t+C2)eω0t|t=0=(C10+C2)eω00=C21=C2C2=x(0)=x0x(0)={C1eω0tω0(C1t+C2)eω0t}|t=0,C2=x0=C1eω00ω0(C10+x0)eω00=C11ω0(x0)1=C1ω0x0C1=x(0)+ω0x0=v0+ω0x0 x(t)=C(t)eω0t=(C1t+C2)eω0t={(v0+ω0x0)t+x0}eω0t={(v0+γx0)t+x0}eγtω0=γ=eγt{x0+(v0+γx0)t}=x0eγt(1+γt)+v0eγt(t)x0,v0

バネマスダンパー系,過減衰(ω0<γ)の場合

バネマスダンパー系

ω0<γ

λ1,2=γ±γ2ω02λ1,2=γ±γ2ω02=γ±ω2ω2=|γ2ω02|,γ2ω02=ω2=γ±ω C1,2=x02±v0+γx02γ2ω02|γ>ω0C1,2=x02±v0+γx02γ2ω02=x02±v0+γx02ω2=x02±v0+γx02ωω0 x(t)=C1eλ1t+C2eλ2t=C1e(γ+ω)t+C2e(γω)t=C1eγteωt+C2eγteωt=eγt{C1eωt+C2eωt}=eγt[{x02+v0+γx02ω}eωt+{x02v0+γx02ω}eωt]=eγt{x02eωt+v0+γx02ωeωt+x02eωtv0+γx02ωeωt}=eγt{x02(eωt+eωt)+v0+γx02ω(eωteωt)}=eγt{x0eωt+eωt2+v0+γx0ωeωteωt2}=eγt{x0cosh(ωt)+v0+γx0ωsinh(ωt)}cosh(ωt)=eωt+eωt2,sinh(ωt)=eωteωt2=x0eγt{cosh(ωt)+γωsinh(ωt)}+v0eγt{1ωsinh(ωt)}x0,v0=x0eγt{cosh(ωt)+ζζ21sinh(ωt)}+v0eγt{1ωsinh(ωt)}γω=γ|γ2ω02|=γω02ω02|γ2ω02|=γω02|γ2ω02ω02ω02|=γω0|ζ21|=ζ|ζ21|ω=|γ2ω02|,ζ=γω0

バネマスダンパー系,減衰振動(0<γ<ω0)の場合

バネマスダンパー系

0<γ<ω0

λ1,2=γ±γ2ω02|0<γ<ω0λ1,2=γ±γ2ω02=γ±γ2ω02=γ±ω2ω2=|γ2ω02|,γ2ω02=ω2=γ±ω1=γ±ωi C1,2=x02±v0+γx02γ2ω02|0<γ<ω0C1,2=x02±v0+γx02γ2ω02=x02±v0+γx02ω1=x02±v0+γx02ωiii=x02±v0+γx02ω(1)i=x02v0+γx02ωiω0 x(t)=C1eλ1t+C2eλ2t=C1e(γ+ωi)t+C2e(γωi)t=C1eγteωit+C2eγteωit=eγt{C1eωit+C2eωit}{C1eωit+C2eωit}={x02v0+γx02ωi}eωit+{x02+v0+γx02ωi}eωit={x02v0+γx02ωi}{cos(ωt)+isin(ωt)}+{x02+v0+γx02ωi}{cos(ωt)+isin(ωt)}eix=cos(x)+isin(x)={x02v0+γx02ωi}{cos(ωt)+isin(ωt)}+{x02+v0+γx02ωi}{cos(ωt)isin(ωt)}cos(θ)=cos(θ),sin(θ)=sin(θ)={x02v0+γx02ωi}cos(ωt)+i{x02v0+γx02ωi}sin(ωt)+{x02+v0+γx02ωi}cos(ωt)i{x02+v0+γx02ωi}sin(ωt)={x02v0+γx02ωi}cos(ωt)+{x02i+v0+γx02ω}sin(ωt)+{x02+v0+γx02ωi}cos(ωt){x02iv0+γx02ω}sin(ωt)ii=1={x02v0+γx02ωi+x02+v0+γx02ωi}cos(ωt)+{x02i+v0+γx02ωx02i+v0+γx02ω}sin(ωt)={x02+x02}cos(ωt)+{v0+γx02ω+v0+γx02ω}sin(ωt)=x0cos(ωt)+v0+γx0ωsin(ωt) x(t)=eγt{C1eωit+C2eωit}=eγt{x0cos(ωt)+v0+γx0ωsin(ωt)}=x0eγt{cos(ωt)+γωsin(ωt)}+v0eγt{1ωsin(ωt)},=x0eγt{cos(ωt)+ζζ21sin(ωt)}+v0eγt{1ωsin(ωt)}γω=γ|γ2ω02|=γω02ω02|γ2ω02|=γω02|γ2ω02ω02ω02|=γω0|ζ21|=ζ|ζ21|ω=|γ2ω02|,ζ=γω0

バネマスダンパー系,単振動(γ=0 / バネマス系)の場合

バネマスダンパー系

γ=0

λ1,2=γ±γ2ω02|γ=0λ1,2=γ±γ2ω02=0±02ω02=±ω02=±ω01=±ω0i C1,2=x02±v0+γx02γ2ω02|γ=0C1,2=x02±v0+γx02γ2ω02=x02±v0+0x0202ω02=x02±v02ω01=x02±v02ω0iii=x02±v02ω0(1)i=x02v02ω0iω00 x(t)=C1eλ1t+C2eλ2t={x02v02ω0i}eω0it+{x02+v02ω0i}eω0it={x02v02ω0i}{cos(ω0t)+isin(ω0t)}+{x02+v02ω0i}{cos(ω0t)+isin(ω0t)}eix=cos(x)+isin(x)={x02v02ω0i}{cos(ω0t)+isin(ω0t)}+{x02+v02ω0i}{cos(ω0t)isin(ω0t)}cos(θ)=cos(θ),sin(θ)=sin(θ)={x02v02ω0i}cos(ω0t)+i{x02v02ω0i}sin(ω0t)+{x02+v02ω0i}cos(ω0t)i{x02+v02ω0i}sin(ω0t)={x02v02ω0i}cos(ω0t)+{x02i+v02ω0}sin(ω0t)+{x02+v02ω0i}cos(ω0t){x02iv02ω0}sin(ω0t)ii=1={x02v02ω0i+x02+v02ω0i}cos(ω0t)+{x02i+v02ω0x02i+v02ω0}sin(ω0t)={x02+x02}cos(ωt)+{v02ω0+v02ω0}sin(ωt)=x0cos(ω0t)+v0ω0sin(ω0t)=x0cos(ω0t)+v01ω0sin(ω0t),

バネマスダンパー系,C1, C2 を求める

バネマスダンパー系

C1,C2

sx0+v0+cmx0s2+scm+km=(C1+C2)s(C1λ2+C2λ1)(sλ1)(sλ2) sx0+v0+cmx0=(C1+C2)s(C1λ2+C2λ1) {C1+C2=x0(C1λ2+C2λ1)=v0+cmx0=v0+2γx0s C1=x0C2(C1λ2+C2λ1)=v0+2γx0C1λ2+C2λ1=(v0+2γx0)(x0C2)λ2+C2λ1=v02γx0x0λ2C2λ2+C2λ1=v02γx0C2(λ1λ2)=v02γx0λ2x0C2(λ1λ2)=v0x0(2γ+λ2)C2=v0x0(2γ+λ2)λ1λ2=v0+x0{2γ+(γγ2ω02)}(γ+γ2ω02)(γγ2ω02)=v0+x0(γγ2ω02)2γ2ω02=v0+γx0x0γ2ω022γ2ω02=v0+γx02γ2ω02+x0γ2ω022γ2ω02=v0+γx02γ2ω02+x02=x02v0+γx02γ2ω02C1=x0C2=x0{x02v0+γx02γ2ω02}=x02+v0+γx02γ2ω02 C1,2=x02±v0+γx02γ2ω02γω0

バネマスダンパー系,λ1,λ2を求める

バネマスダンパー系

λ1,2

sx0+v0+cmx0s2+scm+km=(C1+C2)s(C1λ2+C2λ1)(sλ1)(sλ2) s2+scm+km=(sλ1)(sλ2)λ1,2=cm±(cm)241km21=c2m±12(cm)24km=c2m±14((cm)24km)=c2m±(c2m)2km=γ±γ2ω02γ=c2m,ω02=km

バネマスダンパー系,運動方程式,ラプラス変換,逆ラプラス変換,入力なし

バネマスダンパー系

運動方程式

x¨+cmx˙+kmx=0d2xd2t+cmdxdt+kmx=0

ラプラス変換と計算

L[d2xd2t+cmdxdt+kmx]=0L[d2xd2t]+L[cmdxdt]+L[kmx]=0L[ag(t)+bh(t)]=aL[g(t)]+bL[h(t)]s2Xsx0v0+cm(sXx0)+kmX=0L[x]=XL[dxdt]=s2Xx0,x0=x(0)L[d2xd2t]=s2Xsx0v0,v0=x(0)s2Xsx0v0+scmXcmx0+kmX=0 s2X+scmX+kmX=sx0+v0+cmx0(s2+scm+km)X=sx0+v0+cmx0X=sx0+v0+cmx0s2+scm+km=sx0+v0+cmx0(sλ1)(sλ2)λ1,2=s2+scm+km=C1sλ1+C2sλ2=C1(sλ2)+C2(sλ1)(sλ1)(sλ2)=(C1+C2)s(C1λ2+C2λ1)(sλ1)(sλ2)

逆ラプラス変換

L1[X]=L1[C11sλ1+C21sλ2]=C1L1[1sλ1]+C2L1[1sλ2]x(t)=C1eλ1t+C2eλ2tL1[X]=x(t)L1[1s+a]=eat

eを底とした指数凾数(at乗)のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

eatのラプラス変換

f(t)=eatL[f(t)]=0f(t)estdt=0eatestdt=0e(sa)tdtu=(sa)t,t:0,u:0dudt=(sa),dt=1sadu=0eu1sadu=1sa0eudu=1sa[eu]0=1sa[ee0]=1sa[01]=1sa(1)=1sa

sin凾数のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

sin(ωt)のラプラス変換

f(t)=sin(ωt)L[f(t)]=0sin(ωt)estdt=[sin(ωt)1sest]00ωcos(ωt)1sestdtabf(t)g(t)dt=[f(t)g(t)]ababf(t)g(t)dt(f:f,g:g)eatdt=1aeat+C,ddtsin(ωt)=ωcos(ωt)=[sin(ω)1sessin(ω0)1ses0](ωs)0cos(ωt)estdt=[00]+ωs0cos(ωt)estdte=0,sin(0)=0=ωs0cos(ωt)estdt=ωs[[cos(ωt)1sest]00ωsin(ωt)1sestdt]abf(t)g(t)dt=[f(t)g(t)]ababf(t)g(t)dt(f:f,g:g)eatdt=1aeat+C,ddtcos(ωt)=ωsin(ωt)=ωs[[cos(ω)1sescos(ω0)1ses0]ωs0sin(ωt)estdt]=ωs[[011s1]ωs0sin(ωt)estdt]e=0,cos(0)=1,e0=1=ωs[1sωs0sin(ωt)estdt]=ωs2ω2s20sin(ωt)estdt=ωs2ω2s2L[f(t)]L[f(t)]+ω2s2L[f(t)]=ωs2L[f(t)](1+ω2s2)=ωs2L[f(t)]=ωs211+ω2s2=ωs2(1+ω2s2)=ωs2+ω2

冪凾数のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

tnのラプラス変換

f(t)=tnL[f(t)]=0tnestdt=[tn1sest]00ntn11sestdtabf(t)g(t)dt=[f(t)g(t)]ababf(t)g(t)dt(f:f,g:g)eatdt=1aeat+Cddttn=ntn1=[n1ses0n1ses0](ns)0tn1estdt=[00]+ns0tn1estdte=0,sin(0)=0=ns0tn1estdt=nsn1s0tn2estdt=nsn1s2s0testdt=nsn1s2s1s0estdt=nsn1s2s1s[1sest]0=n!sn1s[eses0]n(n1)21=n!=n!sn1s[01]e=0,e0=1=n!sn1s(1)=n!sn+1

積分のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

0tg(u)duのラプラス変換

f(t)=0tg(u)duL[f(t)]=0{0tg(u)du}estdt=[{0tg(u)du}1sest]00g(t)1sestdtabf(t)g(t)dt=[f(t)g(t)]ababf(t)g(t)dt(f:f,g:g)eatdt=1aeat+Cddt0tg(u)dt=ddt[G(t)G(0)]=ddtG(t)ddtG(0)=g(t)0=g(t)=[{0g(u)du}1ses{00g(u)du}1ses0]+1s0g(t)estdte=0,ttf(x)dx=[F(t)F(t)]=0=[00]+1s0g(t)estdt=1sL[g(t)]

微分のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

dfdtのラプラス変換

f(t)=dg(t)dtL[g(t)]=0dg(t)dtestdt=[g(t)est]00dg(t)dt(s)estdtabf(t)g(t)dt=[f(t)g(t)]ababf(t)g(t)dt(f:f,g:g)ddteat=aeat=[g()esg(0)es0](s)0dg(t)dtestdt=[0g(0)1](s)0dg(t)dtestdte=0,e0=1=g(0)+s0dg(t)dtestdt=g(0)+sL[g(t)]=sL[g(t)]g(0)

線形結合のラプラス変換

ラプラス変換

L[f(t)]=0f(t)estdt

ag(t)+bh(t)のラプラス変換

f(t)=ag(t)+bh(t)a,b:t()L[f(t)]=0{ag(t)+bh(t)}estdt=0{ag(t)est+bh(t)est}dt=0ag(t)estdt+0bh(t)estdt=a0g(t)estdt+b0h(t)estdt=aL[g(t)]+bL[h(t)]

定積分の微分

定積分の微分

ddt0tg(x)dx=ddt[G(x)]0t=ddt[G(t)G(0)]=ddtG(t)ddtG(0)=g(t)00G(0)0=g(t)