cycle_shift_delay17a.m

是自己自创的好程序

close all;
clear all;
tic
Tb1=[0,2,5];
Pwd1=[0,-5,-10];
Tb2=[0,4,9];
Pwd2=[0,-8,-15];

Tb3=[0,10,19,23];
Pwd3=[0,-12,-28,-30];

Tb4=[0,1,8,10];
Pwd4=[0,-8,-12,-25];
% Tb=[0];
% Pwd=[0];
Nfft =128;
fs =1e6;
ht1 = zeros(128,1);
ht2=zeros(128,1);
ht3=zeros(128,1);
ht4=zeros(128,1);
for trial=1:100,
ht1(1+Tb1)=10.^(Pwd1/20).*exp(sqrt(-1)*randn(1,length(Tb1)));
ht2(1+Tb2)=10.^(Pwd2/20).*exp(sqrt(-1)*randn(1,length(Tb2)));
ht3(1+Tb3)=10.^(Pwd3/20).*exp(sqrt(-1)*randn(1,length(Tb3)));
ht4(1+Tb4)=10.^(Pwd4/20).*exp(sqrt(-1)*randn(1,length(Tb4)));
%hh=h(1)+[zeros(1,2),h(2)(1:end-2)]+[zeros(1,5),h(3)(1:end-5)];
Ft = ifft(eye(Nfft))*sqrt(Nfft);
F = ifftshift(Ft,2);
hf1= F'*ht1;
hf2= F'*ht2;
hf3= F'*ht3;
hf4= F'*ht4;
sum_hf=hf1+hf2+hf3+hf4;

shift_ht1=ht1;
shift_ht2=[ht2(128-32:128);ht2(1:(128-32-1))];
shift_ht3=[ht3(128-64:128);ht3(1:(128-64-1))];
shift_ht4=[ht4(128-96:128);ht4(1:(128-96-1))];
shift_hf1=F'*shift_ht1;
shift_hf2=F'*shift_ht2;
shift_hf3=F'*shift_ht3;
shift_hf4=F'*shift_ht4;
sum_shift_hf=shift_hf1+shift_hf2+shift_hf3+shift_hf4;
% figure(trial);
% plot(abs(sum_hf));
% hold on;
% plot(abs(sum_shift_hf),'r');
% hold off;
% pause;
no_ration((trial-1)*128+1:trial*128)=sum_hf;
ration((trial-1)*128+1:trial*128)=sum_shift_hf;
end;
toc;