如何在Fortran中基于kx、ky向量构建含重复元素的x_p向量？

你的代码存在几个关键问题，导致无法生成目标数组：

1. 内层循环覆盖赋值：在给x_p赋值的循环中，每个行号i对应的第一列会被内层循环的所有kx(a)反复覆盖，最终所有行的第一列都变成kx(150)，第二列同理全是ky(150)，完全不符合每个kx元素重复150次的要求。
2. 输出变量错误：最后输出循环中使用了未更新的变量i，而非当前循环变量m，导致输出内容全部相同。

正确实现代码 #

以下是符合需求的完整代码：

program build_xp
    implicit none
    integer, parameter :: n = 150
    real :: kx(n), ky(n)
    real :: x_p(n*n, 2)
    integer :: i, j, row_idx

    ! 生成kx和ky向量
    do i = 1, n
        kx(i) = -2.0 + 4.0 * real(i-1)/real(n-1)
        ky(i) = -2.0 + 4.0 * real(i-1)/real(n-1)
    end do

    ! 构建二维数组x_p
    do i = 1, n  ! 遍历kx的每个元素
        do j = 1, n  ! 每个kx元素对应ky的所有元素
            row_idx = (i-1)*n + j
            x_p(row_idx, 1) = kx(i)  ! kx元素重复n次
            x_p(row_idx, 2) = ky(j)  ! ky序列重复n次
        end do
    end do

    ! 输出结果（可选，可根据需求调整输出格式）
    do row_idx = 1, n*n
        write(*, '(2F10.4)') x_p(row_idx, 1), x_p(row_idx, 2)
    end do

end program build_xp

代码说明 #

• 用n作为参数定义向量长度，方便后续修改；
• 通过两层循环计算行号row_idx = (i-1)*n + j，确保每个kx元素对应n行（搭配所有ky元素），第一列自然实现kx元素重复n次，第二列则是ky序列完整重复n次；
• 输出时使用当前循环变量row_idx，避免变量混淆；
• 使用real()将整数转换为实数，避免整数除法导致的精度丢失。

内容的提问来源于stack exchange，提问作者Celso Júnior