求非线性方程 y=ax^b,已知x,y, 求a和b,需要java实现方法

2020-09-14 02:58:36 587 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

$．．天空\.．$ ．．天空\.．墨衍会员 · AI 创作全网分发关注 + 加入社区

给定n，求1/x + 1/y = 1/n （x<=y）的解数。（x、y、n均为正整数）基本的思路首先化简可得 x+yx∗y=1n\frac{x+y}{x*y}=\frac{1}{n}关键一步令令 x=n+ax=n+a y=n+by=n+b这一步相当关键，整个题目的精髓所在。因为x和y和n都是正整数，要想等式成立，x和y作为分母肯定都要大于n，不然等式成立不了代入原公式化简得代入原公式化简得 n2=ab n^2=ab 然后就等价于求解n2n^2的约数阅读详情

求非线性方程 y=ax^b,已知x,y, 求a和b,需要java实现方法
x y
0.05 0.01
0.1 0.04
0.15 0.08
0.2 0.137
0.25 0.207
0.3 0.288
0.35 0.386
0.4 0.502
0.45 0.631
0.5 0.781

我知道结果, 但是不知道计算过程
a=2.7972, b=1.8712

求教..............

...全文

笔试题回顾：求抛物线y=Ax^2+x+B与直线x=C，直线x=D以及x轴所围成的封闭图形的面积求抛物线y=Ax^2+x+B与直线x=C，直线x=D以及x轴所围成的封闭图形的面积。题目描述求抛物线y=Ax^2+x+B与直线x=C，直线x=D以及x轴所围成的封闭图形的面积。输入: 第一行为一个整数T，表示测试数据的组数接下来每行输入四个整数A, B, C, D 1 <= T <= 1000, 1 <= A, B <= 100, -100 <= C <= 100 输出: 每组测试数据输出一个答案，为封闭图形的面积样例输入: 1 2 3 1 2 样例输出: 9 阅读详情

$．．天空\.．$

．．天空\.．

Java 51,407 成员 86,132 内容

+ 关注 + 加入社区

切换为时间正序

浪费七年时间 2020-09-14

打赏
举报

数学的东西都忘了（其实是不会），不过这种题网上搜一下就有答案了用java怎么拟合形如y=ax^b的幂函数

引用

对 y=ax^b 两边同时取 ln，变换成线性方程 ln(y)=ln(a)+b*ln(x)，最小二乘法拟合

ln(x)看作x'，ln(y)看作y'，b看作a'，ln(a)看作b'，线性方程是 y' = a'x' + b'，然后再找下Java最小二乘直线拟合的代码最小二乘法的Java实现然后自己再处理下数据


public class MyLineRegression {

    public static void main(String[] args) {
        double[] x = {0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5};
        double[] y = {0.01, 0.04, 0.08, 0.137, 0.207, 0.288, 0.386, 0.502, 0.631, 0.781};

        double[] lnx = Arrays.stream(x).map(Math::log).toArray();
        double[] lny = Arrays.stream(y).map(Math::log).toArray();

        Map<String, Double> result = new MyLineRegression().lineRegression(lnx, lny);

        System.out.println("a = " + Math.pow(Math.E, result.get("b").doubleValue()));
        System.out.println("b = " + result.get("a").doubleValue());
        System.out.println("r = " + result.get("r").doubleValue());
    }

    /**
     * 最小二乘法
     *
     * @param X
     * @param Y
     * @return y = ax + b, r
     */
    public Map<String, Double> lineRegression(double[] X, double[] Y) {
        if (null == X || null == Y || 0 == X.length
                || 0 == Y.length || X.length != Y.length) {
            throw new RuntimeException();
        }

        // x平方差和
        double Sxx = varianceSum(X);
        // y平方差和
        double Syy = varianceSum(Y);
        // xy协方差和
        double Sxy = covarianceSum(X, Y);

        double xAvg = arraySum(X) / X.length;
        double yAvg = arraySum(Y) / Y.length;

        double a = Sxy / Sxx;
        double b = yAvg - a * xAvg;

        // 相关系数
        double r = Sxy / Math.sqrt(Sxx * Syy);
        Map<String, Double> result = new HashMap<String, Double>();
        result.put("a", a);
        result.put("b", b);
        result.put("r", r);

        return result;
    }

    /**
     * 计算方差和
     *
     * @param X
     * @return
     */
    private double varianceSum(double[] X) {
        double xAvg = arraySum(X) / X.length;
        return arraySqSum(arrayMinus(X, xAvg));
    }

    /**
     * 计算协方差和
     *
     * @param X
     * @param Y
     * @return
     */
    private double covarianceSum(double[] X, double[] Y) {
        double xAvg = arraySum(X) / X.length;
        double yAvg = arraySum(Y) / Y.length;
        return arrayMulSum(arrayMinus(X, xAvg), arrayMinus(Y, yAvg));
    }

    /**
     * 数组减常数
     *
     * @param X
     * @param x
     * @return
     */
    private double[] arrayMinus(double[] X, double x) {
        int n = X.length;
        double[] result = new double[n];
        for (int i = 0; i < n; i++) {
            result[i] = X[i] - x;
        }

        return result;
    }

    /**
     * 数组求和
     *
     * @param X
     * @return
     */
    private double arraySum(double[] X) {
        double s = 0;
        for (double x : X) {
            s = s + x;
        }
        return s;
    }

    /**
     * 数组平方求和
     *
     * @param X
     * @return
     */
    private double arraySqSum(double[] X) {
        double s = 0;
        for (double x : X) {
            s = s + Math.pow(x, 2);
            ;
        }
        return s;
    }

    /**
     * 数组对应元素相乘求和
     *
     * @param X
     * @return
     */
    private double arrayMulSum(double[] X, double[] Y) {
        double s = 0;
        for (int i = 0; i < X.length; i++) {
            s = s + X[i] * Y[i];
        }
        return s;
    }

}

得到 a = 2.7971554678111508 b = 1.8711560608493885 r = 0.9998278572384243

求非线性方程 y=ax^b,已知x,y, 求a和b,需要java实现方法

相关推荐