vec.m

本程序是QAM的调制程序

%向量化算子
%就是把M行N列的矩阵转化为2M*N维的列向量
function P=vec(A);
% real=zeros(1,2*M*N)
% imag=zeros(1,2*M*N)
% real(:,1:end)=(reshape(A)+reshape(conj(A)))/2%提取矩阵A的实部，并把它转化为一维数组,有M*N个元素
% imag(:,1:end)=(reshape(A)-reshape(conj(A)))/2%提取矩阵A的虚部，并把它转化为一维数组，有M*N个元素
% C=zeros(1,2*M*N)
% C(:,1:2:end) = reshape(,2,N/2) 
% C(:,2:2:end) = (1/sqrt(2))*(kron(ones(1,N/2),[-1;1]).*flipud(reshape(conj(s),2,N/2)))
% C()=reshape(real,1,1:2:2*M*N)
% C()=reshape(imag,1,2:2:2*M*N)
B=zeros(1,2);
B=size(A);
M=B(1);
N=B(2);
%real=zeros(1,M*N)
% imag=zeros(1,M*N)
real=(reshape(A,1,M*N)+reshape(conj(A),1,M*N))./2;
imag=((reshape(A,1,M*N)-reshape(conj(A),1,M*N))./2)*(-j);
%real()=(reshape(A,1,M*N)+reshape(conj(A),1,M*N)).\2 %提取矩阵A的实部，并把它转化为一维数组,有M*N个元素
%imag()=(reshape(A,1,M*N)-reshape(conj(A),1,M*N)).\2 %提取矩阵A的虚部，并把它转化为一维数组，有M*N个元素
C=zeros(1,2*M*N);
C(:,1:2:end)=reshape(real,1,M*N);
C(:,2:2:end)=reshape(imag,1,M*N);
% C()=reshape(real,1,1:2:2*M*N)
% C()=reshape(imag,1,2:2:2*M*N)
C=C';
P=C;