對(duì)于一些微分方程來(lái)說(shuō)，數(shù)值解法對(duì)于求解具有很好的幫助，因?yàn)殡y以求得其原方程。

比如方程：

但是我們知道了它的初始條件，這對(duì)于我們疊代求解很有幫助，也是必須的。

那么現(xiàn)在我們也用Python去解決這一些問(wèn)題，一般的數(shù)值解法有歐拉法、隱式梯形法等，我們也來(lái)看看這些算法對(duì)疊代的精度有什么區(qū)別？

```python```pythonimport numpy as npfrom scipy.integrate import odeintfrom matplotlib import pyplot as pltimport os#先從odeint函數(shù)直接求解微分方程#創(chuàng)建歐拉法的類class Euler: #構(gòu)造方法，當(dāng)創(chuàng)建對(duì)象的時(shí)候，自動(dòng)執(zhí)行的函數(shù) def __init__(self,h,y0):#將對(duì)象與對(duì)象的屬性綁在一起self.h = hself.y0 = y0self.y = y0self.n = 1/self.hself.x = 0self.list = [1]#歐拉法用list列表,其x用y疊加儲(chǔ)存self.list2 = [1]self.y1 = y0#改進(jìn)歐拉法用list2列表，其x用y1疊加儲(chǔ)存self.list3 = [1]self.y2 = y0#隱式梯形法用list3列表，其x用y2疊加儲(chǔ)存 #歐拉法的算法，算法返回t,x def countall(self):for i in range(int(self.n)): y_dere = -20*self.list[i] #歐拉法疊加量y_dere = -20 * x y_dere2 = -20*self.list2[i] + 0.5*400*self.h*self.list2[i] #改進(jìn)歐拉法疊加量 y_dere2 = -20*x(k) + 0.5*400*delta_t*x(k) y_dere3 = (1-10*self.h)*self.list3[i]/(1+10*self.h) #隱式梯形法計(jì)算 y_dere3 = (1-10*delta_t)*x(k)/(1+10*delta_t) self.y += self.h*y_dere self.y1 += self.h*y_dere2 self.y2 =y_dere3 self.list.append(float('%.10f' %self.y)) self.list2.append(float('%.10f'%self.y1)) self.list3.append(float('%.10f'%self.y2))return np.linspace(0,1,int(self.n+1)), self.list,self.list2,self.list3step = input('請(qǐng)輸入你需要求解的步長(zhǎng):')step = float(step)work1 = Euler(step,1)ax1,ay1,ay2,ay3 = work1.countall()#畫(huà)圖工具pltplt.figure(1)plt.subplot(1,3,1)plt.plot(ax1,ay1,’s-.’,MarkerFaceColor = ’g’)plt.xlabel(’橫坐標(biāo)t’,fontproperties = ’simHei’,fontsize =20)plt.ylabel(’縱坐標(biāo)x’,fontproperties = ’simHei’,fontsize =20)plt.title(’歐拉法求解微分線性方程步長(zhǎng)為’+str(step),fontproperties = ’simHei’,fontsize =20)plt.subplot(1,3,2)plt.plot(ax1,ay2,’s-.’,MarkerFaceColor = ’r’)plt.xlabel(’橫坐標(biāo)t’,fontproperties = ’simHei’,fontsize =20)plt.ylabel(’縱坐標(biāo)x’,fontproperties = ’simHei’,fontsize =20)plt.title(’改進(jìn)歐拉法求解微分線性方程步長(zhǎng)為’+str(step),fontproperties = ’simHei’,fontsize =20)plt.subplot(1,3,3)plt.plot(ax1,ay3,’s-.’,MarkerFaceColor = ’b’)plt.xlabel(’橫坐標(biāo)t’,fontproperties = ’simHei’,fontsize =20)plt.ylabel(’縱坐標(biāo)x’,fontproperties = ’simHei’,fontsize =20)plt.title(’隱式梯形法求解微分線性方程步長(zhǎng)為’+str(step),fontproperties = ’simHei’,fontsize =20)plt.figure(2)plt.plot(ax1,ay1,ax1,ay2,ax1,ay3,’s-.’,MarkerSize = 3)plt.xlabel(’橫坐標(biāo)t’,fontproperties = ’simHei’,fontsize =20)plt.ylabel(’縱坐標(biāo)x’,fontproperties = ’simHei’,fontsize =20)plt.title(’三合一圖像步長(zhǎng)為’+str(step),fontproperties = ’simHei’,fontsize =20)ax = plt.gca()ax.legend((’$Eular$’,’$fixed Eular$’,’$trapezoid$’),loc = ’lower right’,title = ’legend’)plt.show()os.system('pause')

對(duì)于歐拉法，它的疊代方法是：

改進(jìn)歐拉法的疊代方法：

隱式梯形法：

對(duì)于不同的步長(zhǎng)，其求解的精度也會(huì)有很大的不同，我先放一幾張結(jié)果圖：

補(bǔ)充：基于python的微分方程數(shù)值解法求解電路模型

安裝numpy（用于調(diào)節(jié)range） 和 matplotlib（用于繪圖）

在命令行輸入

pip install numpy pip install matplotlib電路模型和微分方程模型1

無(wú)損害，電容電壓為5V，電容為0.01F，電感為0.01H的并聯(lián)諧振電路

電路模型1

微分方程1

微分方程2