matlab解微分方程的问题
题目
matlab解微分方程的问题
4x”+35.78x’+2000x=F(t)
F(t)=100 0
答案
似乎还缺初始条件啊
function hh
x0=[10 5];
[T1,X1] = ode45(@odefun1,0:0.01:0.1,x0)
[T2,X2] = ode45(@odefun2,0.1:0.01:0.2,[X1(end,1),X1(end,2)])
plot(T1,X1,T2,X2)
function dx=odefun1(t,x)
x1=x(1);x2=x(2);
dx1=x1;
dx2=(1000-35.78*x1-2000*x2)/4;
dx=[dx1;dx2];
function dx=odefun2(t,x)
x1=x(1);x2=x(2);
dx1=x1;
dx2=(1000*(t-0.2)-35.78*x1-2000*x2)/4;
dx=[dx1;dx2];
结果:
T1 =
0
0.0100
0.0200
0.0300
0.0400
0.0500
0.0600
0.0700
0.0800
0.0900
0.1000
X1 =
10.0000 5.0000
10.1005 0.3512
10.2020 0.3181
10.3045 0.3160
10.4081 0.3138
10.5127 0.3124
10.6184 0.3104
10.7251 0.3087
10.8329 0.3066
10.9417 0.3047
11.0517 0.3027
T2 =
0.1000
0.1100
0.1200
0.1300
0.1400
0.1500
0.1600
0.1700
0.1800
0.1900
0.2000
X2 =
11.0517 0.3027
11.1628 -0.2416
11.2750 -0.2423
11.3883 -0.2394
11.5027 -0.2363
11.6183 -0.2335
11.7351 -0.2305
11.8530 -0.2277
11.9722 -0.2247
12.0925 -0.2220
12.2140 -0.2190
举一反三
我想写一篇关于奥巴马的演讲的文章,写哪一篇好呢?为什么好
最新试题
热门考点
- 仿照例句写一个句子.例句:海上的夜是柔和的,是寂静的,是梦幻的.仿句:秋天的雨是
- When I said students are lazy ,I ______to you
- 英语种1.23.4.5的写法
- 请提供一个有关成功之后因自满而失败的例子(外国名人),最好英文,中文也行
- 一个三角形的三个内角度数的比是1:2:3 .这个三角行的三个内角分别是多少度?它是什麽三角形?
- 如图,在平面直角坐标系中点c(-3,0),点A,B分别在x轴,y轴的正半轴上,且满足√OB²-3 +|OA―1|=0.
- 解释一下下面的句子
- 仿句:母亲是一把伞,为我们遮挡人生的风雨.(仿两句)
- Designer,i miss you,I have missed you ,I regret
- 世界上国土面积跨两大洲的是哪些国家