Code development platform for open source projects from the European Union institutions

Skip to content
Snippets Groups Projects
Commit a666dc0f authored by Markus Quaritsch's avatar Markus Quaritsch
Browse files

require at least 2 points for extrapolation

parent 713f73a2
No related branches found
No related tags found
No related merge requests found
......@@ -87,10 +87,16 @@ namespace TUGraz.VectoCore.InputData.Reader.ComponentData
return new TransmissionLossMap(entries, gearRatio, gearName);
}
var orig = "";
entries.ForEach(x =>orig +=string.Format("{0},{1},{2}" + Environment.NewLine, x.InputSpeed.AsRPM, x.InputTorque.Value(), x.TorqueLoss.Value()));
entries.ForEach(
x =>
orig +=
string.Format("{0},{1},{2}" + Environment.NewLine, x.InputSpeed.AsRPM, x.InputTorque.Value(), x.TorqueLoss.Value()));
entries = ExtendLossMap(entries);
var extended = "";
entries.ForEach(x => extended += string.Format("{0},{1},{2}" + Environment.NewLine, x.InputSpeed.AsRPM, x.InputTorque.Value(), x.TorqueLoss.Value()));
entries.ForEach(
x =>
extended +=
string.Format("{0},{1},{2}" + Environment.NewLine, x.InputSpeed.AsRPM, x.InputTorque.Value(), x.TorqueLoss.Value()));
return new TransmissionLossMap(entries, gearRatio, gearName);
}
......@@ -121,12 +127,13 @@ namespace TUGraz.VectoCore.InputData.Reader.ComponentData
}
var torqueStep = 500.SI<NewtonMeter>();
foreach (var speedBucket in speedBuckets) {
if (speedBucket.Value.Count == 0)
if (speedBucket.Value.Count < 2) {
continue;
}
double k, d, r;
VectoMath.LeastSquaresFitting(speedBucket.Value, x => x.InputTorque.Value(), x => x.TorqueLoss.Value(), out k, out d,
out r);
for (var inTq = speedBucket.Value.Max(x => x.InputTorque) + torqueStep; inTq <= 2 * maxTorque; inTq += torqueStep) {
entries.Add(new TransmissionLossMap.GearLossMapEntry(speedBucket.Key, inTq, k * inTq + d.SI<NewtonMeter>()));
entries.Add(new TransmissionLossMap.GearLossMapEntry(speedBucket.Key, -inTq, k * inTq + d.SI<NewtonMeter>()));
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment