交通管理部門如何掌握一座橋梁的通行情況?
1、交通管理部門為了掌握一座橋梁的通行情況,在橋梁的一端每隔一段不等的時(shí)間連續(xù)記錄1min內(nèi)通過(guò)橋梁的車輛數(shù)量,連續(xù)觀測(cè)一天24h的通過(guò)車輛車輛數(shù)據(jù)如下表所示。試建立模型分析估計(jì)這一天中總共有多少車輛通過(guò)這座橋梁。
代碼如下:
import numpy as np
import matplotlib.pyplot as plt
from scipy.interpolate import interp1d
plt.grid(linestyle='--') #設(shè)置網(wǎng)格線---這里是虛線
X = np.array([0,2,4,5,6,7,8,9,10.5,11.5,12.5,14,16,17,18,19,20,21,22,23,24])
Y = np.array([2,2,0,2,5,8,25,12,5,10,12,7,9,28,22,10,9,11,8,9,3])
xnew= np.linspace(0,24,24*60)
# 引入插值函數(shù),這里的點(diǎn)較散亂,所以我選取高階的樣條插值
f=interp1d(X,Y,'quadratic'); ynew=f(xnew)
plt.plot(xnew,ynew,c='b')
plt.show()
# 計(jì)算一天的車輛總數(shù)
c=sum(ynew)
print('一天總共:',c)
2、某年美國(guó)舊車價(jià)格的調(diào)查資料如下表所示,其中xi表示轎車的使用年數(shù),yi表示相應(yīng)的平均價(jià)格.試分析用什么形式的曲線擬合表中所給的數(shù)據(jù),并預(yù)測(cè)使用4.5年后轎車的平均價(jià)格大致為多少?
代碼如下:
import matplotlib.pyplot as plt
import numpy as np
from numpy import polyfit,polyval
x=np.linspace(0,10,11)[1:11]
y=np.array([2615,1943,1494,1087,765,538,484,290,226,204])
p=polyfit(x,y,2)
print('p:',p)
ynew=polyval(p,x)
plt.plot(x,ynew)
plt.show()
ypred=polyval(p,4.5)
print('使用4.5年后轎車的平均價(jià)格大致為:',ypred)
3、求下列微分方程組(豎直加熱板的自然對(duì)流)的數(shù)值解
代碼如下:
import matplotlib.pyplot as plt
import numpy as np
from scipy.integrate import odeint
def fT(y,x):
f,df,ddf,T,dT=y; #拆分、降階
return np.array([df,ddf,2*(df*df)-T-3*f*ddf,dT,-2.1*f*dT])
x=np.arange(0,10,0.1)
ynew=odeint(fT,[0.0,0.0,0.68,1.0,-0.5],x)
#plt.rc('font',size=16);plt.rc('font',family='Times New Roman')
plt.plot(x,ynew[:,0],label='f')
plt.plot(x,ynew[:,3],label='T')
plt.legend()
plt.show()
4、某地區(qū)野兔的數(shù)量連續(xù)9年的統(tǒng)計(jì)數(shù)量(單位:十萬(wàn))如下表所示。預(yù)測(cè)t=9,10時(shí)野兔的數(shù)量。
代碼如下:
import matplotlib.pyplot as plt
import numpy as np
x=np.linspace(0,8,9)
y=np.array([5,5.9945,7.0932,8.2744,9.5073,10.7555,11.9804,13.1465,14.2247])
plt.scatter(x,y,c='r',label='Data point')
p=np.polyfit(x,y,1)
print('p為:',p)
ynew1=np.polyval(p,x)
plt.plot(x,ynew1)
plt.legend()
plt.show()
#從圖中可以看出一次多項(xiàng)式的擬合效果就很好,所以接下來(lái)預(yù)測(cè)時(shí)采用一次多項(xiàng)式即可
y9=np.polyval(p,9)
y10=np.polyval(p,10)
print('t=9野兔的數(shù)量:',y9)
print('t=10野兔的數(shù)量:',y10)

發(fā)表評(píng)論
請(qǐng)輸入評(píng)論內(nèi)容...
請(qǐng)輸入評(píng)論/評(píng)論長(zhǎng)度6~500個(gè)字
最新活動(dòng)更多
-
7月22-29日立即報(bào)名>> 【線下論壇】第三屆安富利汽車生態(tài)圈峰會(huì)
-
7.30-8.1火熱報(bào)名中>> 全數(shù)會(huì)2025(第六屆)機(jī)器人及智能工廠展
-
7月31日免費(fèi)預(yù)約>> OFweek 2025具身智能機(jī)器人產(chǎn)業(yè)技術(shù)創(chuàng)新應(yīng)用論壇
-
免費(fèi)參會(huì)立即報(bào)名>> 7月30日- 8月1日 2025全數(shù)會(huì)工業(yè)芯片與傳感儀表展
-
即日-2025.8.1立即下載>> 《2024智能制造產(chǎn)業(yè)高端化、智能化、綠色化發(fā)展藍(lán)皮書》
-
8月5日立即報(bào)名>> 【在線會(huì)議】CAE優(yōu)化設(shè)計(jì):醫(yī)療器械設(shè)計(jì)的應(yīng)用案例與方案解析
推薦專題
- 1 AI 眼鏡讓百萬(wàn) APP「集體失業(yè)」?
- 2 豆包前負(fù)責(zé)人喬木出軌BP后續(xù):均被辭退
- 3 一文看懂視覺(jué)語(yǔ)言動(dòng)作模型(VLA)及其應(yīng)用
- 4 “支付+”時(shí)代,支付即生態(tài) | 2025中國(guó)跨境支付十大趨勢(shì)
- 5 中國(guó)最具實(shí)力AI公司TOP10
- 6 深圳跑出40億超級(jí)隱形冠軍:賣機(jī)器人年入6.1億,港股上市
- 7 特斯拉Robotaxi上路,馬斯克端上畫了十年的餅
- 8 “AI六小虎”到了下一個(gè)賽點(diǎn)
- 9 張勇等人退出阿里合伙人
- 10 AI視頻,攪動(dòng)1.5萬(wàn)億市場(chǎng)