本帖最后由 娇娇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