matlab中求复数的相位,怎么求复数相位

本帖最后由 娇娇9228 于 2017-12-11 20:03 编辑

宋老师，非常感谢您的耐心解答。我还有一个问题想要请教您。OFDM脉冲在发送时，会在后面加上一个循环后缀码，这样在接收端估计信号开始位置时，当有E个采样点的误差时，能通过接收端数据的傅里叶变换然后求相位算出E.

没有估计误差的公式和有估计误差的公式如下，这是我从论文里截图的，已经查过好多文件，公式应该没有错误。  接收端没有估计误差

1.png (10.26 KB, 下载次数: 3)

接收端没有估计误差

2017-12-11 20:01 上传

接收端有E个估计误差

2.png (14.03 KB, 下载次数: 1)

2017-12-11 20:01 上传

之前的那个问题就是有E个估计误差，求E。根据公式算出来是对的。但是我想要实现OFDM 这个特性就不对了，总求不出E来。代码如下：

OFDM脉冲的生成：将0-24kHZ的带宽分成64个子载波，每个子载波的频率是375HZ，我只想要18-20kHZ发送数据(1或者-1)，其他子载波上的数据都是零。每个子载波只发送一个数据，载波是cos波，因为cos0 = 1，所以18-20khz(48-54个子载波)上还是1或者-1. 不知道这样生成的算OFDM脉冲吗？ 将组成的64个数据进行逆傅里叶变换，并将每个脉冲的前20个数据作为循环后缀，组成发送序列。在接收端经过傅里叶变换算误差E，我现在的问题是这个E总算不出来。

%发送序列生成

data_initial = zeros([1 64]);

num=7;            %总共7个数

rand('seed',0);   %使每次产生相同的随机数

a=rand(1,num);   %生成0~1随机数序列

a(a>0.5)=1;

a(a<=0.5)=-1;

data_initial(48:54) = a;   %发送的序列

data_ifft = ifft(data_initial); %逆傅里叶变化为复数。

add = zeros(1,200);   %200个采样的间隔

cs = data_ifft(1:20);  %前20个值为循环序列

data_send = [data_ifft cs add]; %发送序列

data_receive1 = fft(data_send(1:64),64);   %模拟接收端没有误差

E =  3;

data_receive = fft(data_send(1+E : 64+E),64);  %模拟接收端有E个误差

phase_after = angle(data_receive);

phase_before = angle(data_receive1);

phase_dif = phase_after - phase_before;

phase_dif1 = unwrap(phase_dif);

for n = 1:64

E (n) = phase_dif1(n)/n /(2*pi)*64;

end