大学物理实验6 简谐振动的合成

云南大学软件学院 实验报告
课程： 大学物理实验 学期： 任课教师： 班级： = 学号： == 序号： = 姓名： = 成绩：
实验6 简谐振动的合成
内容一 同方向同频率简谐振动的合成
两个简谐振动的方程为
使用matlab 编写程序，求x1,x2,合振动的波形，讨论相位差对合成振动的影响。

相位差至少讨论4中情况（1.为0°;2.为180°；3.小于180°；4.大于180°），要求所有波形画在同一个figure 中。

()()
⎩⎨
⎧+=+=222111cos cos ϕωϕωt A x t A x
内容二 相互垂直方向同频率简谐振动的合成
两个简谐振动的方程为
使用matlab 编写程序，求x,y,合振动的波形，讨论相位差对合成振动的影响。

相位差至少讨论4中情况（1.为0°;2.为180°；3.小于180°；4.大于180°），要求所有波形画在同一个figure 中。

()()
⎩⎨
⎧+=+=y y x x t A y t A x ϕωϕωcos cos
内容三相互垂直方向不同频率简谐振动的合成(李萨如图形) 使用matlab编写程序，画李萨如图形，要求：
1.至少4种频率比
2.至少8种相位差
3.所有图形画在同一个figure中，添加标注。

