-
Notifications
You must be signed in to change notification settings - Fork 35
/
Just_for_test.m
102 lines (81 loc) · 1.79 KB
/
Just_for_test.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
% just for test
% Note:
% Release date: 2019/06/10
% Author: xiaoli
% Email: xiaoli644@qq.com
% Copyright (C) 2019 xiaoli
% Usage:
% This program is used to test some ideas.
% 说明:
% 本文主要用来进行一些测试。
% 注意:
% 这个不是主程序,而是用于调试。
%% 用于测试
% 方法二
% 直角扫描并作图
result2 = my_srp_phat002(gcc_all, mic, v, fs);
% 方法三
% 球坐标扫描并做图
result3 = my_srp_phat003(gcc_all, mic, v, fs);
% 其它
% 输入球坐标并进行第一次收缩迭代
result4 = my_srp_phat004(gcc_all, mic2, v, fs);
a1 = [0.1,10.1;60,90;0,90];
b1 = [1,10,10];
a2 = my_spaceShrink(a1, b1, gcc_all, mic2, v, fs)
b2 = [1,2,2];
a3 = my_spaceShrink(a2, b2, gcc_all, mic2, v, fs)
b3 = [1,0.5,0.5];
a4 = my_spaceShrink(a3, b3, gcc_all, mic2, v, fs)
% 最终结果如下,感觉还是很不错的
% a4 =
% 0.1000 10.1000
% 78.0787 79.0787
% 39.6756 40.6756
% my_rthetaToXYZ(7.9525,78.5787,40.1756)
% 5.9559 5.0288 1.5748
%% 测试sinc插值
z = [];
tt_center = 1195;
tt = (tt_center-5):0.01:(tt_center+5);
for kk = tt
zcac = my_sinc(xcorr34((tt_center-10):(tt_center+10)),(tt_center-10):(tt_center+10),kk);
z = [z, zcac];
end
figure;
stem(tt,z);
tt_center-5+0.01*(find(z==max(z))-1)
%% 测试用
clc;
close all;
for kk = 1:441
aaa(kk) = find(test_cac(:,kk)==max(test_cac(:,kk)));
end
plot(aaa,'o-');
hold on;
stem(aaa);
grid on;
%% 用于sinc画图3*3
figure;
n = 1200:1240;
x = gcc_all2(2,n);
num = 10;
[y, m] = my_sinc_vector(x, n, num);
subplot(211);
plot(n,x,'-o','linewidth',2);
hold on;
stem(n,x);
legend('原始GCC13');
xlabel('delay/points');
grid on;
subplot(212);
plot(m,y,'linewidth',2);
hold on;
stem(m,y);
legend('sinc插值后GCC13');
xlabel('delay/points');
grid on;
xlim([1215,1225]);
%% 空域收缩绘图
my_srp_phat002(gcc_all, mic, v, fs);
my_srp_phat003(gcc_all, mic, v, fs);