You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
34 lines
861 B
34 lines
861 B
clear;close all;
|
|
|
|
Nfft = 7000;
|
|
|
|
f1 =0.2e9;f2=7.8e9;
|
|
|
|
freq_array = linspace(f1,f2,Nfft);
|
|
|
|
freq = freq_array;
|
|
|
|
%% begin of SpaceTime data
|
|
folder = './TimeDomain/SpaceTime';
|
|
|
|
INCNAME = sprintf('%s/%s.TD_Vinc',folder,'PBGDDMST_src');
|
|
REFNAME1 = sprintf('%s/%s.TD_Vref',folder,'PBGDDMST_recv');
|
|
REFNAME2 = sprintf('%s/%s.TD_Vref',folder,'PBGDDMST_src');
|
|
|
|
v_in = load(INCNAME);v_in = v_in(:,[1,2]);
|
|
t = v_in(:,1);dt = t(2)-t(1);N = ceil(length(t));
|
|
v_in = v_in(1:N,2);
|
|
|
|
v_ref1 = load(REFNAME1);v_ref2 = load(REFNAME2);
|
|
v_ref1 = v_ref1(:,[1,2]);v_ref1 = v_ref1(1:N,2);
|
|
v_ref2 = v_ref2(:,[1,2]);v_ref2 = v_ref2(1:N,2);
|
|
|
|
|
|
[X_ref1] = time_to_frequency_domain([v_ref1(:,1);zeros(Nfft-length(v_ref1),1)],dt,freq_array,0);
|
|
|
|
[v_ref1_test] = frequency_to_time_domain(X_ref1,freq_array,t);
|
|
plot(t,v_ref1(:,1));
|
|
hold on;
|
|
plot(t,v_ref1_test,':o');
|
|
legend('Original','IFFT');
|
|
|
|
|