MATALB 信号时延_matlab_信号时延_时延_

% 针对程序delaytime.m函数进行测试 clear all; close all; clc Fs = 15e9; % 采样频率 dt = 1/Fs; % 时间间隔 T = 2e-6; % 啁啾信号持续时间长度 【us】 delay = T/4; % 回波信号相对发射信号延迟时间 t = 0:dt:T; % 原始信号时间向量 % 啁啾信号 f0 = 2e9; % 起始频率 【Hz】 B = 500e6; % 带宽 【Hz】s k = B/T; % 啁啾率 y = cos(2*pi*(f0.*t + 0.5*k*t.^2)); N = length(t); % 原始信号信号长度 f = (-N/2:N/2-1)*(Fs/N); % 以零为中心的频率轴 [t1,f1,y1,y2] = delaytime(t,y,Fs,delay); % 原始信号补零及时延 Y = fftshift( fft(y) ); % 发射信号fft Y1 = fftshift( fft(y1) ); % 参考信号fft Y2 = fftshift( fft(y2) ); % 回波信号fft %--------------------绘制发射信号时域和频域波形-----------------------------% fig1 = figure('name','发射信号','color','w'); % 定义图像句柄、图像名称、背景颜色 subplot(2,1,1); plot(t,y); % 绘制原始信号时域波形(发射信号) hold on; title('发射信号');xlabel('Time (s)');ylabel('Amplitude'); subplot(2,1,2) plot(f,abs(Y).^2/N); % 绘制原始信号频域波形(发射信号) title('Amplitude Spectrum');xlabel('Frequency (Hz)');ylabel('Amplitude'); % 原始信号时域频域波形——保存 filepath = pwd; % 保存当前工作目录 cd('.\figure'); % 跳到下一级figure文件夹下 print(gcf,'.\发射信号20201105_500M.png','-dpng','-r300'); % 输出格式png cd(filepath); % 切回原工作目录 %--------------绘制发射信号时域和频域波形 end------------------------------% %--------------------绘制参考信号时域和频域波形-----------------------------% fig2 = figure('name','参考信号','color','w'); % 定义图像句柄、图像名称、背景颜色 subplot(2,1,1); plot(t1,y1); % 绘制原始信号时域波形 hold on; title('参考信号');xlabel('Time (s)');ylabel('Amplitude'); subplot(2,1,2) plot(f1,abs(Y1).^2/length(y1)); % 绘制原始信号频域波形 title('Amplitude Spectrum');xlabel('Frequency (Hz)');ylabel('Amplitude'); % 原始信号时域频域波形——保存 filepath = pwd; % 保存当前工作目录 cd('.\figure'); % 跳到下一级figure文件夹下 print(gcf,'.\参考信号20201105_500M.png','-dpng','-r300'); % 输出格式png cd(filepath); % 切回原工作目录 %--------------绘制回波信号时域和频域波形 end------------------------------% %--------------------绘制参考信号时域和频域波形-----------------------------% fig3 = figure('name','回波信号','color','w'); % 定义图像句柄、图像名称、背景颜色 subplot(2,1,1); plot(t1,y2); % 绘制原始信号时域波形 hold on; title('参考信号');xlabel('Time (s)');ylabel('Amplitude'); subplot(2,1,2) plot(f1,abs(Y2).^2/length(y2)); % 绘制原始信号频域波形 title('Amplitude Spectrum');xlabel('Frequency (Hz)');ylabel('Amplitude'); % 原始信号时域频域波形——保存 filepath = pwd; % 保存当前工作目录 cd('.\figure'); % 跳到下一级figure文件夹下 print(gcf,'.\回波信号20201105_500M.png','-dpng','-r300'); % 输出格式png cd(filepath); % 切回原工作目录 %--------------绘制回波信号时域和频域波形 end------------------------------%