资源描述：
《自由边界条件下克拉尼图形的模拟》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、课程设计题目：自由边界条件下克拉尼图形的模拟学院：物理工程学院姓名：文美乐专业：应用物理1203班学号：121410020322一维行波t=0:pi/20:4*pi;%0到4pi，每隔Pi/20取一个数x=0:0.1*pi:4*pi;%0到4pi，每次增加0.1pim=moviein(40);%定义一个40列的矩阵m,以存储40帧画面，每列对应一帧画面fori=1:40;y=sin(x+t(i));plot(x,y,'o');m(:,i)=getframe;%依次将40帧画面存入矩阵mendmovie(m,10,20)%将40帧画面演示10次，每秒2
2、0帧画面图形：一维驻波t=0:pi/20:4*pi;%0到4pi，每隔Pi/20取一个数x=0:0.1*pi:4*pi;%0到4pi，每次增加0.1pim=moviein(40);%定义一个40列的矩阵m,以存储40帧画面，每列对应一帧画面fori=1:40;y1=sin(x+t(i));y2=sin(x-t(i));y=y1+y2;%两列行波叠加plot(x,y1,x,y2,x,y);m(:,i)=getframe;%依次将40帧画面存入矩阵mendmovie(m,10,20)%将40帧画面演示10次，每秒20帧画面图形：三维驻波t=[0:pi/2
3、0:4*pi];%0到4*pi，每隔pi/20取一个数x=[0:0.1*pi:4*pi];%0到4*pi,每次增加0.1*piy=x;[X,Y]=meshgrid(x,y);%定义数组x,y，m=moviein(40);%定义一个40列的矩阵m,以存储40帧画面，每列对应一帧画面fori=1:40;Z=sin(X+t(i))+sin(X-t(i));%定义z函数mesh(X,Y,Z);%画三维图m(:,i)=getframe;%依次将40帧画面存入矩阵mendmovie(m,10,20)%将40帧画面演示10次，每秒20帧画面图形：双原子clear;
4、clc%清除窗口A=[1,-1;-1,1];[w,v]=eig(A)%矩阵A的特征值w，特征向量vt=[0:pi/20:2*pi];%0到2pi，每隔Pi/20取一个数z1=w(1)*cos(v(4)*t);%第一个小球的位移z2=w(4)*cos(v(4)*t);%第二个小球的位移figure;holdon;boxon;r=0.1;%小球半径xr1=r*cos(t)+0.2;yr1=r*sin(t)+z1(1);%第一个小球xr2=r*cos(t)-0.2;yr2=r*sin(t)+z2(1);%第二个小球hb1=fill(xr1,yr1,'r')
5、;%绘制第一个小球hb2=fill(xr2,yr2,'y');%绘制第二个小球axis([-17-1.21.2]);%坐标轴范围设置axisequal;%坐标轴属性设置fori=1:length(t)set(hb1,'YData',r*sin(t)+z1(i));%移动第一个小球set(hb2,'YData',r*sin(t)+z2(i));%移动第二个小球pause(0.05);%延时播放end图形：十原子运动clc;clearall;A=[1-100000000;-12-10000000;0-12-1000000;00-12-100000;000
6、-12-10000;0000-12-1000;00000-12-100;000000-12-10;0000000-12-1;00000000-11];[W,V]=eig(A)%矩阵A的特征值w，特征向量vlin=5;w1=(W(1,lin));w2=(W(2,lin));w3=(W(3,lin));w4=(W(4,lin));w5=(W(5,lin));w6=(W(6,lin));w7=(W(7,lin));w8=(W(8,lin));w9=(W(9,lin));w10=(W(10,lin));%振幅wv1=sqrt(V(1,1));v2=sqrt(
7、V(2,2));v3=sqrt(V(3,3));v4=sqrt(V(4,4));v5=sqrt(V(5,5));v6=sqrt(V(6,6));v7=sqrt(V(7,7));v8=sqrt(V(8,8));v9=sqrt(V(9,9));v10=sqrt(V(10,10));%频率vt=linspace(0,2*pi);z1=w1*cos(v10*t);z2=w2*cos(v10*t);z3=w3*cos(v10*t);z4=w4*cos(v10*t);z5=w5*cos(v10*t);z6=w6*cos(v10*t);z7=w7*cos(v10*
8、t);z8=w8*cos(v10*t);z9=w9*cos(v10*t);z10=w10*cos(v10*t
