【船舶运动】船舶运动仿真【含Matlab源码 3494期】

最新推荐文章于 2024-08-10 22:15:42 发布

Matlab领域

最新推荐文章于 2024-08-10 22:15:42 发布

阅读量817

点赞数 17

分类专栏： Matlab物理应用（初级版）文章标签： matlab

本文链接：https://blog.csdn.net/weixin_63266434/article/details/135549699

版权

Matlab物理应用（初级版）专栏收录该内容

278 篇文章 10 订阅

订阅专栏

在这里插入图片描述

⛄一、获取代码方式

获取代码方式1：
完整代码已上传我的资源：【船舶运动】基于matlab船舶运动仿真【含Matlab源码 3494期】
点击上面蓝色字体，直接付费下载，即可。

获取代码方式2：
付费专栏Matlab物理应用（初级版）

备注：
点击上面蓝色字体付费专栏Matlab物理应用（初级版），扫描上面二维码，付费29.9元订阅海神之光博客付费专栏Matlab物理应用（初级版），凭支付凭证，私信博主，可免费获得1份本博客上传CSDN资源代码（有效期为订阅日起，三天内有效）；
点击CSDN资源下载链接：1份本博客上传CSDN资源代码

⛄二、部分源代码

function F=luoxuanjiang(X)
global L B Cp Cb rou V n Dp Cpp wp;

%% 需要输入的变量量：u，v,r,δ
u=X(1);
v=X(2);
r=X(3)pi/180;
rn=rL/V;%无因次r
xpn=-0.4;%无因次xp
beta=-atan(v/u);%漂角β
delta=X(4)*pi/180;%转化为弧度

%% --推力系数的计算–
kp(1)=(6Cb+0.9L/B-5.7)10^(-3);
kp(2)=(14Cb-0.5*L/B-0.3)*10^(-3);

wp0=0.55Cb-0.2;
for i=1:2
wp(i)=wp0-kp(i)(beta-xpnrn);
ua(i)=(1-wp(i))u;
JP(i)=ua(i)/(nDp);
kT(i)=-0.4315JP(i)+0.6179;
T(i)=roun^2*Dp4kT(i);
end

%% --推力减额系数的计算–
tp0=0.5Cp-0.18;
BR=beta+0.95L*(r/V);
xB=-0.068;
lcb=xB100/L;
A=(B/L)/(1.3(1-Cb)-3.1lcb);
Kt=0.00023(AL/Dp)-0.028;
tp=tp0-KtBR;
%% --桨力和力矩计算–
if delta>=0
a=1;
else
a=-1;
end
XP=(1-tp)(T(1)+T(2));
YP=0;
NP=a(1-tp)Cpp(T(1)-T(2));
F=[XP YP NP];
%XPn=XP/(0.5rouV^2*L2);
%YPn=YP/(0.5rouV^2*L2);
%NPn=NP/(0.5rouV^2*L3);
%F=[XPn YPn NPn];