加频偏相偏

公式表示

向信号中添加频偏相偏，如下公式所示：

$t_o(n)=(I+jQ)e^{j(2\pi f_{0}t+\theta )}$

其中$f_0$为频偏，$\theta$为相偏。上式中$t=n{T}_s$,$T_s$为采样间隔（采样率$F_s$的倒数），$(I+jQ)$为iq两路信号。

实现

Fs = 200000  #200khz
f0 = 100000  #假设频偏100khz
theta = 90 /360*2*np.pi   #90度 
y = x[:,0]+x[:,1]*1j        #复值化
y = y * np.exp(1j * (2 * np.pi * f0 / Fs * np.arange(1024)+theta))  #频偏+相偏

注：
1、最好使用np中的exp，不要使用math中的。因为np中的更好结合矩阵进行运算。

2、range()返回的是range对象，主要在for循环中使用，其他情况不建议使用。
np.arange在需要模拟长串规律数列参与计算时更加好用。
range与arange 学习链接