channel.m

来自「Abstract：Noise frequency modulation(FM)j」· M 代码 · 共 31 行

function channel(a1,a2,d1,d2,fc)

load  td_sets td_sets;
load bc_chunks bc_chunks;

% 多径信道
copy1 = zeros(size(bc_chunks));
for i=1+d1:length(bc_chunks)
    copy1(i) = a1*bc_chunks(i-d1);
end
copy2 = zeros(size(bc_chunks));
for i=1+d2:length(bc_chunks)
    copy2(i) = a2*bc_chunks(i-d2);
end
bc_chunks = bc_chunks+copy1+copy2;

% 产生随机噪声
Gaussian_y=wgn(fc,1,1)';
bc_chunks = bc_chunks+Gaussian_y;
save Gaussian_y Gaussian_y;
subplot(211);
plot(bc_chunks);

p=fft(bc_chunks);
pp=p.*conj(p);
cp=10*log10(pp);
subplot(212)
plot(cp);