matlab实习八.docVIP

任 务：_ 0 （任务号_组号）___ 课程名称： MATLAB程序设计及应用 　 课程代码： 999345 实习周数：　　　　　　　一周　　　　 　　　 学　　分： 实习单位：　　　物理科学与技术学院　　　　　　 实习地点：图书馆　　　　　　　　　　 实习时间：　　2012.7.21.13:00　 至 15:00 实习的目的和任务： 动画制作： 1. 通过Help学习moviein, getframe, movie函数。 2. 阅读有关水波干涉的分析的文档“水波干涉的动画演示.doc”,并调试wave.m。请认真阅 读wave.m的代码并对每行给出注释，并保存。 3. 根据所学内容编写一个程序，演示多个振动源(2个)的水波干涉现象。 4．2个文件（2个m文件）连同实习报告压缩打包，文件名为东08.rar。 实习主要内容：（仿真结果及分析，含源程序＋注释＋图） %水波干涉的动画演示 a=4;w=4;A0=0.3;v=1; %定义参量 [x,y]=meshgrid([-4*pi:pi/20:4*pi],[-4*pi:pi/20:4*pi]); %将给定区域划分为平面网格 r1=sqrt((x-a).^2+y.^2)+5; %计算第一个波源的波动函数 c1=A0./r1; %定义c1 r2=sqrt((x+a).^2+y.^2)+5; %计算第二个波源的震动函数 c2=A0./r2; %定义c2 t=0:0.1:10; %确定时间范围 for j=1:41 %for循环 z=c1.*cos(w*(t(j)-r1/v))+c2.*cos(w*(t(j)-r2/v)); %干涉函数 surf(x,y,z); %在给定区域画出三维曲面 colormap(gray); %给曲面定色为灰色 axis equal; %设定各坐标轴刻度相同 shading interp; %以插值形式为图形的像点着色 view([180 90]); %定义视角为方位角180，仰角90 axis([-15 15 -15 15 -0.1 0.1]) %定义坐标轴取值范围 axis off %取消对坐标轴的限制 m(j)=getframe; %返回一个movie frame, end %循环结束 movie(m,10) %播放m文件中图形10次 %??2¨?éé?μ??ˉ?-?Yê? a=4;w=4;A0=0.3;v=1; [x,y]=meshgrid([-4*pi:pi/20:4*pi],[-4*pi:pi/20:4*pi]); r1=sqrt((x-a).^2+y.^2)+5; c1=A0./r1; r2=sqrt((x+a).^2+y.^2)+5; c2=A0./r2; r3=sqrt((x-2*a).^2+y.^2)+5; c3=A0./r3; r4=sqrt((x+2*a).^2+y.^2)+5; c4=A0./r4; t=0:0.1:20; for j=1:41 z=c1.*cos(w*(t(j)-r1/v))+c2.*cos(w*(t(j)-r2/v))... +c3.*cos(w*(t(j)-r3/v))+c4.*cos(w*(t(j)-r4/v)); surf(x,y,z); colormap(gray);