-
Notifications
You must be signed in to change notification settings - Fork 0
/
Oppgave3.asv
100 lines (81 loc) · 2.56 KB
/
Oppgave3.asv
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
98
99
100
%FIR 1
h1 = [0.0002, 0.0001, -0.0001, -0.0005, -0.0011, -0.0017, -0.0019, ...
-0.0016, -0.0005, 0.0015, 0.0040, 0.0064, 0.0079, 0.0075, 0.0046, ...
-0.0009, -0.0084, -0.0164, -0.0227, -0.0248, -0.0203, -0.0079, ...
0.0127, 0.0400, 0.0712, 0.1021, 0.1284, 0.1461, 0.1523, 0.1461, ...
0.1284, 0.1021, 0.0712, 0.0400, 0.0127, -0.0079, -0.0203, -0.0248, ...
-0.0227, -0.0164, -0.0084, -0.0009, 0.0046, 0.0075, 0.0079, 0.0064, ...
0.0040, 0.0015, -0.0005, -0.0016, -0.0019, -0.0017, -0.0011, ...
-0.0005, -0.0001, 0.0001, 0.0002];
%FIR 2
h2 = [-0.0002, -0.0001, 0.0003, 0.0005, -0.0001, -0.0009, -0.0007, ...
0.0007, 0.0018, 0.0005, -0.0021, -0.0027, 0.0004, 0.0042, 0.0031, ...
-0.0028, -0.0067, -0.0023, 0.0069, 0.0091, -0.0010, -0.0127, ...
-0.0100, 0.0077, 0.0198, 0.0075, -0.0193, -0.0272, 0.0014, 0.0386, ...
0.0338, -0.0246, -0.0771, -0.0384, 0.1128, 0.2929, 0.3734, 0.2929, ...
0.1128, -0.0384, -0.0771, -0.0246, 0.0338, 0.0386, 0.0014, -0.0272, ...
-0.0193, 0.0075, 0.0198, 0.0077, -0.0100, -0.0127, -0.0010, 0.0091, ...
0.0069, -0.0023, -0.0067, -0.0028, 0.0031, 0.0042, 0.0004, -0.0027, ...
-0.0021, 0.0005, 0.0018, 0.0007, -0.0007, -0.0009, -0.0001, 0.0005, ...
0.0003, -0.0001, -0.0002];
fs = 100;
%Frekvensspekter av de to filtrene
[X1,f1] = frekspekin3190(h1,1000,100);
[X2,f2] = frekspekin3190(h2,1000,100);
%figure
%imagesc(y1(1:200,1:150))
%colormap gray
%colorbar
%plot(t_array, y1(t > 0.5 & t < 1.5, 1))
[Xseis1,fseis1] = frekspekin3190(y1(t > 0.5 & t < 1.5,1),1000,fs);
%figure
%plot(fseis1,20*log(abs(Xseis1)))
L = 249;
t1 = tukeywin(L,1);
tukey_y1= y1(t > 0.5 & t < 1.5, 1).*t1;
[Xseis2,fseis2] = frekspekin3190(tukey_y1,1000,fs);
%figure
%plot(fseis1,20*log(abs(Xseis2)))
figure
imagesc(offset1,t*1000,y1)
figure
plot(t_array, y1(t > 0.9 & t < 1.5, 1))
%{
n = 10;
%[Xse,fse] = frekspekin3190(seismogram1,1000,fs)
seismog1 = seismogram1(1:n,1:n);
figure
plot(t(1:n),seismog1)
%plot(fse,abs(Xse))
title('Plott av h1(t) og h2(t)')
xlabel('Tid')
ylabel('Amplitude')
[Xseis1,fseis1] = frekspekin3190(y2,1000,fs);
figure
plot(fseis1,20*log10((abs(Xseis1))))
y1 = konvin3190(h1,seismogram1,0);
title('Offset vs Tid')
xlabel('Offset')
ylabel('Tid')
%plot(t,y1)
figure
imagesc(seismogram1(1:150,1:100 ))
colormap gray
colorbar
%set(gca, 'YDir','reverse')
%plot(t,y2)
%plot(t,seismogram1)
%{
figure
subplot(2,1,1);
plot(f1,20*log10(abs(X1)))
title('Frekvensspekter av h1')
xlabel('Frekvens')
ylabel('Amplitude')
subplot(2,1,2);
plot(f2,20*log10(abs(X2)))
title('Frekvensspekter av h2')
xlabel('Frekvens')
ylabel('Amplitude')
%}
%}