Question

我正在尝试将matlab / gnu倍频程代码转换为gfortran语言代码。

我已经完成了struc的声明，以形成一个用matlab编写的数组。但是，由于该数组问题，我不理解，编译器不会处理SIZE命令。以下是我要翻译的八度代码

realloc

八度代码的输出为

function test1
clc;
close all;
start= 5;
function out = function1(start)
myarray = struct(...
    'list1', {'A', 'Cc', 'B', 'E', 'F', 'G', 'H', 'I'}, ...
    'list2', {1, 2,  3,  4, 5, 6, 7, 8},...
    'list3', {3, 3, 3, 3, 3, 2, 2, 2});
function list1 = function2(list2)
    done = false;
    myarray    
    numel(myarray)
    ra = randi([1 numel(myarray)], 1, 1)
    myarray(ra).list2;

    list1 = myarray(ra);
    if list1.list1 == 'E'
        list1.list1 = round(randn(1,1) * 5);
    end
   end
  x = {}; 
  y = [];
  list2 = 0;    
    list1 = function2(list2);
    list2 = list2 + list1.list3;
  out = struct('x', {x}, 'y', {y});
 end
function1(5)
end

我的gfortran代码是

test1

myarray =

1x8 struct array containing the fields:

list1
list2
list3

ans =  8
ra =  1
ans =

scalar structure containing the fields:

x = {}(0x0)
y = [](0x0)

我从build命令得到的错误消息是：

function function1(start) result(out1)

 type myarray
 character, dimension(8)::list1=(/'A','C','B','E','F','G','H','I'/); 
 integer, dimension(8)::list2 =(/1,2,3,4,5,6,7,8/); 
 integer, dimension(8)::list3 =(/3,3,3,3,3,2,2/);
  end type 

  integer :: start;
  integer :: out1;
  integer :: sd;

 contains

  function function2(list2) result(list1)

   integer, intent(in) :: list2; ! input
   integer :: list1;             ! output  
   logical :: done;  
   integer, dimension(1,1) :: ra;
   real    :: rnd1;

    done = .FALSE.


    call random_number(rnd1);
    ra = nint(rnd1*5);
    sd= size(myarray);
   ! ra = nint(rnd1*size(myarray));
       print*, "random number is ", ra;
     myarray(ra)

 end function function2

 end function function1

program test1
use iso_fortran_env
implicit none

integer                 :: start, xout;
real                    :: ra;
integer                 :: function1;

start =5;
xout=function1(5);
end program test1

我相信问题是我对数组结构的声明。我在那里丢失了一些东西。有什么建议么？请注意，在list1值的gfortran代码中，我必须将字符“ Cc”更改为“ C”以使其起作用。

Answer 1

与Matlab情况不同，struct创建具有给定值的对象，而Fortran中的type语句则没有。

相反，type myarray定义对象的外观。没有创建对象，我们需要做类似的事情

type myarray
   ...   ! The definition of components
end type myarray
type(myarray) obj ! Declare an object of the defined type.

此后，obj是您的关注对象。

但是，还有更多注意事项。与

type myarray
  character, dimension(8)::list1=(/'A','C','B','E','F','G','H','I'/); 
  integer, dimension(8)::list2 =(/1,2,3,4,5,6,7,8/); 
  integer, dimension(8)::list3 =(/3,3,3,3,3,2,2/);
end type

您不会（再次）创建对象。当您做创建一个对象（使用type(myarray) obj时，该对象很可能以指定的值开头。这与Matlab struct所期望的不太一样您可以阅读有关默认初始化和构造函数的详细信息。

来到size(myarray)，仅声明为type(myarray) obj的对象是 scalar 对象。这是“结构的数组”和“数组的结构”之间的区别。 “结构” myarray包含三个长度为8的数组。要改为具有结构数组，则：

type myarray
  character list1
  integer list2, list3
end type myarray
type(myarray), dimension(8) :: array

但是，您将必须构造数组。也许

array = [myarray('A',1,3), myarray('C',2,3), ...]

其他问题的答案也涉及构建此类结构数组的细节。

正确的struc数组声明

1 个答案: