ambfn7.m

matlab中实现雷达信号处理。 matlab中实现雷达信号处理。

% "ambfn7.m" - builds an MMI and calls dedicated m-files to calcualte the signal 
% complex envelope, plot signals phase, amplitude and frequency structure, and
% ambiguity function and periodic ampbiguity function
%
% Written by Nadav Levanon and Eli Mozeson, Dept. of EE-Systems, Tel Aviv University
%
% ambfn7.m includes a "User defined" line to the preset values list. The line  
% is activated each time the user changes one of the signal parameters
% ambfn7.m also allows the user to save/load the signal and plot parameters 
%
clear all
close all

% create the parameters input figure
inputfig=figure;
set(inputfig,'Position',[6 51 520 690],'MenuBar','none','visible','off');

% create results figure (ambiguity function plot)
ambfig=figure;
set(ambfig,'Position',[274 52 749 670],'Name','ambiguity function plot','visible','off');

% create results figure (signal parameters plot)
sigfig=figure;
set(sigfig,'Position',[274 52 749 670],'Name','signal parameters plot','visible','off');

% create results figure (Autocorelation and Spectrum plot)
acffig=figure;
set(acffig,'Position',[274 52 749 670],'Name','Autocorrelation and Spectrum plot','visible','off');

% create results figure (Periodic ambiguity plot)
pambfig=figure;
set(pambfig,'Position',[274 52 749 670],'Name','Periodic ambiguity plot','visible','off');

% set default values for the signal and plot - single carrier parameters
acode=1;										% signal amplitude modulation flag
pcode=1;										% signal phase modultaion flag
u_amp=ones(1,13);							% signal amplitude vector
u_phase=zeros(1,13);						% signal phase vector
u_basic=u_amp.*exp(j*u_phase*pi);	% signal complex envelope (no frequency modulation)

fcode=1;										% signal frequency modulation flag
f_basic=zeros(1,13);						% signal frequncy modulation vector

% set default values for the signal and plot - ambiguity grid parameters
F=5;											% maximal displayed Doppler is F/Mtb
r=10;											% signal is sampled r times in each tb (r/tb sampling rate)	
N=100;										% number of points in the delay axis for ambiguity plot is 2N
K=50;df=F/K/length(u_amp);				% number of points in the Doppler axis for ambiguity plot is K

T=1;											% maximal delay displayed for the ambiguity plot is T*Mtb

% initilaize flag of signal calculation
sigflag=0;									% signal was not calculated yet

% initilaize strings for title of the plots
titlest='';

% create amplitude modulation input row
y_u_amp=530;

u_amp_legend = uicontrol(inputfig,'Style','text',...
   'String','Amplitude','Pos',[25 y_u_amp 100 20 ]);

u_amp_text = uicontrol(inputfig,'Style','edit',...
   'String',num2str(u_basic),'Position',[125 y_u_amp 300 20],...
   'Callback',[...
      'sigflag=0;'...
      'set(preset,''value'',12);',...
      'u_amp=str2num(get(u_amp_text,''String''));'...
      'set(u_amp_text,''String'',num2str(u_amp));']);

u_amp_on = uicontrol(inputfig,'Style','radio',...
   'Pos',[425 y_u_amp 20 20 ],'Value',1,'Callback',[...
      'sigflag=0;'...
      'set(preset,''value'',12);',...
      'set(u_amp_on,''value'',get(u_amp_on,''value'')),'...
      'acode=get(u_amp_on,''value'');']);

% create phase modulation input row
y_u_phase=505;

u_phase_legend = uicontrol(inputfig,'Style','text',...
   'String','Phase/pi','Pos',[25 y_u_phase 100 20 ]);

u_phase_text = uicontrol(inputfig,'Style','edit',...
   'String',num2str(u_phase),'Position',[125 y_u_phase 300 20],...
   'Callback',[...
      'sigflag=0;'...
      'set(preset,''value'',12);',...
   	'u_phase=str2num(get(u_phase_text,''String''));'...
      'set(u_phase_text,''String'',num2str(u_phase));']);

u_phase_on = uicontrol(inputfig,'Style','radio',...
   'Pos',[ 425 y_u_phase 20 20 ],'Value',1,'Callback',[...
      'sigflag=0;'...
      'set(preset,''value'',12);',...
   	'set(u_phase_on,''value'',get(u_phase_on,''value'')),'...
      'pcode=get(u_phase_on,''value'');']);

