1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
| type: custom:plotly-graph
title: COP vs Tlift – QUADRATIC fit, last 96h, COP_Tlift resample 60secs
hours_to_show: 96
refresh_interval: 0
no_default_entities: true
raw_plotly_config: true
entities:
- entity: sensor.heatpump_1_temperature_water_in
internal: true
filters:
- resample: 60s
- force_numeric: true
- store_var: tout
- entity: sensor.heatpump_1_temperature_outside
internal: true
filters:
- resample: 60s
- force_numeric: true
- store_var: tamb
- entity: sensor.heatpump_1_quatt_cop
internal: true
filters:
- resample: 60s
- force_numeric: true
- store_var: cop_raw
- entity: sensor.heatpump_1_compressor_frequency
internal: true
filters:
- resample: 60s
- force_numeric: true
- store_var: freq
- entity: ""
x: |-
$fn(({vars}) => {
const n = Math.min(vars.tout?.xs?.length||0, vars.tamb?.xs?.length||0, vars.cop_raw?.xs?.length||0, vars.freq?.xs?.length||0);
const x = [], y = [], colors = [], custom = [];
for (let i = 0; i < n; i++) {
const timestamp = (vars.tout.xs[i] || 0);
const d = new Date(timestamp);
// Bulletproof Dutch date/time – always correct
const day = String(d.getDate()).padStart(2,'0');
const month = String(d.getMonth()+1).padStart(2,'0');
const hour = String(d.getHours()).padStart(2,'0');
const min = String(d.getMinutes()).padStart(2,'0');
const datetime = `${day}-${month} ${hour}:${min}`;
const tlift = (vars.tout.ys[i] ?? NaN) - (vars.tamb.ys[i] ?? NaN);
const cop = vars.cop_raw.ys[i] ?? NaN;
const freq = vars.freq.ys[i] ?? 0;
const tamb = (vars.tamb.ys[i] ?? NaN).toFixed(1);
const tout = (vars.tout.ys[i] ?? NaN).toFixed(1);
if (!isNaN(tlift) && !isNaN(cop) && cop > 1.5 && cop < 15 && freq > 20) {
x.push(tlift);
y.push(cop);
colors.push(cop);
custom.push([datetime, tlift.toFixed(1), cop.toFixed(2), freq.toFixed(0), tamb, tout]);
}
}
this.xdata = x;
this.ydata = y;
this.colors = colors;
this.custom = custom;
return x;
})
"y": $fn(() => this.ydata || [])
type: scattergl
mode: markers
name: COP_Tlift:96h/60s
marker:
color: $fn(() => this.colors || [])
colorscale: Rainbow
size: 5
opacity: 0.82
line:
width: 1
color: black
symbol: circle-open
colorbar:
title: COP
titleside: right
x: 1.02
thickness: 18
hovertemplate: >
<b>%{customdata[0]}</b><br> Tlift: <b>%{customdata[1]} °C</b><br>
COP: <b>%{customdata[2]}</b><br> Hz:
<b>%{customdata[3]}</b><br> T_ambient: <b>%{customdata[4]} °C</b><br>
T_water: <b>%{customdata[5]} °C</b><extra></extra>
customdata: $fn(() => this.custom || [])
- entity: ""
x: |-
$fn(() => {
if (!window.qfit) return Array.from({length:100},(_,i)=>i*0.7);
const {xmin,xmax} = window.qfit;
return Array.from({length:100},(_,i)=>xmin + i*(xmax-xmin)/99);
})
"y": |-
$fn(({vars}) => {
const n = Math.min(vars.tout?.xs?.length||0, vars.tamb?.xs?.length||0, vars.cop_raw?.xs?.length||0, vars.freq?.xs?.length||0);
let sx=0, sx2=0, sx3=0, sx4=0, sy=0, sxy=0, sx2y=0, nvalid=0, sf=0, high=0;
const points = [];
for(let i=0;i<n;i++){
const x = (vars.tout.ys[i]??NaN) - (vars.tamb.ys[i]??NaN);
const y = vars.cop_raw.ys[i]??NaN;
const f = vars.freq.ys[i]??0;
if(!isNaN(x)&&!isNaN(y)&&y>1.5&&y<15&&f>20){
sx+=x; sx2+=x*x; sx3+=x*x*x; sx4+=x*x*x*x;
sy+=y; sxy+=x*y; sx2y+=x*x*y;
sf+=f;
if(f>60) high++;
nvalid++;
points.push({x,y});
}
}
if(nvalid<30) {
window.qfit = null;
return Array(100).fill(5);
}
// Normal equations for ax² + bx + c
const N = nvalid;
const A = [
[sx4, sx3, sx2],
[sx3, sx2, sx ],
[sx2, sx , N ]
];
const B = [sx2y, sxy, sy];
// Tiny 3×3 Gauss elimination
const mat = A.map((row,i)=>[...row, B[i]]);
for(let i=0;i<3;i++){
let pivot = mat[i][i];
for(let j=i+1;j<3;j++){
const f = mat[j][i]/pivot;
for(let k=i;k<4;k++) mat[j][k] -= f*mat[i][k];
}
}
const sol = [0,0,0];
for(let i=2;i>=0;i--){
let sum = mat[i].slice(i+1,3).reduce((a,v,j)=>a+v*sol[i+1+j],0);
sol[i] = (mat[i][3] - sum)/mat[i][i];
}
const [a,b,c] = sol;
// R²
const meanY = sy/N;
let ss_res=0, ss_tot=0;
points.forEach(p=>{
const pred = a*p.x*p.x + b*p.x + c;
ss_res += (p.y-pred)**2;
ss_tot += (p.y-meanY)**2;
});
const r2 = ss_tot>0 ? 1 - ss_res/ss_tot : 0;
const tlifts = points.map(p=>p.x);
const xmin = Math.max(0, Math.min(...tlifts));
const xmax = Math.min(70, Math.max(...tlifts));
const ymin = Math.max(1.5, Math.min(...points.map(p=>p.y)));
const ymax = Math.min(12, Math.max(...points.map(p=>p.y)));
window.qfit = {a,b,c,r2,meanF:sf/N,xmin,xmax,ymin,ymax};
console.clear();
console.log('%cQUADRATIC FIT – WINTER 2025','font-weight:bold;color:#c00');
console.log(`Points: ${nvalid} | Tlift ${xmin.toFixed(1)}–${xmax.toFixed(1)} °C`);
console.log(`COP = ${a<0?a.toFixed(6):a.toFixed(6)}·Tlift² ${b>=0?'+' : ''}${b.toFixed(4)}·Tlift + ${c.toFixed(3)}`);
console.log(`R² = ${r2.toFixed(4)} | Avg Hz: ${(sf/N).toFixed(1)}`);
// Return smooth curve
return Array.from({length:100},(_,i)=>{
const x = xmin + i*(xmax-xmin)/99;
return a*x*x + b*x + c;
});
})
type: scattergl
mode: lines
name: Quadratic Fit
line:
color: red
width: 4
- entity: ""
x: "$fn(() => window.qfit ? [window.qfit.xmax * 0.8] : [50])"
"y": "$fn(() => window.qfit ? [window.qfit.ymax * 0.95] : [9])"
text: |-
$fn(() => {
if (!window.qfit) return "calculating…";
const {a,b,c,r2,meanF} = window.qfit;
return `COP = ${a<0?'':''}${a.toFixed(6)}·Tlift² ${b>=0?'+':''}${b.toFixed(4)}·Tlift + ${c.toFixed(2)}
R² = ${r2.toFixed(4)}
Avg Hz: ${meanF.toFixed(0)}`;
})
type: scattergl
mode: text
textfont:
size: 15
color: red
family: Arial Black
showlegend: false
layout:
height: 650
width: 1050
plot_bgcolor: white
paper_bgcolor: white
xaxis:
title: Tlift (°C)
range: "$fn(() => window.qfit ? [window.qfit.xmin, window.qfit.xmax] : [0,60])"
yaxis:
title: COP
range: "$fn(() => window.qfit ? [window.qfit.ymin, window.qfit.ymax] : [1.5,10])"
hovermode: closest
legend:
x: 0.01
"y": 0.99
bgcolor: rgba(255,255,255,0.8)
uirevision: quadratic-winter-final-v2 |
/u/843743/crop653fae5274538_cropped.png?f=community)
:strip_exif()/f/image/5u5ze6uC8RZvJegikCxcyuDW.jpg?f=fotoalbum_large)
:no_upscale():strip_icc():strip_exif()/f/image/5u5ze6uC8RZvJegikCxcyuDW.jpg?f=user_large)
:no_upscale():strip_icc():strip_exif()/f/image/78rU2mPO8eOL8ZCVJMmIrq0b.jpg?f=user_large)
En gaat de "besparing" van de isolatie weer op aan warmte in andere vertrekken. 
:no_upscale():strip_icc():strip_exif()/f/image/u5MvCgJs4PimrS36yHpnQPj1.jpg?f=user_large)
:no_upscale():strip_icc():strip_exif()/f/image/rMPYtGYWCFlNW9vXsONE7TW3.jpg?f=user_large)
:strip_icc():strip_exif()/u/1872142/crop6658870902644_cropped.jpg?f=community)
:strip_icc():strip_exif()/u/1353620/crop68ffebed6da13_cropped.jpg?f=community)
:strip_exif()/f/image/omNm80k9QVoHBG2LDZeYSYQS.png?f=user_large)
:strip_icc():strip_exif()/u/282325/crop6159e3e38b9a8_cropped.jpg?f=community)
:strip_exif()/u/281030/crop601496411ce87.gif?f=community)
:no_upscale():strip_icc():strip_exif()/f/image/tXpTj23ap0nS8kvzenDcTIHl.jpg?f=user_large)
:no_upscale():strip_icc():strip_exif()/f/image/6ZdaRv0IG3PURd1RblO3J6BV.jpg?f=user_large)
/u/1205612/crop5f90439f46a56_cropped.png?f=community)
):no_upscale():strip_icc():strip_exif()/f/image/aHjIZCjW5o7kkkGWoGcH7ioY.jpg?f=user_large)
:no_upscale():strip_icc():strip_exif()/f/image/3udOLwyNIgQPatmUjVjJqY7h.jpg?f=user_large)
:strip_exif()/f/image/r3P7joCVeTfGhqiJAXjRkrxR.png?f=user_large)
:strip_exif()/f/image/xr3WKCS11hk8sPPbNkwuhyXx.png?f=user_large)
:no_upscale():strip_icc():strip_exif()/f/image/aH1T11BqhDy7QSigcXH3s4g5.jpg?f=user_large)
:strip_exif()/u/33665/GTicon.gif?f=community)
:strip_icc():strip_exif()/u/559112/crop648d83b854457_cropped.jpg?f=community)
:strip_exif()/f/image/nKnCzREnnCGHvJpbbL7hPBR9.png?f=user_large)
:strip_icc():strip_exif()/u/389729/crop589c305016610_cropped.jpeg?f=community)
:strip_exif()/f/image/avDWVnz2wFDojImlUPJ598Ot.png?f=user_large)
:strip_exif()/f/image/XAq01592m0NHdRe52dG8d6uR.png?f=user_large)
:strip_exif()/f/image/FnPpXr9Z2RorxnZNVAdrLYw9.png?f=user_large)
:strip_exif()/f/image/aH1sWp5SUWpim9yxCBCHgIu2.png?f=user_large)
:strip_exif()/f/image/osECPBbKHVSKqvJFinl3QB0b.png?f=user_large)
/u/511600/crop685e83d8663bd_cropped.png?f=community)
:strip_exif()/f/image/6DJoswqmqHljtWc4oBq7Zmnf.png?f=user_large)
:strip_exif()/f/image/QH9c7r1bB6ZI9mTU3s0YyE3X.png?f=user_large)
:strip_exif()/u/396930/crop593905cbcd58f_cropped.gif?f=community)
:strip_exif()/f/image/ueGxhV7aVeQXPUX37FIn8CR4.png?f=user_large)
:no_upscale():strip_icc():strip_exif()/f/image/CUWNB0oe8D6gFcAfILeyZx6q.jpg?f=user_large)
:strip_exif()/f/image/JAc0eq5sQWswDsGdNprVFXik.png?f=user_large)
:strip_exif()/f/image/GQwnGMeUJUve0XIzSntCJNrS.png?f=user_large)
:strip_icc():strip_exif()/u/36378/crop5a3f931ecfec0_cropped.jpeg?f=community)
:strip_exif()/f/image/zAxbrSnyP2KdzO4tO9PxviyB.png?f=user_large)