jww100 发表于 2008-10-29 05:02

请问这个卷积程序错在哪?

是用直接积分法计算ht和ft的卷积。可是不知道程序错在哪里,还望各位高人指点。

N=1024;
T=0.1;
delt=T/N;
t1=0:delt:(6*N-1)*delt;
ht=10*exp(-20*t1).*sin(20*pi*t1);
t2=0:delt:(3*N-1)*delt;
ft=5*sin(20*pi*t2);

syms tao;
for i=1:9*N-1
    t=(i-1)*delt;
    hx=10*exp(-20*tao)*sin(20*pi*tao);
    fx=5*sin(20*pi*(t-tao));
    hxfx=hx*fx;
    fconv(i)=int(hxfx,tao,-inf,inf);
end
figure;
plot(double(fconv));

[ 本帖最后由 jww100 于 2008-10-29 05:04 编辑 ]

施少敏 发表于 2008-10-29 10:53

看一下conv2的内部函数就知道了,我也不大懂
页: [1]
查看完整版本: 请问这个卷积程序错在哪?