-
Notifications
You must be signed in to change notification settings - Fork 1
/
gwf_stimulated_echo_pte.m
88 lines (86 loc) · 2.5 KB
/
gwf_stimulated_echo_pte.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
function [gwf, rf, dt] = ctc_wf_stimulated_echo_pte
gwf = [
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.99896 -0.05943
-0.00000 0.99843 -0.06777
0.00000 0.99769 -0.07793
-0.00000 0.99666 -0.09017
-0.00000 0.99522 -0.10482
-0.00000 0.99322 -0.12240
0.00000 0.99037 -0.14364
-0.00000 0.76560 -0.16445
-0.00000 0.47148 -0.18552
0.00000 0.15806 -0.20730
-0.00000 -0.16435 -0.22954
0.00000 -0.48516 -0.25198
0.00000 -0.79386 -0.27438
-0.00000 -0.95929 -0.28498
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 -0.49989 -0.86693
-0.00000 -0.54445 -0.83967
0.00000 -0.58922 -0.80887
0.00000 -0.63522 -0.77327
-0.00000 -0.68353 -0.73092
0.00000 -0.73544 -0.67867
-0.00000 -0.79261 -0.61092
0.00000 -0.85730 -0.51623
0.00000 -0.93161 -0.36547
-0.00000 -0.99810 -0.07255
-0.00000 -0.96475 0.26594
0.00000 -0.89535 0.44700
0.00000 -0.82965 0.55959
-0.00000 -0.76886 0.64055
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
-0.00000 -0.50949 0.86132
0.00000 -0.45115 0.89327
-0.00000 -0.39265 0.92048
0.00000 -0.33322 0.94362
-0.00000 -0.27170 0.96314
0.00000 -0.20641 0.97921
0.00000 -0.13508 0.99157
0.00000 -0.05471 0.99923
0.00000 0.03821 1.00000
0.00000 0.14647 0.98995
-0.00000 0.26917 0.96385
0.00000 0.39931 0.91761
0.00000 0.52660 0.85097
-0.00000 0.64306 0.76677
-0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
0.00000 0.00000 0.00000
-0.00000 0.00000 0.00000
-0.00000 0.37106 0.92940
-0.00000 0.94780 -0.32114
0.00000 0.59884 -0.80178
-0.00000 0.53718 -0.84434
-0.00000 0.49613 -0.86909
0.00000 0.46617 -0.88552
0.00000 0.44320 -0.89724
-0.00000 0.42505 -0.90597
-0.00000 0.41047 -0.91268
0.00000 0.39863 -0.91791
0.00000 0.00000 0.00000
];
rf = ones(size(gwf,1),1);
rf([(19:36)+1 (65+1):end]) = -1;
rf([(37:46)+1]) = 0;
dt = 1e-3;
gwf = gwf.*rf * 0.08;