振动力学课程设计题目采用MATLAB 对所选的问题进行数值计算和作图，采用高于MATLAB7.4（2007）版本所编写的程序需转换为文本（.txt ）文件, 早于MATLAB7.4（2007）版本所编写的程序可直接采用M 文件传送至QQ ：296637844。

题目如下，其中1，2，3题为必做题，4-38选二题(第一轮:一班01号为第4题, 一班02号为第5题…一班28号为第25题, 二班01号为第26题,…二班17号为第38题, 二班18号为第4题,…二班27号为第13题;第二轮:一班01号为第14题…)。

文件名采用自己的姓名。

考核时间暂定于12月30日。

题目： 1. 编写MA TLAB 程序，根据书本公式（3.1-10）、（3.1-10）作出单自由度系统强迫振动的幅频特性曲线、相频特性曲线。

0.1,0.2,0.3,0.5,0.7,1.0,1.2ς=。

2. 根据书本图4.5-3，分析有阻尼动力减振器的特性。

包括在不同的质量比，频率比，阻尼比条件下结构的响应。

3. 对于图2所示体系，用矩阵迭代法计算其固有频率及振型。

1231,2m m m ===，1230c c c ===，1231,5,8k k k ===，1230,0,0F F F ===，1231,1,1ωωω===。

4. 采用中心差分法计算单自由度体系10105sin(/2)x cx x t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

5. 采用Houbolt 法计算单自由度体系10105sin(/2)x cx x t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

6. 采用Wilson-θ法计算单自由度体系10105sin(/2)x cxx t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

7. 采用Newmark-β法计算单自由度体系10105sin(/2)x cx x t ++= ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

8. 采用中心差分法计算10105sin(/2)2sin()sin(2)x cx x t t t ++=++ ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

9. 采用Houbolt 法计算10105sin(/2)2sin()sin(2)x cxx t t t ++=++ ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

10. 采用Wilson-θ法计算10105sin(/2)2sin()sin(2)x cxx t t t ++=++ ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

11. 采用Newmark-β法计算10105s in (/2)2s in ()s in (2x c x x t t t ++=++ ，当c=3和c=20，000,0x x== 前10s 内的位移，作出其时间位移曲线图。

12. 采用卷积积分法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，分别在()5(),5(),5sin(2)(),(02)F t N t N t N t s =≤≤作用下前10s 内的时间位移曲线。

13. 采用中心差分法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，分别在()5(),5(),5sin(2)(),(02)F t N t N t N t s =≤≤作用下前10s 内的时间位移曲线。

14. 采用Houbolt 法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，分别在()5(),5(),5sin(2)(),(02)F t N t N t N t s =≤≤作用下前10s 内的时间位移曲线。

15. 采用Wilson-θ法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，分别在()5(),5(),5sin(2)(),(02)F t N t N t N t s =≤≤作用下前10s 内的时间位移曲线。

16. 采用Newmark-β法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，分别在()5(),5(),5sin(2)(),(02)F t N t N t N t s =≤≤作用下前10s 内的时间位移曲线。

17. 采用卷积积分法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，在5(),(01)()105(),(12)t N t s F t t N t s ≤≤⎧=⎨-≤≤⎩作用下前10s 内的时间位移曲线。

18. 采用中心差分法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，在5(),(01)()105(),(12)t N t s F t t N t s ≤≤⎧=⎨-≤≤⎩作用下前10s 内的时间位移曲线。

19. 采用Houbolt 法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，在5(),(01)()105(),(12)t N t s F t t N t s ≤≤⎧=⎨-≤≤⎩作用下前10s 内的时间位移曲线。

20. 采用Wilson-θ法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，在5(),(01)()105(),(12)t N t s F t t N t s ≤≤⎧=⎨-≤≤⎩作用下前10s 内的时间位移曲线。

21. 采用Newmark-β法计算单自由度体系m=10kg ，c=3Ns/m ，k=10N/s ，在5(),(01)()105(),(12)t N t s F t t N t s ≤≤⎧=⎨-≤≤⎩作用下前10s 内的时间位移曲线。

22. 求1tan ββ=的前4个解，要求误差小于0.0123. 对于图1所示体系，采用中心差分法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，120,1F F ==，121,1ωω==，11220x xx x ==== 。

图124. 对于图1所示体系，采用Houbolt 法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，120,1F F ==，121,1ωω==，11220x xx x ==== 。

25. 对于图1所示体系，采用Wilson-θ法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，120,1F F ==，121,1ωω==，11220x xx x ==== 。

26. 对于图1所示体系，采用Newmark-β法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，120,1F F ==，121,1ωω==，11220x xx x ==== 。

27. 对于图1所示体系，采用中心差分法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，121,1F F ==，121,2ωω==，11220x xx x ==== 。

28. 对于图1所示体系，采用Houbolt 法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，121,1F F ==，121,2ωω==，11220x xx x ==== 。

29. 对于图1所示体系，采用Wilson-θ法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，121,1F F ==，121,2ωω==，11220x xx x ==== 。

30. 对于图1所示体系，采用Newmark-β法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

121m m ==，121c c ==，121,5k k ==，121,1F F ==，121,2ωω==，11220x xx x ==== 。

31. 对于图2所示体系，采用中心差分法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

1231,2m m m ===，1231c c c ===，1231,5,8k k k ===，1230,1,0F F F ===，1231,1,1ωωω===，1122330x xx x x x ====== 。

图232. 对于图2所示体系，采用Houbolt 法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

1231,2m m m ===，1231c c c ===，1231,5,8k k k ===，1230,1,0F F F ===，1231,1,1ωωω===，1122330x xx x x x ====== 。

33. 对于图2所示体系，采用Wilson-θ法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

1231,2m m m ===，1231c c c ===，1231,5,8k k k ===，1230,1,0F F F ===，1231,1,1ωωω===，1122330x xx x x x ====== 。

34. 对于图2所示体系，采用Newmark-β法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

1231,2m m m ===，1231c c c ===，1231,5,8k k k ===，1230,1,0F F F ===，1231,1,1ωωω===，1122330x xx x x x ====== 。

35. 对于图2所示体系，采用中心差分法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

1231,2m m m ===，1231c c c ===，1231,5,8k k k ===，1232,1,1F F F ===，1231,1,2ωωω===，1122330x xx x x x ====== 。

36. 对于图2所示体系，采用Houbolt 法计算其位移响应，并作出前10s 内每个质点的时间位移曲线。

1231,2m m m ===，1231c c c ===，1231,5,8k k k ===，1232,1,1F F F ===，1231,1,2ωωω===，1122330x xx x x x ====== 。