線性調頻信號的數學分析

• 代碼：

%%-----------------------------------------------------------------------------------------------------------%% %%說明：調用此函數可以輸出線性調頻信號余弦表達式下的信號波形圖及頻譜圖；以及復數表達式下的信號波形圖的實部、虛部及頻譜圖%% %%-----------------------------------------------------------------------------------------------------------%% function [st1,st2] = LFM_signal(A,Phi0,T,B,F0) %* st1:線性調頻信號的復數表達式 **% %* st2:線性調頻信號的余弦表達式 **% %*** A:信號的振幅 ***************% % Phi0:信號的隨機初相 ***********% %*** T:信號時寬 *****************% %*** B:信號帶寬 *****************% %** F0:信號的中頻頻率，即載頻頻率 %%%%%% 信號的參數設置 %%%%% K=B/T; %調頻斜率 Fs=2*B; %采樣頻率 Ts=1/Fs; %采樣周期 N=T/Ts; %采樣點數%%%%% 線性調頻信號的兩種表達方式 %%%%% t=linspace(-T/2,T/2,N); st1=A*exp(1j*(2*pi*F0*t+pi*K*t.^2+Phi0)); %線性調頻信號的復數表達式 st2=A*cos(2*pi*F0*t+pi*K*t.^2+Phi0); %線性調頻信號的余弦表達式figure(1); subplot(3,1,1); plot(t*1e6,real(st1)); xlabel('時間/us'); ylabel('實部') title('線性調頻信號的實部'); grid on; axis tight;subplot(3,1,2); plot(t*1e6,imag(st1)); xlabel('時間/us'); ylabel('虛部') title('線性調頻信號的虛部'); grid on; axis tight;subplot(3,1,3); freq=linspace(-Fs/2,Fs/2,N); plot(freq*1e-6,fftshift(abs(fft(st1)))); %先對st做傅里葉變換得到頻譜，然后取幅度值，再將其移動到頻譜中心 xlabel('頻率/MHz'); ylabel('幅度譜') title('線性調頻信號的頻譜'); grid on; axis tight;figure(2); subplot(2,1,1); plot(t*1e6,real(st2)); xlabel('時間/us'); ylabel('實部') title('線性調頻信號的實部'); grid on; axis tight;subplot(2,1,2); freq=linspace(-Fs/2,Fs/2,N); plot(freq*1e-6,fftshift(abs(fft(st2)))); %先對st做傅里葉變換得到頻譜，然后取幅度值，再將其移動到頻譜中心 xlabel('頻率/MHz'); ylabel('幅度譜') title('線性調頻信號的頻譜'); grid on; axis tight; end

• 結果：

當信號的中心頻率為$F_0=0$時

close all; clear all; clc; %%調用LFM_signal函數，觀察結果 A=1; %發射信號的振幅 Phi0=0; %發射信號的隨機初相 T=10e-6; %信號時寬 B=30e6; %信號帶寬 F0=0; %中頻頻率，即載頻頻率 [st1,st2]=LFM_signal(A,Phi0,T,B,F0);

當信號的中心頻率$F_0=5MHz$時

close all; clear all; clc; %%調用LFM_signal函數，觀察結果 A=1; %發射信號的振幅 Phi0=0; %發射信號的隨機初相 T=10e-6; %信號時寬 B=30e6; %信號帶寬 F0=5e6; %中頻頻率，即載頻頻率 [st1,st2]=LFM_signal(A,Phi0,T,B,F0);

總結

以上是生活随笔為你收集整理的雷达原理---线性调频信号的MATLAB仿真的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。