test_sig.m

来自「高斯加性白噪声下采用MPSK调制的系统BER性能分析」· M 代码 · 共 50 行

function [I1,Q1]=test_sig(Fl,Fs,x,delta_fc,dela,B)
% [I1,Q1]=test_sig(Fl,Fs,t1,delta_fc,dela);
% [I1,Q1]=test_sig(3000,30000,1,x,1000,pi/4);
%% s(t)=a(t)cosWct+b(t)sinWct
%% I1(t)=a(t)cos(wt-dela)+b(t)sin(wt-dela);
%% Q1(t)=-a(t)sin(wt-dela)+b(t)cos(wt-dela)
%% w=Wc-W0;
%Fl=3000;Fs=30000;
%delta_fc=100;
%dela=pi/4;
%t1=1;

error(nargchk(3,6,nargin));
if nargin==3
   delta_fc=0;
   dela=0;
   B=0;
end
if nargin==4
   dela=0;
   B=0;
end
if nargin==5
   B=0;
end

if Fs/Fl~=floor(Fs/Fl)
   error('Fs/Fl must be integer!')
   return
end
x1=real(x);
x2=imag(x);
z=ones(1,Fs/Fl);
X1=x1'*z;
X2=x2'*z;
[r,c]=size(X1);
X1=X1.';
X2=X2.';
X1=reshape(X1,1,r*c);
X2=reshape(X2,1,r*c);
t=0:1/Fs:(r*c-1)/Fs;
h=0.5+0.5*cos(2*pi*Fl*(t-0.5/Fl+B/Fl));
X1=X1.*h;
X2=X2.*h;
R=(X1+j*X2).*exp(j*2*pi*delta_fc*t+j*dela);
I1=real(R);
Q1=imag(R);
return