clc,clear,close all
nz=230;
nx=640;
dx=10;
dz=10;
dk2=1/dx/nx;
dk1=1/dz/nz;
k2=-floor(nx/2)*dk2:dk2:floor(nx/2)*dk2;
k1=-floor(nz/2)*dk1:dk1:floor(nz/2)*dk1;
facc=fopen('data.dat','rb');
acc=fread(facc,[nz,nx],'float');
fclose(facc);
Facc=abs(fftshift(fft2(acc)));
Tacc=log(Facc+1);
imagesc(k2,k1,Facc);colormap(jet);
%mx=max(max(Tacc));mn=min(min(Tacc));
caxis([0 1]);
c=colorbar;%('position',[0.91 0.585 0.04 0.295]);
figure_FontSize=16;
set(get(c,'label'),'string','','FontSize',figure_FontSize,'FontName','Times New Roman')
set(c,'YTick',0:0.5:1.0)
set(c,'yTickLabel',num2str(get(c,'yTick')','%.1f'))
set(c,'FontSize',figure_FontSize);
set(gca,'FontName','Times New Roman','FontSize',figure_FontSize);
set(gca,'xaxislocation','top','FontWeight','bold');
xlabel('Wavenumber (1/m)');ylabel('Wavenumber (1/m)');
set(get(gca,'XLabel'),'FontSize',figure_FontSize);
set(get(gca,'YLabel'),'FontSize',figure_FontSize);
set(gca,'XLim',[-floor(nx/2)*dk2 floor(nx/2)*dk2]);
set(gca,'XTick',[-floor(nx/2)*dk2:dk2*300:floor(nx/2)*dk2],'FontSize',figure_FontSize);
set(gca,'xTickLabel',num2str(get(gca,'xTick')','%.2f'),'FontSize',figure_FontSize);
set(gca,'YLim',[-floor(nz/2)*dk1 floor(nz/2)*dk1]);
set(gca,'YTick',[-floor(nx/2)*dk2:dk2*300:floor(nx/2)*dk2],'FontSize',figure_FontSize);
set(gca,'yTickLabel',num2str(get(gca,'yTick')','%.2f'),'FontSize',figure_FontSize)
MATLAB绘制波数谱
于 2023-05-25 16:19:53 首次发布