|
|
@@ -18,7 +18,7 @@ xDecim = 1:1:PointsNum/R;
|
|
|
Fc = 3;
|
|
|
Fs = 50;
|
|
|
|
|
|
-#########################################################################################################################
|
|
|
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
ReadInDataId = fopen('C:/S5243_FFT_REPO/src/src/Sim/InputSignal.txt','r');
|
|
|
InDataSignal = fscanf(ReadInDataId,FormatSpec);
|
|
|
fclose(ReadInDataId);
|
|
|
@@ -36,7 +36,7 @@ FilterCoefs = fscanf(ReadFilterCoefsId,FormatSpec);
|
|
|
fclose(ReadFilterCoefsId);
|
|
|
|
|
|
|
|
|
-#########################################################################################################################
|
|
|
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
|
|
FreqBandOrig = Fs*(1:(PointsNum))/PointsNum;
|
|
|
FreqbandDecim = (Fs/2*R)*(1:(PointsNum))/PointsNum;
|
|
|
@@ -45,7 +45,7 @@ InDataFft = abs(fft(InDataSignal+randn(size(InDataSignal)))/PointsNum);
|
|
|
InDataFftDb = 20*log10(InDataFft);
|
|
|
|
|
|
FilteredDataFft = abs(fft(FilteredData+randn(size(FilteredData)))/PointsNum);
|
|
|
-##FilteredDataFft = abs(fft(FilteredData)/PointsNum);
|
|
|
+% FilteredDataFft = abs(fft(FilteredData)/PointsNum);
|
|
|
FilteredDataFftDb = 20*log10(FilteredDataFft);
|
|
|
FilteredDataFftDb = FilteredDataFftDb-max(FilteredDataFftDb);
|
|
|
|
|
|
@@ -54,36 +54,36 @@ ImpuseRespFftDb = 20*log10(ImpuseRespFft);
|
|
|
ImpuseRespFftDb = ImpuseRespFftDb-max(ImpuseRespFftDb);
|
|
|
|
|
|
Fband = 0:1/PointsNum:1-1/PointsNum;
|
|
|
-##Fband = 0:1/PointsNum:0.5-1/PointsNum;
|
|
|
-
|
|
|
-##InDataFftDb = InDataFftDb(1:length(Fband));
|
|
|
-##FilteredDataFftDb = FilteredDataFftDb(1:length(Fband));
|
|
|
-##ImpuseRespFftDb = ImpuseRespFftDb(1:length(Fband));
|
|
|
-
|
|
|
-#########################################################################################################################
|
|
|
-##figure('name','In Data Time/Freq', 'Numbertitle', 'off')
|
|
|
-##subplot(2,1,1)
|
|
|
-##plot(x, InDataSignal)
|
|
|
-##grid on;
|
|
|
-##grid minor;
|
|
|
-##title("In Signal")
|
|
|
-##xlabel("Time")
|
|
|
-##ylabel("Amp")
|
|
|
-##
|
|
|
-##subplot(2,1,2)
|
|
|
-##plot(Fband,InDataFftDb)
|
|
|
-##grid on;
|
|
|
-##grid minor;
|
|
|
-##title("Single-Sided Amplitude Spectrum of InData")
|
|
|
-##xlabel("f (Hz)")
|
|
|
-##ylabel("Amp")
|
|
|
+%Fband = 0:1/PointsNum:0.5-1/PointsNum;
|
|
|
+
|
|
|
+%InDataFftDb = InDataFftDb(1:length(Fband));
|
|
|
+%FilteredDataFftDb = FilteredDataFftDb(1:length(Fband));
|
|
|
+%ImpuseRespFftDb = ImpuseRespFftDb(1:length(Fband));
|
|
|
+
|
|
|
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
|
+%figure('name','In Data Time/Freq', 'Numbertitle', 'off')
|
|
|
+%subplot(2,1,1)
|
|
|
+%plot(x, InDataSignal)
|
|
|
+%grid on;
|
|
|
+%grid minor;
|
|
|
+%title("In Signal")
|
|
|
+%xlabel("Time")
|
|
|
+%ylabel("Amp")
|
|
|
+%
|
|
|
+%subplot(2,1,2)
|
|
|
+%plot(Fband,InDataFftDb)
|
|
|
+%grid on;
|
|
|
+%grid minor;
|
|
|
+%title("Single-Sided Amplitude Spectrum of InData")
|
|
|
+%xlabel("f (Hz)")
|
|
|
+%ylabel("Amp")
|
|
|
|
|
|
figure('name','Frequency compare', 'Numbertitle', 'off')
|
|
|
subplot(3,1,1)
|
|
|
plot(Fband,InDataFftDb)
|
|
|
grid on;
|
|
|
grid minor;
|
|
|
-##axis([0 25 -70 80]);
|
|
|
+%axis([0 25 -70 80]);
|
|
|
title("Single-Sided Amplitude Spectrum of Data before filter")
|
|
|
xlabel("Time")
|
|
|
ylabel("Amp")
|
|
|
@@ -92,7 +92,7 @@ subplot(3,1,2)
|
|
|
plot(Fband,ImpuseRespFftDb)
|
|
|
grid on;
|
|
|
grid minor;
|
|
|
-##axis([0 50 -60 1]);
|
|
|
+%axis([0 50 -60 1]);
|
|
|
title("Cic Frequency response")
|
|
|
xlabel("f (Hz)")
|
|
|
ylabel("Amp")
|
|
|
@@ -101,27 +101,27 @@ subplot(3,1,3)
|
|
|
plot(Fband,FilteredDataFftDb)
|
|
|
grid on;
|
|
|
grid minor;
|
|
|
-##axis([0 25 -70 80]);
|
|
|
+%axis([0 25 -70 80]);
|
|
|
title("Single-Sided Amplitude Spectrum of Data after filter")
|
|
|
xlabel("f (Hz)")
|
|
|
ylabel("Amp")
|
|
|
|
|
|
-##figure('name','Impulse Response Time/Freq', 'Numbertitle', 'off')
|
|
|
-##subplot(2,1,1)
|
|
|
-##plot(x, ImpulseResp)
|
|
|
-##grid on;
|
|
|
-##grid minor;
|
|
|
-##title("In Signal")
|
|
|
-##xlabel("Time")
|
|
|
-##ylabel("Amp")
|
|
|
-##
|
|
|
-##subplot(2,1,2)
|
|
|
-##plot(Fband,ImpuseRespFftDb)
|
|
|
-##grid on;
|
|
|
-##grid minor;
|
|
|
-##title("Single-Sided Amplitude Spectrum of InData")
|
|
|
-##xlabel("f (Hz)")
|
|
|
-##ylabel("Amp")
|
|
|
+%figure('name','Impulse Response Time/Freq', 'Numbertitle', 'off')
|
|
|
+%subplot(2,1,1)
|
|
|
+%plot(x, ImpulseResp)
|
|
|
+%grid on;
|
|
|
+%grid minor;
|
|
|
+%title("In Signal")
|
|
|
+%xlabel("Time")
|
|
|
+%ylabel("Amp")
|
|
|
+%
|
|
|
+%subplot(2,1,2)
|
|
|
+%plot(Fband,ImpuseRespFftDb)
|
|
|
+%grid on;
|
|
|
+%grid minor;
|
|
|
+%title("Single-Sided Amplitude Spectrum of InData")
|
|
|
+%xlabel("f (Hz)")
|
|
|
+%ylabel("Amp")
|
|
|
|
|
|
|
|
|
|