そのため、いくつかの磁気データから機器の応答 (位相遅れ) を除去しようとしています。私が使用しようとしているコードは以下のとおりです。ファイルを一緒に正しく再構築する方法がわかりません。私が現在試みているアプローチは、位相遅れを差し引いて信号を再構築することです。彼女を使うことができる別の機能があるかもしれないことに気づきました。位相回転を提案してもらいましたが、私の位相遅れは周波数に依存しているため、何らかの方法でベクトルをたたみ込むことができるようにする必要があります...助けてくれてありがとう!
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')