声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1071|回复: 2

[综合讨论] 积分遇到问题了,谁能帮帮我?

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

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

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

x
积分遇到问题了,谁能帮帮我?
程序如下:

clear ;clc;
syms x y z
f=(1-x.^2).^(1/2);
dfdx=diff(f);
A=int(5e13*cosh(dfdx*z/1e6),z,0,1e-7);
B=int(5e13*sinh(dfdx*z/1e6),z,0,1e-7);
C=(0.04*A-B*(A^2-B^2+0.04^2)^0.5)/(A^2-B^2);
D=dfdx*(1e-7)+1e6*sinh(C)
E=int(D,x,(-1e-4),1e-4)
F=int(E,y,(-0.008),0.008)
x1=-1:0.05:1
G=subs(F,x,x1)

运行结果:
D =

-1/10000000/(1-x^2)^(1/2)*x+1000000*sinh((1000000000000000000*(1-x^2)^(1/2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))/x+25000000000000000000*(1-x^2)^(1/2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)/x*(625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))^2/x^2-625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)^2/x^2+1/625)^(1/2))/(625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))^2/x^2-625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)^2/x^2))


Warning: Explicit integral could not be found.
> In sym.int at 58
In yy at 9

E =

int(-1/10000000/(1-x^2)^(1/2)*x+1000000*sinh((1000000000000000000*(1-x^2)^(1/2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))/x+25000000000000000000*(1-x^2)^(1/2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)/x*(625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))^2/x^2-625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)^2/x^2+1/625)^(1/2))/(625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))^2/x^2-625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)^2/x^2)),x = -1/10000 .. 1/10000)



F =

2/125*int(-1/10000000/(1-x^2)^(1/2)*x+1000000*sinh((1000000000000000000*(1-x^2)^(1/2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))/x+25000000000000000000*(1-x^2)^(1/2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)/x*(625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))^2/x^2-625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)^2/x^2+1/625)^(1/2))/(625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)-exp(-1/10000000000000/(1-x^2)^(1/2)*x))^2/x^2-625000000000000000000000000000000000000*(1-x^2)*(exp(1/10000000000000/(1-x^2)^(1/2)*x)+exp(-1/10000000000000/(1-x^2)^(1/2)*x)-2)^2/x^2)),x = -1/10000 .. 1/10000)



x1 =

Columns 1 through 11

-1.0000 -0.9500 -0.9000 -0.8500 -0.8000 -0.7500 -0.7000 -0.6500 -0.6000 -0.5500 -0.5000

Columns 12 through 22

-0.4500 -0.4000 -0.3500 -0.3000 -0.2500 -0.2000 -0.1500 -0.1000 -0.0500 0 0.0500

Columns 23 through 33

0.1000 0.1500 0.2000 0.2500 0.3000 0.3500 0.4000 0.4500 0.5000 0.5500 0.6000

Columns 34 through 41

0.6500 0.7000 0.7500 0.8000 0.8500 0.9000 0.9500 1.0000

??? Error using ==> maple
Error, (in int) wrong number (or type) of arguments

Error in ==> sym.subs at 165
NEWf = double(sym(maple('map','F',char(sym([NEWexpr{:}])))));

Error in ==> yy at 12
G=subs(F,x,x1)


[ 本帖最后由 twomao 于 2007-8-28 16:59 编辑 ]
回复
分享到:

使用道具 举报

发表于 2007-8-29 16:28 | 显示全部楼层
改用数值积分试试.
另: 从你贴的代码看,好象并没有y,“F=int(E,y,(-0.008),0.008)”应该有问题。
     实际上,E已经是一个常数了。
 楼主| 发表于 2007-8-29 17:17 | 显示全部楼层

回复 #2 xjzuo 的帖子

真的太谢谢你啦.
我想用quad做,但是不知道怎样在inline里,将dfdx=diff(f);的结果中的乘法换为数组乘法.
我再考虑下.
再次感谢楼上朋友的指点.
:handshake
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-9 04:59 , Processed in 0.059158 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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