在Matlab中查找FM信號的相位

Question

我想找到FM信號的相位。 消息信號的Matlab代碼是

m= ones(1,Fs); 
m= [m(1:round(.4*Fs))*1, m(round(.4*Fs)+1:round(.7*Fs))*-2 ...
    m(round(.7*Fs)+1:Fs)*0];

如何找到調頻信號的相位？

Answer 1

您需要在代碼中進行以下更正：

Kf=50; %%% increase Kf no notice effect
Fc=400;
Fs=5000; %%% Increase Fs. Nyquist criterion!
t=linspace(0,1,Fs);
m= ones(1,Fs);
m= [m(1:round(.4*Fs))*1 m(round(.4*Fs)+1:round(.7*Fs))*-2 m(round(.7*Fs)+1:Fs)*0];
fi_t =2*pi*Kf*cumtrapz(t,m);%%% use cumtrapz, not cum. Reverse order
u=10*cos(2*pi*Fc*t+fi_t);
plot(t,u);