clear;
clc;
close?all;
fid=fopen（‘filter.txt‘‘r‘）;
[fir_coffr_data]=fscanf（fid‘%f‘[1256]）;


f0=30e6;???%中心頻率
T=1e-6;????%線性調(diào)頻信號(hào)脈寬
B=16e6;?????%線性調(diào)頻信號(hào)帶寬
K=B/T;??????%調(diào)頻斜率
Fs=1024e6;???%采樣率
Ts=1/Fs;?
D=16;
N=1024;
t=linspace（0TN）;
st=cos（2*pi*f0*t+pi*K*t.^2）;?%線性調(diào)頻信號(hào)（實(shí)信號(hào)）
%??st=exp（j*（2*pi*f0*t+pi*K*t.^2））;%線性調(diào)頻信號(hào)（復(fù)信號(hào)）

figure（1）
freq=linspace（0FsN）;
signal_am_freq=abs（fft（st））;
plot（tsignal_am_freq）;
xlabel（‘頻率?Hz‘）;
ylabel（‘幅度‘）;

?%多相濾波器中抽取倍數(shù)為16倍
S=reshape（st1664）???%將總的信號(hào)分解成一個(gè)矩陣
S3?=?flipud（S）;?
for?m=1:64
??S（:m）=?S3（:m）.*exp（-1*j*（pi/2）*（m））;
%?S（:m）=?S3（:m）.*（-1）^m;

end

?
%?%==========================================================================
%?%多相濾波器通道
%?%==========================================================================
?n=256;??????%濾波器系數(shù)個(gè)數(shù)
for?p?=?1:D
????for?m?=?1:n/D
????????h（pm）?=?fir_coff（（m-1）*D+p）;?
????end
end
for?K=1:D
????S4（K:）=filter（h（K:）1S（K:））;
end
for?K=1:D
????S1（K:）=?S4（K:）.*exp（1*j*（pi/（2*D）*（K-1）））;
end



%?以下是對(duì)每一列進(jìn)行FFT轉(zhuǎn)換，可得

[gh]?=?size（S1）;
for?i=1:h
????????output（:i）=D*ifft（S1（:i）16）;
end
?
?[ymyn]?=?size（output）;
?fsy?=?Fs/D;??%抽取后?采樣點(diǎn)數(shù)的變化
?fy?=?linspace（-fsy/2fsy/2yn）;
?
?figure（2）
?%Y1=real（output（16:））
?ys（1:）?=?abs（fftshift（fft（output（1:））））;
?
????subplot（121）;
????plot（fyys（1:））;
%????ylim（[0?6e5]）;
????xlabel（‘頻率?Hz‘）;
????ylabel（‘幅度‘）;
??axis（[-6.4e7?6.4e7?0?12]）;
????grid?on;
????
???
????ys（2:）?=?abs（fftshift（fft（output（2:））））;
????
????subplot（122）;
????plot（fyys（2:））;
%????ylim（[0?6e5]）;
????xlabel（‘頻率?Hz‘）;
????ylabel（‘幅度‘）;
?axis（[-6.4e7?6.4e7?0?12]）;
????grid?on;
????
????figure（3）
for?c=1:16
????ys（c:）?=?abs（fftshift（fft（output（c:））））;
????subplot（116c）;
????plot（fyys（c:））;
%????ylim（[0?6e5]）;
????xlabel（‘頻率?Hz‘）;
????ylabel（‘幅度‘）;
????axis（[-6.4e7?6.4e7?0?20]）;
????grid?on;
end

%?
%?
%?M=output（1:）

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件???????5698??2015-04-07?18:38??filter.txt

?????文件???????2083??2018-06-07?17:57??channel.m

-----------?---------??----------?-----??----

?????????????????7781????????????????????2