2013-10-28 19 views
5

Zaman içinde özel bir ölçümle, x bir gappy zaman serisine sahibim. Ayrıca belirsizliğin bir ölçüsü var sx (belirli bir aralıkta x standart sapma). Ben sadece temelde zaman serilerini boşluk doldurmak istiyorum, ama ben ölçüm belirsizliğini ve umarım enterpolasyon belirsizliğini yaymak istiyorum. Geçmişte, bunu başarmak için Tikhonov regülasyonunu kullandım (diğer bir deyişle, akıcı bir kısıtlama ekleyerek), ama scipy'den bir paket dışı rutini kullanmayı tercih ediyorum. Scipy spline enterpolasyon rutinlerinin bir smoothing parametresi aldığını görebiliyorum (sx için türetilmiş olacaktır), fakat enterpolasyonlu dizinin 'belirsizliğini hesaplamıyor.Belirsizliğin yorumlanması ve hesaplanması

Bunu soruyorum çünkü bunun önemsiz bir hesaplama olduğunu ve bu yeteneğin uygun olup olmadığını bilip bilmediğini görmüyorum.

+1

Bunu yapmanın bir yolu, bir gauss işleminin kullanılmasıdır. Scikits.learn bir göz atın: http://scikit-learn.org/stable/modules/gaussian_process.html (Ayrıca, kriging belirli bir türden bir çingene sürecidir ve enterpolasyon dünyasında karşılaşacağınız bir yöntemdir. çok) –

+0

Örnek kullanarak, tam bir örnek için: http://scikit-learn.org/stable/auto_examples/gaussian_process/plot_gp_regression.html –

+0

Evet, bunun farkındayım, ama scipy spline rutinlerini umuyordum "kutunun dışında" belirsizlik sağlayabilir. – Jose

cevap

1

"Ölçüm belirsizliğini ve umarız, enterpolasyon belirsizliğini yaymak" ile ne kastettiğinizden emin değilim, ancak sizi doğru anlarsam, gürültülü bir ölçüm ve belirsizlikle ilgili ek bilgileriniz olur. Gürültü beyaz gürültü ise veya beyazlatılabilirse, normal ölçümü düzeltmek ve boşluklardaki ölçümleri tahmin etmek için Kalman filter kullanabilirsiniz. Bir Kalman filtresi, gelecekteki ölçümleri tahmin etmek için bir durum modeli kullanır ve bir kovaryans matrisi kullanarak sinyalin beyaz gürültü bileşenini modeller. Belirsizlik ölçümü, filtrenin tahmin ve ölçüm arasındaki "güven" ağırlığını değiştirmek için uygulanabilir. Geçici belirsizlik önlemini belirlemek için kovaryans matrisini kullanabilirsiniz.

+0

Aslında, bir filtreden daha yumuşak bir Kalman olmasını tercih ederim, ama Scipy'nin Spline enterpolasyon rutinlerinin belirsizliği hesaba katıp katmadığını sordum (bunu hesaplayabilirsin, ama önemsiz olmayabilir!) – Jose

+0

Üzgünüz, Sana yardım edemem, hiç kullanmadım. Kalman pürüzsüzlüğü, daha genel Kalman filtresi konseptinin özel bir çalışma koşulu. –

İlgili konular