1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
| clear all;
clc;
L=200;
B=20*10;%Bandwidth
Ts=1/B;%sampling time
a=1;%amplitude
sigmat=125*Ts;
v=exp(-Ts/sigmat);%variance of noise
std=sqrt(v);%standard deviation
sum=0;
num=0;
den=0;
h=[];
nb=5000;
for L=100:100:300
res=[];
for Lr=1:20
for i=1:nb
for j=1:L
h(j)=(std^(j-1))*gaus_noise(1)/sqrt(2);
den=den+abs(h(j))^2;
end
for k=1:Lr
num=num+abs(h(k))^2;
end
rpe=num/den;
sum=sum+rpe;
end
res=[res sum/nb];
end
figure,plot([1:20],res);
xaxis('Lr');
yaxis('E(rp)')
legend(''100','200','300');
end |