Question

我一直在努力解决问题，所以如果有人可以提供任何建议或示例，我们将非常感激。使用Fortran90。

计划目的：

从我选择的数量中删除文件中的随机行。我能想到的最好方法是使用随机数来对应一个行号。

目前的作用：

每次生成新的随机数并将它们输出到单独的文件中。

问题：

（1）它不生成可能对应于行号的整数。（2）我不知道如何使用这些数字来删除文件中的行。

program random1
implicit none
integer :: i, seed, removed
real :: r
open (unit=10,file='random.dat')
removed=5

call init_random_seed() 
do i=1,removed
call random_number(r)
write(10,*) r
end Do

end program random1

subroutine init_random_seed()
        integer :: i,n,clock
        integer, dimension(:),allocatable :: seed

        call random_seed(size=n)
        allocate(seed(n))

        call system_clock(count=clock)

        seed=clock+37*(/(i-1,i=1,n)/)
        call random_seed(put=seed)

        deallocate(seed)
end subroutine

谢谢！

Answer 1

以下是答案的一些片段。首先是一些声明

integer :: num_lines ! number of lines in file
integer :: ix        ! loop index variable
real :: fraction     ! what fraction of lines are to be deleted
logical, dimension(:), allocatable :: lines_index
real, dimension(:), allocatable :: rands

现在有些可执行文件

read(*,*) num_lines  ! or figure it out some other way
read(*,*) fraction   ! likewise 

allocate(rands(num_lines)) ! no error checking
call random_number(rands)
allocate(lines_index(num_lines), source=rands<fraction)      ! no error checking

现在lines_index(ix)为false，您可以删除文件的第ix行。至于实际删除文件中的行，我建议你逐行读取文件，只将那些不被删除的行写出到另一个文件。像这样的东西可能会起作用

do ix = 1, num_lines
    read(infile,*) aline
    if(lines_index(ix)) write(outfile,*) aline
end do

请注意，我采取的方法并不保证将删除20％（或您设置fraction的任何行）行，只会删除最有可能被删除的行数。如果您想保证删除n行，请执行

之类的操作

integer :: num_lines ! number of lines in file
integer :: ix, jx    ! loop index variables
integer :: n         ! number of lines to delete
integer, dimension(:), allocatable :: lines_index    ! line numbers for deletion
real :: rand

read(*,*) n

allocate(del_ix(n))         

do ix = 1,n
    call random_number(rand)
    lines_index(ix) = 1.0+num_lines*rand   ! lines_index(ix) will be between 1 and num_lines
end do

这种方法并不能保证不会多次选择同一行进行删除，您必须编写一些代码来处理这种情况。然后继续：

do ix = 1, num_lines
    read(infile,*) aline
    if(any(lines_index==ix)) then
        ! do not write the line
    else
        write(outfile,*) aline
    end if
end do

使用随机数删除文件的随机行。 FORTRAN90

1 个答案: