华工电信数学实验数值仿真.docxVIP

数学实验报告实验四 数值仿真学 院： 电子与信息学院 班 级： 2010级信工五班 学 生： 陈耿涛 序 号： 3号 实验内容 1．针对单自由度阻尼系统，（1）在无强迫力的情形下，固定阻尼系数，观察自然频率对振动模态的影响；（2）固定阻尼系数和自然频率，选择强迫项的一种形式，观察稳态增益对振动模态的影响。2．针对无强迫力的单自由度阻尼系统，选择适当的阻尼系数和固有频率，动态模拟弹簧的振动情况，观察运动规律。3．针对单自由度阻尼系统，选择适当的阻尼系数和固有频率，创建弹簧振动的Simulink仿真模型；通过选择不同的输入信号，观察振动曲线变化。实验过程（1）实验资料中给出的参考程序是在固定自然频率时，观察阻尼系数对振动模态的影响，现在为了观察自然频率对振动模态的影响，将阻尼系数固定，选取自然频率w的一系列值，对这些自然频率用数值方法求解微分方程得到相应的数值解，将这些解曲线绘制在同一图形上，编制程序如下：编辑exp4_1函数并保存为exp4_1；代码如下：function exp4_1clear;clf; global w u cx0(1)=1;x0(2)=0;u=0;c=0.5;n=3;tspan=linspace(0,4,100);xx=[];hold on ,xlabel(t),ylabel(x),for i=1:2:10; w=i; [t,x]=ode45(exp4_1f,tspan,x0);%求微分方程的数值解。 xx=[x(:,1),xx]; text(t(10),x(10,1),[\leftarrow w=,num2str(w)],fontsize,15) plot(t,x(:,1)),end 编辑微分方程函数exp4_1f(t,x)保存为exp4_1f；代码如下：function dx=exp4_1f(t,x) %建立微分方程global w u cdx=zeros(2,1);dx(1)=x(2);dx(2)=u-c*w*x(2)-w^2*x(1); 在matlab命令运行窗口中运行exp4_1后得出的曲线如下由绘制得出的图像可知，当阻尼系数固定时随着自然频率的增加，曲线振动的衰减速度越来越快，衰减到同一幅度所用的时间越来越少，物体振动到达平衡点的时间久越短。（2）将阻尼系数和自然频率固定，观察有强迫振动项的时候对振动模态的影响。根据要求，现在将输入的强迫振动项 u为10，改变不同的K值，得到稳态增益对振动模态的影响， 编制程序如下：编辑exp4_2函数并保存为exp4_2；代码如下：function exp4_2clear;clf; global u w kx0(1)=1;x0(2)=0;w=10;n=3;u=10;tspan=linspace(0,4,100);kk=[1,3,5,7];hold on,xlabel(t),ylabel(x),for i=1:length(kk); k=kk(i); [t,x]=ode45(exp4_2f,tspan,x0);%求微分方程的数值解。 text(t(10),x(10,1),[\leftarrow k=,num2str(k)],fontsize,15) plot(t,x(:,1));end编辑微分方程函数exp4_2f(t,x)保存为exp4_2f；代码如下：function dx=exp4_2f(t,x)%建立微分方程global u w kc=0.3;dx=zeros(2,1);dx(1)=x(2);dx(2)=k*u-c*w*x(2)-w^2*x(1);在matlab命令运行窗口中运行exp4_2后得出的曲线如下 由上图可见，图像在前半部分几乎重合，在图像的后半部分随着的不同逐渐出现了分化，各自向偏离x=0的位置振动，即振动的前半部分由振动系统的阻尼系数和自然频率决定，在经过一段时间以后它对振动的影响逐渐趋近于0（即时间长后不起作用），而振动的后半部分主要由强迫振动项决定，它不随时间衰减，最后振动的情况与强迫振动项密切相关。=K*u，并且有K值越大，物体运动稳定后，偏离平衡点的距离越远。对由方程刻画的单自由度阻尼系统做动画仿真，设弹簧的平衡点u=0，且弹簧在水平方向运动，以一小正方形代表振动物体，由于弹簧的方程为，为方便期起见，取=0，可以推导出方程的解为其中A、由初始条件决定，现取A=1，= 0，则弹簧的横坐标的运动规律为编辑exp4_3程序代码如下：animinit(onecart1 Animation) axis([-2 6 -10 10]); hold on; u=0;c=0.1;w=10; xy=[ 0 0 0 0 u u u+1 u+1 u u; -1.2 0 1.2 0 0 1.2 1.2 -1.2 -1.2 0];