求助！FFT比值校正方法数据的处理

S0704082

本人编写了一个MATLAB频谱校正的程序，信号有两个频率，但是两个频率下幅值和相位为多组，我要求的是将一个频率下的所有幅值和所有相位分别放在一个矩阵里，该如何处理，循环的单个程序已经运行出来。我的信号为
fs=5120;
N=1024;
t=(0:N-1)/fs;
A1=[5 10 6 8 1 3];A2=[1.5 9 4.6 8.5 7.5 6.1];theta1=[10 20 30 40 50 60];theta2=[70 80 90 100 110 120];
for i=1:6
s=A1(i)*cos(2*pi*250.658*t+theta1(i)*pi/180)+A2(i)*cos(2*pi*600.52*t+theta2(i)*pi/180);
我的要求就是频谱校正后，诊断的结果A11、A22为诊断结果，和A1和A2对应