问答文章1 问答文章501 问答文章1001 问答文章1501 问答文章2001 问答文章2501 问答文章3001 问答文章3501 问答文章4001 问答文章4501 问答文章5001 问答文章5501 问答文章6001 问答文章6501 问答文章7001 问答文章7501 问答文章8001 问答文章8501 问答文章9001 问答文章9501

matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t)是单位阶跃函数,初始状态为0

发布网友 发布时间:2022-04-23 16:51

我来回答

1个回答

热心网友 时间:2023-10-09 19:21

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计算即可。
matlab用三种方法求解二阶微分方程x''+0.2x'=0.4x=0.2u(t),u(t...

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)=...

...解决线性微分方程x''+0.2x'=0.4x=0.2u(t),u(t)某时程函数

y = C2*exp(t*(41^(1/2)/10 - 1/10)) + C3/exp(t*(41^(1/2)/10 + 1/10))

MATLAB微分方程问题 dsolve('D2y-t*Dy=1-t^2','y(0)=0.1,Dy(0)=0.2...

function xdot=exf(t,x)xdot=[x(2);1-t^2+t*x(2)];end 运行结果

取h=0.2,用四阶经典的龙格一库塔方法求解下列初值问题;

以dy/dx=y-2x/y,其中初始条件y(0)=1为例,通过MATLAB编程实现四阶龙格-库塔算法,并将结果与改进的欧拉算法进行对比。这种算法保持了四阶龙格-库塔法精度高的优点,而且数值积分程序计算量小,仿真速度较之一般实时四阶龙格-库塔法可提高约3. 5位。

用R-K方法求解微分方程的初值问题:(取h=0.2)dy/dx=-y+x+1,y(0)=1...

考试可能会要你求x=0.2的值,只要一步迭代就可以了,你必须要记住公式,那里有四个K,一定要搞清楚怎么来的.只有这样你才真正把龙格库塔法学到手,等到哪天没有MATLAB了,要你用C或什么语言来写你才能写得出来,即使不编程,你可以自己手算!R-K是解常微的最经典数值方法,一定要掌握!

设函数f(x,y)=积分上限xy,下限0 sint/(1+t^2)dt求在(0.2)f对x的二

Z'x=【sin(xy)】/【1+(xy)^2】*y★ Z'y=【sin(xy)】/【1+(xy)^2】*x☆ 再由②得到 dz=★dx+☆dy。数学学习方法:1、重视基础。大多数学生都认为学习数学最重要的就是做题,尤其是高年级的同学,天天都埋头写数学卷子,写了好多好多,成绩也没有变好,于是就对数学绝望了。其实,...

1000(1+x)的平方=1440 求x

1000* (1+x)^2=1440 解: (1+x)^2 = 1.44 1+x = 1.2 或 1+x = -1.2 x = 0.2或 x= -2.2

求一篇数值线性代数的课程设计!急!!!(要求含有Matlab程序)

�6�1 线性多步法有四阶阿达姆斯外插公式和内插公式 (三)用Matlab软件求常微分方程的数值解 [t,x]=solver(’f’,ts,x0,options)注意:1、在解n个未知函数的方程组时,x0和x均为n维向量,m-文件中的待解方程组应以x的分量形式写成. 2、使用Matlab软件求数值解时,高阶微分方程必须等价地变换成一阶...

欧姆定律

45.在一段导体两端加2伏特电压时,通过它的电流强度是0.4安培,这段导体的电阻是___欧姆;如果在它两端不加电压,通过它的电流强度是___安培,这段导体的电阻是___欧姆。 46.如图69所示,电源电压U=12伏(保持不变),电阻R1=100欧姆,变阻器R2的最大电阻值是60欧姆,闭合电键K,当滑动片P从R2的左端到右端的过程...

y' +2xy=x为什么不是可分离变量的微分方程?

y' +2xy=x y' = x(1-2y)∫dy/(1-2y)=∫ x dx -(1/2)ln|1-2y| = (1/2)x^2 +C ln|1-2y| = -x^2 -2C 1-2y = C'e^(-x^2)y = [1- C'e^(-x^2) ]/2

声明声明:本网页内容为用户发布,旨在传播知识,不代表本网认同其观点,若有侵权等问题请及时与本网联系,我们将在第一时间删除处理。E-MAIL:11247931@qq.com
幼儿园老师体罚学生会受到什么样的处罚? 幼儿园老师体罚孩子构成犯罪吗?法律对老师体罚孩子怎样规定的? 怎恶把照片做成黑莓8700的壁纸最好是高清的。 为什么苹果手机勿扰模式还有电话打来? 为什么iphone阻止来电还是能打进来? 高分悬赏[北京五中分校]作息时间 有一本bl小说讲一个特种兵穿越,和军人结婚的,其中有一个片段是他在 一台电脑能配2块不同型号的内存条使用吗? MUGEN 谁有会死星的零子人物包?求发. mugen 求一个人物包,宇宙意志大蛇,不是准神的那个,是隔离的。最好... 移动直播和点播业务,对技术的要求有什么区别 MATLAB中我想定义一个单位阶跃函数然后引用绘制一个u(t)构成的函数图,求大神介绍,初学 直播和点播的区别? f(t)=(2-e-2t)u(t)的MATLAB程序怎么写 电脑开机按F8无反应? 电信网络电视直播和点播是什么区别 MATLAB中 function f=u(t) f=(t>0) 是什么意思 笔记本键盘上的F8键怎么找? 用matlab求信号u(t-3) 免费网络电视直播和点播有什么区别 电脑键盘上的F8键的作用是什么? U(t)是单位阶跃函数,请问用matlab怎么画y=t*u(t)的图象。用plot,向量法。求教 键盘上F8是干什么的? MATLAB里怎么画单位阶跃函数图象 输入u(t)怎么不行。。 matlab中,如何定义 阶跃响应u(t)函数,然后又如何调用,初学者求详细点的解释啊 力学之父是谁? U(t)在MATLAB中如何表示 阿基米德:倔强的科学之王生是谁? 阿基米德定律是谁发明的 是不是冬至过完了以后白天的时间就长了? 直播和点播视频有哪些差异? MATLAB 对二阶方程X'-0.3X'+1.4X=0.5U(t)建模仿真,u(t)是单位阶跃函数。 直播系统和点播系统各有什么优劣势? 联想电脑键盘上F8与锁怎么分开使用? 如何用matlab编程使其输出e的-t次方乘以阶跃函数u(t)的信号波形 聚好看上的直播和点播怎么切换? 请问在matlab中怎么写这个程序? 键盘突然F8,F10突然用不了了,怎么样玩奇迹回事? 网络电视机顶盒 网络电视直播机和点播机具体有什么区别 能收多少个台 哪个好用 价格能差多少 u(1-t)如何用u(t)表示 什么是ITV直播、点播和回看? 电热水器内胆坏了可以修吗? 最新PPLIVE界面上的"直播"和"点播"有什么区别 键盘上的F6键和F8键不能正常使用 用matlab求f(t)=e∧(-at)*u(t)的傅里叶变换 pps上直播和点播的区别是什么,为什么点播看不了 怎样利用MATLAB求系统零状态响应 热水器内胆漏水还能修吗? PPLIVE 与PPS里面的点播跟直播有什么不同啊 收费吗? matlab的阶跃函数怎么表示?