[interp1] problème NaN

**Newenda** · 14/09/2013, 19h28

Bonjour,

J'essaye de réaliser une interpolation simple :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

yi = interp1(x,Y,xi)

où

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
Y =
 
    0.0300    0.0600    0.0300
    0.0300    0.0300    0.0300
   -0.0300    0.0300    0.1200
    0.0600    0.1201    0.0300
    0.1259    0.1798    0.0599
    0.1801    0.2402    0.1801
    0.1259    0.1199    0.2398
    0.2402    0.4203    0.3002
    0.2462    0.3002    0.4803
    0.3656    0.6593    0.7792
    0.5524    1.6811    2.1614
    0.9171    2.1579    2.5774
    0.8466    3.6024    4.2028
    0.7805    2.4016    1.9213
    0.2997    0.5994    0.4196
    0.0600    0.1201    0.0300
    0.0300    0.1291    0.1937

et xi un vecteur de dimension 1*760 croissante.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
xi = 
 
  Columns 1 through 6
 
         0    0.0003    0.0007    0.0010    0.0014    0.0017
 
  Columns 7 through 12
 
    0.0021    0.0024    0.0028    0.0031    0.0035    0.0038
 
  Columns 13 through 18
 
    0.0042    0.0045    0.0049    0.0052    0.0056    0.0059
 
  Columns 19 through 24
 
    0.0063    0.0066    0.0070    0.0073    0.0077    0.0080
 
  Columns 25 through 30
 
    0.0084    0.0087    0.0091    0.0094    0.0098    0.0101
 
  Columns 31 through 36
 
    0.0105    0.0108    0.0112    0.0115    0.0119    0.0122
 
  Columns 37 through 42
 
    0.0126    0.0129    0.0133    0.0136    0.0140    0.0143
 
  Columns 43 through 48
 
    0.0147    0.0150    0.0154    0.0157    0.0160    0.0164
 
  Columns 49 through 54
 
    0.0167    0.0171    0.0174    0.0178    0.0181    0.0185
 
  Columns 55 through 60
 
    0.0188    0.0192    0.0195    0.0199    0.0202    0.0206
 
  Columns 61 through 66
 
    0.0209    0.0213    0.0216    0.0220    0.0223    0.0227
 
  Columns 67 through 72
 
    0.0230    0.0234    0.0237    0.0241    0.0244    0.0248
 
  Columns 73 through 78
 
    0.0251    0.0255    0.0258    0.0262    0.0265    0.0269
 
  Columns 79 through 84
 
    0.0272    0.0276    0.0279    0.0283    0.0286    0.0289
 
  Columns 85 through 90
 
    0.0293    0.0296    0.0300    0.0303    0.0307    0.0310
 
  Columns 91 through 96
 
    0.0314    0.0317    0.0321    0.0324    0.0328    0.0331
 
  Columns 97 through 102
 
    0.0335    0.0338    0.0342    0.0345    0.0349    0.0352
 
  Columns 103 through 108
 
    0.0356    0.0359    0.0363    0.0366    0.0370    0.0373
 
  Columns 109 through 114
 
    0.0377    0.0380    0.0384    0.0387    0.0391    0.0394
 
  Columns 115 through 120
 
    0.0398    0.0401    0.0405    0.0408    0.0412    0.0415
 
  Columns 121 through 126
 
    0.0419    0.0422    0.0426    0.0429    0.0433    0.0436
 
  Columns 127 through 132
 
    0.0440    0.0443    0.0447    0.0450    0.0454    0.0457
 
  Columns 133 through 138
 
    0.0461    0.0464    0.0468    0.0471    0.0475    0.0478
 
  Columns 139 through 144
 
    0.0482    0.0485    0.0488    0.0492    0.0496    0.0499
 
  Columns 145 through 150
 
    0.0503    0.0506    0.0510    0.0513    0.0516    0.0520
 
  Columns 151 through 156
 
    0.0523    0.0527    0.0530    0.0534    0.0537    0.0541
 
  Columns 157 through 162
 
    0.0544    0.0548    0.0551    0.0555    0.0558    0.0562
 
  Columns 163 through 168
 
    0.0565    0.0569    0.0572    0.0576    0.0579    0.0583
 
  Columns 169 through 174
 
    0.0586    0.0590    0.0593    0.0597    0.0600    0.0604
 
  Columns 175 through 180
 
    0.0607    0.0611    0.0614    0.0618    0.0621    0.0625
 
  Columns 181 through 186
 
    0.0628    0.0632    0.0635    0.0639    0.0642    0.0646
 
  Columns 187 through 192
 
    0.0649    0.0652    0.0656    0.0659    0.0663    0.0666
 
  Columns 193 through 198
 
    0.0670    0.0673    0.0677    0.0680    0.0684    0.0687
 
  Columns 199 through 204
 
    0.0691    0.0694    0.0698    0.0701    0.0705    0.0708
 
  Columns 205 through 210
 
    0.0712    0.0715    0.0719    0.0722    0.0726    0.0729
 
  Columns 211 through 216
 
    0.0732    0.0736    0.0739    0.0743    0.0746    0.0750
 
  Columns 217 through 222
 
    0.0753    0.0757    0.0760    0.0764    0.0767    0.0771
 
  Columns 223 through 228
 
    0.0774    0.0778    0.0781    0.0785    0.0788    0.0792
 
  Columns 229 through 234
 
    0.0795    0.0799    0.0802    0.0806    0.0809    0.0813
 
  Columns 235 through 240
 
    0.0816    0.0820    0.0823    0.0827    0.0830    0.0834
 
  Columns 241 through 246
 
    0.0837    0.0841    0.0844    0.0848    0.0851    0.0855
 
  Columns 247 through 252
 
    0.0858    0.0862    0.0865    0.0869    0.0872    0.0876
 
  Columns 253 through 258
 
    0.0879    0.0882    0.0886    0.0889    0.0893    0.0896
 
  Columns 259 through 264
 
    0.0900    0.0903    0.0907    0.0910    0.0914    0.0917
 
  Columns 265 through 270
 
    0.0921    0.0924    0.0928    0.0931    0.0935    0.0938
 
  Columns 271 through 276
 
    0.0942    0.0945    0.0949    0.0952    0.0956    0.0959
 
  Columns 277 through 282
 
    0.0963    0.0966    0.0970    0.0973    0.0977    0.0980
 
  Columns 283 through 288
 
    0.0984    0.0987    0.0991    0.0994    0.0997    0.1001
 
  Columns 289 through 294
 
    0.1004    0.1008    0.1011    0.1015    0.1018    0.1022
 
  Columns 295 through 300
 
    0.1025    0.1029    0.1032    0.1036    0.1039    0.1043
 
  Columns 301 through 306
 
    0.1046    0.1050    0.1053    0.1057    0.1060    0.1064
 
  Columns 307 through 312
 
    0.1067    0.1071    0.1074    0.1078    0.1081    0.1085
 
  Columns 313 through 318
 
    0.1088    0.1092    0.1095    0.1098    0.1102    0.1105
 
  Columns 319 through 324
 
    0.1109    0.1112    0.1116    0.1119    0.1123    0.1126
 
  Columns 325 through 330
 
    0.1130    0.1134    0.1137    0.1140    0.1144    0.1147
 
  Columns 331 through 336
 
    0.1151    0.1154    0.1158    0.1161    0.1165    0.1168
 
  Columns 337 through 342
 
    0.1172    0.1175    0.1179    0.1182    0.1198    0.1199
 
  Columns 343 through 348
 
    0.1200    0.1201    0.1202    0.1203    0.1204    0.1206
 
  Columns 349 through 354
 
    0.1207    0.1208    0.1209    0.1210    0.1212    0.1213
 
  Columns 355 through 360
 
    0.1214    0.1215    0.1216    0.1218    0.1219    0.1220
 
  Columns 361 through 366
 
    0.1221    0.1222    0.1224    0.1225    0.1226    0.1227
 
  Columns 367 through 372
 
    0.1228    0.1229    0.1231    0.1232    0.1233    0.1234
 
  Columns 373 through 378
 
    0.1235    0.1237    0.1238    0.1239    0.1240    0.1241
 
  Columns 379 through 384
 
    0.1242    0.1244    0.1245    0.1246    0.1247    0.1248
 
  Columns 385 through 390
 
    0.1250    0.1251    0.1252    0.1253    0.1254    0.1255
 
  Columns 391 through 396
 
    0.1257    0.1258    0.1259    0.1260    0.1261    0.1263
 
  Columns 397 through 402
 
    0.1264    0.1265    0.1266    0.1267    0.1268    0.1270
 
  Columns 403 through 408
 
    0.1271    0.1272    0.1273    0.1274    0.1276    0.1277
 
  Columns 409 through 414
 
    0.1278    0.1279    0.1280    0.1282    0.1283    0.1284
 
  Columns 415 through 420
 
    0.1285    0.1286    0.1288    0.1289    0.1290    0.1291
 
  Columns 421 through 426
 
    0.1292    0.1293    0.1295    0.1296    0.1297    0.1298
 
  Columns 427 through 432
 
    0.1299    0.1301    0.1302    0.1303    0.1304    0.1305
 
  Columns 433 through 438
 
    0.1307    0.1308    0.1309    0.1310    0.1311    0.1312
 
  Columns 439 through 444
 
    0.1314    0.1315    0.1316    0.1317    0.1318    0.1320
 
  Columns 445 through 450
 
    0.1321    0.1322    0.1323    0.1324    0.1325    0.1327
 
  Columns 451 through 456
 
    0.1328    0.1329    0.1330    0.1331    0.1333    0.1334
 
  Columns 457 through 462
 
    0.1335    0.1336    0.1337    0.1338    0.1340    0.1341
 
  Columns 463 through 468
 
    0.1342    0.1343    0.1344    0.1346    0.1347    0.1348
 
  Columns 469 through 474
 
    0.1349    0.1350    0.1351    0.1353    0.1354    0.1355
 
  Columns 475 through 480
 
    0.1356    0.1358    0.1359    0.1360    0.1361    0.1362
 
  Columns 481 through 486
 
    0.1363    0.1365    0.1366    0.1367    0.1368    0.1369
 
  Columns 487 through 492
 
    0.1371    0.1372    0.1373    0.1374    0.1375    0.1376
 
  Columns 493 through 498
 
    0.1378    0.1379    0.1380    0.1381    0.1382    0.1384
 
  Columns 499 through 504
 
    0.1385    0.1386    0.1387    0.1388    0.1389    0.1391
 
  Columns 505 through 510
 
    0.1392    0.1393    0.1394    0.1395    0.1397    0.1398
 
  Columns 511 through 516
 
    0.1399    0.1400    0.1401    0.1402    0.1404    0.1405
 
  Columns 517 through 522
 
    0.1406    0.1407    0.1408    0.1410    0.1411    0.1412
 
  Columns 523 through 528
 
    0.1413    0.1414    0.1416    0.1417    0.1418    0.1419
 
  Columns 529 through 534
 
    0.1420    0.1421    0.1423    0.1424    0.1425    0.1426
 
  Columns 535 through 540
 
    0.1427    0.1429    0.1430    0.1431    0.1432    0.1433
 
  Columns 541 through 546
 
    0.1435    0.1436    0.1437    0.1438    0.1439    0.1440
 
  Columns 547 through 552
 
    0.1442    0.1443    0.1444    0.1445    0.1446    0.1448
 
  Columns 553 through 558
 
    0.1449    0.1450    0.1451    0.1452    0.1453    0.1455
 
  Columns 559 through 564
 
    0.1456    0.1457    0.1458    0.1459    0.1461    0.1462
 
  Columns 565 through 570
 
    0.1463    0.1464    0.1465    0.1466    0.1468    0.1469
 
  Columns 571 through 576
 
    0.1470    0.1471    0.1473    0.1474    0.1475    0.1476
 
  Columns 577 through 582
 
    0.1477    0.1478    0.1480    0.1481    0.1482    0.1483
 
  Columns 583 through 588
 
    0.1484    0.1486    0.1487    0.1488    0.1489    0.1490
 
  Columns 589 through 594
 
    0.1491    0.1493    0.1494    0.1495    0.1496    0.1497
 
  Columns 595 through 600
 
    0.1499    0.1500    0.1501    0.1502    0.1503    0.1505
 
  Columns 601 through 606
 
    0.1506    0.1507    0.1508    0.1509    0.1510    0.1512
 
  Columns 607 through 612
 
    0.1513    0.1514    0.1515    0.1516    0.1518    0.1519
 
  Columns 613 through 618
 
    0.1520    0.1521    0.1522    0.1523    0.1525    0.1526
 
  Columns 619 through 624
 
    0.1527    0.1528    0.1530    0.1531    0.1532    0.1533
 
  Columns 625 through 630
 
    0.1534    0.1535    0.1537    0.1538    0.1539    0.1540
 
  Columns 631 through 636
 
    0.1541    0.1543    0.1544    0.1545    0.1546    0.1547
 
  Columns 637 through 642
 
    0.1548    0.1550    0.1551    0.1552    0.1553    0.1554
 
  Columns 643 through 648
 
    0.1556    0.1557    0.1558    0.1559    0.1560    0.1561
 
  Columns 649 through 654
 
    0.1563    0.1564    0.1565    0.1566    0.1567    0.1569
 
  Columns 655 through 660
 
    0.1570    0.1571    0.1572    0.1573    0.1575    0.1576
 
  Columns 661 through 666
 
    0.1577    0.1578    0.1579    0.1580    0.1582    0.1583
 
  Columns 667 through 672
 
    0.1584    0.1585    0.1586    0.1588    0.1589    0.1590
 
  Columns 673 through 678
 
    0.1591    0.1592    0.1593    0.1595    0.1596    0.1597
 
  Columns 679 through 684
 
    0.1598    0.1599    0.1601    0.1602    0.1603    0.1604
 
  Columns 685 through 690
 
    0.1605    0.1607    0.1608    0.1609    0.1610    0.1611
 
  Columns 691 through 696
 
    0.1612    0.1614    0.1615    0.1616    0.1617    0.1618
 
  Columns 697 through 702
 
    0.1620    0.1621    0.1622    0.1623    0.1624    0.1625
 
  Columns 703 through 708
 
    0.1627    0.1628    0.1629    0.1630    0.1631    0.1633
 
  Columns 709 through 714
 
    0.1634    0.1635    0.1636    0.1637    0.1638    0.1640
 
  Columns 715 through 720
 
    0.1641    0.1642    0.1643    0.1644    0.1646    0.1647
 
  Columns 721 through 726
 
    0.1648    0.1649    0.1650    0.1651    0.1653    0.1654
 
  Columns 727 through 732
 
    0.1655    0.1656    0.1657    0.1659    0.1660    0.1661
 
  Columns 733 through 738
 
    0.1662    0.1663    0.1664    0.1666    0.1667    0.1668
 
  Columns 739 through 744
 
    0.1669    0.1670    0.1672    0.1673    0.1674    0.1675
 
  Columns 745 through 750
 
    0.1676    0.1678    0.1679    0.1680    0.1681    0.1682
 
  Columns 751 through 756
 
    0.1683    0.1685    0.1686    0.1687    0.1688    0.1689
 
  Columns 757 through 760
 
    0.1691    0.1692    0.1693    0.1694

