Генерация Carrierwave, добавление шума и фильтра с помощью фильтра нижних частот в Octave - PullRequest
Новая
       99

Генерация Carrierwave, добавление шума и фильтра с помощью фильтра нижних частот в Octave

0 голосов
/ 16 апреля 2020

У меня есть следующее упражнение из университета:

Битовая последовательность [1 1 0 1 1 1 0 0] со скоростью передачи данных DR = 50 бит / с. Эта битовая последовательность модулирует амплитуду несущего сигнала с частотой F C = 1000 Гц. Во время передачи к модулированному сигналу добавляется белый шум с отклонением сигма = 0,5. Сигнал демодулируется в приемнике и фильтруется фильтром нижних частот с частотой среза FG = F C. Нанесите на график демодулированный и отфильтрованный сигнал приема и сравните его с сигналом передачи. Встроенные функции не допускаются, кроме обсужденных в лекции!

До сих пор я мог генерировать шум, используя следующий код: 

t = linspace(0,0.1,1001);
sigma = 0.5;
avg = 0;
noise = sigma * randn(size(t))+avg;
plot(t,noise)

Я делал что-то подобное в предыдущем упражнении, но с некоторыми отличиями. У меня была прямоугольная волна, к которой я добавил шум, затем я отфильтровал результат с помощью фильтра нижних частот (не встроенной функции) и сравнил прямоугольную волну по умолчанию с той, к которой я добавил шум, и отфильтровал его (результат). Код для этого приведен ниже. 

pkg load signal;
t = linspace(0,0.1,1001);
y = square(2*pi*50*t)*0.5 + 0.5;
mw = 0;
sigma = 0.5;
rauschen = sigma*randn(size(t))+mw;
yr =  y + rauschen;
Y = fft(y);
Yr = fft(yr);
f = linspace(0,1/t(2),length(t));
plot(f,abs(Y),f,abs(Yr))
fg = 400;
Yru = Yr;
Yru(find(f>fg || f<1/t(2)-fg)) = 0;
find(f==fg)
ans =  41
Yru(41:end-41)=0;
plot(f,abs(Yru))
yru = ifft(Yru);
plot(t,real(yru));

На рисунке ниже показан результат предыдущего упражнения.

enter image description here

My main проблема в том, что я не могу генерировать необходимую несущую волну / сигнал ..

1 Ответ

0 голосов
/ 17 апреля 2020

Благодаря этой статье Я думаю, что я понял это. Вот код и результат показан на рисунке ниже. 

pkg load signal;
bits = [ 1 1 0 1 1 1 0 0 ];
bitrate =  50;
n = 1000;
T = length(bits)/bitrate;
N = n*length(bits);
dt = T/N;
t = 0:dt:T;
x = zeros(1,length(t));
lastbit = 1;
for i=1:length(bits)
  if bits(i)==1
    x((i-1)*n+1:i*n) = -lastbit;
    lastbit = -lastbit;
  else x((i-1)*n+1:i*n) = lastbit;
  end
end
x = (x + 1)*0.5;
plot(t, x, 'Linewidth', 3);
figure;
sigma = 0.5;
rauschen = sigma * randn(size(t));
plot(t,rauschen)
figure;
xr = x + rauschen;
plot(t,xr)
figure;
Xr = fft(xr);
f = linspace(0,1/t(2),length(t));
fc = n;
Xru = Xr;
Xru(find(f==fc)) = 0;
plot(f,abs(Xru));
figure;
ans = find(f==fc);
ans
Xru(ans:end-ans)=0;
plot(f,abs(Xru));
figure;
xru = ifft(Xru);
plot(t,real(xru));
hold
plot(t,x)
#axis([0 t(end) -0.5 1.5])

enter image description here

...