b.m

来自「无线移动信道的matlab源码 用于对无线通信信道的编程实现」· M 代码 · 共 35 行

clear all;
% A=zeros(1,800);
% A(1,1:50)=1;A(351:400)=1;
% plot(abs(fftshift((fft(x)))))

% t=0:0.002:4;
% x=zeros(1,length(t));
% for m=1:52
%     x=x+sin(2*pi*(m*312500+5000000000)*t);
% end  
% plot(x);
% X=fft(x,64);
% P=X.*conj(X)/64;
% f=(64*312500+5000000000)*(1:64)/64;
% plot(f,abs(X));


t=0:0.001:0.6;
x=zeros(1,length(t));
for m=1:52
    x=x+sin(2*pi*((m-1)*50+50)*t);
end  
 f_max=400;
  m_s=1;
  t_0=0;
T_s=0.2E-6;
     
    [C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,q_l,T]=F_S_K_p(9,'tu',f_max);
    [y_t,T,t_0]=F_S_K(x,f_max,m_s,T,t_0,q_l,C1,F1,TH1,C2,F2,TH2,F01,F02,RHO,F_RHO,T_s);
    y=y_t+1.5*randn(1,length(t));
Y=fft(y,512);
P=Y.*conj(Y)/512;
f=1000*(0:255)/512;
plot(f,P(1:256))