atlab怎么画常微分向量图
MATLAB解常微分方程组并画图的语句 -
求解常微分方程一般使用ode系列函数（其中最常用的是ode45），基本做法就是把微分方程用函数表示出来，然后调用ode函数求解即可。对于本题而言，使用匿名函数描述微分方程更为方便一些。参考代码k0=8; k1=1; K=1; k2=5; n=2; % 常数定义ds = @(t,s)[k0-k1*(1+(s(2)/K)^n)*s(1)后面会介绍。
首先D2y（0）他不可能等于0，二阶微分方程只可能有两个初始值。你想想Dy（0）0,y（0）1带进去D2y（0）Dy（0）y（0）0 ，D2y（0）自然就出来了先定义微分方程组函数，建立rigid.m文件function dy = rigid(t,y)dy = zeros(2,1); % a column vector dy(1) = y(2) ;还有呢？
求解常微分方程一般使用ode系列函数（其中最常用的是ode45），基本做法就是把微分方程用函数表示出来，然后调用ode函数求解即可。对于本题而言，使用匿名函数描述微分方程更为方便一些。参考代码k0=8; k1=1; K=1; k2=5; n=2; % 常数定义ds = @(t,s)[k0-k1*(1+(s(2)/K)^n)*s(1)后面会介绍。
首先D2y（0）他不可能等于0，二阶微分方程只可能有两个初始值。你想想Dy（0）0,y（0）1带进去D2y（0）Dy（0）y（0）0 ，D2y（0）自然就出来了先定义微分方程组函数，建立rigid.m文件function dy = rigid(t,y)dy = zeros(2,1); % a column vector dy(1) = y(2) ;还有呢？
实现方法：1、根据常微分方程，建立其自定义函数fun=@(rho,z)sqrt((0.5+0.5/sqrt(1-rho^2))^2-1)2、利用matlab的ode函数求解。rho,z]=ode45(fun,tspan,z0);3、根据z和ρ的数值绘制其z（ρ）的函数图plot(rho,z);xlabel('ρ'),ylabel('z(ρ)')编程后运行，可得到如下z（ρ）..
本视频展示如何用matlab绘制散圆状态图，可用于相关科研数据绘图！
matlab用dsolve解常微分方程并画图, 例如y'=1-y^2,y(0)=0. -
>> 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)
function y=fun(t,x)y=zeros(3,1);x1=x(1);x2=x(2);x3=x(3);y(1)=x1*(1-x1/150000-0.5*x2/30000);y(2)=x2*(-1-x2/30000+0.5*x1/150000-2*x3/2200);y(3)=x3*(-1-x3/2200+2*x2/30000);2解方程并画图[t,x]=ode23(@fun,[0,1],[121000,27000,2000])还有呢？
求解常微分方程?作图。matlab -
1、首先，建立自定义常微分方程函数，ode_fun(t,y)2、其二，用ode45函数求解其常微分方程的数值解[t,y]=ode45(@ode_fun,tspan,X0);其中：tspan为时间范围，X0为初始条件3、最后，用plot函数绘制t—y(t)，t—y'(t)关系图，以及y(t)—y'(t)关系图。编程后运行，可得到如下结果。
指定点x1(0)=1,x2(0)=2的相图[t,x]=ode45(@(t,x)[x(1)-x(1)^3-x(2);x(1)],[0,10],[1;2]);plot(x(:,1),x(:,2));axis([-3,3,-3,3])再送你一个完整的相图figure [x1,x2]=meshgrid(linspace(-3,3));streamslice(x1,x2,x1-x1.^3-x2,x1)
matlab画常微分方程的相图dx/dt=-y,dy/dt=x -
你的常微分方程缺初始条件。为了说明问题的解决方法，假定其初始条件为x(0)=1、y(0)=2。则其实现代码为syms x(t) y(t)Dx=diff(x,1);Dy=diff(y,1);[x,y]=dsolve(Dx==-y,Dy==x,x(0)==1,y(0)==2)t=0:0.1:10;x=eval(x);y=eval(y);plot(x,y),grid on xlabel('x还有呢？
M是个常数吧 我设为10了，你可以改为你的，对应参数我改为x，y了，这个程序是正确的 我试过了，能出来图function zode clc;clear all;[t,X] = ode45(@fun,[0:0.01:180],[1 1])m=X(:,1);n=X(:,2);figure plot(m,n);xlabel('x');ylabel('y');figure plot(t,m);xl有帮助请点赞。