heropmj007 发表于 2006-4-13 12:15

happy救我！！！积分问题！！！

function Q=inte(x, y)
wp=93*2*pi/60;wd=95*2*pi/60;
r1=250;r2=48;
for rp=24:16:144
o1=-atan(r2/rp);o2=atan(r2/rp);
Q = dblquad(inline('sqrt((wp*rp)^2+wd^2*(rp^2+y^2)-2*rp^2*wp*wd+(2*wp*rp*wd-2*rp*wd^2)*y*cos(abs(x)))'), o1, o2, rp-48, rp+48);
plot(rp,Q);
hold on
end

发现：
Q = dblquad(inline('sqrt((wp*rp)^2+wd^2*(rp^2+y^2)-2*rp^2*wp*wd+(2*wp*rp*wd-2*rp*wd^2)*y*cos(abs(x)))'), o1, o2, rp-48, rp+48)

这个二次积分有问题，我不知道怎么弄。

我的目的就是对每个rp对应一个积分值，然后以rp为横坐标值，积分值为纵坐标值画出图来，二次积分以x和y为积分变量。我都问了n个人的了，都没弄出来。帮我看看哈，谢谢拉！！

lana2000 发表于 2006-4-13 12:18

好多人指名要找happy哦，他是不是超厉害的阿？？

tomding 发表于 2006-4-13 15:27

好象很厉害的

happy 发表于 2006-4-14 21:58

http://forum.vibunion.com/thread-10318-1-1.html

已回复，换function

tomding 发表于 2006-4-16 10:02

qiang
页: [1]
查看完整版本: happy救我！！！积分问题！！！