main.m

2psk的调制以及误码率等系统性能比较

% 多径型信号产生，分析

clear;

% 仿真基本参数
M=10^4;                 % T内的采样点数，太少则时延量化精度不够。
T=1/(24.3*1000);        % 符号时宽
h=T/M;                  % 采样间隔，步长
fc=1.9*10^9;            % 载波频率
w=2*3.1415926535897932*fc;
Np=10;                 % 路径数
DTR=0.1;               % 最大时延差Dt与符号时宽T的比值  
Dt=T*DTR;              % 最大时延差,绝对值。实际，+，-，超前，滞后
Ts=T-2*Dt;              % 基本信号时宽

% 基本信号g0(t)的产生
Ms=round(Ts/h);         % 信号的采样点数。向最接近的整数取整
g0=sin_wav(Ts,1,0,Ms);  % 正弦波。 持续期，谐波次数，初相位，采样点数。能量归一。
%Vf=var(f)
M_Dt=round(Dt/h);
g0c=shift_r(g0,Ms,M_Dt,M);
pwr=pwr_av(g0c,M);
g0c=g0c/sqrt(pwr*T);
%R_g0c=projct(g0c,g0c,T,M);
%subplot(4,1,1); plot(g0);
%subplot(4,1,2); plot(g0c);

% 多径信号产生
alf0=0.1;
num_g=2;
for k=1:num_g
    tao=unifrnd(-Dt,Dt,1,Np);   % 时延差。1行，Np列，-Dt到Dt 连续均匀分布的随机数。
    M_tao=round((Dt+tao)/h);    % 移动的采样点数
    alf=unifrnd(alf0,1,1,Np);   % 衰减因子。1行，Np列，alf0到1 连续均匀分布的随机数。
    g(k,1:M)=0;
    for i=1:Np
        a(i)=alf(i)*cos(w*tao(i));
        p=shift_r(g0,Ms,M_tao(i),M);
        g(k,:)=g(k,:)+p*a(i);
        pwr=pwr_av(g(k,:),M);
        g(k,:)=g(k,:)/sqrt(pwr*T);   % 能量归一
    end
end    
%subplot(4,1,3); plot(g(1,:));
%subplot(4,1,4); plot(g(2,:));

% 求内积
for k=1:num_g
    R0(k)=projct(g(k,:),g0c,T,M);
end
for k=1:num_g
    for i=1:num_g
        Rg(k,i)=projct(g(k,:),g(i,:),T,M);
    end
end
R0
Rg