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

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)%绘第三阶振型函数曲线
>>


附件给出了各阶的振型曲线

[ 本帖最后由 lxq 于 2006-10-7 21:09 编辑 ]

feifeifool

肯定是振型函数的方程写的不对！自己查查看。

fanghuikeer

振型函数不对,找本结构力学书好好看看,画图的程序没有错误...

tiange

我是按照《结构动力学》（盛宏玉，合肥工业大学出版社 2005）中的公式做的
请大家继续给予指点
谢谢

附 书中的公式

toes

http://forum.vibunion.com/forum/ ... =page%3D1#pid149457