clear all;
realpart=importdata('C:\Users\zxy\Desktop\shiyan\realpart.txt');
imagpart=importdata('C:\Users\zxy\Desktop\shiyan\imagpart.txt');
x=realpart+1j*imagpart;
t=-1.99:0.01:2;
figure;
plot(t,x);
title('Real part of chirp signal');
xlabel('时间/s');
ylabel('幅度/V');
grid on;
Xw=fftshift(fft(x));
figure;
plot(t,abs(Xw));
title('Magnitude spectrum of chirp signal');
xlabel('频率/MHz');
ylabel('幅度/V');
grid on;
Hw=conj(Xw);
figure;
plot(t,Hw);
title('Hw');
xlabel('频率/MHz');
ylabel('幅度/V');
grid on;
ht=abs(ifftshift(ifft(Hw)));
figure;
plot(t,ht);
title('ht');
xlabel('时间/s');
ylabel('幅度/V');
grid on;
Yw=Xw.*Hw;
yt=ifftshift(ifft(Yw));
yt1=20*log10(abs(yt)./max(abs(yt)));
figure;
plot(t,yt1);
title('yt');
xlabel('时间/s');
ylabel('幅度/V');
grid on;
axis([-1 1 -50 0]);