Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2010-11-05, 16:03   #1
zcc8607
初级会员
 
注册日期: 2010-10-29
帖子: 2
声望力: 0
zcc8607 正向着好的方向发展
默认 请教符号积分时遇到的问题

想解一个二重积分,看网上有人说可以先做符号积分int,再做数值积分quadv。可是我int那里就出问题了,不知道怎么解决。

clear;
syms t w
L=0.08;
a=0.1;b=0.08;c=-0.15;uz=0.2;vz=0.2;wz=0.1;
f1=[(-3/20+w)/((2/25-t)*(2/25-conj(t))+(-3/20+w)*(-3/20+conj(w)))^(3/2)+(3/20-w)/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20+w)*(-3/20+conj(w)))^(3/2)+(-3/20-w)/((2/25-t)*(2/25-conj(t))+(-3/20-w)*(-3/20-conj(w)))^(3/2)+(3/20+w)/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20-w)*(-3/20-conj(w)))^(3/2),
0,
1/5/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20+w)*(-3/20+conj(w)))^(3/2)+1/5/(1/25+(2/25-t)*(2/25-conj(t))+(-3/20-w)*(-3/20-conj(w)))^(3/2)];
f2=int(f1,t,-0.1,0.1); %%%%%出错的地方
f2=inline(f2/L);
result=quadv(f2,w,-wz,-wz-L);

f1有点复杂,是个行向量,带有两个变量t和w;int是想对t求积分,quadv想对w求积分。

??? Attempted to access e(1); index out of bounds because numel(e)=0.
Error in ==> sym.sym>char2sym at 518
sk = x(s(k):e(k));
Error in ==> sym.sym at 95
S = char2sym(x);
Error in ==> sym.maple at 92
result = sym(result,'keepijalias');
Error in ==> sym.int at 51
r = reshape(maple('map','int',f(,[x.s '=(' a.s ')..(' b.s ')']),size(f));
Error in ==> tmpy at 10
f2=int(f1,t,-0.1,0.1);

错误信息如上,不知道怎么解决。请各位大虾帮帮忙,谢谢!
zcc8607 当前离线   回复时引用此帖
回复

主题工具
显示模式

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

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



所有时间均为北京时间。现在的时间是 23:43


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