jens/GnuRadio
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

- Xs not involved in calculation. Use X instead

Browse Source 
- added plot legend

git-svn-id: http://moon:8086/svn/software/trunk/projects/GnuRadio@447 b431acfa-c32f-4a4a-93f1-934dc6c82436
This commit is contained in:
jens
2019-05-23 17:20:36 +00:00
parent 124def841f
commit 956ff67a08
1 changed files with 4 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files
matlab/peak_detect.m
+4 -4
View File
@@ -90,7 +90,7 @@ for k=1:numFrames
% Peak detection % Peak detection
Xs = (1-as)*Xs + as*X; Xs = (1-as)*Xs + as*X;
gate_peak = (Xs-Xm) >= Ad; gate_peak = (X-Xm) >= Ad;
peak_count = gate_peak.*peak_count + gate_peak; peak_count = gate_peak.*peak_count + gate_peak;
peaks = peak_count > 0; peaks = peak_count > 0;
@@ -129,7 +129,7 @@ for k=1:numFrames
peak_pos = find(peaks); peak_pos = find(peaks);
% Sort peaks % Sort peaks
[v, n] = sort((Xs-Xm)(peak_pos), 'descend'); [v, n] = sort((X-Xm)(peak_pos), 'descend');
peak_cand_sorted = peak_pos(n); peak_cand_sorted = peak_pos(n);
@@ -216,9 +216,9 @@ for k=1:numFrames
% Output % Output
subplot(2, 1, 1) subplot(2, 1, 1)
plot(1:Nh, X, peak_pos, X(peak_pos), 'ro', 1:Nh, peak_box_values, 'm-'); grid; axis([0, Nh, -60, 40]); plot(1:Nh, X, peak_pos, X(peak_pos), 'ro', 1:Nh, peak_box_values, 'm-'); grid; axis([0, Nh, -60, 40]); legend('X', 'Peaks', 'BBox')
subplot(2, 1, 2) subplot(2, 1, 2)
plot(1:Nh, Xs-Xm, 1:Nh, Xm); grid; axis([0, Nh, -60, 40]); plot(1:Nh, Xs-Xm, 1:Nh, Xm); grid; axis([0, Nh, -60, 40]); legend('X_s - Xm', 'X_m')
pause(1/30); pause(1/30);
end end