91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

  • 大小: 2KB
    文件類型: .m
    金幣: 1
    下載: 0 次
    發布日期: 2021-05-08
  • 語言: Matlab
  • 標簽:

資源簡介

實現用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

評論

共有 條評論

相關資源