signal_base.m

matlab的调制解调

% 这个函数产生基带信号，并画出图形，
% 有一点是一定要注意的：基带信号和调制信号的采样频率应该是一样的；
% 还有采样点数也要一样，这就保证了基带信号和调制信号的持续
% 时间是一样的；
% 输入的x是001001的代码；
function y=signal_base(x,fs,n,plot_flag)
['-------','retaining time of the base symbol',...
    '   ',num2str(n/fs),'s','---']
 if size(x,1)<=size(x,2)
    x=x';
 else
    x=x;
 end;   % 把x变换成列向量；
% 把x变换成双极性码；
   for i=1:size(x,1)
       if x(i)==0
          x(i)=-1;
       elseif x(i)==1
          x(i)=1;
       end;
   end;
   t=(0:n-1)/fs; % 这里不用；
   y=[];
   for i=1:size(x,1)
       temp=x(i)*ones(1,n);
       y=[y,temp];
   end;
   y=y';
   if plot_flag==1
      figure;
      subplot(211);
      t_t=(0:size(y,1)-1)/fs;
      plot(t_t,y); 
      xlabel('时间(s)'); ylabel('基带信号幅度');
      title('变换为双极性码后的基带信号'); grid on;
      subplot(212);
      f_y=abs(fftshift(fft(y)));
      f_f=fs*(-size(y,1)/2:(size(y,1)/2-1))/size(y,1);
      plot(f_f,f_y); 
      xlabel('频率(Hz)'); ylabel('频谱幅度');
      title('基带信号的频谱'); grid on;
   elseif plot_flag==0
   end;