%------------------------------------------------------------------%
%---         程序名稱：Radar正交變換Matlab仿真       ---%
%---         前端帶奇偶抽樣。                                     ---%
%---             作者：王旭東                                        ---%
%---             時間：2003年6月25日                          ---%
%---             版本：v0000.00.02                               ---%
%---         版權所有違者必究！                                  ---%
%-----------------------------------------------------------------%
%-------------------------- 定義常數 --------------------------%
f=28;
fs=100; 
signal_start=10;
signal_end=180;
total_dot=200;
dets=1/fs;
PI=3.1415926;
aa=0.5846832;
bb=0.138025;
Am=1;
%------------------------------------------------------------------%
%------- 頻率范圍：B/2<--->3*B/2 其中 B=fs/4 -----------%
B=fs/4;
if (f>=B/2)&(f<=3*B/2)
    disp('%--- Start Simulation of Radar_Hilbert ---%');
else
    disp('%- ERROR !: f exceed the required scope! -%');
end;
%-------------------------------------------------------%
%---------------------- 生成噪聲 --------------------%
%noise=randn(1,total_dot);  %產生高斯隨機噪聲。%注：每次-
                            %調用產生的噪聲不同。                            
noise=zeros(1,total_dot); %不加噪聲。
MaxNoise=max(abs(noise));
%-------------------------------------------------------%
%--------------------- 生成仿真信號 ----------------%
for i=1:signal_start-1
    data_in(i)=0;
end;

for i=signal_start:signal_end
    data_in(i)=Am*sin(2*PI*f*(i-signal_start)*dets);
end;
for i=signal_end+1:total_dot
    data_in(i)=0;
end;
data_in=(data_in+noise)/(Am+MaxNoise);
%-------------------------------------------------------%
%--------------------- 理論結果 ---------------------%
IQ_real=hilbert(data_in);
I_real=real(IQ_real);
Q_real=imag(IQ_real);
I_reals4_0=zeros(1,total_dot/4);
Q_reals4_0=zeros(1,total_dot/4);
for i=1:4:total_dot-3
    I_reals4_0((i+3)/4)=I_real(i);
    Q_reals4_0((i+3)/4)=Q_real(i);
end;
%-------------------------------------------------------%    
%--------------------- 仿真 I 通道 -------------------%
data_in_odd=zeros(1,total_dot/2);%-- 前端帶奇偶抽樣。
for odd_i=1:2:total_dot-1
    data_in_odd((odd_i+1)/2)=data_in(odd_i);
end;
for n=1:total_dot/2
    if n==1 
        I(n)=-aa*data_in_odd(1);
    else
        I(n)=aa*(I(n-1)-data_in_odd(n))+data_in_odd(n-1);
    end;
end;

%Is2_odd=zeros(1,total_dot/4);%-- Is2 odd even --%
%for i=1:2:total_dot/2-1
%    Is2_odd((i+1)/2)=I(i);
%end;
%Is2_even=zeros(1,total_dot/4);
%for i=2:2:total_dot/2
%    Is2_even(i/2)=I(i);
%end;

Is4_0=zeros(1,total_dot/8);%-- Is4 0,1,2,3 --%
for i=1:4:total_dot/2-3
    Is4_0((i+3)/4)=I(i);
end;
%Is4_1=zeros(1,total_dot/8);
%for i=2:4:total_dot/2-2
%    Is4_1((i+2)/4)=I(i);
%end;   
%Is4_2=zeros(1,total_dot/8);
%for i=3:4:total_dot/2-1
%    Is4_2((i+1)/4)=I(i);
%end;
%Is4_3=zeros(1,total_dot/8);
%for i=4:4:total_dot/2
%    Is4_3(i/4)=I(i);
%end;    
%-------------------------------------------------------%
%--------------------- 仿真 Q 通道 -----------------%
data_in_even=zeros(1,total_dot/2);%-- 前端帶奇偶抽樣。
for even_i=2:2:total_dot
    data_in_even(even_i/2)=data_in(even_i);
end;
for n=1:total_dot/2
    if n==1
        Q(1)=bb*data_in_even(1);
    else
        Q(n)=bb*(Q(n-1)+data_in_even(n))+data_in_even(n-1);
    end;
