matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0
来源:学生作业帮 编辑:大师作文网作业帮 分类:数学作业 时间:2024/11/11 01:31:16
matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0
x''+0.2x'-0.4x=0.2u(t), x(0)=x'(0)=0.
1.你用Laplace Transform 做吧.
L(x'')+0.2L(x')-0.4L(x)=0.2L(u(t))
(p^2)*X(p)-p*x(0)-x'(0)+0.2(p*X(p)-x(0))-0.4*X(p)=0.2/p,
(p^2+0.2*p-0.4)*X(p)=0.2/p,
X(p)=0.2/(p^3+0.2*(p^2)-0.4p)
然后,你查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).
根据上述一阶微分方程组编写M函数文件DyDt.m
function ydot=DyDt(t,y)
if t>0
u(t)=1
else
u(t)=0
end
ydot=[y(2);0.4*y(1)-0.2*y(2)+0.2*u(t)];
tspan=[0,30];%solution interval
y0=[0;0];
[tt,yy]=ode45(@DyDt,tspan,y0)
3.用差分法吧.
h=t/N, N是你给定的正整数,一般取20吧.ti=it/N.如果t>0
x''(ti)=(x(i+1)-2x(i)+x(i-1))/(h^2)
x'(ti)=(x(i+1)-x(i-1))/2h
x(ti)=x(i)
u(ti)=u(i)=1
然后你写出相对应的差分方程组:
(x(i+1)-2x(i)+x(i-1))/(h^2)+0.2*(x(i+1)-x(i-1))/2h-0.4*x(i)=0.2
化成矩阵式用MATLAB计算即可.
1.你用Laplace Transform 做吧.
L(x'')+0.2L(x')-0.4L(x)=0.2L(u(t))
(p^2)*X(p)-p*x(0)-x'(0)+0.2(p*X(p)-x(0))-0.4*X(p)=0.2/p,
(p^2+0.2*p-0.4)*X(p)=0.2/p,
X(p)=0.2/(p^3+0.2*(p^2)-0.4p)
然后,你查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).
根据上述一阶微分方程组编写M函数文件DyDt.m
function ydot=DyDt(t,y)
if t>0
u(t)=1
else
u(t)=0
end
ydot=[y(2);0.4*y(1)-0.2*y(2)+0.2*u(t)];
tspan=[0,30];%solution interval
y0=[0;0];
[tt,yy]=ode45(@DyDt,tspan,y0)
3.用差分法吧.
h=t/N, N是你给定的正整数,一般取20吧.ti=it/N.如果t>0
x''(ti)=(x(i+1)-2x(i)+x(i-1))/(h^2)
x'(ti)=(x(i+1)-x(i-1))/2h
x(ti)=x(i)
u(ti)=u(i)=1
然后你写出相对应的差分方程组:
(x(i+1)-2x(i)+x(i-1))/(h^2)+0.2*(x(i+1)-x(i-1))/2h-0.4*x(i)=0.2
化成矩阵式用MATLAB计算即可.
MATLAB 对二阶方程X'-0.3X'+1.4X=0.5U(t)建模仿真,u(t)是单位阶跃函数.
u(x,y)为二元函数,x、y为自变量,a(x),b(y)为一元函数,求解微分方程:du(x,y)=a(x)u(x,y)
用Matlab编程求解 二阶微分方程:4*d^2y(t)/dt^2+y(t)=dx(t)/d(t)-0.5x(t)
如何在MATLAB中对二阶方程x ’-0.3x ’+1.4x=0.5u(t)建模仿真,
MATLAB里怎么画单位阶跃函数图象 输入u(t)怎么不行.
matlab du/dt=d(du)/dx^2 x属于(0,1),t属于(0,T]u(0,t)=u(1,t)=0u(x,
MATLAB阶跃函数K=exp^(-0.05*(60*(n-1)+x))u(50-(60*(n-1))),这个式子里面n
用matlab解微分方程,题目是:ax''''''+bx''''+cx''+dx=0,其中x=esinωt,求解ω
matlab已知连续系统函数,求在x=sin(t)u(t)、x=e^(-t)*x(t)的输入下的输出
单位阶跃函数u(t)是功率信号,功率是多少呢?
在matlab中如何用laplace变换求解微分方程u''(t)-u=exp(t),y(0)=0,y'(0)=0
du=2/u^2+1 u(0)=5 自变量为t 用matlab 解次微分方程,且画出u随t的变化,