nathan.beaud
/
SignalLab2
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

almost the end

This commit is contained in:
Arnaud REY DU BOISSIEU 2023-03-21 17:44:22 +01:00
parent 6ecd441fd5
commit 5c2bdd7c93
3 changed files with 19 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

BIN
overview_butter(red)andfir.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 283 KiB

BIN
overview_zoomedcomparison.png Normal file
View File

Binary file not shown.
Side by Side

After

Width:  |  Height:  |  Size: 127 KiB

25
speech_analysis.m
View File

@ -15,12 +15,25 @@ y_down = downsample(y,Factor_of_reduction);
y_decimate = decimate(y,Factor_of_reduction);
audiowrite("modulator_new.wav", y_down, Fs/Factor_of_reduction);
plot (t,y,'b')
hold on
t2 = (0:length(y_down)-1)/(Fs/Factor_of_reduction);
plot (t2,y_down,'r')
t3 = (0:length(y_decimate)-1)/(Fs/Factor_of_reduction);
%plot (t,y,'b')
%hold on
%t2 = (0:length(y_down)-1)/(Fs/Factor_of_reduction);
%plot (t2,y_down,'r')
%t3 = (0:length(y_decimate)-1)/(Fs/Factor_of_reduction);
plot (t3,y_decimate,'g');
%plot (t3,y_decimate,'g');
Fc = 1000;
[b] = fir1 (30,2*Fc/Fs);
y_b = filter(b,1,y);
[a,c] = butter(8,2*Fc/Fs);
y_ac = filter(a,c,y);
plot (t,y_b);
hold on;
plot (t,y_ac,'r');
freqz(a,c);
freqz(b,1);