Dear Yanhua,
I have not looked at your attachment, but I would not use the cce.m script. It is far better to apply MBC3 and then azimuth-average the MBC-transformed matrices. Please see the MBC3 documentation for information on what it is doing. As discussed in my post dated May 07, 2010 in the following forum topic: Fast Linearized Models - #4 by Jason.Jonkman, MBC does not eliminate the periodicity, it simply transforms the rotating states of a model into the nonrotating frame and further azimuth-averaging of the lightly-periodic MBC-transformed matrices (i.e. averaging the 36 matrices into 1 single matrix) has been shown to be quite adequate for controls design.
Best regards,