回 帖 发 新 帖 刷新版面

主题:[讨论]FORTRAN能通过求二导确定曲线拐点吗?

下面这组数共60行,为x,y两变量互相对应的数       
          x           y

        2868         507
        2863         623
        2858         878
        2853        1002
        2848        1166
        2843        1102
        2838        1183
        2833        1086
        2828        1191
        2823        1146
        2818        1367
        2813        1328
        2808        1535
        2803        1409
        2798        1600
        2793        1461
        2788        1644
        2783        1503
        2778        1671
        2773        1534
        2768        1716
        2763        1597
        2758        1775
        2753        1647
        2748        1795
        2743        1662
        2738        1807
        2733        1698
        2728        1844
        2723        1737
        2718        1870
        2713        1756
        2708        1877
        2703        1756
        2698        1896
        2693        1829
        2688        1969
        2683        1874
        2678        1981
        2673        1896
        2668        2018
        2663        1956
        2658        2084
        2653        2004
        2648        2121
        2643        2052
        2638        2154
        2633        2074
        2628        2184
        2623        2144
        2618        2269
        2613        2220
        2608        2318
        2603        2268
        2598        2370
        2593        2340
        2588        2438
        2583        2406
        2578        2480
        2573        2482

    想确定x-y对应曲线的拐点位置(即第几行)?试过MATLAB,拟合后圆滑再求二导,结果很不理想.特请教各位大侠,应该FORTRAN就能办了。
    有些人说可调用库函数,我不清楚,还有上面的数只是我文件中得一组,这样的数有9750组,所以还要读取和数组。谢谢!  [em18][em18]

回复列表 (共2个回复)

沙发

既然是离散的,咱就离散地处理呗。
读入时求相邻点的比值(割线斜率,近似一导)
读入完后再求一遍比值(二导)
然后就能知道啦:)

板凳

直接离散处理,不过只能根据你设置的误差,近似的求得拐点结果。一般来说,真正的拐点就在你的数值结果里头,不过数值结果不一定是正确的。

我来回复

您尚未登录,请登录后再回复。点此登录或注册