|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
syms fs ls r m0 q k x L s0
[email=f=@(x)sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2]f=@(x)sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2[/email]);
Q=quad(f,0,L)
结果为:??? Undefined function or method 'isfinite' for input arguments of type 'sym'.
Error in ==> quad at 75
if ~isfinite(y(1))
或改为:
syms fs ls r m0 q k x L s0
[email=f=@sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2]f=sqrt(1+(((m0+q./(k*k))*(k*sinh(k*x)-k*tanh(0.5*k*L)*cosh(k*x))-q*x+0.5*q*L)./s0).^2[/email]);
Q=int(f,x,0,L)
结果为:Warning: Explicit integral could not be found.
> In sym.int at 58
Q =
int((1+((m0+q/k^2)*(k*sinh(k*x)-k*tanh(1/2*k*L)*cosh(k*x))-q*x+1/2*q*L)^2/s0^2)^(1/2),x = 0 .. L)
请教高人,这是为什么啊?怎么不能求解呢?
谢谢!
[ 本帖最后由 simonwoo 于 2009-6-10 16:39 编辑 ] |
|