zzc7.m

基于FDMA背景下的斜投影算子区分多用户信息

%基于FDMA背景下的斜投影算子区分多用户信息，在此定义用户1和用户2

close all
clc
clear all

t=0.001:0.001:1;   %抽样率的设置为1000

%------信号系统及干扰系统，将用户1和用户2信息载波分别定义为H空间和S空间------%

H1=cos(2*pi*100*t);    %用户1系统的实部，载波的频率为100HZ
H2=sin(2*pi*100*t)*i;  %用户1系统的虚部
H=[H1.' H2.'];            
S1=cos(2*pi*150*t);    %用户2系统的实部，载波的频率为150HZ
S2=sin(2*pi*150*t)*i;  %用户2系统的虚部
S=[S1.' S2.'];

%---------斜投影算子--------%

Ps=S*inv((S'*S))*S';             %S空间的投影算子
Ps1=eye(1000,1000)-Ps;           %S空间的正交投影算子
Ehs=H*inv(H'*Ps1*H)*H'*Ps1;     %沿着S空间到H空间的投影算子，即斜投影算子
Ehs1=inv(H'*Ps1*H)*H'*Ps1;

%--------信号及干扰---------%

m=cos(2*pi*5*t)+1;       %用户1信息，在此定义为余弦，频率为5HZ
b=sin(2*pi*5*t);       %用户2信息，在此定义为正弦，频率为5HZ


%--------实现斜投影算子对有用信号的提取，-------%
for n=1:1000
    
    %----以下小段语句实现的是多用户信号参量与载波参量对应时间的乘积---%
    f1=m(n); 
    f2=b(n);
    g1=H.*f1;
    g2=S.*f2;
 
               
    x1=awgn(g1,10);   %给信号加高斯白噪声           
    y1=x1+g2;         %综合信号 

    y2(n)=y1(n,1)+y1(n,2);
   
    g=Ehs1*y1;        %解调信息
    
    j(n)=g(1);
   
end


%-------仿真结果的观察图表--------%
subplot(3,2,1);plot(t,m);title('有用信号');
subplot(3,2,2);plot(angle(m));title('有用信号相位');
subplot(3,2,3);plot(t,real(y2));title('综合信号实部');
subplot(3,2,4);plot(angle(y2));title('综合信号相位');
subplot(3,2,5);plot(t,real(j));title('解调信息实部');
subplot(3,2,6);plot(angle(j));title('解调信息相位');