Mercurial > dive4elements > river
annotate flys-artifacts/src/main/java/de/intevation/flys/artifacts/math/fitting/SQPow.java @ 3133:6f6a8e5fb7a4
Theme updates for sector average wq curves/points
flys-artifacts/trunk@4740 c6561f87-3c4e-4783-a992-168aeb5c3f6f
author | Christian Lins <christian.lins@intevation.de> |
---|---|
date | Thu, 21 Jun 2012 12:26:08 +0000 |
parents | 6391ca82141c |
children | 585c0b49e696 |
rev | line source |
---|---|
3100
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
1 package de.intevation.flys.artifacts.math.fitting; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
2 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
3 public class SQPow |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
4 extends Function |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
5 { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
6 public static final Derivative DERIVATIVE = |
3104
6391ca82141c
SQ: Name derivative S'(Q).
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
3100
diff
changeset
|
7 new Derivative("S'(Q) = a*b*Q^(b-1)") { |
3100
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
8 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
9 @Override |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
10 public de.intevation.flys.artifacts.math.Function |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
11 instantiate(double [] parameters) |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
12 { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
13 final double a = parameters[0]; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
14 final double b = parameters[1]; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
15 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
16 return new de.intevation.flys.artifacts.math.Function() { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
17 @Override |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
18 public double value(double Q) { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
19 return a*b*Math.pow(Q, b-1); |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
20 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
21 }; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
22 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
23 }; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
24 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
25 public SQPow() { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
26 super( |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
27 "sq-pow", |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
28 "S(Q) = a*Q^b", |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
29 new String [] { "a", "b" }); |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
30 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
31 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
32 @Override |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
33 public double value(double x, double [] parameters) { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
34 return parameters[0]*Math.pow(x, parameters[1]); |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
35 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
36 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
37 @Override |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
38 public double [] gradient(double q, double [] parameters) { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
39 double a = parameters[0]; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
40 double b = parameters[1]; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
41 double q_b = Math.pow(q, b); |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
42 return new double [] { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
43 q_b, |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
44 a*q_b*Math.log(q), |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
45 }; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
46 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
47 |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
48 @Override |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
49 public Derivative getDerivative() { |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
50 return DERIVATIVE; |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
51 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
52 } |
c0abfb695332
SQ: Add new power function.
Sascha L. Teichmann <sascha.teichmann@intevation.de>
parents:
diff
changeset
|
53 // vim:set ts=4 sw=4 si et sta sts=4 fenc=utf8 : |