轴系扭转振动仿真（SIMULINK）

操作步骤
1、创建模型文件 点击MATLAB中的SIMULINK

单击Blank Mode,创建空白模块，进去Simulink编辑环境。
2、打开库文件
选择菜单栏中的View—Library Browser命令，弹出如图所示窗口

3、放置模块
在模块库中，选择Simulink–Commonly Used Block 中的1个常熟模块C公司谈谈、2个增益模块Gain\2个积分模块Intergrator,将其拖动到模型中。
选择Dsp System Toolbox—Sinks库中的滤波器模块Scope，将其拖动到模型中。
在Simulink—Math 中选择Add加法模块，拖动到模型中。
4、仿真模型中参数的设定
设置Gain模块中增益值为400，Gain1中增益值为4000；Constant设置为240、默认积分模块Intergrator参数，Add加法模块设置三个减法。连接模块，结果如图

5、仿真分析
单击工具栏中的Run按钮，弹出Scope对话框，在滤波器中显示分析结果，如图所示

6、转化方程组

7、时间响应曲线与平面曲线
（1）创建函数文件verderpol.m

function[ xn ] = verderpol( t , x )
global mu;
xn=[x(2);0.2*mu-400*x(1)-2600*x(2)];
end`

（2）在命令行窗口输入下面程序`（前三行输入一行点一下空格，最后四行一块输入点空格）

global mu;
mu=1200;y0=[1200;0]
[t,x]=ode45(@verderpol,[0,1200],y0);
subplot(1,2,1);plot(t,x);
title('时间响应曲线')
subplot(1,2,2);plot(x(:,1),x(:,2))
title('平面曲线')

（3）执行程序后，弹出图形界面