EDA365电子论坛网

标题: 求y（t）数值解并绘制其图像，与解析解对比精度。 [打印本页]

作者: YNhhU 时间: 2020-5-15 13:36
标题: 求y（t）数值解并绘制其图像，与解析解对比精度。

微分方程y^'' (x)-(2-1/x) y^' (x)+(1-1/x)y(x)=x^2 e^(-5x)，满足边界条件y(1)=π，y(π)=1，求y（t）数值解并绘制其图像，与解析解对比精度。

作者: Zedd 时间: 2020-5-15 14:46
你这是带边界条件的微分方程，采用有限差分法求解数值解，采用dsolve求解解析解。

作者: sharkN 时间: 2020-5-15 16:04
典型的边值微分方程问题（BVP ODE），Matlab可用bvp4c命令试试，下面使用更简单的1stOpt求解：

Variable x=[1,pi],y=[pi,1];
Plot x[x],y,y'[y2];
ODEFunction y''=x^2*exp(-5*x)-(1-1/x)*y+(2-1/x)*y';

欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)