我很難找到資料點的插值。這條線應該稍微像乙個負的反二次曲線(比如倒「c」)。在
因為這不是乙個函式(x可以有多個y值),所以我不確定要使用什麼樣的插值。在
我在想,也許我應該翻轉軸來建立插值點/線,使用單變數樣條線,然後在繪製時將其翻轉回去?在
這是乙個單獨點的圖表:
這是我的**:import datetime as dt
import matplotlib.pyplot as plt
from scipy import interpolate
file = open_file("010217.hdf5", mode = "a", title = 'sondrestrom1')
all_data = file.getnode('/data/table layout').read()
file.close()
time = all_data['ut1_unix'] #time in seconds since 1/1/1970
alt = all_data['gdalt'] #all altitude points
electrondens = all_data['nel'] #all electron density points
x =
y =
positions =
for t in range(len(time)): #looking at this specific time, find all the respective altitude and electron density points
if time[t] == 982376726:
#finding the date
datetime1970 = dt.datetime(1970,1,1,0,0,0)
seconds = long(time[t])
newdatetime = datetime1970 + dt.timedelta(0, seconds)
time1 = newdatetime.strftime('%y-%m-%d %h:%m:%s')
title = "electron density vs. altitude at "
title += time1
plt.plot(x,y,"o")
plt.title(title)
plt.xlabel('electron density (log_10[ne])')
plt.ylabel('altitude (km)')
plt.show()
MATLAB一維插值和二維插值
插值問題描述 已知 乙個函式上的若干點,但函式具體表示式未知,現在要利用已知的若干點求在其他點處的函式值,這個過程就是插值的過程.1.一維插值 一維插值就是給出y f x 上的點 x1,y1 x2,y2 xn,yn 由此求出y f x 在點xa處的值ya的值.實現一維插值使用interp1命令,使用...
插值查詢 插值搜尋
這是一種和二分比較相似的查詢的演算法,不過不同的是,對於分布比較均勻的較大的陣列,插值查詢有時能夠一次就搜尋到位.為什麼能夠這麼快呢 看網上沒有什麼關於這種演算法的描述,我就來描述一下吧.首先要知道一點,這種搜尋方式只能夠針對順序表進行,再乙個要理解順序表中的乙個特點,在順序表中查詢是否存在乙個值,...
MATLAB二維插值和三維插值
插值問題描述 已知 乙個函式上的若干點,但函式具體表示式未知,現在要利用已知的若干點求在其他點處的函式值,這個過程就是插值的過程.1.一維插值 一維插值就是給出y f x 上的點 x1,y1 x2,y2 xn,yn 由此求出y f x 在點xa處的值ya的值.實現一維插值使用interp1命令,使用...