MATLAB：其元素取决于其索引的矩阵

计算科学 matlab 矩阵 线性系统
2021-12-13 08:45:20

我正在尝试将功能

f(μ,ν)=iν−μ∑019Hμ−ν(7j)+δμ,ν ,

μ,ν=−3,−2,...2,3转换成一个 7x7 矩阵，其中是第一类 Hankel 函数，表示 Kronecker Delta。Hδ

在 MATLAB 上如何做到这一点？

1个回答

一个简单的方法：

F = eye(7); % gives you the Kronecker delta part
for mu = -3:3
  for nu = -3:3
    isum = 0; % temporary variable for the sum
    for j = 0:19
      isum = isum + besselh(mu-nu,7.0*j);
    end
    F(4+mu,4+nu) = F(4+mu,4+nu) + (1i)^(nu-mu) * isum;
  end
end

对于大于 7x7 的系统，明智的做法是预先计算的每个可能值的贝塞尔函数和并查找所需的值。μ−ν

其它你可能感兴趣的问题
上一篇高效的特征求解器 下一篇并行 FORTRAN 90 编程初学者的最佳书籍和笔记