雷达信号检测matlab,calplotsig7.m 源代码在线查看 - matlab中实现雷达信号处理.rar 资源下载 虫虫电子下载站...

% calplotsig7.m - Written by Eli Mozeson and Nadav Levanon% used by ambfn7 for calculation and plot of the signal when the signal is defined% by u_amp, u_phase and u_freq (single carrier signal)% output variables include u(t) amd t on which u is definedtb=1; if (acode==1)*(pcode==1), u_basic=u_amp.*exp(j*u_phase*pi);elseif (acode==1)*(pcode==0), u_basic=u_amp;elseif (acode==0)*(pcode==1), u_basic=exp(j*u_phase*pi);else u_basic=ones(size(u_amp));end m_basic=length(u_basic);if r==1 dt=tb; m=m_basic; uamp=abs(u_basic); phas=uamp*0; phas=angle(u_basic); if fcode==1 phas=phas+2*pi*dt*cumsum(f_basic); end uexp=exp(j*phas); u=uamp.*uexp;else % i.e., several samples within a bit dt=tb/r; % interval between samples ud=diag(u_basic); ao=ones(r,m_basic); m=m_basic*r; u_basic=reshape(ao*ud,1,m); % u_basic with each eleement repeated r times uamp=abs(u_basic); phas=angle(u_basic); u=u_basic; if fcode==1 ff=diag(f_basic); phas=2*pi*dt*cumsum(reshape(ao*ff,1,m))+phas; uexp=exp(j*phas); u=uamp.*uexp; endendtscale=[0:length(uamp)-1]/r;tscale1=[0 0:length(uamp)-1 length(uamp)-1]/r;dphas=[NaN diff(phas)]*r/2/pi;figure(sigfig), clf, hold off % plot the signal parameterssubplot(3,1,1)plot(tscale1,[0 abs(uamp) 0],'k','linewidth',1.5)ylabel(' Amplitude ')titlest=presetvalues(get(preset,'value')).Name;title(titlest);axis([-inf inf 0 1.2*max(abs(uamp))])subplot(3,1,2)plot(tscale, phas,'k','linewidth',1.5)axis([-inf inf -inf inf])ylabel(' phase [rad] ')subplot(3,1,3)plot(tscale,dphas*ceil(max(tscale)),'k','linewidth',1)axis([-inf inf -inf inf])xlabel(' \itt / t_b ')ylabel(' \itf * Mt_b ')% variables for ambigity calculationst=tscale;u=u;