計算機仿真報告word版(編輯修改稿)
2025-04-19 03:35 本頁面
【文章內容簡介】 解：用MATLAB編程求解如下：%改用Euler法求初值問題的數(shù)值解function FZSYZ1t0=0。tf=3。x0=1。h1=。t=[t0:h1:tf]。n=length(t)。u=x0。uu(1)=u。for i=2:n du=u(t(i1))^2。 u=du*h1+u。 uu(i)=u。enduuplot(t,uu)得到的實驗結果如下圖所示 試將（21）方程改為用ode23算法調用函數(shù)求解，并試比較結果。解：用MATLAB編程求解如下：function xdot=fun21(t,x)xdot=x2*t./x。function fzsy22t0=0。tf=1。tol=1e6。x0=1。trace=1。[t,x]=ode23(39。fun2139。,t0,tf,x0,tol,trace)。plot(t,x)得到的實驗結果如下圖所示：利用ode23或ode45求解線性時不變系統(tǒng)微分方程，并繪制出曲線，式中 解：用MATLAB編程求解如下：function xdot=fun21(t,x)A=[ 1。1 ]。xdot=A*x。function fzsy22t0=0。tf=4。tol=1e6。x0=[0。1]。trace=1。[t,x]=ode23(39。fun2139。,t0,tf,x0,tol,trace)。plot(t,x)得到的實驗結果如下圖所示：求出與的單位階躍響應，并分別求出狀態(tài)空間模型。解：（1）的狀態(tài)空間模型求解如下：function shiyan2b1=[2]。a1=[1 2 1]。sys=t
范文總結相關推薦
鄂ICP備17016276號-1