普通磨削、超声辅助磨削表面形貌代码及表面轮廓提取

进一步了解+qq1937565361

clear
clc

N=1000;%生成大小
delta=0.4;%表面均方根粗糙度
betax=200;
betay=10;%
C=6;%功率谱密度,

L=10;%
dx=L/N;dy=dx;
NN=-N/2:N/2-1;
[Nx,Ny]=meshgrid(NN,NN);%生成平面网格点
taux=dx.*Nx;tauy=dy.*Ny;

%%生成具有指定自相关函数的粗糙表面
eta=randn(N,N);%高斯分布白噪声
A=fft2(eta);%傅里叶变换
R=zeros(N,N);%生成N*N的零矩阵
R=delta^2*exp(-2.3*((taux/betax).^2+(tauy/betay).^2).^0.5);%自相关函数
Gz=1/(2*pi^2).*fft2(R);%功率谱密度函数
H=(Gz/C).^0.5;%传递函数
Z=H.*A;%表面高度的傅里叶变换
z=ifft2(Z);%表面高度分布
z = abs(z) * 100;
figure(1);
surf(abs(z));% surf2stl('surf_roughness.stl',1,1,z) % 生成模型
shading interp
lighting phong