0

所以我有一个简单的程序，我想在其中移动一个简单的正弦函数 pi/2。现在我知道这将非常容易，只需插入一个偏差（即 A*sin(2*pi*frequency + bias)。但这个程序是测试理论的简单方法。我需要移动复杂的磁数据，但是shift 与频率有关。所以要弄清楚如何做到这一点，我只想将这个正弦波移动一个设定的移位，但我想在频域中进行。虽然下面的代码没有显示任何错误，但它确实无法正确移动数据并影响幅度。代码如下。谢谢！

clear all
time = 1:0.01:2*pi; %Create time vector
mag = sin(2*pi*time);
Y = fft(mag); %transform
Pnewr = pi/2;
%t = angle(mag);
t = imag(Y);
%t_fin = t-Pnewr; %Subtract the pahse delay from the original phase vector
R=real(Y);
I=t;
k = I./R
Phi = tan(k);
PhiFinal = Phi-Pnewr;
PhiFinal = PhiFinal'
IFinal = R * atan(PhiFinal);
spec=complex(R,IFinal);
Finalspec = ifft(spec); %Invert the transform
Final = Finalspec;
plot(time,mag);
hold on
plot(time,Final,'r')
grid on
4

1 回答 1

0

一方面，您没有正确地重新组合虚构和真实的组件，因为

PhiFinal = PhiFinal'
IFinal = R * atan(PhiFinal);

实际上是一个点积，而不是逐个元素的积。一般来说，您似乎没有正确使用复杂的关系。下面产生一个干净的相移：

time = 1:0.01:2*pi; 
mag = sin(2*pi*time);
Y = fft(mag); %transform

Pnewr = pi/2;

R = real(Y);
I = imag(Y);

It = abs(Y);
% It = sqrt(I.^2 + R.^2);

Phi= angle(Y); % <-- using matlab function, equivalent to `Phi = atan2(I, R);`

k = I./R;
Phi0 = atan(k);

figure, subplot(121), plot(Phi0,Phi,'.'), xlabel('Phi from tan'), ylabel('Phi from ''angle'''), grid on, axis('tight')

PhiFinal = Phi-Pnewr; % <-- phase shift

IFinal = It .* sin(PhiFinal);
RFinal = It .* cos(PhiFinal);
spec= RFinal + 1i*IFinal;
Final = ifft(spec); %Invert the transform

subplot(122)
plot(time,mag);
hold on
plot(time,real(Final),'r')
plot(time,imag(Final),'r:')
grid on
axis('tight')
legend('Initial','Final'),xlabel('t'), ylabel('I')

mag*mag'   % <-- check that total power is conserved
Final*Final'

这些图显示了使用tanvs matlab angle（使用atan2）计算的相位，以及相移的结果（右图）：

在此处输入图像描述

于 2013-08-14T09:59:25.553 回答