very slow any other way to do it?

2 次查看（过去 30 天）

sharif 2014-7-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/139984-very-slow-any-other-way-to-do-it

编辑： Cedric 2014-7-3

 clear, clc, close all;
N=1024; Cp=64; zf=[]; NTx=8; Npath=100; tap=randi([3 TAP]); ds=Cp; s=N+N-1; N1=Npath*Npath; 
N2=Npath*NTx;
for aa=1:Npath
    for bb=1:tap
        p(aa,bb)=randi([1 ds]);
        p(aa,:)=sort(p(aa,:),'descend');
        xx(aa,:)=(abs(randn(1,tap)+(randn(1,tap))*1i));
        xx(aa,:)=sort(xx(aa,:));
    end
end
% Loop for adding random positions & normalized impulse response (h_N)
for cc=1:Npath
    for dd=1:tap
        h(cc,p(cc,dd))=xx(cc,dd);
        U(cc)=sqrt(sum((h(cc,:).^2)));
        h_N(cc,:)=h(cc,:)./U(cc);
        h_N(cc,:)=circshift(h_N(cc,:),[0 -min(p(cc,:))+1]);
    end
end
bk= randi ([0,15],NTx,N);
Dk=C(bk+1);
d=ifft((Dk).*sqrt(N));
tx=[d(:,end-Cp+1:end) d];
* | _for ii=1:Npath
    for hh=1:N+Cp
        for jj=1:NTx           
            [Txx(jj,:),zf]=filter(h_N(ii,:),1,tx(jj,:),zf);
           Txxx(ii,hh)=sum(Txx(:,hh));
        end
    end    
end_ | *