24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1176  |  回复: 0

qianbianzhu

新虫 (小有名气)

[求助] 如何计算这个二重积分(有图有程序)

如何计算这个二重积分(有图有程序)
tz.JPG
下面是我对这个公式中的H1写的一段程序,但是怎么计算也算不出来

a1=1/1000;
a2=11/1000;
b1=5/1000;
b2=30/1000;
h1=0.5/1000;
h2=0.7/1000;
z=0;
x=(1:30)/1000;
%y=(1:30)/1000;
y=1/1000;
syms a b h
R11=sqrt(((b2-b1)*a/(a2-a1)+y).^2+(z-h).^2+(x-a).^2);
R12=sqrt(((b2-b1)*a/(a2-a1)-y).^2+(z-h).^2+(x-a).^2);
HB11=(((b2-b1)*a/(a2-a1)-y)./R11+((b2-b1)*a/(a2-a1)+y)./R12).*(z-h)./((z-h).^2+(x-a).^2);
F=vpa(int(HB11,h,h1,h2));
H=vpa(int(F,a,a1,a2));


然后尝试了下面用泰勒展开
k=taylor(HB11,h);
kk=taylor(k,a);
r=int(kk,a,a1,a2);
rr=int(r,h,h1,h2);

计算出来的结果感觉不对
请大家帮我看看应该怎么改才好呢?谢谢了
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 qianbianzhu 的主题更新
信息提示
请填处理意见