% 多徑型信號產生，分析

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