zxqk.net
当前位置:首页 >> mAtlAB画微分方程相图 >>

mAtlAB画微分方程相图

clear;clc [x,y]=meshgrid(linspace(-5,5)); streamslice(x,y,y.*(4+x.^2-y.^2),-x.*(-2-x.^2+y.^2)); xlabel('x');ylabel('y');

你的常微分方程缺初始条件.为了说明问题的解决方法,假定其初始条件为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 onxlabel('x'),ylabel('y')运行后得到x,y的相平面图

题主的图形属于相平面图,即x(t)与y(t)的关系曲线图.得到该关系曲线图可以按下列方法来做.1、根据给出的常微分方程组,自定义其函数ode_fun(t,x);2、设定参数a、b、k(应该是已知的);3、利用ode45函数,求出t、x、y[t,x]=ode45(@ode_fun,tspan,x0);tspan时间范围,x0x、y的初值4、利用plot函数,绘出x、y的关系曲线图plot(x(:,1),x(:,2),'r-'); %作图;xlabel('x(t)'),ylabel('y(t)'); %描述x轴,y轴title('x(t)随y(t)变化的相位图');5、按要求编程并调整运行得到如下结果.

把下面的2113程5261序放到一4102个m文件里1653,直内接运行容 function [T,Y] = stest tspan = [0,5]; Dy0 = 0; y0 = 30*pi/180; [T,Y] = ode45(@f,tspan,[y0,Dy0]); plot(T,Y); xlabel('t/s'); ylabel('y,dy') legend('y','dy') function dy = f(t,y) dy = zeros(2,1); dy(1) = y(2); dy(2) = -3*sin(y(1));

你没给初值,我随便设一个[1 1]新建一个脚本文件,输入:clear allclcf=@(t,x)([x(2);-2.6*x(2)^3+0.6*x(2)-0.1*x(1)]);[t,x]=ode45(f,[0 1],[1 1]);plot(t,x(:,1),t,x(:,2))xlabel('t')legend('x','y')

使用ode系列的函数求解微分方程组.d = @(t,x)[x(2); -x(1).^2-2*x(1)-x(2)/2]; [t x] = ode45(d, [0 5.1], [0; 4]); % 这个tfinal自己调出来的一个近似值 plot(x(:,1), x(:,2))

初始条件不完整, 函数图像是无法确定的.以下仅针对y在t趋于负无穷时有限,给出函数图像,

把下面的程序放到一个m文件里,直接运行 function [t,y] = stest tspan = [0,5]; dy0 = 0; y0 = 30*pi/180; [t,y] = ode45(@f,tspan,[y0,dy0]); plot(t,y); xlabel('t/s'); ylabel('y,dy') legend('y','dy') function dy = f(t,y) dy = zeros(2,1); dy(1) = y(2); dy(2) = -3*sin(y(1));

图片:程序: tspan=[0 20]; y0=[0.5 0 0.5 ]; y1=[0 0 1]; [t0,Y0]=ode45('dif_dsdedi',tspan,y0) [t1,Y1]=ode45('dif_dsdedi',tspan,y1) figure subplot(2,1,1) plot(t0,Y0(:,1),'-b',t0,Y0(:,2),'-g',t0,Y0(:,3),'-k') grid on legend('s','e','i') xlabel('t-time') title('y0=[0.5 0 0.

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

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