实验计算机辅助船舶设计new

《实验计算机辅助船舶设计new》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、小挠度曲线算例编程代码：clearx=[0400800120016002000240028003060];y=[0130232.6320.4348389329.41650];n=length(x);dy(1)=0.361;dy(n)=-0.673;%求解系数c（i）fori=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)));fori=2:n-1;fff=2*(x(i+1)-x(i-1))-a(i

2、-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;endfori=(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);forj=1:1:

3、3059x1(j+1)=x1(j)+1;endfori=1:1:n-1;i=1;forj=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;endendplot(x1,y1,'*')a,b运行结果：c=

4、1.0e-003*Columns1through7-0.0245-0.04100.0174-0.12110.0906-0.1576-0.0889Columns8through9-0.1419-0.0029a=Columns1through70.50000.28570.26920.26800.26800.26790.2679Columns8through90.21440b=1.0e-003*Columns1through7-0.0450-0.0361-0.0152-0.09680.0484-0.1814-0.12

5、69Columns8through9-0.1425-0.0029对运行结果的分析：小挠度插值曲线：clearx=[8.1258.499.4859.69.95910.16610.2];y=[0.07740.0990.280.60.7081.21.82.177];n=length(x);k(1)=0.01087;k(n)=100;X=x(1);Y=y(1);z=atan((y(n)-y(1))/(x(n)-x(1)));%进行坐标转换fori=1:naaa(i)=x(i)-X;bbb(i)=y(i)-Y;x1=aa

6、a*cos(z)+bbb*sin(z);y1=bbb*cos(z)-aaa*sin(z);enddy(1)=tan(atan(k(1))-z);dy(n)=tan(atan(k(n))-z);x=x1;y=y1;plot(x,y,x,y,'co')pauseholdoff%求解系数c（i）fori=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)));fori=2:n-1;fff=2*(x(i+1)

7、-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;endfori=(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);endcab=0.01;j=1;xx(j)=x(1);whilexx(j)<

8、=x(n)j=j+1;xx(j)=xx(j-1)+ab;endk=j;xx(k)=x(n);fx(1)=y(1);fx(k)=y(n);fori=1:n-1i=1;forj=2:k-1if(xx(j)-x(i))*(xx(j)-x(i+1))>0;i=i+1;endfff1=m(i)*(xx(j)-x(i));fff2=(x(i+1)-xx(j))*(xx(j)-x(i));fff3=