%%时域FIR宽带波束形成
%%
clear all;clc;close all；
c=1500;
d=0.375;
N=32;
sam_fre=5000;
bandwidth=[1900 2100];
f_step=10;
FilterStep=64;
signalenth=16384;
%%
hwait=waitbar(0,'please wait');
for i1=1:1:length(ScanRange);
k=i1/lenth(ScanRange);
FilterParameters_1=zeros(N,FilterStep);
TimeDelay_1=(0:N1-1)*d*cosd(ScanRange(i1))/c;
FilterSampDelays=-fix(TimeDelay_1*sam_fre);
TimeDelay_2=TimeDelay_1+FilterSampDelays*(1/sam_fre);
for i2=1:1:N
dn=zeros(signalenth,1);
xn=zeros(signalenth,1);
for i3=1:1:signalenth
for i4=1:1:(sam_fre/2/f_step+1)
xn(i3)=xn(i3)+cos(2*pi*(i4-1)*f_step)*(i3/sam_fre);
if ((i4-1)*f_step)>=min(bandwidth) && (i4-1)*f_step<=max(bandwidth)
dn(i3)=dn(i3)+cos(2*pi*((i4-1)*f_step)*(i3/sam_fre-TimeDelay_2(i2))+2*pi*((i4-1)*f_step)/sam_fre*(FilterStep-1)/2);
end
end
end
Rxx=zeros(FilterStep,FilterStep);
Rxd=zeros(FilterStep,1);
for i3=1:1:(signalenth-FilterStep+1)
Rxx=Rxx+fliplr(xn(i3:(i3+FilterStep-1)))*fliplr(xn(i3:(i3+FilterStep-1)))'/(signalenth-FilterStep+1);
Rxd=Rxd+fliplr(xn(i3:(i3+FilterStep-1)))*dn(i3)/(signalenth-FilterStep+1);
end
FilterParameters_1(i2,:)=(pinv(Rxx)*Rxd)';
end
FilterParameters=[FilterParameters_1 FilterSampDelays'];
str=['正在运行',num2str(100*k),'%'];
waitbar(k,hwait,str);
pause(0.05);
save(['' num2str(ScanRange(i1)) '.mat'],'FilterParameters','-v7.3');
clear FilterParameters;
clear FilterParameters_1;
end
close(hwait);
评论6