- 相關推薦
實驗計算機輔助船舶設計方案
實驗計算機輔助船舶設計方案
小撓度曲線算例
編程代碼:
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;
%求解系數(shù)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
c
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;
【實驗計算機輔助船舶設計方案】相關文章:
制備實驗的設計方案03-10
計算機輔助與設計學生求職信12-08
《夜晚的實驗》教學設計方案11-26
實驗設計方案8篇12-04
簡單生物實驗設計方案12-01
教育實驗設計方案參考12-02
實驗設計方案九篇03-10
實驗設計方案3篇12-09
實驗設計方案(精選25篇)12-04