ntjm.net
当前位置:首页 >> mAtlAB求解微分方程 >>

mAtlAB求解微分方程

那y1,y2呢?那个在y(3)是y的三次导的意思.上面所有的y1,y2我认为也是求导的意思(你有一些没有加括号) Dsolve('D3y+t*Dy*D2y+t^2*Dy*y^2=exp(-t*y)','y(0)=2','Dy(0)=0','D2y(0)=0') 但是没有精确解.你可以做数值解,下面做数值解:在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*diff(u,t,2)+4*diff(u,t)+2*u; syms t y; y=dsolve(['D4y+10*D3y+35*D2y+50*Dy

dsolve()求解微分方程后,使用ezplot()作图,参考代码:>> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间

syms u;dsolve('D2y+Dy+3*y^3-u=0');不过此答案没有显示解换个简单的比喻y'=2x可以用 dsolve('Dy=2*x',x);

用matlab求解微分方程组并画图,其实现过程:1、首先应根据已知微分方程组,编写其函数 dy(1)=-0.10982*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3)); dy(2)=1441.5074*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3));2、用ode45函数求解其微分方程组的数值解 [t,y]=ode45(@ode_fun,tspan,y0)3、然后用plot绘图函数,绘出t-Cg和t-θ关系曲线图4、编程后运行结果

syms A B C x y;>> dsolve('Dy=A-B*y-C*x*y','x')ans =-1/2*A*erf(1/2*(-2*C)^(1/2)*x-B/(-2*C)^(1/2))*exp(-1/2*x*(2*B+C*x)-1/2*B^2/C)/C*2^(1/2)*(-pi*C)^(1/2)+exp(-1/2*x*(2*B+C*x))*C1

这说明此微分方程没有显示解,应用数值方法求解,例如ode45等函数就一个例子(选自matlab帮助文档)解入下微分方程组1,建立方程组函数function dy = rigid(t,y)dy = zeros(3,1); % a column vectordy(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),'.')结果

syms x y; y=dsolve('D2y+2*Dy+2*y=0','y(0)=1','Dy(0)=0') 解得 y = exp(-t)*cos(t)+exp(-t)*sin(t)

>> syms r k positive>> dsolve('Dy = r*x*(1-x/k) - x^2/(1+x^2)')ans =C2 - (t*(k*x^2 + r*x^2 + r*x^4 - k*r*x^3 - k*r*x))/(k*x^2 + k)

[x,yz]=ode23(@(x,yz)[f(x,yz(1),yz(2));g(x,yz(1),yz(2))],[x0,xend],[y0,z0])

jingxinwu.net | rjps.net | 9371.net | tuchengsm.com | nmmz.net | 网站首页 | 网站地图
All rights reserved Powered by www.ntjm.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com