MATLAB FFT ILA

%测试频谱
close all;clear all;tic;%关闭图像以及清空数据，不让影响数据的初始化

%aa = load('1M_26.prn');%将数据导入matlab，注意matlab文件需要与CSV文件在同一目录下，否则需要指定位置aa = csvread('C://iladata.csv');
aa = csvread('iladata.csv');
% used=aa(:,3);
used_i=aa(:,1);%提取某一列，这个数据对应ILA中的一列
used_q=aa(:,2);%提取某一列，这个数据对应ILA中的一列
used=used_i+ j*used_q;%j*used_q;
% A（i，：） 提取矩阵A的第 i行
% A（：，i） 提取矩阵A的第 i列
Fs=110e6;   %采样频率
P=11e6;
N=10;
L=2^N;     %数据点数
A=1;
t=1:1:Fs/10;
% theta=90;
% sig=A*exp(j*2*pi*t/Fs + theta);

times=1:1:length(used);
% carrier=exp(-j*2*pi*times/Fs*8e6);
%carrier=exp(-j*2*pi*times/Fs);

NFFT = 2^nextpow2(L); % Next power of 2 from length of y           
% Y = fft(carrier.'.*used,NFFT)/L;    %  FFT   
 Y = fft(used.^2,NFFT)/L;    %  FFT   
AMP = fftshift(abs(Y));
WIDTH = linspace(-Fs/2, Fs/2, L);%频率坐标，单位Hz
figure(1)
plot(WIDTH,10*log10(abs(AMP)));
title('信号的频谱');
xlabel('频率(MHz)');
grid on;
% hold on
% stem(max(10*log(abs(AMP))))
abcd=atan(used_q/used_i);
figure(2)
plot(abcd);
title('角度');
xlabel('角度值');

scatterplot(used);

toc;