高阶微分方程matlab

---

东野林19278593483问： matlab求解高阶微分方程 - ？
龙泉驿区和畅回答： clear syms t y=dsolve('D2y=( -3 *cos(2*t)) +2*sin(t)+t-3.8','Dy(0)=0','y(0)=0','t'); ezplot(y,[0 10]); hold on f=@(t,x)[x(2); -3 *cos(2*t)+2*sin(t)+t-3.8]; [t x1]=ode45(f,[0,10],[0 0]); plot(t,x1(:,1),'*');grid on title('连续的是解析解,星号的是数值解');

东野林19278593483问： 如何用matlab 高阶微分方程求解 - ？
龙泉驿区和畅回答： f=@(t,y)[y(2);y(3);-.5*y(1)*y(3)]; r=@(a,b)[a(1);a(2);b(2)-1]; ini=bvpinit(linspace(0,10,5),[0 0 0]); sol=bvp4c(f,r,ini);%边值问题用bvp4c求解,ode45只能求初值问题 x=linspace(0,10); y=deval(sol,x); plot(x,y(1,:))

东野林19278593483问： MATLAB求解高阶常微分方程组 - ？
龙泉驿区和畅回答： function yy = odex3(t,x) yy=zeros(4,1); yy=[x(2);sin(1.5*t)-4*x(2)+2*x(4)-5*x(1)+2*x(3);x(4);0.5*x(2)-x(4)+x(1)-1.5*x(3)];

东野林19278593483问： matlab解高阶常微分方程 - ？
龙泉驿区和畅回答： dY(2) = 2*m/h^2*[U - En]*Y(1); 这个语句的m h等变量 看看都是几维的 只能是一个数值才行

东野林19278593483问： matlab软件能求解高阶线性微分方程吗 - ？
龙泉驿区和畅回答： 可以的 有ode函数 一般可以用ode45函数,这个是四阶,五阶Runge-Kutta单步算法.一般的线性微分方程都能解的.需要按照要求格式来编写方程函数

东野林19278593483问： 我论文需要用MATLAB求解微分方程,Dx= - 0.85*x+0.25*y+0.35*z - 0.8*x^2*y,Dy=0.8*x^2*y - 1.2*y,Dz=0.9*x+0.8*y - 0.3*z用MATLAB画出该微分方程的图像 - ？
龙泉驿区和畅回答：[答案] 用dsolve求不出解析解,用数值解法吧,最好能提供初始条件,这里假设t=0时,x0=1,y0=2,z0=3.function hh[t,x]=ode23t(@fun,0:1,[1 2 3])plot(t,x)function yhat=fun(t,x)dxdt=-0.85*x(1)+0.25*x...

东野林19278593483问： 如何用matlab 求解高阶方程? - ？
龙泉驿区和畅回答： 用roots(p)函数求解 其中p是方程的各次方的系数组成的矩阵 比如 你这道题的求根程序为: roots(p)matlab运行后的结果为:ans =10.3180 4.0487 1.8166 + 1.3141i 1.8166 - 1.3141i

东野林19278593483问： 高分求助,用MATLAB求高阶微分方程的数值解作图 - ？
龙泉驿区和畅回答： [t,y] = ode45(@fun,[0 100],[0; 0; 10600; 0]); plot(y(:,1),y(:,3)); function dy = fun(t,x) dy = [ x(2); -0.20811*(0.0587*x(4)+0.9983*x(2)); x(4); -10.01476+0.20811*(0.0587*x(2)-0.9983*x(4)); ]; end 初始条件 x'(0),y'(0)未给,这里给0

东野林19278593483问： 用matlab,怎么解高阶方程?????? - ？
龙泉驿区和畅回答： ode方法是解微分方程的.ode45也不是说解决4,5阶的方程,而是一种45阶的龙哥库塔方法.solve解决的是代数方程组问题.

---