%初始化
clear all
close all;
clc;
N=512; %训练符号总长度或数据符号长度
N_symbol=10;
L=128; %训环前缀长度
ebn0=15;
sr=250000;
ml=2;
br=sr*ml;
deltad=0
deltaf=2.45
%产生train
pn=rand(1,N/2)>0.5;
pn=reshape(pn,N/4,2);
[ich,qch]=qpskmod(pn,N/4,1,ml);
kmod=sqrt(2);
x=ich*kmod+qch*kmod*i;
y=ifft(x);
y=reshape(y,1,N/4);
train=[y,y,-y,-y];
%生成信号
input_stream=rand(1,N*N_symbol*ml)>0.5;
%串并变换
para_stream=reshape(input_stream,N,N_symbol*ml);
%qpsk调制
[ich,qch]=qpskmod(para_stream,N,N_symbol,ml);
x=ich*kmod+qch*kmod*i;
%ifft变换
y=ifft(x);
train=reshape(train,N,1);
y(:,1)=train;
ich=real(y);
qch=imag(y);
%加入训环前缀
[ich,qch]=giins(ich,qch,N,L,N_symbol);
%Attenuation Calculation 衰减计算
spow=sum(ich.^2+qch.^2)/N/N_symbol;
attn=0.5*spow*sr/br*10.^(-ebn0/10);
attn=sqrt(attn);
%加时延
[ich,qch]=delay(ich,qch,length(ich),deltad);
y=ich+i*qch;
for k=1:length(ich)
y(k)=y(k)*exp(j*2*pi*deltaf*k/N); %加频移
end
ich=real(y);
qch=imag(y);
%加人AWNG
[ich,qch]=comb(ich,qch,attn);
y=ich+i*qch;
figure();
plot(y,'.');
%计算公式
for d=1:N
for n=1:N/4
z10(n)=conj(y(n+d))*y(n+d+N/4);
end
for n=1:N/4
z11(n)=conj(y(n+d+N/2))*y(n+d+N/4+N/2);
end
p1(d)=sum(z10)+sum(z11);
for n=1:N/2
z20(n)=abs(y(n+d+N/4)).^2;
end
for n=1:N/2
z21(n)=abs(y(n+d+N/4+N/2)).^2;
end
p2(d)=sum(z20)+sum(z21);
end
figure();
plot(p1)
figure();
plot(p2)
p=abs(p1).^2./(p2.^2);
figure();
plot(p);
grid;
[a,b]=max(p)
基于Matlab实现Minn定时同步.zip
版权申诉
176 浏览量
2023-06-13
22:27:13
上传
评论
收藏 10KB ZIP 举报
基于Matlab实现Minn定时同步.zip （7个子文件） 
基于Matlab实现Minn定时同步 
giins.m 846B
1.png 6KB delay.m 741B Park_Cheon.asv 1KB qpskmod.m 893B Minn.m 2KB comb.m 258B资源评论
Matlab科研辅导帮
- 粉丝: 1w+
- 资源: 7474
