读取相关谱并画出频谱图

numCh = 4
corFile = 'E:\TEST\MATLAB\fengchao\corData\11.bin';
fd = fopen(corFile, 'r');
corData = fread(fd, [512, numCh*2000], 'float32');%读出数据
ncol = 512/2 - 10;% 去掉前面5个和后面5个数据
sumCorr = zeros(2000, ncol, numCh); % 2000*ncol*4
corData = corData(1:2:end, :) + 1j*corData(2:2:end, :);% 相邻两个数组成复数
nn = 5;
for i = 1:numCh
    sumCorr(:, :, i) = corData(nn+1:256-nn, i:numCh:end)';
end
fclose(fd);
%..................数据读取完毕...................%

for chn = 1:numCh
    amp(:,:,chn) = 10*log10(abs(sumCorr(:,:,chn)));
    amp(:,1:2,chn) = 55;
    figure(201)
    subplot(2,3,chn)
    image(amp(:,:,chn),'cdatamapping','scaled')
    titleName = strcat('幅度谱',num2str(chn));
    title(titleName);
    phase(:,:,chn) = angle(sumCorr(:,:,chn));
    figure(202)
    subplot(2,3,chn)
    image(phase(:,:,chn),'cdatamapping','scaled')
    titleName = strcat('相位图',num2str(chn));
    title(titleName);
end