GEFs/alkhairy2022.m at main · AnalyticModeling/GEFs · GitHub

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
addpath matlab_wrappers
pyenv(ExecutionMode="OutOfProcess")

function fig2_2022
    tiledlayout(2, 2, TileSpacing='tight')

    % first two cochleas
    c1 = Cochlea(Ap=[0.11], bp=[1], Bu=[7], CF0=1);
    [betas1, reals1, imags1] = c1.plot_wavenumber(betas=logspace(log10(0.5), log10(2), 10000), show=false);
    c2 = Cochlea(Ap=[0.055], bp=[1], Bu=[7], CF0=10);
    [betas2, reals2, imags2] = c2.plot_wavenumber(betas=logspace(log10(0.5), log10(2), 10000), show=false);

    nexttile
    semilogx(betas1, reals1)
    hold on
    semilogx(betas2, reals2, '--')
    ylabel('Re\{k\} (mm^{-1})')
    xlabel('\beta')
    yline(0, ':k')
    xline(1, ':k')

    nexttile
    semilogx(betas1, imags1)
    hold on
    semilogx(betas2, imags2, '--')
    ylabel('Im\{k\} (mm^{-1})')
    xlabel('\beta')
    yline(0, ':k')
    xline(1, ':k')

    % last two cochleas
    c1 = Cochlea(Ap=[0.11], bp=[1], Bu=[7], CF0=1);
    [betas1, reals1, imags1] = c1.plot_impedance(betas=logspace(log10(0.5), log10(2), 10000), show=false);
    c2 = Cochlea(Ap=[0.055], bp=[1], Bu=[7], CF0=10);
    [betas2, reals2, imags2] = c2.plot_impedance(betas=logspace(log10(0.5), log10(2), 10000), show=false);

    nexttile
    semilogx(betas1, reals1)
    hold on
    semilogx(betas2, reals2, '--')
    ylabel('Re\{Z\}/(2\pi\rho l)')
    xlabel('\beta')
    yline(0, ':k')
    xline(1, ':k')

    nexttile
    semilogx(betas1, imags1)
    hold on
    semilogx(betas2, imags2, '--')
    ylabel('Im\{Z\}/(2\pi\rho l)')
    xlabel('\beta')
    yline(0, ':k')
    xline(1, ':k')
end

fig2_2022