matlab求解多元微分方程组
答：用matlab求解微分方程组并画图，其实现过程：1、首先应根据已知微分方程组，编写其函数 dy(1)=-0.10982*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3));dy(2)=1441.5074*y(1)*(1-y(2))^(1/3)/(1-(1-y(2))^(1/3));2、用ode45函数求解其微分方程组的数值解 [t,y]...
答：用数值方法吧，把t=0时的y值输入，即可解出来。function hh [t,Y]=ode45(@myfun,[0:0.1:1],[0 0 0 0])function dy=myfun(t,y)y1=y(1);y2=y(2);y3=y(3);y4=y(4);dy1=y2;dy2=0.7143*(y1*y4*y4-9.81*sin(y3));dy3=y4;dy4=(280.2913*y4*cos(y3)+420.437...
答：)/(4*exp(a*t/4));dy(2)=-(a/4)*(exp(a*t/4))*(y(1)+0.5)+(a/4)*y(2)-y(3)*((exp(a*t/4))^2);dy(3)=4*y(2);end end 但微分方程组似乎是刚性的，不过换用ode15s、ode23s等适合刚性系统的算法效果也不理想（可以调用ode*函数时不返回参数，观察求解的过程）。
答：首先，自定义求解微分方程组的函数，即 function dy=func(t,y)dy=zeros(3,1);dy(1)=10*(y(1)+y(2)); %dx/dt dy(2)=20*y(1)-y(2)-y(1)*y(3); %dy/dt dy(3)=y(1)*y(2)-8*y(3); %dz/dt 其二，确定初值，由于当x(0)=0，y(0)=0，z(0)=0 该系统是无法触发...
答：如何用MATLAB写求解微分方程组的代码？这个代码可以按下列几个方面来写：第一个，根据题主提供的微分方程组，创建其自定义函数，即内容为 dydx=ode45_fun (t,z)x=z(1);y=z(2);dydx1=r1*x*(1-x/K1).*x/(x+e)-q*x*y/(a+x)-m1*x;dydx2=r2*y*(1-y/K2)+e*q*x*y/(a+x)...
答：在Matlab下输入：edit zhidao_feiying.m，然后将下面两行百分号之间的内容，复制进去，保存 function y=zhidao_feiying(t,x)a=1;b=2;c=3;d=4;比方说 f=a*x+y-z;g=b*sin(x*y)-c*cos(z);h=d*y-a*x;注意x用x(1)代，y用x(2)代,z用x(3)代 f=a*x(1)+x(2)-x(3);g...
答：这说明此微分方程没有显示解，应用数值方法求解，例如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)...
答：Simulink是matlab下基于框图的用于复杂动态系统建模的工具。今以 y"+4y'+4y=cos2t（0≤t≤10），y(0)=1，y'(0)=-2 二阶微分方程为例，用Simulink求解的步骤如下：第一步、分析微分方程的特点确定模块 我们将y"+4y'+4y=cos2t方程改写成y"=cos2t-4y'-4y，从右边可知，式子是由cos2t外加载...
答：如何用matlab求解微分方程组的解析解?为了说明问题。特列举一个例子。例：x"(t)+2x'(t)=x(t)+2y(t)-exp(-t)，y'(t)=4x(t)+3y(t)+4exp(-t)对于线性微分方程组，一般可以用dsolve()函数求解直接得到其解析解，其计算过程如下：syms t x(t) y(t) %声明变量 D2x=diff(x,2); %求...
答：dsolve()求解微分方程后，使用ezplot()作图，参考代码： >> y = dsolve('Dy=y*cos(x)/(1+2*y^2)','y(0)=1','x');>> ezplot(y,[-10,10]) % [-10,10]为x区间
网友评论：
柏才13779264187：
用matlab怎么求解多元一次微分方程组 -
14451爱新觉罗卿
： 如果是常微分方程,可以用dsolve函数.该函数可以解单变量常微分方程或者多变量常微分方程组,所以5个变量也不在话下.调用格式如下:[y1,,yN]=dsolve(eqns)...
柏才13779264187：
如何用matlab中的ode45求解多元微分方程 -
14451爱新觉罗卿
： 示例如下: odefun=@(t,y) (y+3*t)/t^2; %定义函数tspan=[0 2.3]; %求解区间 !!!!!!!!!!!!y0=-2; %初值[t,y]=ode45(odefun,tspan,y0);其中求解区间为0~2.3,那么最后结果的y的最后一个y(length(y))就是y(2.3)
柏才13779264187：
我论文需要用MATLAB求解微分方程,Dx= - 0.85*x+0.25*y+0.35*z - 0.8*x^2*y,Dy=0.8*x^2*y - 1.2*y,Dz=0.9*x+0.8*y - 0.3*z用MATLAB画出该微分方程的图像 -
14451爱新觉罗卿
：[答案] 用dsolve求不出解析解,用数值解法吧,最好能提供初始条件,这里假设t=0时,x0=1,y0=2,z0=3.function hh[t,x]=ode23t(@fun,0:1,[1 2 3])plot(t,x)function yhat=fun(t,x)dxdt=-0.85*x(1)+0.25*x...
柏才13779264187：
用matlab命令求微分方程d^2y/dx^2+2*dy/dx+2y=o,满足初始条件y(0)=1,dy/dx(0)=0的解,并绘制出方程解y(t)的时间曲线图 -
14451爱新觉罗卿
：[答案] 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) 祝你学习愉快!
柏才13779264187：
如何使用MATLAB求解微分方程 -
14451爱新觉罗卿
： 可以用desolve方程:r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var').eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t.
柏才13779264187：
matlab里怎么求解多元微分方程,急! -
14451爱新觉罗卿
： ode45函数,看帮助文档.很容易的.
柏才13779264187：
怎么用matlab解微分方程 dy/dx=A - By - Cxy -
14451爱新觉罗卿
：[答案] syms A B C x y; >> dsolve('Dy=A-B*y-C*x*y','x') ans = -1/2*A*erf(1/2*(-2*C)^(1/2)*x-B/(-2*C)^(1/2))*exp(-1/2*x*(2*B+C*x)-1/2*B^2/C)/C*2^(1/2)*(-pi*C)^(1/2)+exp(-1/2*x*(2*B+C*x))*C1
柏才13779264187：
matlab求解微分方程的通解问题解某微分方程:x^2*D2y+x*Dy+(x^2 - 1/2)*y=0,初值:y(pi/2)=2,Dy(pi/2)= - 2/pi,书本标准答案为:ans=2^(1/2)*pi*^(1/2)/x^(1/2)*sin... -
14451爱新觉罗卿
：[答案] syms x y y=sqrt(2*pi/x)*sin(x);y1=diff(y,'x');y2=diff(y,'x',2); h=x^2*y2+x*y1+(x^2-1/2)*y;h=simplify(h) 结果: h = -1/4/x*2^(1/2)*sin(x)*pi^(1/2)/(1/x)^(1/2) 可见不等于0,=y,说明原答案有问题
柏才13779264187：
matlab怎么求解多元偏微分方程,例如有个四个元 -
14451爱新觉罗卿
： 如果是四个一维的可以考虑pdepe
