crt.m

用MATLAB实现数字信号处理的程序

M1=5*7*11;
M2=3*7*11;
M3=3*5*11;
M4=3*5*7;
M=3*5*7*11;

i1=0;i2=0;i3=0;i4=0;
while(k1~=1)
    k1=mod(M1*i1,3);
    i1=i1+1;
end
    C1=k1;
while(k2~=1)
    k2=mod(M1*i2,5);
    i2=i2+1;
end
    C2=i2;
while(k3~=1)
    k3=mod(M1*i3,7);
    i3=i3+1;
end
    C3=i3;
while(k4~=1)
    k4=mod(M4*i4,11);
    i4=i4+1;
end
    C4=i4;
C={C4,C3,C2,C1}

%x=M1*C1*u1+M2*C2*u2+M3*C3*u3+M4*C4*u4;
x=M1*C1*1+M2*C2*4+M3*C3*2+M4*C4*9;
X=mod(x,M)

u1=mod(394,3);
u2=mod(394,5);
u3=mod(394,7);
u4=mod(394,11);
U={u4,u3,u2,u1}

t1=mod(214,3);
t2=mod(214,5);
t3=mod(214,7);
t4=mod(214,11);
t={t4,t3,t2,t1}