声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 929|回复: 1

[综合讨论] 请教该积分为何不能求解

[复制链接]
发表于 2009-6-10 16:37 | 显示全部楼层 |阅读模式

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

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

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 编辑 ]

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2009-6-10 16:39 | 显示全部楼层
结果说明没有解析解,可以采用数值积分的办法求解。

评分

1

查看全部评分

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

本版积分规则

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

GMT+8, 2024-11-14 18:04 , Processed in 0.069523 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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