亲测可用
- %%线性调频信号经过匹配滤波
- T=10e-6; %脉冲持续时间10us
- B=30e6; %线性调频信号的带宽30MHz
- K=B/T; %调频斜率
- Fs=10*B;Ts=1/Fs; %采样频率和采样间隔
- N=T/Ts;
- t=linspace(-T/2,T/2,N);
- St=exp(j*pi*K*t.^2); %线性调频信号
- Ht=exp(-j*pi*K*t.^2); %匹配滤波
- Sot=conv(St,Ht); %线性调频信号经过匹配滤波
- subplot(411)
- L=2*N-1;
- t1=linspace(-T,T,L);
- Z=abs(Sot);
- Z=Z/max(Z); %归一化
- Z=20*log10(Z+1e-6);
- Z1=abs(sinc(B.*t1)); %辛格函数
- Z1=20*log10(Z1+1e-6);
- t1=t1*B;
- plot(t1,Z,t1,Z1,'r.');
- axis([-15,15,-50,inf]);grid on;
- legend('匹配滤波','辛格');
- xlabel('时间/s ');
- ylabel('幅值/dB');
- title('线性调频信号经过匹配滤波');
- subplot(412) %展开
- N0=3*Fs/B;
- t2=-N0*Ts:Ts:N0*Ts;
- t2=B*t2;
- plot(t2,Z(N-N0:N+N0),t2,Z1(N-N0:N+N0),'r.');
- axis([-inf,inf,-50,inf]);grid on;
- set(gca,'Ytick',[-13.4,-4,0],'Xtick',[-3,-2,-1,-0.5,0,0.5,1,2,3]);
- xlabel('时间/s');
- ylabel('幅度/dB');
- title('线性调频信号经过匹配滤波（补零展开之后）');
- Z_value = Z(2950:3050)
- subplot(413)
- plot(Z_value);
- grid on;
- [value, locmax ] = max( Z_value(1:40));
- locmax
- value
- subplot(414)
- plot(Z);
- sum = 0;
- for mm = 1:length(Z)
- sum = sum +Z(mm)
-
- end
- sum = -sum;
- sum
|
QQ好友和群
QQ空间
腾讯微博
腾讯朋友
收藏
淘帖
顶
踩
