fortran将txt文件中的一列数据赋值给一维数组

yangpeinian

问题已经解决了，谢谢各位！

yangpeinian

necrohan 发表于 2020-5-4 08:53
传数据文件

-3.0326800e+00
  -3.0303801e+00
  -3.0281079e+00
  -3.0258633e+00
  -3.0236461e+00
  -3.0214563e+00
  -3.0192937e+00
  -3.0171582e+00
  -3.0150497e+00
  -3.0129681e+00
  -3.0109133e+00
  -3.0088851e+00
  -3.0068835e+00
  -3.0049084e+00
  -3.0029595e+00
  -3.0010369e+00
  -2.9991404e+00
  -2.9972699e+00
  -2.9954252e+00
  -2.9936063e+00
  -2.9918131e+00
  -2.9900453e+00
  -2.9883031e+00
  -2.9865861e+00
  -2.9848943e+00
  -2.9832276e+00
  -2.9815859e+00
  -2.9799691e+00
  -2.9783770e+00
  -2.9768095e+00
  -2.9752666e+00
  -2.9737480e+00
  -2.9722538e+00
  -2.9707838e+00
  -2.9693378e+00
  -2.9679158e+00
  -2.9665176e+00
  -2.9651431e+00
  -2.9637923e+00
  -2.9624650e+00
  -2.9611611e+00
  -2.9598805e+00
  -2.9586230e+00
  -2.9573886e+00
  -2.9561771e+00
  -2.9549885e+00
  -2.9538226e+00
  -2.9526793e+00
  -2.9515585e+00
  -2.9504601e+00
  -2.9493839e+00
  -2.9483299e+00
  -2.9472979e+00
  -2.9462879e+00
  -2.9452997e+00
  -2.9443332e+00
  -2.9433883e+00
  -2.9424648e+00
  -2.9415628e+00
  -2.9406819e+00
  -2.9398223e+00
  -2.9389836e+00
  -2.9381659e+00
  -2.9373690e+00
  -2.9365927e+00
  -2.9358370e+00
  -2.9351018e+00
  -2.9343870e+00
  -2.9336924e+00
  -2.9330179e+00
  -2.9323634e+00
  -2.9317289e+00
  -2.9311141e+00
  -2.9305190e+00
  -2.9299434e+00
  -2.9293873e+00
  -2.9288506e+00
  -2.9283331e+00
  -2.9278346e+00
  -2.9273552e+00
  -2.9268947e+00
  -2.9264529e+00
.
.
.
数据文件是这样一列7001行的数据

yangpeinian

vvt 发表于 2020-5-4 10:02
语法和逻辑上都是正常的。

fortran运行后结果显示的波形是这样的

yangpeinian

但是我想得到的数据波形要是这样的，不知道哪里出了问题

yangpeinian

vvt 发表于 2020-5-4 10:02
语法和逻辑上都是正常的。

好的，谢谢您！

yangpeinian

vvt 发表于 2020-5-4 11:17
你可以把数据文件打包上传上来~

感谢大哥，文件已上传，请指教。

yangpeinian

vvt 发表于 2020-5-4 15:09
我试了一下，读出来的数据画图如你给出的图片2一样。

那么这张图，是如何得到的呢？

首先表示感谢您的有效回复。这个是我把文件中的数据赋值给一个一维数组，然后让这个数组是声场计算中的声源。计算后提取声源处的波形，就得到是第一幅图的波形。图1和图2的起跳点大概一致，有没有可能是数据计算的有效位数问题？但是我定义的是real(8),精度应该够吧

yangpeinian

vvt 发表于 2020-5-4 15:09
我试了一下，读出来的数据画图如你给出的图片2一样。

那么这张图，是如何得到的呢？

所以可以肯定的是，数据赋值应该没有问题。但是在计算过程中可能有什么有效位数的问题。

yangpeinian

vvt 发表于 2020-5-4 17:09
这需要你提供计算过程的代码，相关逻辑对应的公式，才能判断。

麻烦请指教！不胜感激！

yangpeinian

vvt 发表于 2020-5-4 17:09
这需要你提供计算过程的代码，相关逻辑对应的公式，才能判断。

程序是在原来有限差分的基础上改的。之前的程序声源直接是一个声源函数，现在的是把声源波形的txt文件数据赋值，设为新的声源。程序可能不好看懂，模型就是包括轴向和径向两个方向的声场计算，需要得到不同时刻不同位置的跟声场相关的物理量。总之就是时间循环里嵌套了一个轴向网格循环，里面再嵌套了一个径向网格循环。