jmfs.net
当前位置:首页>>关于MATLAB微分方程问题 dsolve('D2y-t*Dy=1-t^2','y(0...的资料>>

MATLAB微分方程问题 dsolve('D2y-t*Dy=1-t^2','y(0...

问题出在你定义exf(t,x)有错误. 要改成下列形式才对. function xdot=exf(t,x) xdot=[x(2);1-t^2+t*x(2)]; end 运行结果

这是我运行的结果:>> syms x y;>> y=dsolve('D2y+Dy-2*y=2*x','y(0)=0,Dy(0)=1','x')y =-1/2*exp(-2*x)+exp(x)-1/2-x你的2y与2x中间没有打乘号*

y=dsolve('D3y=D2y-Dy-y+t^2','y(0)=0','Dy(0)=1','D2y(0)=-1');tt=linspace(0,10,1001); %表示t从0到10for i=1:1001t=tt(i);yy(i)=eval(y);endplot(tt,yy)

>> y=dsolve('Dy=1-y^2','y(0)=0') y = tanh(t)>> t=-2*pi:0.1*pi:2*pi;>> y=tanh(t);>> plot(t,y)

y=dsolve('D2y-7*(1-y^2)*Dy+y=0','y(0)=1,Dy(0)=0','t')命令窗口输出:警告: Explicit solution could not be found.所以采用上述方法找不到解析解.数值解:function dx=a(t,x) dx=zeros(2,1); dx(1)=x(2); dx(2)=7*(1-x(1)^2)*x(2)-x(1); [T,Y]=ode15s('a',[0 50],[1 0]) plot(T,Y(:,1))上半段用m文件写成函数形式并保存文件名为a.m,下半段在命令窗口输入,得到图像:

y(n)代表了上面的微分方程的第n个解.比如你这个的y(2)=C3*x - C3^2,也是上面的微分方程的一个解.

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)祝你学习愉快!

ydot=2/sqrt(1-t^2); % 1-t^2 s=dsolve('dy=2/sqrt(1+t^2)','y(0)=1') % 1+t^2 ,能一样吗

dsolve('D2y/Dx2=Dy/Dx-3*x^2+6*x') dsolve('(1+x)*D2y/Dx2+Dy/Dx=0') matlab里 一个一个复制进去,回车就行 第一个解 ans = [ x(t)=-1+(1-2*Int(diff(y(t),t)*diff(x2(t),t)/diff(y(t),$(t,2)),t)-2*C1)^(1/2), x(t)=-1-(1-2*Int(diff(y(t),t)*diff(x2(t),t)/diff(y(t),$(t,2)),t)-2*C1)^(1/2)] 第二个有些问题,自己看看吧

>> y=dsolve('d2y-dy+y==1','y(0)==1','y(1)==0') y =1 - (exp(-1/2)*exp(t/2)*sin((3^(1/2)*t)/2))/sin(3^(1/2)/2)

相关文档
so1008.com | fpbl.net | zxqs.net | zxpr.net | qwfc.net | 网站首页 | 网站地图
All rights reserved Powered by www.jmfs.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com