end;

%Qs2_odd=zeros(1,total_dot/4);%-- Qs2 odd even --%
%for i=1:2:total_dot/2-1
%    Qs2_odd((i+1)/2)=Q(i);
%end;
%Qs2_even=zeros(1,total_dot/4);
%for i=2:2:total_dot/2
%    Qs2_even(i/2)=Q(i);
%end;

Qs4_0=zeros(1,total_dot/8);%-- Qs4 0,1,2,3 --%
for i=1:4:total_dot/2-3
    Qs4_0((i+3)/4)=Q(i);
end;
%Qs4_1=zeros(1,total_dot/8);
%for i=2:4:total_dot/2-2
%    Qs4_1((i+2)/4)=Q(i);
%end; 
%Qs4_2=zeros(1,total_dot/8);
%for i=3:4:total_dot/2-1
%    Qs4_2((i+1)/4)=Q(i);
%end; 
%Qs4_3=zeros(1,total_dot/8);
%for i=4:4:total_dot/2
%    Qs4_3((i)/4)=Q(i);
%end; 
%-------------------------------------------------------%
%-------------------- 圖形化輸出 -------------------%
%figure(1);
%plot(data_in);
%title '原始輸入數據';

%n=1:total_dot/2;
%figure(2);
%plot(n,I,n,Q);
%title '正交變換原始輸出';

%i=1:total_dot/4;%----------Is2,Qs2----------%
%figure(3);
%plot(i,Is2_odd,i,Qs2_odd);
%title '正交變換2采樣輸出（odd-odd）';

%figure(4);
%plot(i,Is2_odd,i,Qs2_even);
%title '正交變換2采樣輸出（odd-even）';

%figure(5);
%plot(i,Is2_even,i,Qs2_odd);
%title '正交變換2采樣輸出（even-odd）';

%figure(6);
%plot(i,Is2_even,i,Qs2_even);
%title '正交變換2采樣輸出（even-even）';

j=1:total_dot/8;%--------Is4,Qs4-----------%
figure(7);
plot(j,Is4_0,'-b',j,Qs4_0,'-r');
title '正交變換4采樣輸出（0-0）';

%figure(8);
%plot(j,Is4_0,j,Qs4_1);
%title '正交變換4采樣輸出（0-1）';

%figure(9);
%plot(j,Is4_0,j,Qs4_2);
%title '正交變換4采樣輸出（0-2）';

%figure(10);
%plot(j,Is4_0,j,Qs4_3);
%title '正交變換4采樣輸出（0-3）';

%figure(11);
%plot(j,Is4_1,j,Qs4_0);
%title '正交變換4采樣輸出（1-0）';

%figure(12);
%plot(j,Is4_1,j,Qs4_1);
%title '正交變換4采樣輸出（1-1）';

%figure(13);
%plot(j,Is4_1,j,Qs4_2);
%title '正交變換4采樣輸出（1-2）';

%figure(14);
%plot(j,Is4_1,j,Qs4_3);
%title '正交變換4采樣輸出（1-3）';

%figure(15);
%plot(j,Is4_2,j,Qs4_0);
%title '正交變換4采樣輸出（2-0）';

%figure(16);
%plot(j,Is4_2,j,Qs4_1);
%title '正交變換4采樣輸出（2-1）';

%figure(17);
%plot(j,Is4_2,j,Qs4_2);
%title '正交變換4采樣輸出（2-2）';

%figure(18);
%plot(j,Is4_2,j,Qs4_3);
%title '正交變換4采樣輸出（2-3）';

%figure(19);
%plot(j,Is4_3,j,Qs4_0);
%title '正交變換4采樣輸出（3-0）';

%figure(20);
%plot(j,Is4_3,j,Qs4_1);
%title '正交變換4采樣輸出（3-1）';

%figure(21);
%plot(j,Is4_3,j,Qs4_2);
%title '正交變換4采樣輸出（3-2）';

%figure(22);
%plot(j,Is4_3,j,Qs4_3);
%title '正交變換4采樣輸出（3-3）';
k=1:total_dot/4;
figure(23);
plot(k,I_reals4_0,'-b',k,Q_reals4_0,'-r');
axis([0,25,-1,1])
title '正交變換的理論值';