帮忙看一下带参数积分的程序，有什么错误

maigicku

我想画一条D-Q的曲线，表达式如下图所示：

其中，


Q的取值为Q=[8 9 13 18 26 36 52 72]，u的值为7，最后我画出来的曲线感觉不太对，请高手帮忙看一下，程序有什么错误..谢谢~~
程序如下：

clear;clc
u=7;
%%%%%%%%%%%%%%%%%%%i>0时%%%%%%%%%%%%%%%%%%%%%%%%%
fch1=@(Q) arrayfun(@(i) quadl(@(x) abs(x-i.*Q).*1/pi*u./(u^2+x.^2),(i-1/2).*Q,(i+1/2).*Q),1:100);
Q=[8 9 13 18 26 36 52 72];
ych1=zeros(size(Q));
for j=1:length(Q)
    ych1(j)=sum(fch1(Q(j)));
end
%%%%%%%%%%%%%%%%%%i=0时%%%%%%%%%%%%%%%%%%%%%%%%%%
fch2=@(Q) quadl(@(x) abs(x).*1/pi*u./(u^2+x.^2),0,0.5.*Q);
Q=[8 9 13 18 26 36 52 72];
ych2=zeros(size(Q));
for j=1:length(Q)
    ych2(j)=sum(fch2(Q(j)));
end
dch=2*ych1+2*ych2;
plot(Q,dch,'-sg')

画出的曲线如下：


[ 本帖最后由 maigicku 于 2009-10-26 17:27 编辑 ]