所以我试图从一些磁性数据中删除仪器响应(相位滞后)。我试图使用的代码如下。我无法想象如何正确地重建文件。我现在尝试的方法是减去相位滞后并重建信号。我意识到可能还有其他功能可以使用她。我有相位旋转建议,但我的相位滞后是频率依赖的,所以我真的需要能够以某种方式卷积矢量...谢谢你的帮助!
load('coil') % Load coil response data
time = 1/40:1/40:288000/40; %Create time vector for later plotting
%Removing Coil Response
T2=load(mag); %Load data
y=abs(real(T2.y.data));
t = angle(T2.y.data); %Create a phase vector for the original data
Fs = 40; %Sampling frequency
NFFT = size(y);
NFFT=NFFT(1);
j = [y,t]; %Create a matrix of amplitude and phase
Y = fft2(j); %2D transform
F = ((0:1./NFFT:1-1./NFFT)*Fs);
F1=coil(:,1); %Load coil phase delay
P=coil(:,3);
Pnew=interp1(F1,P,F); %Create the phase delay vector of the same length
Pnew(1) = Pnew(2);
Pnew=Pnew';
t_fin = t-Pnew %Subtract the pahse delay from the original phase vector
%Mag_Final = ifft(abs(Y));
%Phase_Final = ifft(phaseF);
j_fin = [y,t_fin]; %Create matrix of amplitude versus correct phase
Final = ifft(j_fin); %Invert the transform
y = struct('data',Final);
save('y')