L=1; %悬臂梁的度
C1=10; %假设 C1
a1=1.875/L;
a2=4.694/L;
a3=7.855/L;
a4=7*pi/2*L;
x=0:0.1:1;
y1=C1*( cos(a1*x)-cosh(a1*x) + (cos(a1*L)+cosh(a1*L))*(sinh(a1*x)-sin(a1*x))/(sin(a1*L)+sinh(a1*L)) );
%第一阶振型函数
y2=C1*( cos(a2*x)-cosh(a2*x) + (cos(a2*L)+cosh(a2*L))*(sinh(a2*x)-sin(a2*x))/(sin(a2*L)+sinh(a2*L)) );
%第三阶振型函数
y3=C1*( cos(a3*x)-cosh(a3*x) + (cos(a3*L)+cosh(a3*L))*(sinh(a3*x)-sin(a3*x))/(sin(a3*L)+sinh(a3*L)) );
%第三阶振型函数
y4=C1*( cos(a4*x)-cosh(a4*x) + (cos(a4*L)+cosh(a4*L))*(sinh(a4*x)-sin(a4*x))/(sin(a4*L)+sinh(a4*L)) );
%第四阶振型函数
m1=diff(y1,2);
m2=diff(y2,2);
m3=diff(y3,2);
m4=diff(y4,2);
x1=0.01:0.01:1;
subplot(4,1,1)
plot(x,m1)%绘第一阶振型函数曲线
subplot(4,1,2)
plot(x,m2)%绘第二阶振型函数曲线
subplot(4,1,3)
plot(x,m3)%绘第三阶振型函数曲线
subplot(4,1,4)
plot(x,m4)%绘第四阶振型函数曲线 |