Mais la matrice yi ne contient que des NaN.

Sauriez vous pourquoi? j'ai essayé avec toutes les options d'interp1 possibles.

Merci d'avance

**Jean Dumoncel** · 14/09/2013, 20h18

Bonjour,

c'est normal d'obtenir des nan pour les valeurs de xi définies en dehors de l'intervalle des valeurs de x, puisque dans ce cas, c'est une extrapolation et non plus une interpolation. Tu peux essayer d'ajouter l'option 'extrap' pour faire l'approximation de ces valeurs. Mais à moins que la méthode d'interpolation utilisée colle bien à tes valeurs, plus tu t'éloignes des valeurs de x, moins les résultats seront "justes".

**FLB** · 15/09/2013, 03h20

Salut,
pour mieux gérer la méthode utilisée pour effectuer ton extrapolation, je te conseil de passer par un polyfit ou un spline.
Une fois que tu as clairement identifié ta courbe grâce à l'échantillon tu peux chercher les valeurs estimées aux abscisses qui t'intéressent.

**Newenda** · 15/09/2013, 11h14

Merci pour vos réponses,

Je suis d'accord d'obtenir des NaN en yi pour des valeurs de xi correspondant en dehors de x. Mais il n'y a pas de raison que le calcul donne des NaN pour l'intervalle qui va bien...alors que dans ce cas, il y a des NaN partout.

