forked from davidjuliancaldwell/MATLAB_ECoG_code
-
Notifications
You must be signed in to change notification settings - Fork 0
/
shannon.m
40 lines (34 loc) · 875 Bytes
/
shannon.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
area=pi*0.23^2/4;
current=1e-3;
time=1e-3;
charge=current*time;
charge_density=charge/area;
figure;plot(log10(1e6*charge), log10(1e6*charge_density),'ro')
hold on
time=1e-4;
charge=current*time;
charge_density=charge/area;
plot(log10(1e6*charge), log10(1e6*charge_density),'ro')
time=1e-5;
charge=current*time;
charge_density=charge/area;
plot(log10(1e6*charge), log10(1e6*charge_density),'ro')
current=1e-2;
time=1e-3;
charge=current*time;
charge_density=charge/area;
plot(log10(1e6*charge), log10(1e6*charge_density),'ro')
time=1e-4;
charge=current*time;
charge_density=charge/area;
plot(log10(1e6*charge), log10(1e6*charge_density),'ro')
time=1e-5;
charge=current*time;
charge_density=charge/area;
plot(log10(1e6*charge), log10(1e6*charge_density),'ro')
k = 1.85;
d = [-1:0.1:1.5];
q = k - d;
plot(d,q);
xlabel('log 10 charge/phase')
ylabel('log 10 charge density')