source_coding1.m

该代码实现了am

e=zeros(8,161)
b=zeros(8,161)
e(1,2:161)=s
b(1,1:160)=s
h(1,1)=2*sum(e(1,:).*b(1,:))
g(1,1)=sum(e(1,:).*e(1,:))+sum(b(1,:).*b(1,:))
k(1,1)=h(1,1)/g(1,1)                                                  %k1,e0(n)=e(1,:),b0(n)=b(1,:)
    
e(2,:)=e(1,:)-k(1,1)*b(1,:)
b(2,:)=b(1,:)-k(1,1)*e(1,:)
h(1,2)=2*sum(e(2,:)+b(2,:))
g(1,2)=sum(e(2,:).*e(2,:))+sum(b(2,:).*b(2,:))
k(1,2)=h(1,2)/g(1,2)                                                  %k2,e1(n)=e(2,:),b1(n)=b(2,:)
    
e(3,:)=e(2,:)-k(1,2)*b(2,:)
b(3,:)=b(2,:)-k(1,2)*e(2,:)
h(1,3)=2*sum(e(3,:)+b(3,:))
g(1,3)=sum(e(3,:).*e(3,:))+sum(b(3,:).*b(3,:))
k(1,3)=h(1,3)/g(1,3)                                                 %k3,e2(n),b2(n)
    
e(4,:)=e(3,:)-k(1,3)*b(3,:)
b(4,:)=b(3,:)-k(1,3)*e(3,:)
h(1,4)=2*sum(e(4,:)+b(4,:))
g(1,4)=sum(e(4,:).*e(4,:))+sum(b(4,:).*b(4,:))
k(1,4)=h(1,4)/g(1,4)                                                 %k4,e3(n),b3(n)
    
e(5,:)=e(4,:)-k(1,4)*b(4,:)
b(5,:)=b(4,:)-k(1,4)*e(4,:)
h(1,5)=2*sum(e(4,:)+b(4,:))
g(1,5)=sum(e(4,:).*e(4,:))+sum(b(4,:).*b(4,:))
k(1,5)=h(1,5)/g(1,5)                                                 %k5,e4(n),b4(n)

e(6,:)=e(5,:)-k(1,5)*b(5,:)
b(6,:)=b(5,:)-k(1,5)*e(5,:)
h(1,6)=2*sum(e(5,:)+b(5,:))
g(1,6)=sum(e(5,:).*e(5,:))+sum(b(5,:).*b(5,:))
k(1,6)=h(1,6)/g(1,6)                                                 %k6,e5(n),b5(n)

e(7,:)=e(6,:)-k(1,6)*b(6,:)
b(7,:)=b(6,:)-k(1,6)*e(6,:)
h(1,7)=2*sum(e(6,:)+b(6,:))
g(1,7)=sum(e(6,:).*e(6,:))+sum(b(6,:).*b(6,:))
k(1,7)=h(1,7)/g(1,7)                                                 %k7,e6(n),b6(n)

e(8,:)=e(7,:)-k(1,7)*b(7,:)
b(8,:)=b(7,:)-k(1,7)*e(7,:)
h(1,8)=2*sum(e(7,:)+b(7,:))
g(1,8)=sum(e(7,:).*e(7,:))+sum(b(7,:).*b(7,:))
k(1,8)=h(1,8)/g(1,8)                                                 %k8,e7(n),b7(n)