我按照书本所作出的振型曲线，怎么就是不正确，请大家继续给予指导

tiange

大家看看我的这个求悬臂梁 振型函数 的程序 错在那里？

若一个悬臂梁，长度L=30;
那么，按照《结构动力学》里所给的振型函数公式，编程序求前三阶振型图像，但是怎么就和书中的不一样呢？？

程序（设C1=10）：
>> L=30;  %悬臂梁的度
>> C1=10;    %假设 C1
>> a1=1.875/L;  
>> a2=4.694/L;
>> a3=7.855/L;
>> x=0:1:500;
>> y1=C1*(cos(a1*x)-cosh(a1*x)+(cos(a1*L)+cos(a1*L))*(sin(a1*L)-sin(a1*x))/(sin(a1*L)+sinh(a1*L))); %第一阶振型函数
>> y2=C1*(cos(a2*x)-cosh(a2*x)+(cos(a2*L)+cos(a2*L))*(sin(a2*L)-sin(a2*x))/(sin(a2*L)+sinh(a2*L)));%第三阶振型函数
>> y3=C1*(cos(a3*x)-cosh(a3*x)+(cos(a3*L)+cos(a3*L))*(sin(a3*L)-sin(a3*x))/(sin(a3*L)+sinh(a3*L)));%第三阶振型函数
>> subplot(3,1,1)
>> plot(x,y1)%绘第一阶振型函数曲线
>> subplot(3,1,2)
>> plot(x,y2)%绘第三阶振型函数曲线
>> subplot(3,1,3)
>> plot(x,y3)%绘第三阶振型函数曲线
>>

附：1.前三阶的振型曲线
        2.书中的公式


许多朋友说我的公式可能有问题，但是我是根据书中（《结构动力学》，盛宏玉，合肥工业大学出版社，2005）中所给的公式做出来的。怎么就是不正确呢
请大家继续给予指点
谢谢

无水1324

问题可能出在解超越方程上！
对于非线性的超越方程准确求解是很困难的，需要多种方法对比看解的误差！

Hermes

我原来参照王光远老师的<建筑结构的振动>一书中的公式做过类似的振型图，给你我的程序，你参考一下：
sl=1.875;    %第一振型
%sl=4.7;    %第二振型
%sl=7.87;  %第三振型
l=10;
s=sl/l;
x=0:0.01:10;
y=((sin(s*x)-sinh(s*x))-(sin(sl)+sinh(sl))/(cos(sl)+cosh(sl))*(cos(s*x)-cosh(s*x)))
plot(x,y)

由于当时是借的图书馆的书，现在我这里没有这本书，也没有看过你说的那本书，不过从程序上看，可能是公式上的问题，具体的你可以参照书上的思路自己推一下。