test.m

rake接收机matlab实现 完成相关接受的信道相关估计算法

clear;
prompt={'Simulink Time:','First Path Attenuation:','Second Path Attenuation:','Channel Delay Time:','Simple Time of realmultipath'};
name=['Multipath Test Dimension'];
line=1;
defaultanswer={'10','0.8','0.6', '0.007','0.310'};
glabel=inputdlg(prompt,name,line,defaultanswer);
time=str2num(char(glabel(1,1)));
am1=str2num(char(glabel(2,1)));
am2=str2num(char(glabel(3,1)));
de=str2num(char(glabel(4,1)));
samtime=str2num(char(glabel(5,1)));
hold off;
subplot(1,2,1);
sim('Rakesimulink.mdl') ;
mit=samtime;
x=0:mit:(length(attunation1)-1)*mit;
plot(x,attunation1);
hold on;
plot(x,attunation2,'-.');
xlabel('time/s');
ylabel('attunation/%');
title('attunation');
axis([0 time 0 1]);
subplot(1,2,2);
plot(x,delay);
xlabel('time/s');
ylabel('delay/s');
title('delay');
axis([0 time 0 5*de]);
u=0;
for i=1:length(delay)
    if delay(i)>0
        u=u+1;
        if u==1
            position=i;
        end
    else u=u;
    end
end
Psum1=sum(attunation1);
Ave_att1=Psum1/u;
Psum2=sum(attunation2);
Ave_att2=Psum2/u;
Psum3=sum(delay);
Ave_delay=Psum3/u;
Time_discover=(position-2)*mit;
Rate_search=u/length(delay);
disp('第一径衰减:');
disp(Ave_att1);
disp('第二径衰减：');
disp(Ave_att2);
disp('信道时延：');
disp(Ave_delay);
disp('搜索时间：');
disp(Time_discover);
disp('搜索效率：');
disp(Rate_search);