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

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

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

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

x

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

账号		自动登录	找回密码
密码			我要加入

[编程技巧] happy救我！！！积分问题！！！

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