交通管理部門如何掌握一座橋梁的通行情況?
1、交通管理部門為了掌握一座橋梁的通行情況,在橋梁的一端每隔一段不等的時(shí)間連續(xù)記錄1min內(nèi)通過橋梁的車輛數(shù)量,連續(xù)觀測(cè)一天24h的通過車輛車輛數(shù)據(jù)如下表所示。試建立模型分析估計(jì)這一天中總共有多少車輛通過這座橋梁。
代碼如下:
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ù)量(單位:十萬)如下表所示。預(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)式的擬合效果就很好,所以接下來預(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)更多
-
即日-10.29立即報(bào)名>> 2024德州儀器嵌入式技術(shù)創(chuàng)新發(fā)展研討會(huì)
-
10月31日立即下載>> 【限時(shí)免費(fèi)下載】TE暖通空調(diào)系統(tǒng)高效可靠的組件解決方案
-
即日-11.13立即報(bào)名>>> 【在線會(huì)議】多物理場(chǎng)仿真助跑新能源汽車
-
11月14日立即報(bào)名>> 2024工程師系列—工業(yè)電子技術(shù)在線會(huì)議
-
12月19日立即報(bào)名>> 【線下會(huì)議】OFweek 2024(第九屆)物聯(lián)網(wǎng)產(chǎn)業(yè)大會(huì)
-
即日-12.26火熱報(bào)名中>> OFweek2024中國(guó)智造CIO在線峰會(huì)
推薦專題
- 1 Intel宣布40年來最重大轉(zhuǎn)型:年底前裁員15000人、拋掉2/3房產(chǎn)
- 2 因美封殺TikTok,字節(jié)股價(jià)骨折!估值僅Meta1/5
- 3 宏山激光重磅發(fā)布行業(yè)解決方案,助力智能制造產(chǎn)業(yè)新飛躍
- 4 國(guó)產(chǎn)AI芯片公司破產(chǎn)!白菜價(jià)拍賣
- 5 具身智能火了,但規(guī)模落地還需時(shí)間
- 6 國(guó)產(chǎn)英偉達(dá)們,抓緊沖刺A股
- 7 三次錯(cuò)失風(fēng)口!OpenAI前員工殺回AI編程賽道,老東家捧金相助
- 8 英特爾賦能智慧醫(yī)療,共創(chuàng)數(shù)字化未來
- 9 英偉達(dá)的麻煩在后頭?
- 10 將“網(wǎng)紅”變成“商品”,AI“爆改”實(shí)力拉滿
- 高級(jí)軟件工程師 廣東省/深圳市
- 自動(dòng)化高級(jí)工程師 廣東省/深圳市
- 光器件研發(fā)工程師 福建省/福州市
- 銷售總監(jiān)(光器件) 北京市/海淀區(qū)
- 激光器高級(jí)銷售經(jīng)理 上海市/虹口區(qū)
- 光器件物理工程師 北京市/海淀區(qū)
- 激光研發(fā)工程師 北京市/昌平區(qū)
- 技術(shù)專家 廣東省/江門市
- 封裝工程師 北京市/海淀區(qū)
- 結(jié)構(gòu)工程師 廣東省/深圳市