声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1738|回复: 2

[编程技巧] 实函数积分得到虚数,怎么就不能积分了呢

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

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

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

x
其中,a=(-625/144*cos(25/12*x)-625/144*cosh(25/12*x)+4132590367507980625/1297036692682702848*sin(25/12*x)+4132590367507980625/1297036692682702848*sinh(25/12*x))*(-5508409/202500*cos(2347/450*x)-5508409/202500*cosh(2347/450*x)+2526579984188513246993/91197892454252544000*sin(2347/450*x)+2526579984188513246993/91197892454252544000*sinh(2347/450*x)),对它积分,程序出错,说是reshape出错,一项一项地积分,得到的却是虚变函数,在0-0.9积分,得到虚数,而a只是我的程序中的一个矩阵的元素,其它元素都能积分,它就不能,请高手指点,怎么样才能不用一项一项积分,也能计算出来,而我试过
sinx+cosx+sinhx+coshx的积分是个实数,是不是我计算错了而得到了虚数呢
回复
分享到:

使用道具 举报

发表于 2007-5-8 17:57 | 显示全部楼层
请将你的处理贴出来看看.
%%%%%%%%%%%%
%%% 我没有看到你说的问题:
定义inline函数a,积分: Q = quadl(a,0,2)
%%% 结果为实数:
Q =
-4.1434e+004

评分

1

查看全部评分

 楼主| 发表于 2007-5-9 19:04 | 显示全部楼层
xiexie你的回复,我对k2=int(diff(ph2,2)).'*(diff(ph2,2),0,l)求积分,ph2如下所示
ph2 =[cos(25/12*x)-cosh(25/12*x)-64527948014827/87901297004728*sin(25/12*x)+64527948014827/87901297004728*sinh(25/12*x), cos(2347/450*x)-cosh(2347/450*x)-7688618371577599/7549204801302143*sin(2347/450*x)+7688618371577599/7549204801302143*sinh(2347/450*x)]
总是积不出来,我才先算出diff(ph2,2)).'*(diff(ph2,2)),然后对各个元素积分,我先前的发的那个帖子里面的a就是diff(ph2,2)).'*(diff(ph2,2))的第(2,2)个元素

望指点,太感谢你了

[ 本帖最后由 eight 于 2007-5-9 22:59 编辑 ]
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-25 00:42 , Processed in 0.052329 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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