matlab常微分求解

plot等函数，初始条件为 ， 。求 在区间[5,10]的数值解，并画出 、 的图形。[提示：请使用ode45已知描述某动态电路的微分方程为 并画出 x(t)的图形。[提示已知描述某动态电路的微分方程为10x"-4x'+7x=10sin(t) ，初始条件为x(o)=4,x'(0)=9 。求 在区间[5,10]的数值解

9931 7.8092 9.8193
14:10],x)
dx1=x(2);
dx2=(10*sin(t)-7*x(1)+4*x(2))/10.6301
13.0356 -12.9294
15.4885
-2;
f=[dx1;dx2],'x(0)=4.7416 -15.3169
15.1502
14.8617
11.4256
13.3185 -9.1938 3.0449 -4.7628 5.8506
4.0000 9.2128
10.1156
15:,1))

%【2】解析解验证
X=dsolve('10*D2x-4*Dx+7*x=10*sin(t)'.1063 -16.8962
8.3521 8.0571
6.9032 9.7907 4;

结果：（图形运行程序后可见）
t =

0
0.1000
0.2000
0.3000
0.4000
0.5000
0.6000
0.7000
0.8000
0.9000
1.0000
1.1000
1.2000
1.3000
1.4000
1.5000
1.6000
1.7000
1.8000
1.9000
2.0000
2.1000
2.2000
2.3000
2.4000
2.5000
2.6000
2.7000
2.8000
2.9000
3.0000
3.1000
3.2000
3.3000
3.4000
3.5000
3.6000
3.7000
3.8000
3.9000
4.0000
4.1000
4.2000
4.3000
4.4000
4.5000
4.6000
4.7000
4.8000
4.9000
5.0000
5.1000
5.2000
5.3000
5.4000
5.5000
5.6000
5.7000
5.8000
5.9000
6.0000
6.1000
6.2000
6.3000
6.4000
6.5000
6.6000
6.7000
6.8000
6.9000
7.0000
7.1000
7.2000
7.3000
7.4000
7.5000
7.6000
7.7000
7.8000
7.9000
8.0000
8.1000
8.2000
8.3000
8.4000
8.5000
8.6000
8.7000
8.8000
8.9000
9.0000
9.1000
9.2000
9.3000
9.4000
9.5000
9.6000
9.7000
9.8000
9.9000
10.0000

x =

4.1993
15:0.1918
11.4298 -18,[4 9])
plot(t,x(.9684
13.6083 8.1964
-4.2645 -19.8049
-6.2700 -20.3045
-8.3199 -20.6865
-10.4022 -20.9427
-12.5041 -21.0655
-14.6118 -21.0480
-16.7108 -20.8839
-18.7855 -20.5673
-20.8201 -20.0936
-22.7998 -19.4593
-24.7076 -18.6606
-26.5267 -17.6947
-28.2403 -16.5604
-29.8319 -15.2580
-31.2850 -13.7890
-32.5835 -12.1566
-33.7112 -10.3652
-34.6522 -8.4207
-35.3915 -6.3305
-35.9163 -4.1028
-36.2118 -1.7466
-36.2650 0.7274
-36.0642 3.3072
-35.5997 5.9791
-34.8634 8.7283
-33.8489 11.5381
-32.5514 14.3907
-30.9681 17.2665
-29.0980 20.1454
-26.9424 23.0079
-24.5029 25.8321
-21.7833 28.5957
-18.7900 31.2759
-15.5324 33.8501
-12.0223 36.2951
-8.2742 38.5875
-4.3053 40.7040
-0.1355 42.6205
4.2131 44.3144
8.7183 45.7653
13.3546 46.9504
18.0950 47.8484
22.9103 48.4398
27.7696 48.7070
32.6401 48.6344
37.4871 48.2082
42.2741 47.4166
46.9627 46.2495
51.5142 44.6996
55.8926 42.7631
60.0572 40.4351
63.9681 37.7142
67.5860 34.6026
70.8727 31.1059
73.7911 27.2332
76.3048 22.9971
78.3785 18.4139
79.9779 13.5032
81.0701 8.2876
81.6264 2.7902
81.6199 -2.9613
81.0265 -8.9359
79.8261 -15.0997.0050
7.7288
9.7105 -2.6799
7.0712
15.6301
15;)
figure
ezplot(X.4198
15.0000
4.7591 7.5392
9.8624
3.2464 -11.3798 -17.7118 -13.4471
12.4900 8.3413 -10.8763 -1.1502 6.1886 8.1.8124
1.9177 -6.3144 -19.7832
15.4626
14.0019
12.8123
13.7496 1.1772 -8.7867
6.8820
15.8914 0.9352 -0,Dx(0)=9'.5014
10.3737
10.5151 2.6910
-0.4344 -3,[0 10])

function f=fun(t.0547
5.3116 5.2788 -14.7128 9function hhh
clc;clear;
%【1】数值解
[t,x]=ode45(@fun,[0.5399 -5.4802 6