|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
我的程序段是这样的:
clear; g=0.01; n=24; cei=0.01; omga=1; bata=4.5; afa=0.1;
v=0.8; b=0.005*v^2; lame=[1.875 4.694 7.855];
syms x
for t=1:3
fai(t)=cosh(lame(t)*x)-cos(lame(t)*x)-(cosh(lame(t))+cos(lame(t)))*(sinh(lame(t)*x)-sin(lame(t)*x))/(sinh(lame(t))+sin(lame(t)));
end
for i=1:3, for j=1:3
mm(i,j)=int(fai(j)*fai(i), x, 0, 1);
cc(i,j)=cei*int(fai(j)*diff(fai(i), 4), x, 0, 1)/omga;
kk(i,j)=(int(fai(j)*diff(fai(i), 4), x, 0, 1)-bata^2*afa*int((1-x)*fai(j)*diff(fai(i), 2), x, 0, 1)-1/2*bata^2*int((1-x^2)*fai(j)*diff(fai(i), 2), x, 0, 1))/omga^2;
m(i,j)=double(mm(i,j)); c(i,j)=double(cc(i,j)); k(i,j)=double(kk(i,j));
end; end
运行之后一直提示这样的错误:
??? Error using ==> reshape
To RESHAPE the number of elements must not change.
Error in ==> sym.maple at 94
result = reshape(result,size(varargin{3}));
Error in ==> sym.int at 51
r = reshape(maple('map','int',f(:),[x.s '=(' a.s ')..(' b.s ')']),size(f));
怎么回事呢,我改了好久也不行。。太着急了,请大家帮帮忙~~
[ 本帖最后由 ChaChing 于 2009-10-27 20:02 编辑 ] |
|