MATLAB解dy dx=x y,y(0)=1,步长h=0.1
来源:学生作业帮助网 编辑:作业帮 时间:2024/10/01 06:24:14
clear>>symsxy>>z=x.*y;>>ezsurf(x,y,z)我就已经实现好了你也试试吧
方程两边求关x的导数ddx(xy)=(y+xdydx); ddxex+y=ex+y(1+dydx);所以有 (y+xdy
由微分方程dydx=2xy,得dyy=2xdx(y≠0)两边积分得:ln|y|=x2+C1即y=Cex2(C为任意常数)
看不懂问题估计我还没学.
1.建立M文件vdp1001a.m如下:functiondy=vdp1001a(x,y)dy=zeros(2,1);dy(1)=y(2);dy(2)=2*x*y(2)/(1+x^2);2.取x0=0,
[x,y]=meshgrid(1:10,1:10);z=x.*y;surf(x,y,z)不知道lz是不是这个意思?
上边两个图应该是不对吧,因为固定y=0或者固定z=0的时候应该是个抛物线吧~clear;closeall;z=solve('z^2=y^2+x*y','z')[x,y]
方程两边对x求导得2x+y′x2+y=3x2y+x3y′+cosxy′=2x−(x2+y)(3x2y+cosx)x5+x3y−1由原方程知,x=0时y=1,代入上式得y′|x=0=dydx|x=0=1
1:数值求解出对应的值2:按照需求进行作图(x-y;y-y';x-y')3:没了其实你的问题问的就很模糊~
x=[-pi:pi/100:pi];[x,y]=meshgrid(x,x);z=x.*y;surf(x,y,z)shadingflat
x=[0:0.01:2*pi];ezplot('x*y+sin(x)-y');关于ezplot的用法.请在matlab窗口中输入docezplot简单来说就是绘制隐函数定义的图像.
这是一阶线性微分方程,其中P(x)=1,Q(x)=e-x∴通解y=e−∫dx(∫e−x•e∫dxdx+C)=e−x(∫e−x•exdx+C)=e−x(x+C).
ezplot('y*(0.005*(sin(x))^2)+y^2*(1e-6*(cos(x))^4+5e-5*(sin(x))^4+6.1793e-4*(cos(x))^2*(sin(x))^2)-1
dydx要是等式才行吧.如果是的话,这句话就是求这个等式的根,用r表示x.
x=linspace(-1,1);y=linspace(-1,1);[X,Y]=meshgrid(x,y); Z=X.*Y;mesh(X,Y,Z);
x=-2:0.1:2;y=x;[x,y]=meshgrid(x,y);z=x.*y;surf(x,y,z);grid on;xlabel('x.axis');ylabel(&
x1=linspace(-5,-.2,1000);x2=linspace(0.2,5,1000);y1=1./x1;y2=1./x2;plot(x1,y1);holdon;plot(x2,y2);再问
x=linspace(1,5,1000);y=linspace(2,15,1000);f=x.^2+x.*y;fmin=min(f)fmax=max(f)fmin=3fmax=100再问:我对matl
在方程ex+y+cos(xy)=0左右两边同时对x求导,得:ex+y(1+y′)-sin(xy)•(y+xy′)=0,化简求得:y′=dydx=ysin(xy)−ex+yex+y−xsin(xy).