如：
clear
t = 0:0.01:4;
Ax = 1;Ay = 3;
w1 = 1; w2 = 1./2;
w3 = 2./3;w4 = 3./4;w5 = 2./5;
m0 = 0;m1 = 0;
m2 = pi./4;m3 = pi./2;
m4 = 3.*pi./4;m5 = pi;
m6 = 5.*pi./4; m7 =3.*pi./2;m8 = 7.*pi./4; x0 = Ax.*cos(2.*pi*t+m0);
y11 = Ay.*cos(2.*w1.*pi*t+m1);
y12 = Ay.*cos(2.*w1.*pi*t+m2);
y13 = Ay.*cos(2.*w1.*pi*t+m3);
y14 = Ay.*cos(2.*w1.*pi*t+m4);
y15 = Ay.*cos(2.*w1.*pi*t+m5);
y16 = Ay.*cos(2.*w1.*pi*t+m6);
y17 = Ay.*cos(2.*w1.*pi*t+m7);
y18 = Ay.*cos(2.*w1.*pi*t+m8);
y21 = Ay.*cos(2.*w2.*pi*t+m1);
y22 = Ay.*cos(2.*w2.*pi*t+m2);
y23 = Ay.*cos(2.*w2.*pi*t+m3);
y24 = Ay.*cos(2.*w2.*pi*t+m4);
y25 = Ay.*cos(2.*w2.*pi*t+m5);
y26 = Ay.*cos(2.*w2.*pi*t+m6);
y27 = Ay.*cos(2.*w2.*pi*t+m7);
y28 = Ay.*cos(2.*w2.*pi*t+m8);
y31 = Ay.*cos(2.*w3.*pi*t+m1);
y32 = Ay.*cos(2.*w3.*pi*t+m2);
y33 = Ay.*cos(2.*w3.*pi*t+m3);
y34 = Ay.*cos(2.*w3.*pi*t+m4);
y35 = Ay.*cos(2.*w3.*pi*t+m5);
y36 = Ay.*cos(2.*w3.*pi*t+m6);
y37 = Ay.*cos(2.*w3.*pi*t+m7);
y38 = Ay.*cos(2.*w3.*pi*t+m8);
y41 = Ay.*cos(2.*w4.*pi*t+m1);
y42 = Ay.*cos(2.*w4.*pi*t+m2);
y43 = Ay.*cos(2.*w4.*pi*t+m3);
y44 = Ay.*cos(2.*w4.*pi*t+m4);
y45 = Ay.*cos(2.*w4.*pi*t+m5);
y46 = Ay.*cos(2.*w4.*pi*t+m6);
y47 = Ay.*cos(2.*w4.*pi*t+m7);
y48 = Ay.*cos(2.*w4.*pi*t+m8);
y51 = Ay.*cos(2.*w5.*pi*t+m1);
y52 = Ay.*cos(2.*w5.*pi*t+m2);
y53 = Ay.*cos(2.*w5.*pi*t+m3);
y54 = Ay.*cos(2.*w5.*pi*t+m4);
y55 = Ay.*cos(2.*w5.*pi*t+m5);
y56 = Ay.*cos(2.*w5.*pi*t+m6);
y57 = Ay.*cos(2.*w5.*pi*t+m7);
y58 = Ay.*cos(2.*w5.*pi*t+m8);
subplot(5,8,1);plot(x0,y11);Axis([-4 4 -4 4]);text(-16,0,'ω =
1');text(-2,7,'0');text(-16,7,'相位差');
subplot(5,8,2);plot(x0,y12);Axis([-4 4 -4 4]);text(-2,7,'π/4'); subplot(5,8,3);plot(x0,y13);Axis([-4 4 -4 4]);text(-2,7,'π/2'); subplot(5,8,4);plot(x0,y14);Axis([-4 4 -4 4]);text(-2,7,'3π/4'); subplot(5,8,5);plot(x0,y15);Axis([-4 4 -4 4]);text(-2,7,'π');
subplot(5,8,6);plot(x0,y16);Axis([-4 4 -4 4]);text(-2,7,'5π/4'); subplot(5,8,7);plot(x0,y17);Axis([-4 4 -4 4]);text(-2,7,'3π/2'); subplot(5,8,8);plot(x0,y18);Axis([-4 4 -4 4]);text(-2,7,'7π/4'); subplot(5,8,9);plot(x0,y21);Axis([-4 4 -4 4]);text(-16,0,'ω = 1/2'); subplot(5,8,10);plot(x0,y22);Axis([-4 4 -4 4]);
subplot(5,8,11);plot(x0,y23);Axis([-4 4 -4 4]);
subplot(5,8,12);plot(x0,y24);Axis([-4 4 -4 4]);
subplot(5,8,13);plot(x0,y25);Axis([-4 4 -4 4]);
subplot(5,8,14);plot(x0,y26);Axis([-4 4 -4 4]);
subplot(5,8,15);plot(x0,y27);Axis([-4 4 -4 4]);
subplot(5,8,16);plot(x0,y28);Axis([-4 4 -4 4]);
subplot(5,8,17);plot(x0,y31);Axis([-4 4 -4 4]);text(-16,0,'ω = 2/3'); subplot(5,8,18);plot(x0,y32);Axis([-4 4 -4 4]);
subplot(5,8,19);plot(x0,y33);Axis([-4 4 -4 4]);
subplot(5,8,20);plot(x0,y34);Axis([-4 4 -4 4]);
subplot(5,8,21);plot(x0,y35);Axis([-4 4 -4 4]);
subplot(5,8,22);plot(x0,y36);Axis([-4 4 -4 4]);
subplot(5,8,23);plot(x0,y37);Axis([-4 4 -4 4]);
subplot(5,8,24);plot(x0,y38);Axis([-4 4 -4 4]);
subplot(5,8,25);plot(x0,y41);Axis([-4 4 -4 4]);text(-16,0,'ω = 3/4'); subplot(5,8,26);plot(x0,y42);Axis([-4 4 -4 4]);
subplot(5,8,27);plot(x0,y43);Axis([-4 4 -4 4]);
subplot(5,8,28);plot(x0,y44);Axis([-4 4 -4 4]);
subplot(5,8,29);plot(x0,y45);Axis([-4 4 -4 4]);
subplot(5,8,30);plot(x0,y46);Axis([-4 4 -4 4]);
subplot(5,8,31);plot(x0,y47);Axis([-4 4 -4 4]);
subplot(5,8,32);plot(x0,y48);Axis([-4 4 -4 4]);
subplot(5,8,33);plot(x0,y51);Axis([-4 4 -4 4]);text(-16,0,'ω = 2/5'); subplot(5,8,34);plot(x0,y52);Axis([-4 4 -4 4]);
subplot(5,8,35);plot(x0,y53);Axis([-4 4 -4 4]);
subplot(5,8,36);plot(x0,y54);Axis([-4 4 -4 4]);
subplot(5,8,37);plot(x0,y55);Axis([-4 4 -4 4]);
subplot(5,8,38);plot(x0,y56);Axis([-4 4 -4 4]);
subplot(5,8,39);plot(x0,y57);Axis([-4 4 -4 4]);
subplot(5,8,40);plot(x0,y58);Axis([-4 4 -4 4]);。