clear;
clc;
subplot(2,2,1);
%% 风力大小变化
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=5;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'k',v,Bt_2,'k-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=10;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'r',v,Bt_2,'r-.')
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=15;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'b',v,Bt_2,'b-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=20;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'y',v,Bt_2,'y-.')
legend('Vw=5m/s,B0','Vw=5m/s,Bt','Vw=10m/s,B0','Vw=10m/s,Bt','Vw=15m/s,B0','Vw=15m/s,Bt','Vw=20m/s,B0','Vw=20m/s,Bt')
%% 风向大小变化
subplot(2,2,2);
Qw=pi/6;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=10;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'k',v,Bt_2,'k-.');
hold on;
clear;
Qw=pi/4;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=10;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'r',v,Bt_2,'r-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=10;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'b',v,Bt_2,'b-.');
hold on;
legend('Qw=pi/6,B0','Qw=pi/6,Bt','Qw=pi/4,B0','Qw=pi/4,Bt','Qw=pi/2,B0','Qw=pi/2,Bt');
subplot(2,2,3);
%% 流力大小变化
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=5;%%风速风向
Vc=1;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'k',v,Bt_2,'k-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=5;%%风速风向
Vc=2;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'r',v,Bt_2,'r-.')
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=5;%%风速风向
Vc=3;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'b',v,Bt_2,'b-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=5;%%风速风向
Vc=4;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'y',v,Bt_2,'y-.')
legend('Vc=1kn,B0','Vc=1kn,Bt','Vc=2kn,B0','Vc=2kn,Bt','Vc=3kn,B0','Vc=3kn,Bt','Vc=4kn,B0','Vc=4kn,Bt');
%% 流向大小变化
subplot(2,2,4);
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=5;%%风速风向
Vc=4;Qc=pi/6;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'k',v,Bt_2,'k-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=10;%%风速风向
Vc=4;Qc=pi/4;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v;Bt_2-Bt_1];
plot(v,Bt_1,'r',v,Bt_2,'r-.');
hold on;
clear;
Qw=pi/2;
v=2:0.1:20;
L=274;
B=32;
d=13;
Cb=0.8142;
K=2.*v/L;%%旋回性指数
AR=70;%%舵面积
Br=420;
Au=L*d/2;%%受风面积
Vw=10;%%风速风向
Vc=4;Qc=pi/2;%%流向流速
Ad=L*d;%%受流面积
thea=K.*(Vw./v).^2*sin(Qw)+asin(Vc*sin(Qc)./v);%%偏航角：一般指船舶轴线方向与航道方向之间的夹角
Bt_1=(0.041*sqrt(Au/Ad).*(exp(-0.14.*v))*Vw+Vc.*sin(thea))*(26.464+0.408*Cb*L*d/AR-0.033*L/B*L*d/AR-79.114*Cb+0.757*L/B+46.129*Cb^2)*L./v;
Bt_2=Br/2-B/2-(B+L.*sin(thea))/4;
%%Bt=Bt_2-Bt_1;
%%[v
