原理介紹

船體數據來源來自于C++代碼。在C++中，將船體數據導入到EXCEL中，然后再MATLAB中調用EXCEL的船體數據，進行三維建模。在建模過程中，因為船體數據不可能每個點的距離都是無限小，所以我們采用了擬合的辦法，調用了MATLB的函數? griddata 。對船體數據進行三角形插值擬合。最終得到三維船體圖形。如下：

遺傳算法代碼如下

#include <iostream> #include<math.h> #include<time.h> #include<stdlib.h> #include<fstream> #define tot 122*18*16 using namespace std; int main(int argc, const char * argv[]) {double suf[90100];double vol[90100];double fin[122][60];double fin1[122][60];double ran1[122][60];double head[60],head2[60];int min=99999,flag=0;ofstream out("grade");if(!out){cout << "Cannot open grade file.\n";return 1;} for(int t=0;t<90100;t++){suf[t]=0;vol[t]=0;}int i=900; while(i>=0){double ran[122][60];double area[60];for(int t=0;t<60;t++){head[t]=3*0.1*0.1*t*t-32*t*0.1;head2[t]=0.5*t*t*0.01-8*t*0.1;}for(int t=0;t<600;t++){area[t]=0;}for(int z=0;z<60;z++){for(int x=1;x<122;x++){if(x<100){srand((unsigned)time(NULL)+x*z*i*7);ran[x][z]=max((11.22925*pow(x, 0.1)+0.01*(abs(rand()%100)))*(((-z*0.1+12.00)*(double)(z*0.1))/36.00),ran[x-1][z]+(((-z*0.1+12.00)*(double)(z*0.1))/36.00)*0.01*(abs(rand()%100)));ran1[x][z]=max((14.2171*pow(x, 0.05)+0.01*(abs(rand()%100)))*(((-z*0.1+12.00)*(double)(z*0.1))/36.00),ran[x-1][z]+(((-z*0.1+12.00)*(double)(z*0.1))/36.00)*0.01*(abs(rand()%100)));}else{ran[x][z]=ran[x-1][z];ran1[x][z]=ran1[x-1][z];}// cout<<ran[x][z]<<" ";}} for(int z=0;z<60;z++){for(int x=0;x<122;x++){area[z]+=4*(ran[x][z]+ran[x+1][z]);}//cout<<area[z]<<" ";}for(int z=1;z<60;z++){for(int x=1;x<122;x++){double a[5];double p1=0,p2=0;a[0]=sqrt(abs((ran[x][z]-ran[x-1][z])*(ran[x][z]-ran[x-1][z])+((122+head[z])/122)*((122+head[z])/122)));a[1]=sqrt(abs((ran[x][z]-ran[x][z-1])*(ran[x][z]-ran[x][z-1])+1+(((122-x)*(head[z]-head[z-1]))/122)*(((122-x)*(head[z]-head[z-1]))/122)));a[2]=sqrt(abs((ran[x-1][z]-ran[x][z-1])*(ran[x-1][z]-ran[x][z-1])+1+((((122-x+1)*(head[z-1]-head[z]))/122)+(122+head[z-1])/122)*((((122-x+1)*(head[z-1]-head[z]))/122)+(122+head[z-1])/122)));a[3]=sqrt(abs((ran[x][z-1]-ran[x-1][z-1])*(ran[x][z-1]-ran[x-1][z-1])+((122+head[z-1])/122)*((122+head[z-1])/122)));a[4]=sqrt(abs((ran[x-1][z]-ran[x-1][z-1])*(ran[x-1][z]-ran[x-1][z-1])+(((122-x+1)*(head[z]-head[z-1]))/122)*(((122-x+1)*(head[z]-head[z-1]))/122)));p1=(a[0]+a[1]+a[2])/2;p2=(a[2]+a[3]+a[4])/2;suf[i]+=sqrt(abs(p1*(p1-a[0])*(p1-a[1])*(p1-a[2])));suf[i]+=sqrt(abs(p2*(p2-a[2])*(p2-a[3])*(p2-a[4])));}}vol[i]=(area[0]+area[5])*0.5+area[1]+area[2]+area[3]+area[4];if(suf[i]<min){min=suf[i];flag=i;for(int t=0;t<122;t++){for(int q=0;q<60;q++){fin[t][q]=ran[t][q];fin1[t][q]=ran1[t][q];}}}//cout<<suf[i]<<" ";i--;}cout<<flag<<" "<<min<<" "<<vol[flag]<<" "<<suf[flag]<<endl;for(int t=0;t<60;t++){for(int i=0;i<122;i++){out<<0.5*((double)i+((double)(122-i)*head[t]/(double)122))<<" "<<fin[i][t]<<" "<<(double)t<<endl;out<<0.5*((double)i+((double)(122-i)*head[t]/(double)122))<<" "<<-fin[i][t]<<" "<<(double)t<<endl;out<<0.5*((double)244-(double)i-((double)(122-i)*head2[t]/(double)122))<<" "<<fin1[i][t]<<" "<<(double)t<<endl;out<<0.5*((double)244-(double)i-((double)(122-i)*head2[t]/(double)122))<<" "<<-fin1[i][t]<<" "<<(double)t<<endl;}cout<<endl;}out.close(); }

matlab仿真代碼如下：

function?createfigure(xdata1,?ydata1,?zdata1,?cdata1,?X1,?Y1,?Z1)%CREATEFIGURE(xdata1,?ydata1,?zdata1,?cdata1,?X1,?Y1,?Z1)%??XDATA1:??surface?xdata%??YDATA1:??surface?ydata%??ZDATA1:??surface?zdata%??CDATA1:??surface?cdata%??X1:??x?數據的向量%??Y1:??y?數據的向量%??Z1:??z?數據的向量%??由?MATLAB?于?31-Oct-2020?15:47:44?自動生成%?創建?figurefigure1?=?figure;%?創建?axesaxes1?=?axes(&apos;Parent&apos;,figure1);hold(axes1,&apos;on&apos;);%?創建?meshmesh(xdata1,ydata1,zdata1,cdata1,&apos;Parent&apos;,axes1);%?創建?plot3plot3(X1,Y1,Z1,&apos;MarkerSize&apos;,15);%?創建?titletitle(&apos;船體模型&apos;);view(axes1,[-37.5?30]);grid(axes1,&apos;on&apos;);axis(axes1,&apos;tight&apos;);hold(axes1,&apos;off&apos;);%?設置其余坐標區屬性set(axes1,&apos;DataAspectRatio&apos;,[1?1?1]); clc;x=Column1(:,1);%導入數據y=Column2(:,1);%導入數據z=Column3(:,1);%導入數據xlin?=?linspace(min(x),max(x),1000);%對X插值ylin?=?linspace(min(y),max(y),1000);%對Y插值[X,Y]?=?meshgrid(xlin,ylin);%對X，Y進行曲面擬合Z?=?griddata(x,y,z,X,Y,&apos;linear&apos;);%對Z進行擬合mesh(X,Y,Z)?%interpolatedaxis?equal;axis?tight;?hold?on%設置坐標軸plot3(x,y,z,&apos;MarkerSize&apos;,15)?%nonuniform%進行繪制三維圖形title(&apos;船體模型&apos;);%?plot3(x,y,z)

?

總結

以上是生活随笔為你收集整理的基于matlab的船体三维模型仿真的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。