comprcv.m

来自「此程序给出了不同情况的雷达距离测量仿真平台」· M 代码 · 共 59 行

% COMPRESSIVE RECEIVER USING FFT
% ----------------------------------

clear;clc;clf;

% Input Chirp Bandwidth and Pulsewidth

bw=1e6;                          % Receiver bandwidth - Hz
T= 100e-6;                       % Analysis window - secs
fo=1e5;                          % Input signal #1 frequency - Hz
f1=1.165e5;                      % Input signal #2 frequency - Hz

% Input FFT Points

lfft=1024;                       % FFT points

% Set Time Base and Compute Parameters

smp= 1/(2*bw);
pts= bw*T;
t=-T/2:T/(lfft-1):T/2;

% Sample Complex Chirp

s= exp(-j*pi*bw*t.*t/T);

% Sample Input Signals

vt=exp(-j*2*pi*fo*t)+exp(-j*2*pi*f1*t);
svt=s.*vt;SV=fft(svt,lfft);

%Find Fourier Spectrum of Chirp

Y= fft(s,lfft);

% Find Power Spectral Density

Psd=SV.*conj(Y)/lfft;
A=ceil(max(abs(Psd)));
l=length(Y);
Psd1=fftshift(Psd);
f=(1/T)*(-(l-1)/2:(l-1)/2);

% Find Matched Filter Time Waveform

zh=ifft(Psd,lfft);
h=fftshift(zh);
h=h/max(h);

% Plot Analysis Window- Waveform

B=max(abs(h));

plot(t,abs(h));grid
axis([-(T/2) 0 0 B])
title(['Compressive Receiver Output']);
xlabel('Time - sec');ylabel('Amplitude');