|
a |
|
b/Code/datasetgeneration.m |
|
|
1 |
%TRAINING DATASET GENERATION |
|
|
2 |
x=0:0.01:7; |
|
|
3 |
|
|
|
4 |
ph = 2*pi*rand; |
|
|
5 |
y = sin(3*x + ph); |
|
|
6 |
outpure = y; |
|
|
7 |
outnoisy = awgn(y,15); |
|
|
8 |
|
|
|
9 |
%frequency=3, low noise |
|
|
10 |
for i=1:299 |
|
|
11 |
ph = 2*pi*rand; |
|
|
12 |
y = sin(3*x + ph); |
|
|
13 |
outpure = [outpure;y]; |
|
|
14 |
outnoisy = [outnoisy;awgn(y,15)]; |
|
|
15 |
end |
|
|
16 |
|
|
|
17 |
%frequency=3, medium noise |
|
|
18 |
for i=1:300 |
|
|
19 |
ph = 2*pi*rand; |
|
|
20 |
y = sin(3*x + ph); |
|
|
21 |
outpure = [outpure;y]; |
|
|
22 |
outnoisy = [outnoisy;awgn(y,7)]; |
|
|
23 |
end |
|
|
24 |
|
|
|
25 |
%frequency=3, harsh noise |
|
|
26 |
for i=1:300 |
|
|
27 |
ph = 2*pi*rand; |
|
|
28 |
y = sin(3*x + ph); |
|
|
29 |
outpure = [outpure;y]; |
|
|
30 |
outnoisy = [outnoisy;awgn(y,1)]; |
|
|
31 |
end |
|
|
32 |
|
|
|
33 |
%f=5 |
|
|
34 |
for i=1:300 |
|
|
35 |
ph = 2*pi*rand; |
|
|
36 |
y = sin(5*x + ph); |
|
|
37 |
outpure = [outpure;y]; |
|
|
38 |
outnoisy = [outnoisy;awgn(y,13)]; |
|
|
39 |
end |
|
|
40 |
|
|
|
41 |
%frequency=5, medium noise |
|
|
42 |
for i=1:300 |
|
|
43 |
ph = 2*pi*rand; |
|
|
44 |
y = sin(5*x + ph); |
|
|
45 |
outpure = [outpure;y]; |
|
|
46 |
outnoisy = [outnoisy;awgn(y,5)]; |
|
|
47 |
end |
|
|
48 |
|
|
|
49 |
%frequency=5, harsh noise |
|
|
50 |
for i=1:300 |
|
|
51 |
ph = 2*pi*rand; |
|
|
52 |
y = sin(5*x + ph); |
|
|
53 |
outpure = [outpure;y]; |
|
|
54 |
outnoisy = [outnoisy;awgn(y,0.5)]; |
|
|
55 |
end |
|
|
56 |
|
|
|
57 |
%frequency=8, low noise |
|
|
58 |
for i=1:300 |
|
|
59 |
ph = 2*pi*rand; |
|
|
60 |
y = sin(8*x + ph); |
|
|
61 |
outpure = [outpure;y]; |
|
|
62 |
outnoisy = [outnoisy;awgn(y,13)]; |
|
|
63 |
end |
|
|
64 |
|
|
|
65 |
%frequency=8, medium noise |
|
|
66 |
for i=1:300 |
|
|
67 |
ph = 2*pi*rand; |
|
|
68 |
y = sin(8*x + ph); |
|
|
69 |
outpure = [outpure;y]; |
|
|
70 |
outnoisy = [outnoisy;awgn(y,5)]; |
|
|
71 |
end |
|
|
72 |
|
|
|
73 |
%frequency=8, harsh noise |
|
|
74 |
for i=1:300 |
|
|
75 |
ph = 2*pi*rand; |
|
|
76 |
y = sin(8*x + ph); |
|
|
77 |
outpure = [outpure;y]; |
|
|
78 |
outnoisy = [outnoisy;awgn(y,0.5)]; |
|
|
79 |
end |
|
|
80 |
|
|
|
81 |
%frequency=12, low noise |
|
|
82 |
for i=1:300 |
|
|
83 |
ph = 2*pi*rand; |
|
|
84 |
y = sin(12*x + ph); |
|
|
85 |
outpure = [outpure;y]; |
|
|
86 |
outnoisy = [outnoisy;awgn(y,13)]; |
|
|
87 |
end |
|
|
88 |
|
|
|
89 |
%frequency=12, medium noise |
|
|
90 |
for i=1:300 |
|
|
91 |
ph = 2*pi*rand; |
|
|
92 |
y = sin(12*x + ph); |
|
|
93 |
outpure = [outpure;y]; |
|
|
94 |
outnoisy = [outnoisy;awgn(y,5)]; |
|
|
95 |
end |
|
|
96 |
|
|
|
97 |
%frequency=12, harsh noise |
|
|
98 |
for i=1:300 |
|
|
99 |
ph = 2*pi*rand; |
|
|
100 |
y = sin(12*x + ph); |
|
|
101 |
outpure = [outpure;y]; |
|
|
102 |
outnoisy = [outnoisy;awgn(y,0.5)]; |
|
|
103 |
end |
|
|
104 |
|
|
|
105 |
|
|
|
106 |
xlswrite('pure.xlsx',outpure); |
|
|
107 |
xlswrite('noisy.xlsx',outnoisy); |