% create frequency modulation input row
y_freq=480;

freq_legend = uicontrol(inputfig,'Style','text',...
   'String','Frequency*tb','Pos',[25 y_freq 100 20 ]);

freq_text = uicontrol(inputfig,'Style','edit',...
   'String',num2str(f_basic),'Position',[125 y_freq 300 20],...
   'Callback',[...
      'sigflag=0;'...
      'set(preset,''value'',12);',...
      'f_basic=str2num(get(freq_text,''String''));'...
      'set(freq_text,''String'',num2str(f_basic));']);

freq_on = uicontrol(inputfig,'Style','radio',...
   'Pos',[ 425 y_freq 20 20 ],'Value',1,'Callback',[...
      'sigflag=0;'...
      'set(preset,''value'',12);',...
      'set(freq_on,''value'',get(freq_on,''value'')),'...
      'fcode=get(freq_on,''value'');']);

% create F input slider
y_F = 375;

sli_F = uicontrol(inputfig,'Style','Slider','sliderstep',[0.01 0.1],...
   'Position',[125 y_F 300 20],'Min',0,'Max',60,'Value',F,'Callback',[...
      'set(F_cur,''String'',',...
      'num2str(get(sli_F,''Val''))),',...
      'F=get(sli_F,''Val'');',...
      'df=F/K/length(u_amp);']);

F_cur = uicontrol(inputfig,'style','edit','Pos',[ 425 y_F 50 20 ],...
   'String',num2str(get(sli_F,'Value')),'Callback',[...
      'set(sli_F,''Val'',',...
      'str2num(get(F_cur,''String''))),',...
      'F=get(sli_F,''Val'');',...
      'df=F/K/length(u_amp);']);

F_lbl = uicontrol(inputfig,'style','Text',...
   'Pos',[ 25 y_F 100 20 ],'String','F*Mtb');

% create r input slider
y_r = 400;

sli_r = uicontrol(inputfig,'Style','Slider','Position',[125 y_r 300 20],...
   'Min',1,'Max',300,'Value',r,'Callback',[...
      'sigflag=0;'...
      'set(r_cur,''String'',num2str(floor(get(sli_r,''Val'')))),',...
      'r=floor(get(sli_r,''Val''));']);

r_cur = uicontrol(inputfig,'style','edit','Pos',[425 y_r 50 20],...
   'String',num2str(floor(get(sli_r,'Value'))),'Callback',[...
      'sigflag=0;'...
      'set(sli_r,''Val'',str2num(get(r_cur,''String''))),',...
      'r=floor(get(sli_r,''Val''));']);

r_lbl = uicontrol(inputfig,'style','Text',...
   'Pos',[25 y_r 100 20],'String','r');

% create N input slider
y_N = 325;

sli_N = uicontrol(inputfig,'Style','Slider','Position',[125 y_N 300 20],...
   'Min',1,'Max',300,'Value',N,'Callback',[...
      'set(N_cur,''String'',',...
      'num2str(floor(get(sli_N,''Val'')))),',...
      'N=floor(get(sli_N,''Val''));']);

N_cur = uicontrol(inputfig,'style','edit',...
   'Pos',[425 y_N 50 20],'String',num2str(floor(get(sli_N,'Value'))),'Callback',[...
      'set(sli_N,''Val'',',...
      'str2num(get(N_cur,''String''))),',...
      'N=floor(get(sli_N,''Val''));']);

N_lbl = uicontrol(inputfig,'style','Text',...
   'Pos',[25 y_N 100 20],'String','N');

% create T input slider
y_T = 350;

sli_T = uicontrol(inputfig,'Style','Slider','Position',[125 y_T 300 20],...
   'Min',0.01,'Max',1.1,'Value',T,'sliderstep',[0.01 0.1],'Callback',[...
      'set(T_cur,''String'',',...
      'num2str(get(sli_T,''Val''))),',...
      'T=get(sli_T,''Val'');']);

T_cur = uicontrol(inputfig,'style','edit',...
   'Pos',[425 y_T 50 20],'String',num2str(get(sli_T,'Value')),'Callback',[...
      'set(sli_T,''Val'',',...
      'str2num(get(T_cur,''String''))),',...
      'T=get(sli_T,''Val'');']);