8psk_oceangdm_might4ep_QPSK_8psk_8psk调制解调

function [minout] = minilength(bn)%最小距离判决法 %bn为增加了噪声后的信号 N=length(bn); minout=zeros(1,N); d=zeros(8,N); md=zeros(1,N); for i=1:N d(1,i)=(bn(1,i)-cos(1/8*pi))^2+(bn(2,i)-sin(1/8*pi))^2; %计算距离 d(2,i)=(bn(1,i)-cos(3/8*pi))^2+(bn(2,i)-sin(3/8*pi))^2; d(3,i)=(bn(1,i)-cos(5/8*pi))^2+(bn(2,i)-sin(5/8*pi))^2; d(4,i)=(bn(1,i)-cos(7/8*pi))^2+(bn(2,i)-sin(7/8*pi))^2; d(5,i)=(bn(1,i)-cos(9/8*pi))^2+(bn(2,i)-sin(9/8*pi))^2; d(6,i)=(bn(1,i)-cos(11/8*pi))^2+(bn(2,i)-sin(11/8*pi))^2; d(7,i)=(bn(1,i)-cos(13/8*pi))^2+(bn(2,i)-sin(13/8*pi))^2; d(8,i)=(bn(1,i)-cos(15/8*pi))^2+(bn(2,i)-sin(15/8*pi))^2; md(i)=min(d(:,i)); %取最小距离 if (md(i)==d(1,i)) minout(i)=0; end if (md(i)==d(2,i)) minout(i)=1; end if (md(i)==d(3,i)) minout(i)=2; end if (md(i)==d(4,i)) minout(i)=3; end if (md(i)==d(5,i)) minout(i)=4; end if (md(i)==d(6,i)) minout(i)=5; end if (md(i)==d(7,i)) minout(i)=6; end if (md(i)==d(8,i)) minout(i)=7; end end