2楼: Originally posted by csgt0 at 2012-07-06 09:49:44
不太懂卷积，这样？
=deconv(Eout,Ein);
Eout'=conv(Ein,q)+r;
a=Eout'-Eout;
any(a);

3楼: Originally posted by shuangqingx at 2012-07-06 10:25:34
这是要在matlab里做吗？我不会用matlab，现在我用的是origin做的，结果总是呈现震荡。。。...

4楼: Originally posted by csgt0 at 2012-07-06 10:35:00
是matlab，但是算出来的q=0,r=Ein......

for i=1:length(u)   
    k=u(1:i-1)'*flipud(v(1:i-1));
    p(i)=(v(i)-k)/u(1);
end

6楼: Originally posted by csgt0 at 2012-07-06 12:15:47
我用定义算的,不知道对不对
u=Ein;
v=Eout;

for i=1:length(u)   
    k=u(1:i-1)'*flipud(v(1:i-1));
    p(i)=(v(i)-k)/u(1);
end

ec/90/449046_1341548139_596.jpg

8楼: Originally posted by csgt0 at 2012-07-06 13:36:48
好像是这样，刚有点错
clc
clear
load a.mat
p(1)=v(1)/u(1);
for i=2:length(u)   
    k=p(1:i-1)*flipud(v(2:i));
    p(i)=(v(i)-k)/u(1);
end
plot(1:length(u),p)

我按这个的计算的http://wenku ...