**Jean Dumoncel** · 15/09/2013, 14h14

J'ai fait un test avec ton code, et j'obtiens des résultats pour les valeurs de xi qui sont entre 0.1249 et 0.1463.
Es-tu sûr de bien avoir regardé dans yi?
Si all(isnan(yi)) renvoie 0, c'est que yi contient des valeurs qui ne sont pas NaN.

**Newenda** · 15/09/2013, 18h04

Moi rien, y'a bien que des NaN et je me suis pas trompé de yi.

Les valeurs trouvées viennent bien d'une interp1 simple? sans 'method'?

**Jean Dumoncel** · 15/09/2013, 22h01

Voilà le code utilisé:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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
x = [0.1249
    0.1272
    0.1295
    0.1313
    0.1324
    0.1336
    0.1347
    0.1359
    0.1371
    0.1382
    0.1394
    0.1405
    0.1417
    0.1428
    0.1440
    0.1452
    0.1463];
 
Y = [ 0.0300    0.0600    0.0300
    0.0300    0.0300    0.0300
   -0.0300    0.0300    0.1200
    0.0600    0.1201    0.0300
    0.1259    0.1798    0.0599
    0.1801    0.2402    0.1801
    0.1259    0.1199    0.2398
    0.2402    0.4203    0.3002
    0.2462    0.3002    0.4803
    0.3656    0.6593    0.7792
    0.5524    1.6811    2.1614
    0.9171    2.1579    2.5774
    0.8466    3.6024    4.2028
    0.7805    2.4016    1.9213
    0.2997    0.5994    0.4196
    0.0600    0.1201    0.0300
    0.0300    0.1291    0.1937];
 
