用matlab求解矩阵微分方程
答:微分方程K = [k1 -k1 0; -k1 k1+k2 -k2; 0 -k2 k2+k3];C = [c1 -c1 0; -c1 c1+c2 -c2; 0 -c2 c2+c3]dy = @(t,y)[y(4:6); (-K*y(1:3)-C*y(4:6)+[Qy;0;0])./[m1;m2;m3]];初始条件(y4~y6 为 y1'~y3')y0 = zeros(6,1);求解,直接绘图[t,...
答:第一种方法:利用dsolve函数求微分方程的符号解(通解):对于一些不是很难,要求出通解的微分方程,用dsolve函数求解。1、 打开Matlab软件-->点击新建脚本菜单,新建一个脚本文件用于编写微分方程求解程序。2、 输入微分方程求解程序-->点击保存-->点击运行。3、在matlab的命令窗口即可看到求解结果,是...
答:用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(...
答:可以用数值解法 设置初值 步长 进行迭代 直到达到要求
答:先分析,先把矩阵算一下,化简整理下 x‘-4t^2x=0 其中x是t的函数 编程如下;dsolve('dy-4*t^2*y=0')因为matlab中默认的dy表示y对t的一阶导,d2y表示y对t的二阶导,没有dx表示x对t的一阶导 所以会这样编程 运行结果:c2*exp((4*t^3)/3)+c1 ...
答: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外加载...
答:ode求解微分方程的求解函数odefun_wffc这个函数的要求是输入变量是列向量,输出变量也是列向量;初值y0也为列向量;上述说的所有列向量的行数都必须等于你需要求解微分方程未知量的个数。你把微分方程重新发出来,帮你看看。
答:如何用matlab求解二阶微分方程,对于一般的微分方程(组)可以用dsolve()函数求得其解析解或数值解,对于较复杂的微分方程(组)可以用ode45()函数求得其数值解。例如:微分方程,用dsolve和ode45计算t=0.1,0.2时y的值。y = dsolve('D2y+0.5*Dy+2*y =0','Dy(0)=0,y(0)=1')t=0....
答:用dsolve()函数,可以求得其微分方程(组)的解析解。>> syms M C K F y(x)>> D2x=diff(y,2);Dx=diff(y,1);>> y=dsolve(M*D2x+C*Dx+K*x==F)y =C1 + C2*exp(-(C*x)/M) - (K*x^2)/(2*C) + (x*(C*F + K*M))/C^2 式中:C、C1 、 C2为微分方程(组...
网友评论:
通柄17869313318:
Matlab如何求解矩阵微分方程组的符号解例如:dX/dt=HX,X为列向量,H为矩阵H=[2*c,A,0;0,g,0;A,0,a],c,a,g,A为符号常量,初值X(0)=[0,1,0]',求方程组的解X? -
34818甫师
:[答案] 使用dsolve函数 dsolve('Dx=2*c*x+A*y','Dy=g*y','Dz=A*x+A*z','x(0)=0','y(0)=1','z(0)=0') 就是把方程组写出三个微分方程进行计算 祝你学习愉快!
通柄17869313318:
matlab 矩阵间的微分如何求? -
34818甫师
: 这实际上是求Jacobian矩阵 设有向量 v(t)=[v1(t),v2(t),……,vm(t)]' 则有向量函数 f(v)=[f1(v),f2(v),……,fn(v)]' 则有f(v)的obian矩阵: Jacobian=[df1/dv1…df1/dvm ……dfn/dv1…dfn/dvm] 在MATLAB中用函数命令jacobian()来实现.
通柄17869313318:
矩阵微分方程如何用MATLAB求解 -
34818甫师
: 可以用数值解法 设置初值 步长 进行迭代 直到达到要求
通柄17869313318:
Matlab如何求解矩阵微分方程组的符号解 -
34818甫师
: 使用dsolve函数 dsolve('Dx=2*c*x+A*y','Dy=g*y','Dz=A*x+A*z','x(0)=0','y(0)=1','z(0)=0') 就是把方程组写出三个微分方程进行计算 祝你学习愉快!
通柄17869313318:
如何使用MATLAB求解微分方程 -
34818甫师
: 可以用desolve方程:r=dsolve('eqn1','eqn2',...,'cond1','cond2',...,'var').eqn1表示第1个微分方程,cond1表示第1个初始条件,var表示微分方程中的自变量,默认为t.
通柄17869313318:
我论文需要用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画出该微分方程的图像 -
34818甫师
:[答案] 用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...
通柄17869313318:
如何用matlab求解微分方程组 -
34818甫师
: 这说明此微分方程没有显示解,应用数值方法求解,例如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),'.') 结果
通柄17869313318:
matlab求解微分方程 急用 -
34818甫师
: matlab中[]和{}有特殊用法,所以上面的有点问题.可改为:z1=dsolve('Dy=1-a*(1+(-b*y)^c)^(e/2)/f*(1-(-b*y)^(c-1)*(1+(-b*y)^c)^(-e))^2') 另外,给你一个例子看吧:a=1; dsolve('Dy=a*t','t') 得到的结果是:ans = 1/2*a*t^2+C1 从这个例子中,可以看出虽然上面a=1,但是dsolve中不会将其用a的值代入.
通柄17869313318:
matlab 求解微分方程 -
34818甫师
: 先分析,先把矩阵算一下,化简整理下 x'-4t^2x=0 其中x是t的函数 编程如下; dsolve('Dy-4*t^2*y=0') 因为matlab中默认的Dy表示y对t的一阶导,D2y表示y对t的二阶导,没有Dx表示x对t的一阶导 所以会这样编程 运行结果:C2*exp((4*t^3)/3)+C1
通柄17869313318:
MATLAB中怎么解微分方程 -
34818甫师
: 好像没什么特别好的办法: >> syms t q y; q=4*t^3;>> syms Dy>> eq=Dy-(q-2*t);>> y=dsolve(char(subs(eq)),'y(0)=5','t') y = t^4-t^2+5 也就是要先移项把微分方程写成等号右端为零的形式(需额外定义Dy),然后subs代入,再求解.