資源簡介
實現用matelb編程的三維形貌圖程序,利用傅里葉變換自動生成頻譜圖
代碼片段和文件信息
clear
%function?f=roughs3(NrLsigmalc);
%M=N??%total?number?of?sample?points?%為了作fft,N取2的n次方個值
f=roughs3(6420.20.1);%64為取得全部數據點,2為rl-粗糙長度,0.2為方均根的高度-rms?height0.1為相關長度lc;rl=lc*20;
Skl=zeros(64);%Skl=zeros(M)產生M*M的矩陣
for?m=0:64/2-1;%for?m=0:N/2-1
????kx?=?m*2*pi/2;%kx?=?m*2*pi/rL;
????ky=(0:64/2-1)*2*pi/2;%ky=(0:M/2-1)*2*pi/rL;
????Skl(m+1(0:64/2-1)+1)=(0.2^2*0.1^2/(4*pi))*exp(-(kx^2+ky.^2)*(0.1/2)^2);%Skl(m+1(0:M/2-1)+1)=(sigma^2*lc^2/(4*pi))*exp(-(kx^2+ky.^2)*(lc/2)^2);%2D-GAUSS譜密度;
end;
Skl(64-(1:64/2-1)+164-(1:64/2-1)+1)=Skl((1:64/2-1)+1(1:64/2-1)+1);%Skl(N-(1:N/2-1)+1M-(1:M/2-1)+1)=Skl((1:N/2-1)+1(1:M/2-1)+1);
Skl((0:64/2)+164-(1:64/2-1)+1)=Skl((0:64/2)+1(1:64/2-1)+1);%??Skl((0:N/2)+1M-(1:M/2-1)+1)=Skl((0:N/2)+1(1:M/2-1)+1);
Skl(64-(1:64/2-1)+1(0:64/2)+1)=Skl((1:64/2-1)+1(0:64/2)+1);%Skl(N-(1:N/2-1)+1(0:M/2)+1)=Skl((1:N/2-1
- 上一篇:LT碼的編譯碼程序MATLAB
- 下一篇:CNN特征圖可視化
評論
共有 條評論