xi = 0:0.1694/759:0.1694;
 
yi = interp1(x,Y,xi);

La différence avec ton code, c'est que j'ai régénéré le vecteur xi car je ne pouvais pas faire de copier coller. Et là, les lignes 561 à 656 sont différentes de NaN.

**Newenda** · 16/09/2013, 03h52

ha en faite, xi est composé de deux pentes linéaires différentes..

Mais je vois pas en quoi ça pourrait influencer l'apparition de NaN ou non.

Voici pour un test si tu as le temps..

Merci d'avance

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

**Dam2227** · 16/09/2013, 18h20

Je ne sais pas quel genre d'interpolation tu veux faire, mais si tu spécifies une méthode, ça marche. Si tu veux un résultat qui colle au mieux les données, la méthode 'linear' semble être la mieux adaptée. Plusieurs autres méthodes ne marchent pas ton cas ( les valeurs explosent), comme par exemple avec la méthode 'cubic'. Mais ça c'est parce que tu essayes d'extrapoler avec un code pour faire de l'interpolation et dans le cas du 'cubic', le code utilises des polynômes. Les polynômes sont terribles pour faire de l'extrapolation... Mais si tu zoomes dans la région de tes données, après interpolation avec 'cubic', tu vas voir que l'interpolation est bien faites, dans la région où tu as des données. Maintenant, tu vois aussi que avec la méthode 'cubic', le résultat n'est pas super. Et ça, c'est parce que tes données ne sont pas lisses.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
Ys=interp1(x,y,xi,'spline');
 
plot(x,y(:,1),'b');
hold on;
plot(xi,Ys(:,1),'r');

[interp1] problème NaN

MATLAB

Discussions similaires

Partager

Partager