生活随笔
收集整理的這篇文章主要介紹了
Matlab：船舶航向舵数学模型构建
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個參考.
%matlab船舶數(shù)學(xué)模型
%可參考學(xué)習(xí)，不要抄襲。望學(xué)業(yè)進(jìn)步
%船舶數(shù)學(xué)模型構(gòu)建
clc
,clear
,close all
;
%% 船舶設(shè)計參數(shù)
L
=100;B
=20.00;m
=15000;T
=10;
V
=15;Xc
=0.30;Cb
=0.700;A
=20.0;%% 參數(shù)無量綱化
V
=(1852/3600)*V
;%速度國際單位
Rho
=1025;%海水密度
m
=m
*Rho
;
m1
=m
/((1/2)*Rho
*L
^3);
Xc1
=Xc
/L
;
Izz
=m
*(L
^2)/16;
Izz1
=Izz
/((1/2)*Rho
*L
^5);%% 流體動力導(dǎo)數(shù)公式
dYdv
=-(1+0.16*Cb
*B
/T
-5.1*((B
/L
)^2))*pi
*(T
/L
)^2;
dYdr
=-(0.67*B
/L
-0.0033*(B
/T
)^2)*pi
*(T
/L
)^2;
dNdv
=-(1.1*B
/L
-0.041*B
/T
)*pi
*(T
/L
)^2;
dNdr
=-(1/12+0.017*Cb
*B
/T
-0.33*B
/L
)*pi
*(T
/L
)^2;
dYv
=-(1+0.40*Cb
*B
/T
)*pi
*(T
/L
)^2;
dYr
=-(-1/2+2.2*B
/L
-0.080*B
/T
)*pi
*(T
/L
)^2;
dNv
=-(1/2+2.4*T
/L
)*pi
*(T
/L
)^2;
dNr
=-(1/4+0.039*B
/T
-0.56*B
/L
)*pi
*(T
/L
)^2;
Ydt
=3.0*A
/(L
^2);
Ndt
=-(1/2)*Ydt
;%% 流體動力導(dǎo)數(shù)修正
ddYv
=-0.3*Ydt
;
dYv
=dYv
+ddYv
;
dYr
=dYr
+((-1/2)*ddYv
);
dNv
=dNv
+((-1/2)*ddYv
);
dNr
=dNr
+((1/4)*ddYv
);%% 構(gòu)造二自由度狀態(tài)空間型船舶線性數(shù)學(xué)模型
I1
=[m1
-dYdv
,L
*(m1
*Xc1
-dYdr
);m1
*Xc1
-dNdv
,L
*(Izz1
-dNdr
)];
P1
=[V
*dYv
/L
,V
*(dYr
-m1
);V
*dNv
/L
,V
*(dNr
-m1
*Xc1
)];
Q1
=[-(V
^2)*Ydt
/L
;-(V
^2)*Ndt
/L
];
a
=inv
(I1
)*P1
b
=inv
(I1
)*Q1
K
=(b
(1,1)*a
(2,1)-b
(2,1)*a
(1,1))/(a
(1,1)*a
(2,2)-a
(1,2)*a
(2,1))
T
=-((a
(1,1)+a
(2,2))/(a
(1,1)*a
(2,2)-a
(1,2)*a
(2,1))+b
(2,1)/(b
(1,1)*a
(2,1)-b
(2,1)*a
(1,1)))
總結(jié)
以上是生活随笔為你收集整理的Matlab：船舶航向舵数学模型构建的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。
