Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 基础科学 > 数学
数学 A discussion board for Mathematics.
回复
 
主题工具 显示模式
旧 2009-04-13, 09:22   #1
ersansi-1
初级会员
 
注册日期: 2009-04-13
年龄: 40
帖子: 2
声望力: 0
ersansi-1 正向着好的方向发展
默认 求助 如何对复杂的二元符号函数中的一个变量积分

function [intcf1,intcf2]=intcf

v = 0.4*340;
b = 0.09;
omg = 10*pi;

syms xx;
syms t;
intcf1=0;
intcf2=0;



cf1 =exp(omg*(t-xx/v))*(sqrt((xx+b)/(xx-b))-1);

cf2 =exp(omg*(t-xx/v))*(xx-sqrt(xx^2-b^2));
%要对以上两个表达式中xx项进行积分一下是我采用的方式
%直接用int得不到结果,所以我用离散后累加求和的方式
for c=b+0.0001:0.001:20

temp1= subs(cf1,'xx',c);

intcf1=intcf1+temp1*0.001;

end;

for c=b+0.0001:0.001:20

temp2= subs(cf2,'xx',c);

intcf2=intcf2+temp2*0.001;

end;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
结果成了这个样子
Warning: Use "maple restart" to restart Maple or "clear maplemex"to unload the Maple kernel.
Note that you cannot unload the kernel on Macintosh.
> In maple at 59
In intcf at 4
1/1000*exp(10*pi*(t-23251/1360000))*(1/22351*24151^(1/2)*22351(*exp(*exp(10*pi*(t-31261/1360000))*(1/30361*32161^(1/2)*^(1/2)-1)+1/1000*exp(10*pi*(t-4111/1360000))*(1/3211*5011^(1/2)*321... Output truncated. Text exceeds maximum line length for Command Window display.
??? Error using ==> sym.sym>char2sym
Error, integer too large in context

Error in ==> sym.sym at 92
S = char2sym(x);

Error in ==> sym.maple at 92
result = sym(result);

Error in ==> sym.plus at 22
X = maple(A(,'+',B();

Error in ==> intcf at 24
intcf1=intcf1+temp1*0.001;
ersansi-1 当前离线   回复时引用此帖
旧 2009-04-13, 09:23   #2
ersansi-1
初级会员
 
注册日期: 2009-04-13
年龄: 40
帖子: 2
声望力: 0
ersansi-1 正向着好的方向发展
默认 回复: 求助 如何对复杂的二元符号函数中的一个变量积分

希望大家能给与指教,能不能用其他的方法得到积分的结果?
ersansi-1 当前离线   回复时引用此帖
回复

主题工具
显示模式

发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 18:11


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.