用matlab求微分方程组
答:第一种方法:利用dsolve函数求微分方程的符号解(通解):对于一些不是很难,要求出通解的微分方程,用dsolve函数求解。1、 打开Matlab软件-->点击新建脚本菜单,新建一个脚本文件用于编写微分方程求解程序。2、 输入微分方程求解程序-->点击保存-->点击运行。3、在matlab的命令窗口即可看到求解结果,是...
答:第一步、分析微分方程的特点确定模块 我们将y"+4y'+4y=cos2t方程改写成y"=cos2t-4y'-4y,从右边可知,式子是由cos2t外加载信号,一阶导数项和一次变量组成。所以我们,可以用下列功能模块来搭建其仿真模型 cos2t:使用Fcn模块和Clock模块;y':使用lntegrator模块;y':使用lntegrator模块;-4:...
答:1、首先得介绍一下,在matlab中解常微分方程有两种方法,一种是符号解法,另一种是数值解法。在本科阶段的微分数学题,基本上可以通过符号解法解决。2、用matlab解决常微分问题的符号解法的关键命令是dslove命令。该命令中可以用D表示微分符号,其中D2表示二阶微分,D3表示三阶微分,以此类推。值得注意的...
答:用matlab求解微分方程可以按下列思路来求解。1、用dsolve函数,求解该微分方程的通解,T表达式 T=dsolve(diff(T)==-k*(T-21.1));2、求当t=0时的T表达式,eq1=C1/10 + 211/10 3、求当t=1时的T表达式,eq2=(C1*exp(-k))/10 + 211/10 4、用solve函数,求解由2和3组成的方程组的解...
答:如何用matlab求解微分方程并画图,可以先用dsolve()或ode()求出其微分方程(组)的解析解或数值解,然后用plot()绘制其图形。例如:解微分方程 y'=y-2t/y,y(0)=1,0<t<4 1、用dsolve()求解,代码及结果如下 2、用ode45()求解,代码及结果如下 3、当然喽,使用dsolve()或ode(...
答:题主给出的微分方程组,可以这样求解:1、根据题意,该微分方程组可以使用ode45函数来求解 2、完善自定义微分方程组的函数,即 function dy = odefun(t,y)a,b,c,d,e,f都为参数 a=5.985;b=4.709;c=6.959;d=6.998;e=6.385;f=0.336; %为解题而设定 dy(1)=y(3)*c-y(1)*d;...
答:第一步:根据已经微分方程组和相关系数,自定义求解微分方程组的函数,其函数名 odefun,其参数为【t,z】第二步:由于未知初始条件,用随机数初定,即z0=rand(1,6)/1000;第三步:确定时间t的范围,如tspan=[0 50];第四步:使用ode45函数,求其数值解,即 [t,z]=ode45(@(t,z)odefun(t,z...
答:首先定义变量符号 syms s x t;使用字符串定义公式 eq1 = 'diff(x,t)= 0.11*s*x/(0.06*x+s)';eq2 = 'diff(s,t)= (-1)*0.029*x-0.234*s/(0.006*x+s)-s/(30*s^2+300*s+0.03)';接下来是求解 比如现在要求解s,下面的代码可以给出s的表达式 solution = solve(eq1,...
答:1、用function函数体,自定义微分方程组函数,其函数格式为 function f = func(t,x)。。。%微分方程组表达式 end 2、确定t的范围及步长,如t=0~10 t0=0; %起始点 tf=10; %终止点 n=100; %等份 h=(tf-t0)/n; %步长 3、确定x的初值,如x0=[0,0];4、调用欧拉方法函数,求解微分...
答:如何用matlab求解微分方程的数值解和解析解?一、微分方程的数值解可以ode函数来求解。其求解方法:1、自定义微分方程组函数,odefun(x,y)2、定义x【0,10】间若干等份的数值,如n=50 3、定义y的初值,即 y0=[2,7];4、使用ode45函数求y(x),z(x)数值解,即 [x,y]=ode45(@odefun,x,y0)...
网友评论:
花苛19320174642:
matlab解微分方程组 -
65651傅胞
: >> s=dsolve('Dy=-2*2.825678*10^-7*x*Dx','22050000*sin(atan(Dy/Dx))/(7056000-7500*t)-1.624=-5.652*10^-7*Dx^2-6.652*10^-1442537.6*cos(atan(Dy/Dx))/(7056000-7500*t)','t') Warning: Explicit solution could not be found. > In dsolve at 197s ...
花苛19320174642:
怎样用matlab求解复杂的微分方程组 -
65651傅胞
: [x y]=dsolve('D2x=(x*x+3*x)*Dy','D2y=-(x*x+3*x)*Dx','x(0)=0','y(0)=0','Dx(0)=a','Dy(0)=a','D2x(0)=a','D2y(0)=a')%%%%%%%%%%%%%%%%%%%%%%%%%%% 设定的X(1)=dx/dt; x(2)=x; x(3)=dy/dt; x(4)=y;%%%%%%% function dx=shuzhijie(t,x) ...
花苛19320174642:
如何用matlab求解微分方程组 -
65651傅胞
: 这说明此微分方程没有显示解,应用数值方法求解,例如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),'.') 结果
花苛19320174642:
MATLAB 微分方程组求解 -
65651傅胞
: 除了有参数,还要有x,y,x',y'的初始值 看你的方程,类似有阻力的抛体运动 这里随便给了个初值,x(0)=0 y(0)=10 x'(0)=10 y(0)'=0 可以理解为物体从10米高的高度,以10米每秒的速度做平抛运动 m=1;k=0.2;g=10;%参数 ts=0:0.01:2;%时间区...
花苛19320174642:
matlab分段微分方程组求解 -
65651傅胞
: 用matlab求解分段微分方程组,可以这样来考虑. 1、建立分段函数,ode_fun1(t,y);ode_fun2(t,y);ode_fun3(t,y) 2、分三次用ode45求解,其【t,x,y】的数值 [t1,y1]=ode45(@ode_fun1,tspan,X0) [t2,y2]=ode45(@ode_fun2,tspan,X0) [t3,y3]=ode45(@ode_fun3,tspan,X0) 3、最后用plot函数绘出微分方程组解的图形.plot(x1,y1,x2,y2,x3,y3) 运行结果如下图.
花苛19320174642:
如何使用MATLAB求解微分方程 -
65651傅胞
: 可以用desolve方程:r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var').eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t.
花苛19320174642:
matlab中怎样写微分方程组 -
65651傅胞
: 以下列例子,来说明怎样写微分方程组.解 1、建立m-文件rigid.m如下: function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2);2、取t0=0,tf=12,输入命令: [T,Y]=ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'*',T,Y(:,3),'+')3、结果如图 图中,y1的图形为实线,y2的图形为“*”线,y3的图形为“+”线.
花苛19320174642:
matlab微分方程组求解方程组:dx/dt= - 2Ax+uydy/dt=2Ax+auz - (A+u)yx+y+z=1matlab求解如下:>>global A,u>>[x,y,z]=dsolve('Dx=( - 2)*A*x+u*y','Dy=2*A*x+2*u*z - ... -
65651傅胞
:[答案] 把第三个式子也微分. >> [x,y,z]=dsolve('Dx=(-2)*A*x+u*y','Dy=2*A*x+2*u*z-(A+u)*y','Dx+Dy+Dz=0','x(0)=1,y(0)=0,z(0)=0') x = u^2/(A^2+2*u*A+u^2)+2*u*A/(A^2+2*u*A+u^2)*exp(-(A+u)*t)+A^2/(A^2+2*u*A+u^2)*exp(-2*(A+u)*t) y = -(-2*u*A^2/(A^2+2*u*A+u^...
花苛19320174642:
如何用MATLAB软件求微分方程?
65651傅胞
: 有符号法和数值法.一般使用dsolve和ode命令,比如: 符号法 dsolve('Df = f + sin(t)') ans = C4*exp(t) - sin(t)/2 - cos(t)/2 数值法 sol = ode45(@vdp1,[0 20],[2 0]);x = linspace(0,20,100);y = deval(sol,x,1); plot(x,y);
花苛19320174642:
怎样用matlab求解复杂的微分方程组d2x/dt2=(x^2+3x)(dy/dt)d2y/dt2= - (x^2+3x)(dx/dt)当t=0时,x(0)=0,y(0)=0(dx/dt)(0)=a,(dy/dt)(0)=a(d2x/dt2)(0)=a,(d2y/dt2)(0)=a... -
65651傅胞
:[答案] [x y]=dsolve('D2x=(x*x+3*x)*Dy','D2y=-(x*x+3*x)*Dx','x(0)=0','y(0)=0','Dx(0)=a','Dy(0)=a','D2x(0)=a','D2y(0)=a')%%%%%%%%%%%%%%%%%%%%%%%%%%%设定的X(1)=dx/dt;x(2)=x;x(3)=dy/dt;x(4)=y;%%%%%%%function dx=...