Commit 255a6c99 authored by 袁通's avatar 袁通
Browse files

update

parent 2c9e865c
Loading
Loading
Loading
Loading
+1 −1
Original line number Diff line number Diff line

.vscode/*
LAB3/~$b3+袁通+代育津+11810818+11910114.docx
.DS_Store
*.code-workspace

LAB2/P2_10_f.asv

0 → 100644
+74 −0
Original line number Diff line number Diff line
% Problem 2.10(f)

load inputs/lineup.mat

% For y1 ----------------------------------------
ryy1 = conv(y, flip(y));
nryy1 = -(length(y)-1):(length(y)-1);

figure;
plot(nryy1, ryy1);
title('y[n]*y[-n]');
saveas(gcf, "plots/P2_10_f_out1_1.png");
close;
%! Here N is equal to the distance of tow peak.

% N: 1000
% Peak 1: 9967
% Peak 2: 3957
% alpha: 0.4938

N = 1000;
alpha = 0.5;
yo1 = filter(1, [1 zeros(1, N-1) alpha], y);
% sound(yo1, 8192); 

% For y2 ----------------------------------------
ryy2 = conv(y2, flip(y2));
nryy2 = -(length(y2)-1):(length(y2)-1);

figure;
plot(nryy2, ryy2);
title('y_2[n]*y_2[-n]');
saveas(gcf, "plots/P2_10_f_out1_2.png");
close;

% N: 501
% Peak 1: 13160
% Peak 2: 6938
% alpha: 0.1 Too small, the peak is fake !!!!!
for alpha = [0.1 0.2 0.3 0.4 0.5 0.6 0.7]
    N = 500;
    yo2 = filter(1, [1 zeros(1, N-1) alpha], y2);
    sound(yo2, 8192);
end
N = 500;
alpha = 0.2;
yo2 = filter(1, [1 zeros(1, N-1) alpha], y2);
sound(yo2, 8192);

% For y3 ----------------------------------------
ryy3 = conv(y3, flip(y3));
nryy3 = -(length(y3)-1):(length(y3)-1);

figure;
plot(nryy3, ryy3);
title('y_3[n]*y_3[-n]');
saveas(gcf, "plots/P2_10_f_out1_3.png");
close;

% N1: 751
% alpha1: 0.4
% N2: 2252
% alpha2: 0.28
% Peak 1: 14010
% Peak 2: 4884
% Peak 3: 3597

N = 751;
N1 = 2252;
alpha = 0.5;
alpha1 = 0.28;
b = [1 zeros(1, N-1) alpha zeros(1, N1-1)] + [1 zeros(1, N1-1) alpha1 zeros(1, N-1)];
yo3 = filter(1, b, y3);
sound(yo3, 8192);
 No newline at end of file
+13 −0
Original line number Diff line number Diff line
@@ -46,10 +46,17 @@ x2=(-b2+d2^0.5)/(2*a2);
% N: 501
% Peak 1: 13160
% Peak 2: 6938
<<<<<<< Updated upstream
% alpha: 0.7

N = 501;
yo2 = filter(1, [1 zeros(1, N-1) x2], y2);
=======
% alpha: 0.1 Too small, the peak is fake !!!!!
N = 500;
alpha = 0.1;
yo2 = filter(1, [1 zeros(1, N-1) alpha], y2);
>>>>>>> Stashed changes
sound(yo2, 8192);

% For y3 ----------------------------------------
@@ -72,8 +79,14 @@ saveas(gcf, "plots/P2_10_f_out1_3.png");

N = 751;
N1 = 2252;
<<<<<<< Updated upstream
alpha1 = 0.75;
alpha2 = 0.6;
b = [1 zeros(1, 750) alpha1 zeros(1, 1500) alpha2] ;
=======
alpha = 0.5;
alpha1 = 0.28;
b = [1 zeros(1, N-1) alpha zeros(1, N1-1)] + [1 zeros(1, N1-1) alpha1 zeros(1, N-1)];
>>>>>>> Stashed changes
yo3 = filter(1, b, y3);
sound(yo3, 8192);
 No newline at end of file