From 0dec9505c5422a2adbe3b2d1ef5e702608e4983d Mon Sep 17 00:00:00 2001 From: unknown <2779155576@qq.com> Date: Mon, 30 Dec 2024 14:06:53 +0800 Subject: [PATCH] =?UTF-8?q?=E4=BF=AE=E6=94=B9z=5Fdsp.m=E7=9B=B8=E5=85=B3?= =?UTF-8?q?=E5=87=BD=E6=95=B0=E4=BB=A5=E6=89=B9=E9=87=8F=E6=89=AB=E6=8F=8F?= =?UTF-8?q?=E7=BA=BF=E8=B7=AF=E5=8F=82=E6=95=B0=E5=92=8C=E6=B3=A2=E5=BD=A2?= =?UTF-8?q?=20=E4=BF=AE=E6=94=B9TailCorr=5FTest=E7=9A=84=E5=90=8D=E5=AD=97?= =?UTF-8?q?=E4=BE=BF=E4=BA=8E=E5=8C=BA=E5=88=86Verdi=E5=B9=B3=E5=8F=B0?= =?UTF-8?q?=E7=94=A8=E7=9A=84=E8=84=9A=E6=9C=AC=E5=92=8CWindows=E5=B9=B3?= =?UTF-8?q?=E5=8F=B0?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit Signed-off-by: unknown <2779155576@qq.com> --- ...{TailCorr_Test.m => TailCorr_Test_Verdi.m} | 0 script_m/diff_plot_py.m | 17 +- script_m/z_dsp.m | 382 ++++++++++++------ 3 files changed, 268 insertions(+), 131 deletions(-) rename script_m/{TailCorr_Test.m => TailCorr_Test_Verdi.m} (100%) diff --git a/script_m/TailCorr_Test.m b/script_m/TailCorr_Test_Verdi.m similarity index 100% rename from script_m/TailCorr_Test.m rename to script_m/TailCorr_Test_Verdi.m diff --git a/script_m/diff_plot_py.m b/script_m/diff_plot_py.m index 21c5eed..acadcbf 100644 --- a/script_m/diff_plot_py.m +++ b/script_m/diff_plot_py.m @@ -1,5 +1,5 @@ %compare FIL with python script -function diff_plot_py(fs,iir_out, Script_out,title1,title2,a,amp,edge) +function diff_plot_py(fs,iir_out, Script_out,title1,title2,a,amp,edge,fileID) %输入数据长度不等时取其公共部分 N = min(length(iir_out),length(Script_out)); iir_out = iir_out(1:N); @@ -49,13 +49,7 @@ plot_p = @(x)[ text(n(x), diff(x)+diff(x)*0.1, ['(',num2str(n(x)),',',num2str(diff(x)),')'],'color','k'); ]; -%标注出关心的点 -%plot_p(n_edge(1));%下降沿 -%plot_p(n50(1)); %下降沿20ns -%plot_p(n1000(1)); %下降沿1us - ne(1) = 1; -%plot_p(ne(end)); %误差小于万分之一 % [diff_max,R_mpos] = max(abs(diff));%误差最大值 % plot_p(R_mpos); @@ -63,16 +57,15 @@ ne(1) = 1; if a(2) <= 5e-6 plot_p(n_edge(1));%下降沿 % plot_p(R_mpos); -elseif a(2) > 5e-6 +elseif a(2) == 20e-6 plot_p(n50(1)); %下降沿20ns plot_p(n1000(1)); %下降沿1us plot_p(ne(end)); %误差小于万分之一 - fprintf("Falling edge of 20ns~40ns mean :%.4e\t std :%.4e\t",mean(diff(n20_40)),std(diff(n20_40))); - fprintf("Falling edge of 1us~1.1us mean :%.4e\t std :%.4e\t",mean(diff(n1000_1100)),std(diff(n1000_1100))); + fprintf(fileID,"Falling edge of 20ns~40ns mean :%.4e\t std :%.4e\t",mean(diff(n20_40)),std(diff(n20_40))); + fprintf(fileID,"Falling edge of 1us~1.1us mean :%.4e\t std :%.4e\t",mean(diff(n1000_1100)),std(diff(n1000_1100))); % fprintf("The error after falling edge of 1us is:%.4e\t",diff(n1000(1))); % fprintf("The time of erroe less than 1e-4 is :%.4e us\n",(n(ne(end))-n(n_edge(1)))); - fprintf("The mean and std stably less than 1e-4 is :%.4e s\n",(n(n_common)-n(n_edge(1)))); - + fprintf(fileID,"The mean and std stably less than 1e-4 is :%.4e s\n",(n(n_common)-n(n_edge(1)))); end diff --git a/script_m/z_dsp.m b/script_m/z_dsp.m index 8032c08..e69c58f 100644 --- a/script_m/z_dsp.m +++ b/script_m/z_dsp.m @@ -1,7 +1,22 @@ clc;clear;close all % hdlsetuptoolpath('ToolName','Xilinx Vivado','ToolPath','D:\SoftWare\Xilinx\Vivado\2019.2\bin\vivado.bat'); +%%配置参数 +fs_L = 0.75e9; %硬件频率 +fs_H = 12e9; %以高频近似理想信号 +TargetFrequency = 3e9; +G = 1; +DownSample = 2; +simulink_time = 20e-6; %1.5*16e-6;1.5e-3 +intp_mode = 3; %0不内插，1内插2倍，2内插4倍,3内插8倍 +dac_mode_sel = 0; %选择DAC模式，0出八路，1邻近插值，2邻近插值 +route_num = 5; %线路个数 +env_num = 7; %包络个数 +Ideal2Low = fs_H/(fs_L/2); +Ideal2Target = fs_H/TargetFrequency; +%% 添加路径、产生包络、配置S21参数、使用脚本计算 +%%添加路径 % addpath(genpath('D:\Work\EnvData')); % addpath(genpath('D:\Work\EnvData\data-v2')); % addpath(genpath('D:\Work\TailCorr_20241008_NoGit')); @@ -12,43 +27,32 @@ py.importlib.reload(obj1); cd("D:\Work\TailCorr_20241008_NoGit"); obj2 = py.importlib.import_module('wave_calculation'); py.importlib.reload(obj2); -cd("D:\Work\TailCorr"); - -fs_L = 0.75e9; %硬件频率 -fs_H = 12e9; %以高频近似理想信号 -TargetFrequency = 3e9; -Ideal2Low = fs_H/(fs_L/2); -Ideal2Target = fs_H/TargetFrequency; -G = 1; -DownSample = 2; -simulink_time = 20e-6; %1.5*16e-6;1.5e-3 -intp_mode = 3; %0不内插，1内插2倍，2内插4倍,3内插8倍 -dac_mode_sel = 0; %选择DAC模式，0出八路，1邻近插值，2邻近插值 - +cd("D:\Work\TailCorr"); +%%产生包络 %按点数产生理想方波 -amp_rect = 1.5e4; -%单位是ns front是到达时间，flat是持续时间，lagging是后边还有多少个0，会影响脚本的修正时间 -[front(1), flat(1), lagging(1)] = deal(50,100,7400);% 50,100,7400;100ns方波 -[front(2), flat(2), lagging(2)] = deal(50,4000,11500);% 50,4000,11500；4us方波 +% amp_rect = 1.5e4; +% %单位是ns front是到达时间，flat是持续时间，lagging是后边还有多少个0，会影响脚本的修正时间 +% [front(1), flat(1), lagging(1)] = deal(50,100,7400);% 50,100,7400;100ns方波 +% [front(2), flat(2), lagging(2)] = deal(50,4000,11500);% 50,4000,11500；4us方波 +% +% for i = 1:2 +% front_H(i) = front(i)*fs_H/1e9; flat_H(i) = flat(i)*fs_H/1e9; lagging_H(i) = lagging(i)*fs_H/1e9; +% wave_pre{i} = amp_rect*cat(2,zeros(1,front_H(i)),ones(1,flat_H(i)),zeros(1,lagging_H(i)));%脚本的单位是点数 +% end -for i = 1:2 - front_H(i) = front(i)*fs_H/1e9; flat_H(i) = flat(i)*fs_H/1e9; lagging_H(i) = lagging(i)*fs_H/1e9; - wave_pre{i} = amp_rect*cat(2,zeros(1,front_H(i)),ones(1,flat_H(i)),zeros(1,lagging_H(i)));%脚本的单位是点数 -end - -%%% flattop波 +%flattop波 A = 1.5e4; [edge(1), length_flattop(1)] = deal(2,30);%ns，在fsn_L取1时是参数里的length [edge(2), length_flattop(2)] = deal(4,30); [edge(3), length_flattop(3)] = deal(4,50); -[edge(4), length_flattop(4)] = deal(6,50); - -for i = 1:4 +[edge(4), length_flattop(4)] = deal(4,1000); +[edge(5), length_flattop(5)] = deal(100,10000); +for i = 1:5 [edge_H(i), length_H(i)] = deal(edge(i)*fs_H/1e9,length_flattop(i)*fs_H/1e9); - wave_pre{i+2} = flattop(A, edge_H(i), length_H(i), 1); + wave_pre{i} = flattop(A, edge_H(i), length_H(i), 1); end -%%% acz波 +%acz波 amplitude = 1.5e4; carrierFreq = 0.000000; @@ -64,27 +68,47 @@ length_acz(2) = 50; for i = 1:2 length_acz_H(i) = int32(length_acz(i)*fs_H/1e9); - wave_pre{i+6} = real(double(py.acz.aczwave(amplitude, length_acz_H(i), carrierFreq,carrierPhase, dragAlpha,thf, thi, lam2, lam3))); + wave_pre{i+5} = real(double(py.acz.aczwave(amplitude, length_acz_H(i), carrierFreq,carrierPhase, dragAlpha,thf, thi, lam2, lam3))); end -% signalAnalyzer(wave_pre{2},'SampleRate',fs_H); -for i = 1:8 +for i = 1:7 wave_pre{i} = cat(2,wave_pre{i},zeros(1,floor(simulink_time*fs_H))); %校正前的高频信号 wave_preL{i} = wave_pre{i}(1:Ideal2Low:end); %校正前的低频信号 end -% signalAnalyzer(HardwareMeanIntpDataAlign{1},'SampleRate',3e9); +%%S21参数 +amp_real{1}= [0.025 0.015 0.0002 0.2 0 0]; +amp_imag{1}= [0 0 0 0 0 0]; +time_real{1} = [-1/250, -1/650, -1/1600 -1/20 0 0]; +time_imag{1} = [0 0 0 0 0 0]; -%%%python脚本校正结果 -%S21参数 -amp_real = [0.025 0.015 0.0002 0.2 0 0]; -amp_imag = [0 0 0 0 0 0]; -time_real = [-1/250, -1/650, -1/1600 -1/20 0 0]; -time_imag = [0 -1/300 -1/500 0 0 0]; +amp_real{2}= [0.025 0.015 0.0002 0.2 0 0]; +amp_imag{2}= [0 0 0 0 0 0]; +time_real{2} = [-1/250, -1/650, -1/1600 -1/20 0 0]; +time_imag{2} = [0 -1/300 -1/500 0 0 0]; -amp_routing = amp_real + 1j*amp_imag; -time_routing = time_real + 1j*time_imag; -tau = -1./time_routing; +amp_real{3}= [0.025 0.009 0.0002 0.2 0 0]; +amp_imag{3}= [0 0.012 0 0 0 0]; +time_real{3} = [-1/250, -1/650, -1/1600 -1/20 0 0]; +time_imag{3} = [0 -1/300 -1/500 0 0 0]; + +amp_real{4}= [0.025 0.015 0.0002 0.2 0 0]; +amp_imag{4}= [0 0 0 0 0 0]; +time_real{4} = [-1/250, -1/2000, -1/1600 -1/20 0 0]; +time_imag{4} = [0 -1/15 -1/50 0 0 0]; + +amp_real{5}= [0.025 0.009 0.0002 0.2 0 0]; +amp_imag{5}= [0 0.012 0 0 0 0]; +time_real{5} = [-1/250, -1/2000, -1/1600 -1/20 0 0]; +time_imag{5} = [0 -1/15 -1/50 0 0 0]; + +for i = 1:5 + amp_routing{i} = amp_real{1,i} + 1j*amp_imag{1,i}; + time_routing{i} = time_real{1,i} + 1j*time_imag{1,i}; + tau{i} = -1./time_routing{i}; +end + +%%python脚本校正结果 convolve_bound = int8(3); calibration_time = int32(20e3); @@ -93,107 +117,227 @@ sampling_rateL = int64(fs_L/2); sampling_rate = int64(fs_H); %校正后的高频信号 -for i = 1:8 - wave_cal = cell(py.wave_calculation.wave_cal(wave_pre{i}, amp_real, amp_imag, time_real, time_imag, convolve_bound, calibration_time, cal_method, sampling_rate)); - wave_revised{i} = double(wave_cal{1,1}); - wave_calL = cell(py.wave_calculation.wave_cal(wave_preL{i}, amp_real, amp_imag, time_real, time_imag, convolve_bound, calibration_time, cal_method, sampling_rateL)); - wave_revisedL{i} = double(wave_calL{1,1}); +for m = 1:route_num + for n = 1:env_num + wave_cal = cell(py.wave_calculation.wave_cal(wave_pre{1,n}, amp_real{1,m}, amp_imag{1,m}, time_real{1,m}, time_imag{1,m}, convolve_bound, calibration_time, cal_method, sampling_rate)); + wave_revised{m,n} = double(wave_cal{1,1}); + wave_calL = cell(py.wave_calculation.wave_cal(wave_preL{1,n}, amp_real{1,m}, amp_imag{1,m}, time_real{1,m}, time_imag{1,m}, convolve_bound, calibration_time, cal_method, sampling_rateL)); + wave_revisedL{m,n} = double(wave_calL{1,1}); + end + alpha{m} = double(wave_calL{1,2}); + beta{m} = double(wave_calL{1,3}); end +% signalAnalyzer(wave_pre{1,1},'SampleRate',fs_H); + %校正后的低频信号 -alpha = double(wave_calL{1,2}); -beta = double(wave_calL{1,3}); -beta(5:6) = 0; + alpha_wideth=32; beta_width=32; -alphaFixRe = ceil((2^(alpha_wideth-1))*real(alpha)); -alphaFixIm = ceil((2^(alpha_wideth-1))*imag(alpha)); -betaFixRe = ceil((2^(beta_width-1))*real(beta)); -betaFixIm = ceil((2^(beta_width-1))*imag(beta)); - -%%%仿真 -for i = 1:8 - options=simset('SrcWorkspace','current'); - sim('z_dsp',[0,simulink_time]); - sim2m = @(x)reshape(logsout.get(x).Values.Data,[],1); - dout0{i} = sim2m("dout0"); - dout1{i} = sim2m("dout1"); - dout2{i} = sim2m("dout2"); - dout3{i} = sim2m("dout3"); - - N(i) = length(dout0{i}); - cs_wave{i} = zeros(4*N(i),1); - - cs_wave{i}(1:4:4*N) = dout0{i}; - cs_wave{i}(2:4:4*N) = dout1{i}; - cs_wave{i}(3:4:4*N) = dout2{i}; - cs_wave{i}(4:4:4*N) = dout3{i}; - - HardwareMeanIntpData{i} = cs_wave{i};%硬件校正后内插 - DownsamplingBy12GData{i} = wave_revised{i}(1:Ideal2Target:end); - [DownsamplingBy12GDataAlign{i},HardwareMeanIntpDataAlign{i},Delay(i)] = ... - alignsignals(DownsamplingBy12GData{i}(1:round(TargetFrequency*20e-6)),HardwareMeanIntpData{i}(1:round(TargetFrequency*20e-6)),"Method","xcorr"); +%定点化系数 +for i = 1:route_num + alphaFixRe{i} = ceil((2^(alpha_wideth-1))*real(alpha{i})); + alphaFixIm{i} = ceil((2^(alpha_wideth-1))*imag(alpha{i})); + betaFixRe{i} = ceil((2^(beta_width-1))*real(beta{i})); + betaFixIm{i} = ceil((2^(beta_width-1))*imag(beta{i})); end - -% signalAnalyzer(DownsamplingBy12GDataAlign{1},HardwareMeanIntpDataAlign{1},'SampleRate',3e9); - +%% 仿真 +for m = 1:route_num + for n = 1:env_num + optnons=simset('SrcWorkspace','current'); + sim('z_dsp_FIL',[0,simulink_time]); + sim2m = @(x)reshape(logsout.get(x).Values.Data,[],1); + dout0{m,n} = sim2m("dout0"); + dout1{m,n} = sim2m("dout1"); + dout2{m,n} = sim2m("dout2"); + dout3{m,n} = sim2m("dout3"); + + N = length(dout0{m,n}); + cs_wave{m,n} = zeros(4*N,1); + + cs_wave{m,n}(1:4:4*N) = dout0{m,n}; + cs_wave{m,n}(2:4:4*N) = dout1{m,n}; + cs_wave{m,n}(3:4:4*N) = dout2{m,n}; + cs_wave{m,n}(4:4:4*N) = dout3{m,n}; + + HardwareMeanIntpData{m,n} = cs_wave{m,n};%硬件校正后内插 + DownsamplingBy12GData{m,n} = wave_revised{m,n}(1:Ideal2Target:end); + [DownsamplingBy12GDataAlign{m,n},HardwareMeanIntpDataAlign{m,n},Delay(m,n)] = ... + alignsignals(DownsamplingBy12GData{m,n}(1:round(TargetFrequency*20e-6)),HardwareMeanIntpData{m,n}(1:round(TargetFrequency*20e-6)),"Method","xcorr"); + end +end +% signalAnalyzer(wave_revised,'SampleRate',3e9); %% 绘图并保存 close all; + Amp = 1.5e4; FallingEdge = [ - 150e-9,4050e-9,...%矩形波 - 30e-9,30e-9,50e-9,50e-9,...%flattop +% 150e-9,4050e-9,...%矩形波 + 30e-9,30e-9,50e-9,1000e-9,10000e-9,...%flattop 30e-9,50e-9%acz ]; + name = [ - "rect_100ns_校正后的波形_下降沿后10ns.fig","rect_100ns_校正后的波形_下降沿后1us.fig"; - "rect_4us_校正后的波形_下降沿后10ns.fig","rect_4us_校正后的波形_下降沿后1us.fig"; - "flattop_上升沿2ns_持续时间30ns_校正后的波形_下降沿后10ns.fig","flattop_上升沿2ns_持续时间30ns_校正后的波形_下降沿后1us.fig"; - "flattop_上升沿4ns_持续时间30ns_校正后的波形_下降沿后10ns.fig","flattop_上升沿4ns_持续时间30ns_校正后的波形_下降沿后1us.fig"; - "flattop_上升沿4ns_持续时间50ns_校正后的波形_下降沿后10ns.fig","flattop_上升沿4ns_持续时间50ns_校正后的波形_下降沿后1us.fig"; - "flattop_上升沿6ns_持续时间50ns_校正后的波形_下降沿后10ns.fig","flattop_上升沿6ns_持续时间50ns_校正后的波形_下降沿后1us.fig"; - "acz_持续时间30ns_校正后的波形_下降沿后10ns.fig","acz_持续时间30ns_校正后的波形_下降沿后1us.fig"; - "acz_持续时间50ns_校正后的波形_下降沿后10ns.fig","acz_持续时间50ns_校正后的波形_下降沿后1us.fig"; + "第一组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后10ns",... + "第一组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后10ns",... + "第一组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后10ns",... + "第一组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后10ns",... + "第一组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后10ns",... + "第一组S21参数_acz_持续时间30ns_下降沿后10ns",... + "第一组S21参数_acz_持续时间50ns_下降沿后10ns"; + "第二组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后10ns",... + "第二组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后10ns",... + "第二组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后10ns",... + "第二组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后10ns",... + "第二组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后10ns",... + "第二组S21参数_acz_持续时间30ns_下降沿后10ns",... + "第二组S21参数_acz_持续时间50ns_下降沿后10ns"; + "第三组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后10ns",... + "第三组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后10ns",... + "第三组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后10ns",... + "第三组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后10ns",... + "第三组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后10ns",... + "第三组S21参数_acz_持续时间30ns_下降沿后10ns",... + "第三组S21参数_acz_持续时间50ns_下降沿后10ns"; + "第四组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后10ns",... + "第四组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后10ns",... + "第四组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后10ns",... + "第四组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后10ns",... + "第四组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后10ns",... + "第四组S21参数_acz_持续时间30ns_下降沿后10ns",... + "第四组S21参数_acz_持续时间50ns_下降沿后10ns"; + "第五组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后10ns",... + "第五组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后10ns",... + "第五组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后10ns",... + "第五组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后10ns",... + "第五组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后10ns",... + "第五组S21参数_acz_持续时间30ns_下降沿后10ns",... + "第五组S21参数_acz_持续时间50ns_下降沿后10ns"; + "第一组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后1us",... + "第一组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后1us",... + "第一组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后1us",... + "第一组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后1us",... + "第一组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后1us",... + "第一组S21参数_acz_持续时间30ns_下降沿后1us",... + "第一组S21参数_acz_持续时间50ns_下降沿后1us"; + "第二组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后1us",... + "第二组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后1us",... + "第二组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后1us",... + "第二组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后1us",... + "第二组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后1us",... + "第二组S21参数_acz_持续时间30ns_下降沿后1us",... + "第二组S21参数_acz_持续时间50ns_下降沿后1us"; + "第三组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后1us",... + "第三组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后1us",... + "第三组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后1us",... + "第三组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后1us",... + "第三组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后1us",... + "第三组S21参数_acz_持续时间30ns_下降沿后1us",... + "第三组S21参数_acz_持续时间50ns_下降沿后1us"; + "第四组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后1us",... + "第四组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后1us",... + "第四组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后1us",... + "第四组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后1us",... + "第四组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后1us",... + "第四组S21参数_acz_持续时间30ns_下降沿后1us",... + "第四组S21参数_acz_持续时间50ns_下降沿后1us"; + "第五组S21参数_flattop_上升沿2ns_持续时间30ns_下降沿后1us",... + "第五组S21参数_flattop_上升沿4ns_持续时间30ns_下降沿后1us",... + "第五组S21参数_flattop_上升沿4ns_持续时间50ns_下降沿后1us",... + "第五组S21参数_flattop_上升沿4ns_持续时间1000ns_下降沿后1us",... + "第五组S21参数_flattop_上升沿100ns_持续时间10000ns_下降沿后1us",... + "第五组S21参数_acz_持续时间30ns_下降沿后1us",... + "第五组S21参数_acz_持续时间50ns_下降沿后1us"; ]; -Delay_mode = mode(Delay); -for i = 1:8 - start_time(i) = abs(Delay_mode)/(TargetFrequency/1e9)*1e-9;%由于相位修正后会有偏移的点数，所以需要考虑上这个偏移的时间，采样率为3GHz，3个点对应1ns - edge_Align(i) = FallingEdge(i) + start_time(i); - tmp(i) = edge_Align(i) + 10e-9; - a{i} = [start_time(i)-5e-9 tmp(i)];%[1/fs_H 50e-9];[50e-9 1.5e-6],[500e-9+10e-9 tmp-20e-9] - b{i} = [tmp(i) 10e-6]; - fig1 = figure('Units','normalized','Position',[0.000390625,0.517361111111111,0.49921875,0.422916666666667]); - diff_plot_py(TargetFrequency,HardwareMeanIntpDataAlign{i}', DownsamplingBy12GDataAlign{i}(1:floor(TargetFrequency*20e-6)),'HardwareRevised','ScriptRevised',a{i},Amp,edge_Align(i)); - title(name(i,1),Interpreter="none"); - % savefig(name(i,1)); - fig2 = figure('Units','normalized','Position',[0.000390625,0.034027777777778,0.49921875,0.422916666666667]); - diff_plot_py(TargetFrequency,HardwareMeanIntpDataAlign{i}', DownsamplingBy12GDataAlign{i}(1:floor(TargetFrequency*20e-6)),'HardwareRevised','ScriptRevised',b{i},Amp,edge_Align(i)); - title(name(i,2),Interpreter="none"); - % savefig(name(i,2)); + +Delay_mode = mode(Delay,'all'); +fileID = fopen('20241226_output.txt', 'w'); +if fileID == -1 + disp('文件打开失败'); +else + end +for m = 1:route_num + for n = 1:env_num + start_time(n) = abs(Delay_mode)/(TargetFrequency/1e9)*1e-9;%由于相位修正后会有偏移的点数，所以需要考虑上这个偏移的时间，采样率为3GHz，3个点对应1ns + edge_Align(n) = FallingEdge(n) + start_time(n); + tmp(n) = edge_Align(n) + 10e-9; + a{n} = [start_time(n)-5e-9 tmp(n)];%[1/fs_H 50e-9];[50e-9 1.5e-6],[500e-9+10e-9 tmp-20e-9] + b{n} = [tmp(n) 20e-6]; + fig1 = figure('Units','normalized','Position',[0.000390625,0.517361111111111,0.49921875,0.422916666666667]); + diff_plot_py(TargetFrequency,HardwareMeanIntpDataAlign{m,n}', DownsamplingBy12GDataAlign{m,n}(1:floor(TargetFrequency*20e-6)),'HardwareRevised','ScriptRevised',a{n},Amp,edge_Align(n),fileID); + title(name(m,n),Interpreter="none"); + savefig(name(m,n)); + fig2 = figure('Units','normalized','Position',[0.000390625,0.034027777777778,0.49921875,0.422916666666667]); + diff_plot_py(TargetFrequency,HardwareMeanIntpDataAlign{m,n}', DownsamplingBy12GDataAlign{m,n}(1:floor(TargetFrequency*20e-6)),'HardwareRevised','ScriptRevised',b{n},Amp,edge_Align(n),fileID); + title(name(m,n),Interpreter="none"); + savefig(name(m+5,n)); + end +end +fclose(fileID); %% 可视化S21参数 t = 0:1/(1e2):10000; - -for i = 1:1:length(amp_routing) - S21_time(:,i) = amp_routing(i)*exp(time_routing(i)*t); +for m = 1:route_num + for n = 1:1:length(amp_routing{1,m}) + S21_time{m}(:,n) = amp_routing{1,m}(n)*exp(time_routing{1,m}(n)*t); + end + figure + plot(t*1e-9,real(sum(S21_time{m},2))); + grid on + title("s(t)"); end -figure -plot(t*1e-9,real(sum(S21_time,2))); -grid on -title("s(t)"); -savefig("S21参数"); +% savefig("S21参数"); % signalAnalyzer(real(sum(S21_time,2)),'SampleRate',1e11);%时间是1ns，还得加上采样率 % rmpath(genpath('D:\Work\EnvData')); % rmpath(genpath('D:\Work\EnvData\data-v2')); % rmpath(genpath('D:\Work\TailCorr_20241008_NoGit')); -%% 图像可视化 -cd("D:\Work\TailCorr\仿真结果\20241101_125M八倍内插至1G_第1组S21参数") -for i = 1:8 - close all - open(name(i,1)); - open(name(i,2)); - pause() +%% 查看图像 +for m = 1:35 + figure(2*m-1) + figure(2*m) + pause(0.2); end +%% 读图像并查看 +% cd("D:\Work\TailCorr\仿真结果\20241101_125M八倍内插至1G_第1组S21参数"); +close all +for n = 1:route_num + for m = 1:env_num + open(strcat(name(n,m),'.fig')); + open(strcat(name(n+5,m),'.fig')); + pause(); + end +end +%% +close all +fid = fopen('20241226_output.txt','r'); +data = textscan(fid,'Falling edge of 20ns~40ns mean :%s std :%s Falling edge of 1us~1.1us mean :%s std :%s The mean and std stably less than 1e-4 is :%s s'); +fclose(fid); +data{1} = cellfun(@str2num,data{1}); +data{2} = cellfun(@str2num,data{2}); +data{3} = cellfun(@str2num,data{3}); +data{4} = cellfun(@str2num,data{4}); +data{5} = cellfun(@str2num,data{5}); +title_name = ["下降沿后20ns~40ns误差的平均值","下降沿后20ns~40ns误差的标准差","下降沿后1us~1.1us误差的平均值","下降沿后1us~1.1us误差的标准差","加窗参数"]; +err_threshold = [1e-3 1e-3 1e-4 3e-4 5e-5]; +figure() +tiledlayout('flow','TileSpacing','tight') +colors = lines(route_num); + +for m = 1:5 + nexttile + hold on + for i = 1:(route_num) + idx = (i-1)*(length(falling20_mean)/route_num) + 1 : i*(length(falling20_mean)/route_num); + plot(idx,abs(data{m}(idx)),'-o','Color', colors(i, :)); + end + yline(err_threshold(m),'--r'); + title(title_name(m)); + set(gca,'YScale','log'); +end + +%% +figure() +semilogy(abs(falling20_mean),'-o');