方程如下
m1a1=-c1(v1-v0)-k1(z1-z0)+c2(v2-v1)+k2(z2-z1)-Fsin(wt);
m2a2=-c2(v2-v1)-k2(z2-z1)+c3(v3-v2)+k3(z3-z2);
m3a3=-c3(v3-v2)-k1(z3-z2);
数据: m1=15;m2=8.8;m3=43.4;
k1=31000;k2=18000;k3=44130;
c1=830;c2=200;c3=1485;
z0=0;v0=0;z1z2z3初始一阶导数为0,v是z的一阶导数,a是z的二阶导数。
f=1000,w=2
怎么得到z1z2z3的数值解?
如何用matlab求解一个二阶常系数微分方程组
第一种方法:利用dsolve函数求微分方程的符号解(通解):对于一些不是很难,要求出通解的微分方程,用dsolve函数求解。1、 打开Matlab软件-->点击新建脚本菜单,新建一个脚本文件用于编写微分方程求解程序。2、 输入微分方程求解程序-->点击保存-->点击运行。3、在matlab的命令窗口即可看到求解结果,是...
如何用matlab求解一个二阶常系数微分方程组
1. 建立M文件 function dy=vdp(t,y)dy=[y(2);1000*(1-y(1)^2)*y(2)-y(1);];2. 调用MATLAB 函数ODE15S [T,Y]=ode15s('vdp',[0 1000],[2 0]);plot(T,Y(:,1));图像为:
如何用matlab求二阶微分方程?
第一步、分析微分方程的特点确定模块 我们将y"+4y'+4y=cos2t方程改写成y"=cos2t-4y'-4y,从右边可知,式子是由cos2t外加载信号,一阶导数项和一次变量组成。所以我们,可以用下列功能模块来搭建其仿真模型 cos2t:使用Fcn模块和Clock模块;y':使用lntegrator模块;y':使用lntegrator模块;-4:...
二阶微分方程组的MATLAB解法
先降阶为一阶微分方程组 y'=z z'=-z-y+sin(t)编制如下函数m文件 function dy=weifen(t,x)dy=zeros(2,1);y=x(1)z=x(2)dy(1)=x(2);dy(2)=sin(t)-x(2)-x(1);然后用ode45解方程 [t,y]=ode45(@weifen,[0 20],[5 6])plot(t,y)就可以画出y和y'的图像 ...
如何用matlab求解一个二阶常系数微分方程组
用命令:dsolve('S','s1','s2',…,'x')其中S 为方程s1,s1,s3,…为初始条件x 为自变量方程S 中用D 表示求导 数D2,D3,…表示二阶三阶等高阶导数初始条件缺省时给出带任意常数 C1,C2,..的通解自变量缺省值为t 也可求解微分方程组 例 1、dsolve('Dy=1+y^2')结果ans =tan(t+C1)2...
用matlab怎样解二阶符号微分方程
联立(1),(3)与(5)得到三个一阶微分方程:x'= - y*exp(-t\/5)+y' * exp(-t\/5)+1; (1)z=y'; (3)z'= -2*sin(t) (5)既然 z=y' ,用z代替等式(1)中的y' 。而且,因为MATLAB要求所有的导数项在左边,改写等式(3)。得到如下的方程组:x'= - y*exp(-...
如何用matlab求解二阶微分方程,以及程序实例
1、首先看一下matlab求解方程的方法,首先指明所解方程的变量,然后指明方程,未知数和限制条件,最后求解方程。2、在matlab命令行窗口中输入symsx[x,params,conds]=solve(sin(x)==1,'ReturnConditions', true),按回车键可以得到方程解。3、转换一下,可以看到sin(x)=1方程的解是如下图所示 。4、...
matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t...
然后,你查Laplace Transform 表即可。或者用MAtlab 求解,代码如下:F=0.2\/(s^3+0.2*s^2-0.4*s);x(t)=ilaplace(F,s,t)然后,x(t)的表达式就出来了。2.令x1=x,x2=x1'=x',把原方程转化成微分方程组:x1'=x2 x2'=0.4*x1-0.2*x2+0.2u(t),(x1(0),x2(0))=(0,0)...
请问老师这个题用matlab怎么做?
题主给出的二阶微分方程组,用matlab求解其精确解和数值解,可以使用dsolve()函数得到其解析值,使用ode()函数得到其数值解。一、使用dsolve()函数求解 第一步,对x(t)、y(t)进行声明变量,即 syms x(t) y(t)第二步,使用diff()函数分别对x(t)、y(t)进行导数计算 Dy=diff(y,t,1...
如何用matlab求解二阶微分方程,以及程序实例
如何用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....