第十二讲 matlab在理工课程中的应用.ppt

Matlab应用在中文文献检索情况 题名 “matlab” “应用” 检索1417篇 题名：“ matlab” “ 应用 ” “机械” 24篇 建模： 两个振动相加： 两个同方向的振动分别为： 用三角函数关系，可求出振动方程 t=0:0.001:10;%10秒分10000个点 a1=input(振幅1＝);w1=input(相位1＝); a2=input(振幅2＝);w2=input(相位2＝); y1=a1*sin(w1*t); y2=a2*sin(w2*t); y=y1+y2; subplot(3,1,1),plot(t,y1),ylabel(y1) subplot(3,1,2),plot(t,y2),ylabel(y2) subplot(3,1,3),plot(t,y),ylabel(y),xlabel(t) pause,sound(y1); pause(2),sound(y2); pause(2),sound(y),pause subplot(1,1,1) 振幅1＝1.2 相位1＝300 振幅2＝1.5 相位2＝290 6.7 光学 [例6-7-1]两点（双缝）光干涉图案 单色光通过两个窄缝射向屏幕，相当于位置不同的两个同频同相光源向屏幕照射的叠和。由于到达屏幕各点的距离不同，引起相位差，有的点加强，有的点抵消，造成干涉现象。 y ys O z L2 L1 S2 S1 光缝 屏幕 双缝干涉示意图 振幅 光强 两光源到屏幕的距离 光程差 相位差 clc Lambda = input(输入光的波长(单位为 nm): (书上取500) ); Lambda = Lambda * 1e-9; % 将nm换 变为 m d = input(输入两个缝的间距 (单位为 mm): (书上取2) ); d = d * 0.001; % 将mm 变换为 m Z = input(输入缝到屏的距离 (单位为 m): (书上取1) ); yMax = 5*Lambda*Z/d; xs = yMax; % 设定图案的y,x 向范围 Ny=101;ys = linspace(-yMax,yMax,Ny); % y方向分成101点 for i=1:Ny % 对屏上全部点进行循环计算 % 计算第一和第二个光源到屏上各点的距离 L1 = sqrt((ys(i)-d/2).^2 + Z^2 ); L2 = sqrt((ys(i)+d/2).^2 + Z^2 ); Phi = 2*pi*(L2-L1)/Lambda; % 从距离差计算相位差 B(i,:) = 4*cos(Phi/2).^2; % 计算该点光强（设两束光强相同) end % 在屏上画出图象 % clf; figure(gcf); % 清图形窗，将它移到前面 NCLevels = 255; % 确定用的灰度等级 % 定标：使最大光强(4.0)对应于最大灰度级（白色） Br = (B/4.0) * NCLevels; subplot(1,4,1),image(xs,ys,Br); % 画图象 colormap(gray(NCLevels)); % 用灰度级颜色图 subplot(1,4,2), plot(B(:),ys) % 画出沿y 向的光强变化曲线 波长(m)＝0.0000005 光缝距离(m)＝0.002 光栅到屏幕距离(m)＝1 双缝干涉条纹及光强分布 第七章 在力学、机械中的应用 例7-1-4四连杆机构如图所示： 输入杆L1的转角 求输出杆L3转角随时间的变化规律，并求其角速度和角加速度。 建模： 四连杆机构的运动方程X和Y方向的长度关系确定为 （1） （2） 从上述两个方程中削去 ，便可划成一个只包括 和 的方程，给定 ，可求出满足此方程的 由（2）式得 将（1）中的 代以 得出 求得 、 和 后，就可以根据杆1的角速度求出杆3的角速度 为了求能使f=0的 值，利用fzero函数，把f=f（ ） 定义为函数文件ex714f.m 利用global命令定义了杆件的长度参数 function y=ex714f(x) global L0 L1 L2 L3 th