声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1322|回复: 1

[结构振动] 这个程序错在哪里了,我怎么弄不明白

[复制链接]
发表于 2007-5-13 15:58 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
我要对一个悬臂梁振动前四阶振型进行求导,然后用matlab做出其图像,但是总是出现这个错误
的提示:??? Error using ==> plot
Vectors must be the same lengths.
请大家给予指点。


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.01: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);
%对第一阶振型函数求导
m2=diff(y2);
%对第二阶振型函数求导
m3=diff(y3);
%对第三阶振型函数求导
m4=diff(y4);
%对第四阶振型函数求导

plot(x,m1,x,m2,x,m3,x,m4)
%绘各阶振型函数曲线
回复
分享到:

使用道具 举报

发表于 2007-5-13 17:18 | 显示全部楼层
Vectors must be the same lengths.维数不一致

DIFF(X), for a vector X, is [X(2)-X(1)  X(3)-X(2) ... X(n)-X(n-1)].
y1是1*101的,m1应当是100的,怎么画图?

评分

1

查看全部评分

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-2 19:32 , Processed in 0.067155 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表