matlab工程应用基础50_449.pptVIP

Matlab 工 程应用基 础 5.0 第 5 章 综合实例 参数 ω n ＝ 10 ， x 0 ＝ 1 ， v 0 ＝ 0 ，计算的终止时间 t =2 。试求 ξ 从 0.1 到 1 运动方程的解，并画出波形 。 2. MATLAB 编程 编写 M 文件 ex1.m % 首先清空 MATLAB 的工作空间 clear; % 给定初值 wn=10; tf=2; x0=1; v0=0; % 计算不同的 ξ 值所对应的振型 for j=1:10; eta(j)=0.1*j; wd(j)=wn*sqrt(1-eta(j)^2); % 求振幅 A a=sqrt((wn*x0*eta(j)+v0)^2+(x0*wd(j))^2)/wd(j); 第 5 章 综合实例 % 求相位角 phi=atan2(wd(j)*x0,v0+eta(j)*wn*x0); % 设定自变量数组 t t=0:tf/1000:tf; % 求过渡过程 x(j,:)=a*exp(-eta(j)*wn*t).*sin(wd(j)*t+phi); end % 在同一个图形窗口中绘制不同的 ξ 值所对应的振型 plot(t,x(1,:),t,x(2,:),t,x(3,:),t,x(4,:),... t,x(5,:),t,x(6,:),t,x(7,:),t,x(8,:),... t,x(9,:),t,x(10,:)) grid on % 新建一个图形窗口，绘制三维网格图 figure mesh(x) 第 5 章 综合实例 如果改变初始条件令 x 0 ＝ 0 ， v 0 ＝ 1 ，其运动曲线实际上就是 系统的脉冲过渡函数 。 第 5 章 综合实例 二、气体分子运动的麦克斯韦分布曲线 通过本例说明如何 用复杂的数学公式绘制曲线 。 利用气体分子运动的麦克斯韦速度分布律，求氯分子运动 的 速度分布曲线 ，并讨论 温度 T 及 分子量 mu 对速度分布曲 线的影响。 1. 建立计算模型 麦克斯韦速度分布律为 : 3 2 2 ? ? ? mv ? m ? 2 ? ? f ? 4 ? p ? ? ? v ? ex ? 2 ? ? k T k T ? 2 ? ? ? 其中， m --- 分子质量 , m = mu / NA , mu --- 分子量 , NA --- 阿伏加德罗数 k --- 波尔茨曼常数 T ---- 气体的绝对温度 v ---- 分子速度 第 5 章 综合实例 为研究单个参数的影响，先把麦克斯韦分布律编为一个 函数子程序 ，以便重复调用，同时将 常数项 也放在子程 序中。 需要强调的是： 子程序不得与主程序放在同一个 M 文件中 ，只能将 子程序单独做成 M 文件，并放在与主程序同一个工作路径中。 2. MATLAB 编程 ? 首先建立计算麦克斯韦分布律的子程序 mxw.m function f=mxw(T,mu,v) % The subfunction mxw.m of ex2 利用麦克斯韦速度分布律求分子的速度分布曲线的 子程序 %mu 、 v 、 T 分别是分子量、分子速度和气体的绝对温度 k=1.381*10^(-23); % 波尔茨曼常数 NA=6.022*10^23; % 阿伏加德罗数 m=mu/NA % 分子质量 f=4*pi*((m/2*pi*k*T)).^(3/2) .*v.*v.*exp(-m*v.^2./(2*k*T)); 第 5 章 综合实例 ? 编写主程序 ex2.m T=300; mu=28e-3; % 给出 T 和 mu 的值 v=0:1500; % 调出自变量数组 y=mxw(T,mu,v); % 调用子程序 plot(v,y, r) % 绘制分布曲线 hold on % 为了看出不同的 T 和 mu 对曲线形状的 影响，再次给定 T 和 mu