Diff of /Code/datasetgeneration.m [000000] .. [a1f973]

Switch to unified view

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);