butfly.m

附录 光盘说明 本书附赠的光盘包括各章节实例的设计工程与源码

%-------------------------------------------------------------------------%
%--     固定缩小比例的蝶形结算法matlab验证程序。                          --%
%--         作者：Xudong                                                --%
%--         时间：2004年5月6日                                          --%
%------------------------------------------------------------------------%
function [xk1,xk2]=butfly(x1,x2,wn);
%--输入：x1,x2,wn都是8位有符号数（-128～+127），复数形式。
%--输出：xk1,xk2也是8位有符号数（-128～+127），复数形式。
DataWidth=8;
PointDot=32;

xr1=real(x1);
xi1=imag(x1);
xr2=real(x2);
xi2=imag(x2);
wr=real(wn);
wi=imag(wn);

mul1out=xr2*wr;%--四个乘法器的输出。
mul2out=xi2*wi;
mul3out=xr2*wi;
mul4out=xi2*wr;

midsubout=mul1out-mul2out;%--中间的加法器和减法器。
midaddout=mul3out+mul4out;

xkr1=floor(xr1/2)+floor(midsubout/2^DataWidth);%--最后输出前的加法和减法器。
xkr2=floor(xr1/2)-floor(midsubout/2^DataWidth);%--为防止数据遗出输出数据固定缩小2倍。
xki1=floor(xi1/2)+floor(midaddout/2^DataWidth);%--输入的wn事先已经扩大了2^(DataWidth-1)倍。
xki2=floor(xi1/2)-floor(midaddout/2^DataWidth);

xk1=xkr1+i*xki1;%--组合成复数后输出。
xk2=xkr2+i*xki2;