惯性技术作业Word格式文档下载.docx
《惯性技术作业Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《惯性技术作业Word格式文档下载.docx(14页珍藏版)》请在冰豆网上搜索。
end
利用龙格库达法计算微分方程:
计算时间从0至2秒,设定初值分别为2和1
[t,x]=ode45('
lisongkun'
[0,2],[2,1])
t=
0
0.0017
0.0035
0.0052
0.0069
0.0156
0.0243
0.0330
0.0417
0.0766
0.1116
0.1466
0.1816
0.2300
0.2784
0.3268
0.3753
0.4253
0.4753
0.5253
0.5753
0.6239
0.6726
0.7213
0.7700
0.8200
0.8700
0.9200
0.9700
1.0177
1.0655
1.1133
1.1610
1.2110
1.2610
1.3110
1.3610
1.4096
1.4582
1.5068
1.5554
1.6054
1.6554
1.7054
1.7554
1.8054
1.8554
1.9054
1.9554
1.9666
1.9777
1.9889
2.0000
x=
2.00001.0000
2.00170.9497
2.00330.8994
2.00480.8491
2.00620.7987
2.01210.5464
2.01570.2934
2.01720.0401
2.0164-0.2132
1.9912-1.2291
1.9307-2.2231
1.8361-3.1778
1.7091-4.0762
1.4839-5.1989
1.2085-6.1454
0.8921-6.8830
0.5455-7.3877
0.1684-7.6474
-0.2148-7.6311
-0.5900-7.3386
-0.9441-6.7812
-1.2563-6.0041
-1.5255-5.0219
-1.7423-3.8676
-1.8995-2.5810
-1.9937-1.1677
-2.01590.2874
-1.96521.7318
-1.84353.1138
-1.66524.3305
-1.43225.4044
-1.15186.2995
-0.83356.9872
-0.47117.4613
-0.09187.6661
0.29087.5934
0.66297.2467
1.00206.6597
1.30685.8457
1.56694.8320
1.77363.6534
1.92342.3115
2.00370.8864
2.0115-0.5704
1.9467-2.0068
1.8119-3.3714
1.6117-4.6142
1.3533-5.6896
1.0460-6.5598
0.9720-6.7225
0.8963-6.8732
0.8189-7.0115
0.7401-7.1371
>
plot(x(:
1),x(:
2))
xlabel('
r'
)
ylabel('
dr/dt'
title('
相轨迹图形'
求取响应曲线:
plot(t,x(:
1))
t'
r(t)'
时间响应曲线'
垂直阻尼分析过程:
垂直阻尼时,陀螺罗经的微分方程为:
达因.厘米H=4*
为了分析时计算方便特别的令Md=4*
达因.厘米。
+
functionxs=lisongkun1(t,x)
-x
(2)-14.472*x
(1)];
lisongkun1'
0
0.0034
0.0050
0.0067
0.0151
0.0235
0.0319
0.0403
0.0759
0.1115
0.1472
0.1828
0.2312
0.2797
0.3281
0.3765
0.4265
0.4765
0.5265
0.5765
0.6217
0.6668
0.7120
0.7571
0.8071
0.8571
0.9071
0.9571
1.0023
1.0475
1.0928
1.1380
1.1880
1.2380
1.2880
1.3380
1.3879
1.4379
1.4878
1.5378
1.5878
1.6378
1.6878
1.7378
1.7866
1.8355
1.8843
1.9332
1.9499
1.9833
2.00160.9498
2.00320.8996
2.00470.8495
2.00600.7994
2.01170.5499
2.01530.3018
2.01680.0556
2.0162-0.1887
1.9913-1.1989
1.9314-2.1521
1.8388-3.0330
1.7163-3.8283
1.5076-4.7527
1.2591-5.4763
0.9807-5.9838
0.6831-6.2696
0.3669-6.3338
0.0532-6.1720
-0.2468-5.8006
-0.5235-5.2431
-0.7462-4.6040
-0.9377-3.8607
-1.0936-3.0395
-1.2111-2.1670
-1.2948-1.1732
-1.3287-0.1869
-1.31420.7572
-1.25411.6287
-1.16422.3320
-1.04462.9369
-0.90033.4299
-0.73633.8021
-0.53884.0671
-0.33204.1760
-0.12384.1323
0.07873.9452
0.26843.6291
0.43953.2009
0.58672.6813
0.70602.0928
0.79501.4585
0.85170.8040
0.87550.1535
0.8674-0.4706
0.8303-1.0356
0.7672-1.5388
0.6813-1.9659
0.5765-2.3064
0.5372-2.4018
0.4963-2.4859
0.4542-2.5587
0.4109-2.6200
心得体会及遇到的困难
完成作业过程中所遇到的困难:
1,MATLAB的编写程序过程中遇到了困难。
例如:
要计算14.472乘以X的值。
错误写法:
14.472X正确写法:
14.472*X
2,M函数的编写过程中要注意参数间的对应关系。
3,用软件画图形时,如果大小或比例失调,需要调节比例重新画图并分析。
4,软件的使用遇到不的函数使用方法,需要查阅相关的资料。
通过作业的完成所得到的心得体会:
1,数学模型需要自己设定参数。
2,通过这次作业,我增强了解决问题的能力。
3,通过改变参数的值,可以得到不同的比例不同图形,但是图形的大致走向很相似。
4,通过这次作业,我也锻炼的自己的实践和动手能力。
5,在和同学们的交流中我增强了理解和交流的能力。
总结:
当没有阻尼作用时,相轨迹为圆,时间响应是等幅震荡的过程。
当有阻尼作用时,相轨迹为逐渐收敛于原点,时间响应为逐渐衰减的。
所以,阻尼对于二阶系统的快速性有一定的影响,加快响应的快速性,加快收敛。
同时,对系统的稳定性也有一定的影响。
当没有阻尼是系统的时间响应是等幅震荡的,这是不稳定的。
加上阻尼后系统时间响应的超调量就会逐渐衰减,并最终保持在相应的误差范围之内。