MATLAB编程与神经网络模拟实践

55、修改一段代码，使其初始条件为在位置 15 和 35 处 phi = 0.5，并使用 mesh() 函数绘图。假设代码用于模拟某个物理量 phi 在二维空间和时间上的演化，phi 是一个二维数组，第一维表示时间，第二维表示空间位置。代码中还涉及到时间步长 dt、空间步长 dx 和扩散系数 D 等参数。

可将代码中的初始条件部分修改为 `phi(1,15) = 0.5; phi(1,35) = 0.5;`，以下是完整修改后的代码：

```matlab
phi = zeros(1000, 1000);
dt = 0.001; % time from 0 to 1 sec
dx = 0.05; % space from 0 to 50
D = 0.5;
phi(1,15) = 0.5;
phi(1,35) = 0.5; % 修改后的初始条件
for t = 1:999
    for x = 2:999
        phi(t + 1, x) = phi(t,x) + D*dt/dx^2*(phi(t,x + 1)-2*phi(t,x)+phi(t,x - 1));
    end
end
mesh(phi)

56、加载同一只鸟的更长录音，并为多个主题创建一个马尔可夫模型。可从指定途径获取文件 zf_y89.wav 作为更长的录音。使用 wavread() 函数加载该文件。

可按以下步骤操作：

1. 从指定途径获取文件 zf_y89.wav ；
2. 使用 wavread() 函数加载该文件；
3. 基于加载的更长录音为多个主题创建马尔可夫模型。

57、给定一系列观测值，如何确定最可能的隐马尔可夫模型？

MATLAB软件中的统计工具箱提供了一些实现解决这些问题算法的例程，可使用 hmmtrain 函数，它使用鲍姆-韦尔奇算法通过最大似然法来估计转移和发射概率，但需要转移和发射概率矩阵的初始猜测值，调用格式为：

[esttr, estemit] = hmmtrain(seqs, trguess, emitguess)

58、在一个项目中，需要创建一个马尔可夫模型来近似电生理记录中的放电情况。具体步骤如下：首先，要估计内部状态序列，并使用它来估计转移矩阵和发射矩阵。找到文件RA - spike - times.mat，该文件包含区域RA电生理记录的一系列放电时间。将放电时间转换为瞬时放电率（即放电间隔的倒数），并确定三种状态（驱动、强直和爆发后）的合适阈值。利用这三种状态的阈值，确定每个时刻的状态，即在每1毫秒的时间切片内，通过将瞬时放电率与为每种状态设置的放电率阈值进行比较，估计系统的状态，计算出每个时间点的状态后，得到一个状态序列states。接着，将放电数据转换为点过程数据，这将构成输