声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1153|回复: 2

[综合讨论] 四重符号积分程序为何提示Explicit integral could not be found.

[复制链接]
发表于 2007-12-31 09:44 | 显示全部楼层 |阅读模式

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

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

x
源程序!!!!
syms a1 a2 z1 z2
n=2;                  %磁极数
br1=1.1;br2=1.1;     %分别为内、外驱动器永磁体的剩余磁感应强度
l1=0.117;l2=0.052;    %分别为外驱动器和内驱动器的长度
g=0.04825;            %内、外驱动器间的径向间隙
h=0.005;              %外驱动器磁环的径向厚度
R3=0.002;R2=R3+g;R1=R2+h;
l=(l1-l2)/2;                 %内驱动器的起始轴端外驱动器起始轴端的距离
k=1.075;a=0*pi/18;x=-0.05;
f1=int(((R2*R3/(R2^2+R3^2+(z2-z1)^2-2*R2*R3*cos(a2-a1))^1.5))*(z1-z2),a1,a,a+pi);
ff1=inline(f1);
fff1=vectorize(ff1);
sss11=triplequad(fff1,0,pi,0,l1,l+x,l+l2+x);
对变量a1,a2,z1,z2积分。原函数:function=((R2*R3/(R2^2+R3^2+(z2-z1)^2-2*R2*R3*cos(a2-a1))^1.5))*(z1-z2)
提示很多错误
于是改为
syms a1 a2 z1 z2
n=2;                  %磁极数
br1=1.1;br2=1.1;     %分别为内、外驱动器永磁体的剩余磁感应强度
l1=0.117;l2=0.052;    %分别为外驱动器和内驱动器的长度
g=0.04825;            %内、外驱动器间的径向间隙
h=0.005;              %外驱动器磁环的径向厚度
R3=0.002;R2=R3+g;R1=R2+h;
l=(l1-l2)/2;                 %内驱动器的起始轴端外驱动器起始轴端的距离
k=1.075;a=0*pi/18;x=-0.05;
f1=int(((R2*R3/(R2^2+R3^2+(z2-z1)^2-2*R2*R3*cos(a2-a1))^1.5))*(z1-z2),a1,a,a+pi);
temp=int(f1,a2,0,pi)
Warning: Explicit integral could not be found.
请高手帮忙看看。

[ 本帖最后由 eight 于 2007-12-31 13:36 编辑 ]

本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2007-12-31 10:50 | 显示全部楼层
我还是觉得数值积分好

顺便把楼主错误提示贴出来:
第1个程序
??? Error using ==> inline.subsref at 14
Not enough inputs to inline function.

Error in ==> triplequad>innerintegral at 65
fcl = intfcn(xmin, y(1), z(1), varargin{:});

Error in ==> dblquad>innerintegral at 80
fcl = intfcn(xmin, y(1), varargin{:}); %evaluate only to get the class below

Error in ==> quad at 77
y = f(x, varargin{:});

Error in ==> dblquad at 60
Q = quadf(@innerintegral, ymin, ymax, tol, trace, intfcn, ...

Error in ==> triplequad at 53
Q = dblquad(@innerintegral, ymin, ymax, zmin, zmax, tol, quadf, intfcn, ...

评分

1

查看全部评分

 楼主| 发表于 2007-12-31 14:50 | 显示全部楼层
什么原因阿
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-14 03:30 , Processed in 0.072467 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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