error_analyse.m

通过visual DSP

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(signal));

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('求模结果误差分析')