实验计算机辅助船舶设计方案

时间：2024-08-16 04:44:04 设计方案我要投稿

相关推荐

实验计算机辅助船舶设计方案

　　小挠度曲线算例

实验计算机辅助船舶设计方案

　　编程代码：

　　clear

　　x=[0 400 800 1200 1600 2000 2400 2800 3060];

　　y=[0 130 232.6 320.4 348 389 329.4 165 0];

　　n=length(x);

　　dy(1)=0.361;

　　dy(n)=-0.673;

　　%求解系数c（i）

　　for i=1:n-1

　　m(i)=(y(i+1)-y(i))/(x(i+1)-x(i));

　　end

　　a(1)=0.5;

　　b(1)=(m(1)-dy(1))/(2*(x(2)-x(1)));

　　for i=2:n-1;

　　fff=2*(x(i+1)-x(i-1))-a(i-1)*(x(i)-x(i-1));

　　a(i)=(x(i+1)-x(i))/fff;

　　b(i)=(m(i)-m(i-1)-(x(i)-x(i-1))*b(i-1))/fff;

　　end

　　for i=(n-1):-1:1

　　a(n)=0;

　　c(n)=(dy(n)-m(n-1)-(x(n)-x(n-1))*b(n-1))/(2-a(n-1))/(x(n)-x(n-1));

　　b(n)=c(n);

　　c(i)=b(i)-a(i)*c(i+1);

　　end

　　x1=0:1:3060;

　　y1(1)=y(1);

　　y1(3061)=y(n);

　　for j=1:1:3059

　　x1(j+1)=x1(j)+1;

　　end

　　for i=1:1:n-1;

　　i=1;

　　for j=2:1:3060

　　if (x1(j)-x(i))*(x1(j)-x(i+1))>0;

　　i=i+1;

　　end

　　t=c(i)*(2*x(i+1)-x1(j)-x(i))+c(i+1)*(x(i+1)+x1(j)-2*x(i));

　　y1(j)=y(i)+m(i)*(x1(j)-x(i))-((x(i+1)-x1(j))*(x1(j)-x(i)))/(x(i+1)-x(

　　i))*t;

【实验计算机辅助船舶设计方案】相关文章：