【雷达通信】雷达脉冲压缩含Matlab源码

function varargout = test1(varargin)%TEST1 M-file for test1.fig%      TEST1, by itself, creates a new TEST1 or raises the existing%      singleton*.%%      H = TEST1 returns the handle to a new TEST1 or the handle to%      the existing singleton*.%%      TEST1('Property','Value',...) creates a new TEST1 using the%      given property value pairs. Unrecognized properties are passed via%      varargin to test1_OpeningFcn.  This calling syntax produces a%      warning when there is an existing singleton*.%%      TEST1('CALLBACK') and TEST1('CALLBACK',hObject,...) call the%      local function named CALLBACK in TEST1.M with the given input%      arguments.%%      *See GUI Options on GUIDE's Tools menu.  Choose "GUI allows only one%      instance to run (singleton)".%% See also: GUIDE, GUIDATA, GUIHANDLES% Edit the above text to modify the response to help test1% Last Modified by GUIDE v2.5 09-Mar-2009 14:24:57% Begin initialization code - DO NOT EDITgui_Singleton = 1;gui_State = struct('gui_Name',       mfilename, ...                   'gui_Singleton',  gui_Singleton, ...                   'gui_OpeningFcn', @test1_OpeningFcn, ...                   'gui_OutputFcn',  @test1_OutputFcn, ...                   'gui_LayoutFcn',  [], ...                   'gui_Callback',   []);if nargin && ischar(varargin{1})   gui_State.gui_Callback = str2func(varargin{1});endif nargout    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});else    gui_mainfcn(gui_State, varargin{:});end% End initialization code - DO NOT EDIT% --- Executes just before test1 is made visible.function test1_OpeningFcn(hObject, eventdata, handles, varargin)% This function has no output args, see OutputFcn.% hObject    handle to figure% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% varargin   unrecognized PropertyName/PropertyValue pairs from the%            command line (see VARARGIN)% Choose default command line output for test1handles.output = hObject;% Update handles structureguidata(hObject, handles);% UIWAIT makes test1 wait for user response (see UIRESUME)% uiwait(handles.figure1);% --- Outputs from this function are returned to the command line.function varargout = test1_OutputFcn(hObject, eventdata, handles)% varargout  cell array for returning output args (see VARARGOUT);% hObject    handle to figure% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)% Get default command line output from handles structurevarargout{1} = handles.output;% --- Executes on button press in pushbutton1.function pushbutton1_Callback(hObject, eventdata, handles)% hObject    handle to pushbutton1 (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)creat=findobj(gcf,'Tag','radiobutton1');processing=findobj(gcf,'Tag','radiobutton2');global TT=str2double(get(handles.edit1,'string'));global BB=str2double(get(handles.edit2,'string'));if(get(creat,'Value')==get(creat,'Max'))K=B/T;                                    %频率调制斜率Fs=5*B;Ts=1/Fs;                          %计算机仿真的采样频率和采样周期N=T/Ts;                                   %采样点数t=linspace(-T/2,T/2,N);St=exp(j*pi*K*t.^2);                         %产生线性调频信号axes(handles.axes1)handles.axes1_handle=gca;plot(t*1e6,real(St));xlabel('\fontsize{9}时间(us)');title('\fontsize{9}LFM脉冲的时域波形');grid on;axis tight;freq=linspace(-Fs/2,Fs/2,N);% freq=linspace(0,Fs,N);axes(handles.axes2)handles.axes2_handle=gca;% plot(freq*1e-6,abs(fft(St)));plot(freq*1e-6,fftshift(abs(fft(St))));xlabel('\fontsize{9}频率(MHz)');title('\fontsize{9}LFM脉冲的幅频特性');grid on;axis tight;elseif(get(processing,'Value')==get(processing,'Max'))        global RminRmin=str2double(get(handles.edit3,'string'));global RmaxRmax=str2double(get(handles.edit4,'string'));global RR=str2num(get(handles.edit5,'string'));global RCSRCS=str2num(get(handles.edit6,'string'));global winidwinid=str2double(get(handles.edit7,'string'));% T=10e-6; %pulse duration 10us% B=30e6; %chirp frequency modulation bandwidth 30MHz% Rmin=10000;Rmax=15000; %range bin% R=[10500,11000,12000,12008,13000,13002]; %position of ideal point targets% RCS=[1 1 1 1 1 1]; %radar cross section% winid=2;%=========================================================%%ParameterC=3e8; %propagation speedK=B/T; %chirp slopeRwid=Rmax-Rmin; %receive window in meterTwid=2*Rwid/C; %receive window in secondFs=5*B;Ts=1/Fs; %sampling frequency and sampling spacingNwid=ceil(Twid/Ts); %receive window in number%=================================================================%determine proper windowif(winid==0.)    win(1:Nwid)=1.;endif(winid==1.)    win=hamming(Nwid)';endif(winid==2.)    win=chebwin(Nwid,60)';endnosewite=findobj(gcf,'Tag','checkbox1');    if(get(nosewite,'Value')==get(nosewite,'Max'))         nose=randn(1,Nwid);      else        nose=zeros(1,Nwid);        end%==================================================================%%Gnerate the echo t=linspace(2*Rmin/C,2*Rmax/C,Nwid); %receive window%open window when t=2*Rmin/C%close window when t=2*Rmax/C M=length(R); %number of targets td=ones(M,1)*t-2*R'/C*ones(1,Nwid);Srt=RCS*(exp(j*pi*K*td.^2).*(abs(td)<T/2))+nose;%radar echo from point targets Srt=Srt.*win;%=========================================================%%Digtal processing of pulse compression radar using FFT and IFFTNchirp=ceil(T/Ts); %pulse duration in numberNfft=2^nextpow2(Nwid+Nchirp-1); %number needed to compute linear %convolution using FFT algorithmSrw=fft(Srt,Nfft); %fft of radar echot0=linspace(0,T,Nchirp); St=exp(j*pi*K*t0.^2); %chirp signal Sw=fft(St,Nfft); %fft of chirp signalSot=ifft(Srw.*conj(Sw)); %signal after pulse compression%=========================================================N0=1;Z=abs(Sot(N0:N0+Nwid-1));Z=Z/max(Z);Z=20*log10(Z+1e-6);%figure% subplot(211)axes(handles.axes1)handles.axes1_handle=gca;plot(t*1e6,real(Srt));axis tight;xlabel('Time in u sec');ylabel('Amplitude')title('Radar echo without compression');% subplot(212)axes(handles.axes2)handles.axes2_handle=gca;plot(t*C/2,Z)%  axis([11900,12100,-60,0]);axis([Rmin,Rmax,-60,0]);xlabel('Range in meters');ylabel('Amplitude in dB')title('Radar echo after compression');end% --- Executes on button press in pushbutton2.function pushbutton2_Callback(hObject, eventdata, handles)% hObject    handle to pushbutton2 (see GCBO)% eventdata  reserved - to be defined in a future version of MATLAB% handles    structure with handles and user data (see GUIDATA)close