MATLAB绘制波数谱

吃最辣的锅

已于 2023-05-25 16:22:42 修改

阅读量571

点赞数

文章标签： matlab

于 2023-05-25 16:19:53 首次发布

本文链接：https://blog.csdn.net/qq_37611077/article/details/130869860

版权

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)