求微分方程 (1+x)D2y＝2y－4 初始條件 y(0)＝0 y(1)＝2Dy(1)

如果想用inline和ode45解決，不用function窗口，該如何做？

解

1. 對于此類邊值微分方程，ode**函數是無力直接求解的，Matlab提供了bvp解算器。 2. 對干你用的Matlab版本，用@(x)函數(匿名函數)比inline更方便。 3. 請參考本文中的其它例題及相關資料理解下面的代碼。

http://hi.baidu.com/nature_and_man/item/6757103c923be981b611db6c

4. 在Matlab7.1版上，可用： dsolve('D2y = (2*y-4)/(1+t)', 'y(0) = 0', 'Dy(1) = y(1)/2') 求出解析解(符號解)。 如不需要解析解與數值解的對比，可不要第一段代碼。

% 解析解 syms t y = -(1+t)^(1/2)*besseli(1,2*2^(1/2)*(1+t)^(1/2))*(4*i*bessely(0,4*i)-2*bessely(1,4*i)+2^(1/2)*bessely(1,2*i*2^(1/2)))/(2*i*bessely(0,4*i)*besseli(1,2*2^(1/2))-bessely(1,4*i)*besseli(1,2*2^(1/2))+besseli(1,4)*bessely(1,2*i*2^(1/2))-2*besseli(0,4)*bessely(1,2*i*2^(1/2)))+(1+t)^(1/2)*bessely(1,2*i*2^(1/2)*(1+t)^(1/2))*(besseli(1,2*2^(1/2))*2^(1/2)-2*besseli(1,4)+4*besseli(0,4))/(2*i*bessely(0,4*i)*besseli(1,2*2^(1/2))-bessely(1,4*i)*besseli(1,2*2^(1/2))+besseli(1,4)*bessely(1,2*i*2^(1/2))-2*besseli(0,4)*bessely(1,2*i*2^(1/2)))+2; ezplot(y,[0 1])

% 數值解 dydx = @(x,y) [y(2);(2*y(1)-4)/(1+x) ]; %邊值微分方程 res = @(ya,yb) [yb(2) - yb(1)/2;ya(1) - 0 ]; %邊界條件 solinit = bvpinit(linspace(0,1,10),[1 0]); sol = bvp4c(dydx,res,solinit); xint = linspace(0,1,50); Sxint = deval_r(sol,xint);

% 畫圖 hold on plot(xint,Sxint(1,:),'*') title('stevenchang041''s equation.') xlabel('x') ylabel('solution y')

總結

以上是生活随笔為你收集整理的bvp解算器是什么_边值问题(BVP)的Matlab解法的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。