The multiscaling function of Alpert multiwavelet system and Chui-Lian system multiwavelet consists of one symmetric and one antisymmetric scaling functions. In its single level decomposition, only one of 16 subband blocks may be considered as a core part because the antisymmetric scaling function may play the high pass filter. We complete the Alpert multiwavelet system by determining the high pass filter using the algorithm of paraunitary extention principle. Taking the advantage of the symmetric/antisymmetric (SA) feature of Alpert and Chui-Lian multiscaling functions, our iteration (SA iteration) of multiwavelet decomposition continues to the core part (1/16 of subbands) for the image denoising problem while the traditional iteration applies to the one fourth of 16 subbands as core parts. The numerical results demonstrate that our method of symmetric/antisymmetric iteration exhibits performance superior to the traditional ways.