概述
clear all;
clc;
%%%%%%%%%%%首先模拟被检测信号产生待采样信号,在本次仿真中,设定被检测信号为BPSK调制,fc=1Mhz
%%%%%%%%%%%基带信号长度为100bit,信号比特率fb=1Mbps,采样频率fs=100Msps
signal_base=randint(1,100,2);
signal_base_bpsk=2*signal_base-1;
T=1/1000000;
t0=0:(T/100):(99*T/100);
for i=1:100,
if signal_base_bpsk(i)>0
signal_base_pure(((i-1)*100+1):((i)*100))=sin(2*pi*1000000*t0);
else
signal_base_pure(((i-1)*100+1):((i)*100))=-sin(2*pi*1000000*t0);
end;
end;
signal_base_awgn = awgn(signal_base_pure,0);
%%%%%%%%%将在fs为100Msps下采集得到的数据,进行三倍内插%%%%%%%%%%%%%%
t1=0:(T/100):(100*T-T/100);
s1=signal_base_awgn ;
t2=0:(T/300):(100*T-T/100);
s2=interp1(t1,s1,t2);%%%%3倍线性内插
%%%%%%%%%%分别对内插后数据序列s2,延迟100Ts数据序列s3,延迟250Ts数据序列s4,
s3_temp=[zeros(1,100) s2];
s3=s3_temp(1:2048);
s4_temp=[zeros(1,250) s2];
s4=s4_temp(1:2048);
s6_temp=[zeros(1,120) s2];
s6=s6_temp(1:2048);
s7_temp=[zeros(1,300) s2];
s7=s7_temp(1:2048);
s5=s2(1:2048);
%%%%%%%%%%分别计算s5与s3、s4之间的互相关函数
N=2048;
r35=zeros(1,500);
for m=1:500,
for n=1:N-m,
r35(m)=r35(m)+s5(n)*s3(n+m);
end;
end;
r35=r35/N;
r45=zeros(1,500);
for m=1:500,
for n=1:N-m,
r45(m)=r45(m)+s5(n)*s4(n+m);
end;
end;
r45=r45/N;
r65=zeros(1,500);
for m=1:500,
for n=1:N-m,
r65(m)=r65(m)+s5(n)*s6(n+m);
end;
end;
r65=r65/N;
r75=zeros(1,500);
for m=1:500,
for n=1:N-m,
r75(m)=r75(m)+s5(n)*s7(n+m);
end;
end;
r75=r75/N;
最后
以上就是慈祥冷风为你收集整理的信号延迟仿真源码的全部内容,希望文章能够帮你解决信号延迟仿真源码所遇到的程序开发问题。
如果觉得靠谱客网站的内容还不错,欢迎将靠谱客网站推荐给程序员好友。
发表评论 取消回复