L=0.60;%弦长L
F=9.8;%外力F
P=0.00055;%弦密度
f=1.00;%频率
b=(f)^-1*(F/P)^0.5;%波长
A=6;w=2pi/b;k=2pi*f;
x=0:0.01:L;
t=0:500;
for f=50:0.01:58
for t=0:500
b=(f)^-1*(F/P)^0.5;%波长
A=6;w=2pi/b;k=2pif;
x=0:0.01:L;
E_1=Acos(wt+kx);
E_2=Acos(wt-kx-4L*pi/b-pi);
E=E_1+E_2;
end
N=max(E);
end
N=max(E);
f;