ray_channel.m

来自「这是个扩频系统中的PN码仿真程序」· M 代码 · 共 29 行

% function ray_channel=();
%clarke model
clear;
fm=1.2;%fdmax
N=2048;
fs=10000;
data_f=2*fm/(N-1);
% for i=1:floor(N/2)
%     [gauss_noise(i),gauss_noise1(i)]=gngauss;
% end
gauss_noise=randn(1,floor(N/2)).*exp(rand(1,floor(N/2))*2*pi);
gauss_noise1=[conj(fliplr(gauss_noise)),gauss_noise];
for i=1:N
    S_ez(i)=sqrt(1.5/(pi*fm*sqrt(1-(i/fm)^2)));
end
gauss_noise1=gauss_noise1.*S_ez;
gauss_noise1=ifft(gauss_noise1);
gauss_noise1=abs(gauss_noise1);
gauss_noise1=gauss_noise1.*gauss_noise1;
gauss_noise=randn(1,floor(N/2)).*exp(rand(1,floor(N/2))*2*pi);
gauss_noise2=[conj(fliplr(gauss_noise)),gauss_noise];
gauss_noise2=gauss_noise2.*S_ez;
gauss_noise2=ifft(gauss_noise2);
gauss_noise2=abs(gauss_noise2);
gauss_noise2=gauss_noise2.*gauss_noise2;
ray_channel=sqrt(gauss_noise1+gauss_noise2);
plot(ray_channel);