《实验计算机辅助船舶设计方案.docx》由会员分享,可在线阅读,更多相关《实验计算机辅助船舶设计方案.docx(3页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、实验计算机辅助船舶设计方案实验计算机辅助船舶设计方案 编程代码: clearx=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-1m(i)=(y(i+1)-y(i)/(x(i+1)-x(i);enda(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
2、-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;endfor i=(n-1):-1:1a(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);endcx1=0:1:3060;y1(1)=y(1);y1(3061)=y(n);for j=1:1:3059x1(j+1)=x1(j)+1;endfor i=1:1:n-1;i=1;for j=2:1:3060if (x1(j)-x(i)*(x1(j)-x(i+1)0;i=i+1;endt=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; 3 / 3