完整word版matlab模拟光的反射折射.doc

(完整word版)matlab模拟光的反射折射 (完整word版)matlab模拟光的反射折射 PAGE / NUMPAGES (完整word版)matlab模拟光的反射折射 (3).我们第一假定光从空气进入水中， 此时 。用一次函数图像来表示光的流传路线， 过改变光的入射点来改变光的入射角和入射方向， 并依据反射光、 折射光与入射光的斜率关 系来控制反射角、折射角与入射角关系，设计程序以下：  通 function ZHEFAN; set(gcf,doublebuffer,on); axis([-2,2,-1,1]);hold on;  % 光的折射 光的反射 %设置图片属性，双缓存防备图闪耀 %设置坐标轴范围（ x 轴 -2 到 2， y 轴 -1 到  1） rectangle(Position,[-2,-1,4,1],... FaceColor,[0.1,0.6,0.6]);  %给图形下半部分设置成水的颜色 text(1.4,0.4,air,color,r,fontsize,14) text(1.4,-0.4,water,color,r,fontsize,14)  %注明空气和水，而且设置字体 n=1.5;t=0; text(0.7,0.6,[n= ,num2str(n)],fontsize,14);  %设置折射率 %在图形上注明折射率 HL=rectangle(Position,[-0.02,0.4,0.04,0.2],... FaceColor,[ ]); Hi=plot([0,0],[0,0.4],b); Ho=plot([0,0],[0,-1]); Hr=plot([0,0],[0,1],r); xlabel( 请按空格停止 ,... fontsize,14,color,r); k=1;dt=0.1; while k pause(0.5); s=get(gcf,currentkey); if strcmp(s,space); clc;k=0;  %在图形上画出光源 %设置入射光 %设置折射光 %设置反射光 %设置结束提示 %初始化 k 方便循环，设置入射光变化量 %循环体 % 设置每次图形变换间隔为 0.5 秒 end if abs(t)0.3*pi dt=-dt;  %输入空格停止演示，是图形停止 end  %设置入射光变化的范围 t=t+dt; A=t; Xi=0.4*tan(A); set(Hi,xdata,[Xi,0]);  % 入射光发生改变 Ao=asin(sin(A)/n); Xo=Xi+tan(Ao); set(Ho,xdata,[Xi,Xo]);  %绘制最后折射光 Xr=Xi+tan(A); set(Hr,xdata,[Xi,Xr]);  % 绘制最后反射光 set(HL,Facecolor,rand(1,3)) end figure(gcf);  %显示图形窗口 运转程序，获得运动的图形，即入射角改变，折射光和反射光的变化。获得截图以下：