[求助]关于带参数积分以及积分上下限也带参数的问题

maigicku

看了论坛上几篇关于带参数积分的帖子后，偶试着去做下面的积分，发现做不出，总是出现错误：

程序如下：
clear;clc; syms x i Q j; l=zeros(1,201); D(1:81)=0;
for Q=linspace(0,80,81)
    for i=-10000:100:10000
        for j=1:81
            q=Q(j);
            for k=1:201
                l=i(k);
               fun=inline(subs('(abs(x-i.*Q)).^2.*((0.0734/2).*exp((-0.0734).*abs(x)))', {Q,i},{q,l}));
               D(Q+1)=D(Q+1)+sum(quadl(fun,(i-1/2).^Q,(i+1/2).^Q));
            end
        end
    end
end
Q=linspace(0,80,81);
plot(Q,D)

错误信息：
??? Error using ==> sym.sym>char2sym at 414
Not a valid symbolic expression.
Error in ==> sym.sym at 95
   S = char2sym(x);
Error in ==> subs at 60
r = subs(sym(f),varargin{:});
Error in ==> DLaplace at 21
        fun=inline(subs('(abs(x-i.*Q)).^2.*((0.0734/2).*exp((-0.0734).*abs(x)))',...

[ 本帖最后由 maigicku 于 2009-3-27 09:19 编辑 ]

maigicku

哦，偶这个自己已经解决了。。就是太拘泥于积分函数了。。其实不用积分函数也可以做出来。。
先建立一个函数d(q)：
function f = d(q)
delta_x = 1;
max_i =100;
min_i = -100;
f = 0;
for i = min_i:1:max_i
    temp = 0;
    for x = (min_i-0.5).*q:delta_x:(max_i+0.5).*q
        temp = temp + (x-i.*q)*(x-i.*q)*(0.0734/2).*exp((-0.0734).*abs(x));
    end
    f = f + temp;
end

然后运行：
clear;clc
q=linspace(0,80,81);
f=zeros(size(q));
for j=1:length(q)
    f(j)=d(q(j));
end
plot(q,f)

就可以了。。。

[ 本帖最后由 maigicku 于 2009-3-31 09:03 编辑 ]