I am not sure, but it seems that sresid function of MatLab 5.3 is not correct. It is supposed to implement a residualization method of order reduction, but when I want to do so (not a balanced residualization), it do not change A or B or C matrices. Have you encountered this problem?