振动分析MATLAB实例

wykdya

请问哪位大侠有个振动分析的MATLAB实例啊，我需要结合自己课题写个MATLAB 程序，但是我还没开始接触课题，但又不想随便找个程序交上去，所以来问下哪位有，不要特别复杂的，最好有具体的流程，题目的背景什么的，比如一般的一个二自由度的无阻尼振动或者一个弹性杆的固有振型之类的，小弟万分感激

ChaChing

试着自行搜索下吧! 如此学的多:@)
或待高人路过

happy

三自由度系统振动分析的实例

1. %Calculate the natural frequencies,the normalized mode shapes and the FRFS
   
2. clear
   
3. close
   
4. a=0.001;b=0.001;
   
5. M=[100,0,0;0,100,0;0,0,100];
   
6. K1=[2,-1,0;-1,2,-1;0,-1,2];
   
7. K=1.0e+5*K1;
   
8. C=a*M+b*K;           %比例阻尼
   
9. [V,D]=eig(K,M);
   
10. for j=1:1:3
    
11.     omega(j)=sqrt(D(j,j));
    
12.     for k=1:1:3
    
13.         Y(k,j)=V(k,j)/V(3,j);
    
14.     end
    
15. end
    
16. omega
    
17. Y
    
18. w=0:0.1:70;                   %频率取值范围和步长
    
19. H=zeros(701,3,3);
    
20. for j=1:1:701
    
21.     Z=K-w(j)^2*M+w(j)*C*i;    %动刚度矩阵
    
22.     HH=inv(Z);                %频响函数矩阵        
    
23.     H(j,:,:)=HH;
    
24. end;
    
25. Hr=real(H);Hi=imag(H);      
    
26. Habs=abs(H);Hang=(180/pi)*angle(H);
    
27. figure(1)
    
28. subplot(2,1,1),plot(w,Hr(:,1,1)),grid on
    
29. title('实频曲线'),xlabel('\omega'),ylabel('Hr')
    
30. subplot(2,1,2),plot(w,Hi(:,1,1)),grid on
    
31. title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
    
32. figure(2)
    
33. subplot(2,1,1),plot(w,Habs(:,1,1)),grid on
    
34. title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
    
35. subplot(2,1,2),plot(w,Hang(:,1,1)),grid on
    
36. title('相频曲线'),xlabel('\omega'),ylabel('\psi')
    
37. figure(3)
    
38. plot(Hr(:,1,1),Hi(:,1,1)),grid on
    
39. title('乃奎斯特圆'),xlabel('Re(H)'),ylabel('Im(H)')
    
40. %画跨点 的频响
    
41. figure(4)
    
42. subplot(2,1,1),plot(w,Hr(:,2,1)),grid on
    
43. title('实频曲线'),xlabel('\omega'),ylabel('Hr')
    
44. subplot(2,1,2),plot(w,Hi(:,2,1)),grid on
    
45. title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
    
46. figure(5)
    
47. subplot(2,1,1),plot(w,Habs(:,2,1)),grid on
    
48. title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
    
49. subplot(2,1,2),plot(w,Hang(:,2,1)),grid on
    
50. title('相频曲线'),xlabel('\omega'),ylabel('\psi')
    
51. figure(6)
    
52. plot(Hr(:,2,1),Hi(:,2,1)),grid on
    
53. title('乃奎斯特圆'),xlabel('Re(H)'),ylabel('Im(H)')
    
54. %画跨点 的频响
    
55. figure(7)
    
56. subplot(2,1,1),plot(w,Hr(:,3,1)),grid on
    
57. title('实频曲线'),xlabel('\omega'),ylabel('Hr')
    
58. subplot(2,1,2),plot(w,Hi(:,3,1)),grid on
    
59. title('虚频曲线'),xlabel('\omega'),ylabel('Hi')
    
60. figure(8)
    
61. subplot(2,1,1),plot(w,Habs(:,3,1)),grid on
    
62. title('幅频曲线'),xlabel('\omega'),ylabel('|H|')
    
63. subplot(2,1,2),plot(w,Hang(:,3,1)),grid on
    
64. title('相频曲线'),xlabel('\omega'),ylabel('\psi')
    
65. figure(9)
    
66. plot(Hr(:,3,1),Hi(:,3,1)),grid on
    
67. title('乃奎斯特圆'),xlabel('Re(H)'),ylabel('Im(H)')

复制代码

happy

单自由度系统阻尼振动实例
http://forum.vibunion.com/thread-15380-1-1.html

想要流程图有点过分了，除了教师上课讲解需要，一般人对于这类问题是不会做流程图的

ChaChing

应该是要交报告吧:@)

wykdya

回复 ChaChing 的帖子

是啊，哎

zhouyang664

我介意以后不要直接要代码或者流程图什么的！

xiwang001

这方面的东西应该是很多的，想做深就很难了，现在就感觉出现瓶颈了

ChaChing

zhouyang664 发表于 2010-10-28 22:44

                               
登录/注册后可看大图

我介意以后不要直接要代码或者流程图什么的！

针对可能是作业或报告的帖, 以前我的想法与作法是, 不马上回覆或仅仅使用提示!:@)
当然原因除了个人水平/时间有限外:@L, 个人不想扼杀了同学们思考空间, 不小心当了帮凶
当然若是有心学习的话, 将尽量讨论, 不论是多麽基础/无厘头的问题
总之, 个人以为如此才能一代超越一代, 不然怎有未来...
不过每个人的看法不同, 依循自然不同:@)

ChaChing

回复 xiwang001 的帖子

什麽问题? 个人不一定懂, 但所知高手如云
建议可以整理贴出, 大家一起讨论!:@)

lmf001

学习一下！

pengcheng5251

回复 ChaChing 的帖子

顶你！！！！！！！！

M一蓑烟雨

学习中，好难啊。。。

右手在忙

9楼说的好  顶

aaronlbx

学习一下