帮我优化这段代码：f=0:0.00001:1; for k=1:100000 P(k)=(4*pi*F*Ps)/(beta*En)*exp(-4/beta*(pi*k*F)^2)*(cos(pi*k*F*detat))^2*dirac(f*10^6-k*F); end

可以尝试使用向量化来优化这段代码，避免使用 for 循环。具体的优化方式如下：

1. 将 f 和 k 转化成行向量和列向量，以便进行矩阵运算：

f = (0:0.00001:1)';
k = (1:100000)';

2. 将 dirac 函数替换为矩阵运算，使用 Kronecker delta 函数进行代替：

dirac = @(x) (abs(x) < eps);

3. 将 exp 函数中的常数项提取出来，避免重复计算：

c = -4 / beta * (pi * F)^2;

4. 将 cos 函数替换为矩阵运算，使用余弦函数进行代替：

cos2 = @(x) 0.5 * (cos(x) + 1);

5. 使用矩阵运算进行计算，避免使用 for 循环：

P = (4*pi*F*Ps)/(beta*En) * exp(c * (k * pi).^2) .* (cos2(pi*k*F*detat)).^2 .* kron(f*10^6-k*F, ones(size(f)));

这样可以大大提高代码的计算效率。