function [sys,x0,str,ts] = s_function(t,x,u,flag)
switch flag,
case 0,
[sys,x0,str,ts]=mdlInitializeSizes;
case 1,
sys=mdlDerivatives(t,x,u);
case 3,
sys=mdlOutputs(t,x,u);
case {2,4,9}
sys=[];
otherwise
error(['Unhandled flag = ',num2str(flag)]);
end
function [sys,x0,str,ts]=mdlInitializeSizes
sizes = simsizes;
sizes.NumContStates = 3;
sizes.NumDiscStates = 0;
sizes.NumOutputs = 3;
sizes.NumInputs = 3;
sizes.DirFeedthrough = 0;
sizes.NumSampleTimes = 1;
sys = simsizes(sizes);
x0 = [0.01;0;0];
str = [];
ts = [0 0];
function sys=mdlDerivatives(t,x,u)
L=1.46; m=45;density=1025.0;
L2=0.5*density*L*L;
L3=L2*L;
L4=L3*L;
L5=L4*L;
Xdu=-1.5777e-3*L3;Xuu=-5.9e-3*L2;Xu=0;
Ydv=-3.0753e-2*L3;Yv=-4.4961e-2*L2;Yvv=-1.6687e-1*L2;
Iz=8.067;Ndr=-1.6012e-3*L5;Nrr=-1.4432e-3*L5;Nr=-1.1679e-2*L4;
M=[m-Xdu 0 0 ;0 m-Ydv 0 ;0 0 Iz-Ndr];
C=[0 0 -(m-Ydv)*x(2);0 0 -(m-Xdu)*x(1);(m-Ydv)*x(2) -(m-Xdu)*x(1) 0];
D=-[Xu+Xuu*norm(x(1)) 0 0 ;0 Yv+Yvv*norm(x(2)) 0 ;0 0 Nr+Nrr*norm(x(3))];
g=[0;0;0];
V=[x(1);x(2);x(3)];
T=[u(1);u(2);u(3)];
f=[0;0;0];
sys=inv(M)*(T-g-C*V-D*V+f);
function sys=mdlOutputs(t,x,u)
sys(1)=x(1);
sys(2)=x(2);
sys(3)=x(3);
没有合适的资源?快使用搜索试试~ 我知道了~
船舶路径跟纵,只限跟踪,contro微控制器,vlm为船模型Matlab代码.rar
共14个文件
m:7个
emf:4个
asv:2个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
0 下载量 24 浏览量
2024-10-10
11:04:13
上传
评论
收藏 29KB RAR 举报
温馨提示
无干扰情况下 船舶的回转性试验验证 二阶NOMOTO
资源推荐
资源详情
资源评论
收起资源包目录
船舶路径跟纵,只限跟踪,contro微控制器,vlm为船模型Matlab代码.rar (14个子文件)
船舶路径跟纵,只限跟踪,contro微控制器,vlm为船模型Matlab代码
路径跟踪(直线)
n.asv 776B
adapt.m 773B
renyier.emf 38KB
X.m 804B
WLM.m 1KB
tao.m 835B
renyiyi.emf 27KB
n.m 785B
lujinggenzong.mdl 32KB
yuaner.emf 31KB
yuanyi.emf 25KB
Y.m 809B
control.m 1KB
tao.asv 828B
共 14 条
- 1
资源评论
天天Matlab代码科研顾问
- 粉丝: 3w+
- 资源: 1873
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
最新资源
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功