matlab 边界拟合,边界拟合坐标系实例(matlab编程).doc

边界拟合坐标系实例(matlab编程).doc

物理平面坐标为圆，其上下左右为边界拟合坐标的边界，将各边界分成30份，则对于拟合坐标有fori=1:Nx(i,1)=R*cos(-3*pi/4+(i-1)*pi/(2*(N-1)));y(i,1)=R*sin(-3*pi/4+(i-1)*pi/(2*(N-1)));x(i,N)=R*cos(3*pi/4-(i-1)*pi/(2*(N-1)));y(i,N)=R*sin(3*pi/4-(i-1)*pi/(2*(N-1)));endforj=1:Nx(1,j)=R*cos(5*pi/4-(j-1)*pi/(2*(N-1)));y(1,j)=R*sin(5*pi/4-(j-1)*pi/(2*(N-1)));x(N,j)=R*cos(-pi/4+(j-1)*pi/(2*(N-1)));y(N,j)=R*sin(-pi/4+(j-1)*pi/(2*(N-1)));end由边界拟合坐标转换方程220 xxJPQxyyy其中2xy2xyJ对于本问题可将物理平面均匀划分，故0PQ故2xxyy设，并采用中心差分格式，有限差分方程为=11,1,,1,1,,,,2220ijijijijijijijijijijijijxxxxyyyy上述有限差分方程可采用显示的迭代格式求解1,,,1,11,,,1,122nijijijijnijnijijijijnijxxxyyy上式中表示迭代次数,其中2xy2xy而其中、、、、、将采用中心差分近似xyxy1,,2ijijx,1,ijij1,1,1,1,4ijijijijxx1,,2ijijy,1,ijij1,1,1,1,4ijijijijyy给出matlab程序如下N=30;%%准备将物理平面上下左右都分为30份x(N,N)=0;y(N,N)=0;%%设初值，将拟合平面上的每个坐标对应x、y值都设为0e0=1e-4;%%设定的求解精度R=1;%%设定的圆的半径%%下面设定拟合坐标上的边值fori=1:Nx(i,1)=R*cos(-3*pi/4+(i-1)*pi/(2*(N-1)));y(i,1)=R*sin(-3*pi/4+(i-1)*pi/(2*(N-1)));x(i,N)=R*cos(3*pi/4-(i-1)*pi/(2*(N-1)));y(i,N)=R*sin(3*pi/4-(i-1)*pi/(2*(N-1)));endforj=1:Nx(1,j)=R*cos(5*pi/4-(j-1)*pi/(2*(N-1)));y(1,j)=R*sin(5*pi/4-(j-1)*pi/(2*(N-1)));x(N,j)=R*cos(-pi/4+(j-1)*pi/(2*(N-1)));y(N,j)=R*sin(-pi/4+(j-1)*pi/(2*(N-1)));endw=1.2;%%设定为1.2的超松弛因子以加速收敛fork=1:10000%%设定迭代步数，希望在10000步内收敛err=0;%%下面根据中心差分法进行坐标拟合fori=2:N-1forj=2:N-1x1=(x(i+1,j)-x(i-1,j))/2;x2=(x(i,j+1)-x(i,j-1))/2;x12=(x(i+1,j+1)-x(i+1,j-1)+x(i-1,j-1)-x(i-1,j+1))/4;y1=(y(i+1,j)-y(i-1,j))/2;y2=(y(i,j+1)-y(i,j-1))/2;y12=(y(i+1,j+1)-y(i+1,j-1)+y(i-1,j-1)-y(i-1,j+1))/4;alpha=x2^2+y2^2;beta=x1*x2+y1*y2;gama=x1^2+y1^2;tempx=(alpha*(x(i-1,j)+x(i+1,j))+gama*(x(i,j+1)+x(i,j-1))-2*beta*x12)/(2*(alpha+gama))-x(i,j);tempy=(alpha*(y(i-1,j)+y(i+1,j))+gama*(y(i,j+1)+y(i,j-1))-2*beta*y12)/(2*(alpha+gama))-y(i,j);x(i,j)=x(i,j)+w*tempx;y(i,j)=y(i,j)+w*tempy;e1=max(abs(tempx),abs(tempy))*w;err=max(e1,err);endendiferr

下载提示(请认真阅读)1.请仔细阅读文档，确保文档完整性，对于不预览、不比对内容而直接下载带来的问题本站不予受理。

2.下载的文档，不会出现我们的网址水印。

下载文档到电脑，查找使用更方便

8 积分

还剩页未读，继续阅读 关 键 词：MATLAB拟合 MATLAB 坐标系 matlab拟合 Matlab 拟合程序matlab程序 matlab Matlab示例程序 示例程序

  蚂蚁文库所有资源均是用户自行上传分享，仅供网友学习交流，未经上传用户书面授权，请勿作他用。

关于本文