【信号处理】基于Matlab实现POSP驻定相位原理推导线性调频信号的频谱

1 简介

2 完整代码

%% Principle of Srationary Phaseclose all, clear all, clc%% the parameters of radarT = 10e-6;    Tsyn = T;    %脉冲持续时间K = 1e12;B = K*T;    %线性调频信号带宽% K = B/Tsyn ;     %调频斜率a = 1.3;Fs = a*B;     %抽样频率Ts = 1/Fs;    %样本间隔TBP = B*Tsyn;N = T/Ts;%% time domaint = [ -N/2 : N/2-1 ]*T/N ;sita = pi*K*t.^2;St = exp( i*sita ).*( abs( t ) <= Tsyn/2 );    %LFM signal%% frequency domainfre = [ -N/2 : N/2-1 ]*Fs/N ;% -----FFT-----% Sf = fft(St);Sf = (fft(fftshift(St)));% -----PSP-----Sf2 = sqrt(Tsyn*Fs*a)*exp(-j*pi*fre.^2/K+j*pi/4).*(abs(fre)<=(K*Tsyn/2));%% figureh1 = figure( 1 ); set( h1, 'name', 'Linear Frequency Modulation Signal', 'Numbertitle', 'off')subplot(121),plot( fre*1e-6, fftshift(abs( Sf ) ), 'r'); hold onplot(fre*1e-6, abs( Sf2 ))legend('FFT', 'PSP');xlabel( 'Frequency/MHz ' ); ylabel( 'Amplitude' ); title( 'Amplitude-Frequency Specttrum' ); grid on; %generate the figure 1.4subplot(122),plot( fre*1e-6, unwrap( angle( fftshift(Sf) ) ) ,'r'), hold on% plot( fre*1e-6, phase( Sf ), 'g'), hold onplot( fre*1e-6, unwrap( angle( (Sf2) ) ));% plot( fre*1e-6, phase(Sf2) ,'r');legend('FFT', 'PSP');xlabel( 'Frequency/MHz' ), ylabel( 'Radian' ), title( 'Radian-Frequency Specttrum' ), grid on;  figureSf = fftshift(Sf);sf = Sf(N/2*(1-1/a)+1 : N/2*(1+1/a));sf2= Sf2(N/2*(1-1/a)+1 : N/2*(1+1/a));fre2 = [N/2*(1-1/a)+1 : N/2*(1+1/a)]*Fs/N-Fs/2;plot( fre2, unwrap( angle( (sf) ) ), 'r.' ), hold onerror = unwrap( angle( (sf) ) ) - unwrap( angle( (sf2) ) );plotyy( fre2, unwrap( angle( (sf2) ) ), fre2, error, 'plot', 'plot'), hold onlegend('FFT', 'PSP', 'ERROR');xlabel( 'Frequency/MHz' ), ylabel( 'Radian' ), title( 'Radian' ), grid on;  %% display informationdisp('Chirp rate:')Kdisp('Time-Bandwidth Product：' )TBPdisp('the mean of error:')mean(error)disp('the variance of error:')var(error)

3 仿真结果

4 参考文献

[1]唐亮. 频率步进SAR/ISAR成像算法研究[D]. 国防科学技术大学, 2014.

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

部分理论引用网络文献，若有侵权联系博主删除。

【信号处理】基于Matlab实现POSP驻定相位原理推导线性调频信号的频谱相关推荐

驻定相位原理（POSP）以及线性调频信号的频谱
最近看论文遇到了驻定相位原理,问老师直接给了我一本书让我看,看半天只有一段-不是这个方向的,半路出家做毕业设计需要用到这个定理,有错误的话请不吝赐教. 一.驻定相位原理在数字信号处理中,经常需要将一 ...
驻定相位原理（POSP）的简单应用
在SAR雷达成像中,POSP是相当基础重要的一个定理,一般在对回波做傅里叶变换时经常用到,一般在论文的开头就会出现. 下面简单复习一下POSP的步骤: 1:列出傅里叶变换表达式 2:对相位在驻 ...
驻定相位原理在求解频域表达式中的应用
通过前面文章的学习,我们知道了距离弯曲量在时域的表达式,但是做距离弯曲校正时,要在多普勒频域完成.那么在多普勒频域下距离弯曲量如何计算呢?那首先是不是得知道回波信号的频域表达式是什么呢? 这就涉及到了 ...
matlab的dft谱分析,数字信号处理基于matlab(用DFT作谱分析,窗函数的设计)
数字信号处理基于matlab(用DFT作谱分析,窗函数的设计) 1实验一用DFT作谱分析X11111X212344321N108X3COSN1PI/4N208X4SINN2PI/8FIGURESUBP ...
NPC_3level_Inverter：基于MATLAB Simulink的中性点钳位三电平逆变器仿真模型
NPC_3level_Inverter:基于MATLAB Simulink的中性点钳位三电平逆变器仿真模型. 仿真条件:MATLAB Simulink R2015b ID:54506503695843 ...
基于MATLAB Simulink的中性点钳位三电平逆变器仿真模型
NPC_3level_Inverter:基于MATLAB Simulink的中性点钳位三电平逆变器仿真模型. 仿真条件:MATLAB Simulink R2015b ID:54506503695843 ...
基于MATLAB Simulink的中性点钳位五电平逆变器仿真模型
NPC_5level_Inverter:基于MATLAB Simulink的中性点钳位五电平逆变器仿真模型. 仿真条件:MATLAB Simulink R2015b ID:32806503698006 ...
基于matlab的有噪声语音信号处理,基于matlab的有噪声语音信号处理毕设
基于matlab的有噪声语音信号处理毕设 - I - I 大学本科毕业设计论文基于MATLAB的有噪声语音信号处理- II - II 摘要滤波器设计在数字信号处理中占有极其重要的地位,FIR 数字 ...
matlab 音乐信号处理,基于matlab的音乐信号处理和分析.doc
基于matlab的音乐信号处理和分析 1音乐信号的音谱和频谱的观察使用windows下的录音机录制一段音乐信号或者采用其他软件截取一段音乐信号(要求:时间不超过5s,文件格式为WAV) 使用wavr ...

【信号处理】基于Matlab实现POSP驻定相位原理推导线性调频信号的频谱

1 简介

2 完整代码

3 仿真结果

4 参考文献

博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。

【信号处理】基于Matlab实现POSP驻定相位原理推导线性调频信号的频谱相关推荐

最新文章

热门文章