matlab求解微分方程 matlab画图线

1.可以解析求解的微分方程.dsolve() 调用格式为:y=dsolve(f1,f2,.,fmO; y=dsolve(f1,f2,.,fm,'x'); 如下面的例子,求解了微分方程 syms t; u=exp(-5*t)*cos(2*t-1)+5; uu=5.

那y1,y2呢? 那个在y(3)是y的三次导的意思. 上面所有的y1,y2我认为也是求导的意思. 你可以做数值解, 下面做数值解: 在Matlab下输入:editzhidao_rk4_4.m,然后将下.

楼主是要解析解吧 dsolve('dy=a/v-b/v*sqrt(2*g*y)','y(0)=h0','t') ans = (2^(1/2)*a + 2^(1/2)*a*lambertw(0, (2^(1/2)*b*g^(1/2)*exp(-(a + b^2*g*(t/v - (a*(log(-(exp(-(a + 2^(1/2)*b*(g*h.

先分析,先把矩阵算一下,化简整理下 x'-4t^2x=0 其中x是t的函数 编程如下; dsolve('Dy-4*t^2*y=0') 因为matlab中默认的Dy表示y对t的一阶导,D2y表示y对t的二阶导,没有Dx表示x对t的一阶导 所以会这样编程 运行结果:C2*exp((4*t^3)/3)+C1

这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数 就一个例子(选自matlab帮助文档) 解入下微分方程组1,建立方程组函数 function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2);2,求解并绘图 [T,Y] = ode45(@rigid,[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') 结果

y=dsolve('D2y+2*Dy+2*y=0','y(0)=1,Dy(0)=0');t=0:0.01:20;y=eval(y);plot(t,y)祝你学习愉快!

直接解就可以了>> S=dsolve('Dx+5*x+y=exp(t)','Dy-x-3*y=0','x(0)=1','y(0)=0') S = x: [1x1 sym] y: [1x1 sym]>> S.x,S.y ans = -4*exp((-1+15^(1/2))*t)*(1/22+13/330*15^(1/2))+exp(.

用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=.

好像没什么特别好的办法: >> syms t q y; q=4*t^3;>> syms Dy>> eq=Dy-(q-2*t);>> y=dsolve(char(subs(eq)),'y(0)=5','t')y = t^4-t^2+5也就是要先移项把微分方程写成等号右端为零的形式(需额外定义Dy),然后subs代入,再求解.

help desolve