close all;
clear all;

load D:\dsp_experiment\cz\chirp\lfmimag.dat  %載入調頻信號的實部和虛部
load D:\dsp_experiment\cz\chirp\lfmreal.dat
signal=lfmreal+j*lfmimag;
result=abs(fft(signal));                      %用malab函數進行fft變換然后求模

figure(10)
plot(angle())

load D:\dsp_experiment\cz\chirp\lfmmod.dat    %載入DSP處理lfm信號結果

figure(1);
subplot(2,1,1);plot(result);title('matlab處理lfm結果');
subplot(2,1,2);plot(lfmmod);title('DSP處理lfm');

%誤差分析
e=(lfmmod/max(lfmmod))-(result/max(result));
figure(2);
plot(abs(e));title('誤差分析')

%DSP的fft結果求模
 load D:\dsp_experiment\cz\chirp\refft.dat
 load D:\dsp_experiment\cz\chirp\imfft.dat
 signalfft=refft+j*imfft;
 result0=abs(signalfft);

%fft結果誤差分析
e1=(fft(signal)/max(fft(signal)))-(signalfft/max(signalfft))
figure(3);
plot(e1);title('fft結果誤差分析')

%求模結果誤差分析
e2=(lfmmod/max(lfmmod))-(result0/max(result0));
figure(4);
plot(e2);title('求模結果誤差分析')