initial

Author: fischeti

FMReciver.m

@@ -0,0 +1,44 @@
+% Request user input from the command-line for application parameters
+userInput = helperFMUserInput;
+
+% Calculate FM system parameters based on the user input
+[fmRxParams,sigSrc] = helperFMConfig(userInput);
+
+% Create FM broadcast receiver object and configure based on user input
+fmBroadcastDemod = comm.FMBroadcastDemodulator(...
+    'SampleRate', fmRxParams.FrontEndSampleRate, ...
+    'FrequencyDeviation', fmRxParams.FrequencyDeviation, ...
+    'FilterTimeConstant', fmRxParams.FilterTimeConstant, ...
+    'AudioSampleRate', fmRxParams.AudioSampleRate, ...
+    'Stereo', true);
+
+% Create audio player
+player = audioDeviceWriter('SampleRate',fmRxParams.AudioSampleRate);
+
+% Initialize radio time
+radioTime = 0;
+
+% Main loop
+while radioTime < userInput.Duration
+  % Receive baseband samples (Signal Source)
+  if fmRxParams.isSourceRadio
+    [rcv,~,lost,late] = sigSrc();
+  else
+    rcv = sigSrc();
+    lost = 0;
+    late = 1;
+  end
+
+  % Demodulate FM broadcast signals and play the decoded audio
+  audioSig = fmBroadcastDemod(rcv);
+  player(audioSig);
+
+  % Update radio time. If there were lost samples, add those too.
+  radioTime = radioTime + fmRxParams.FrontEndFrameTime + ...
+    double(lost)/fmRxParams.FrontEndSampleRate;
+end
+
+% Release the audio and the signal source
+release(sigSrc)
+release(fmBroadcastDemod)
+release(player)

Reciver.m

@@ -0,0 +1,15 @@
+%SDR Reciver
+
+%recive
+inputnoise = recive(); %TODO
+
+%demodulate
+input = findbegin(inputnoise); %TODO
+recive = demodulate(input); %TODO
+bits = derepet(recive); %TODO
+data = prepostremove(bits); %TODO
+
+%save data
+%TODO
+
+

Transmitter.asv

@@ -0,0 +1,24 @@
+% SDR Transmitter
+
+%Recive data input
+
+data = [1,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,1,0,1,1,0,1,0,1,1,0,0];
+% s = fopen('datei.txt','r');
+% bitsequence = fread(s, 'ubit1')'
+% fclose(s);
+
+% Modulate data
+%prefix =  [1,1,1,1,1,1,1,1];
+%postfix = [1,1,1,1,1,1,1,1];
+
+bits = data;
+send = fm4modulate(bits);
+receive = fm4demodulate(send);
+
+subplot(2,1,1)
+plot(receive,'ro')
+subplot()
+%Send data
+%send(send) %TODO
+% plot(send);
+% plot(fmdemodulate(send),'ro');

Transmitter.m

@@ -0,0 +1,25 @@
+% SDR Transmitter
+
+%Recive data input
+
+data = [1,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,1,0,1,1,0,1,0,1,1,0,0];
+% s = fopen('datei.txt','r');
+% bitsequence = fread(s, 'ubit1')'
+% fclose(s);
+
+% Modulate data
+%prefix =  [1,1,1,1,1,1,1,1];
+%postfix = [1,1,1,1,1,1,1,1];
+
+bits = data;
+send = fm4modulate(bits);
+receive = fm4demodulate(send);
+
+subplot(2,1,1)
+plot(receive,'ro')
+subplot(2,1,2)
+plot(data, 'ro')
+%Send data
+%send(send) %TODO
+% plot(send);
+% plot(fmdemodulate(send),'ro');

findbits.asv

@@ -0,0 +1,15 @@
+function x = findbits (bits,tofind)
+
+x = zeros(1,length(bits)/length(tofind));
+
+for i=1:length(tofind):length(bits)-1
+    
+    x = [x,tofind == bits(1,i:i+length(tofind)-1)]
+        
+        
+end
+
+
+
+
+end

findbits.m

@@ -0,0 +1,10 @@
+function x = findbits (bits,tofind)
+
+x = [];
+
+for i=1:length(tofind):length(bits)-1
+    
+    x = [x,isequal(tofind,bits(1,i:i+length(tofind)-1))]
+         
+end
+end

fm4demodulate.asv

@@ -0,0 +1,38 @@
+function z = fm4demodulate (bits)
+
+
+tau00 = 4;
+tau01 = 8;
+tau10 = 16;
+tau11 = 20;
+
+tauS = 80;
+
+t = 0:length(bits)-1;
+
+for i = 1:length(bits)/tauS
+   z00(i) = cos(2*pi/tau00*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+   z01(i) = cos(2*pi/tau01*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+   z10(i) = cos(2*pi/tau10*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+   z11(i) = cos(2*pi/tau11*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+end
+% 
+
+r00 = repdecode(z00,2);
+r01 = repdecode(z01,2);
+r10 = repdecode(z10,2);
+r11 = repdecode(z11,2);
+% subplot(4,1,1)
+% plot(z00,'ro')
+% subplot(4,1,2)
+% plot(z01,'ro')
+% subplot(4,1,3)
+% plot(z10,'ro')
+% subplot(4,1,4)
+% plot(z11,'ro')
+
+
+
+
+
+end

fm4demodulate.m

@@ -0,0 +1,50 @@
+function z = fm4demodulate (bits)
+
+
+tau00 = 4;
+tau01 = 8;
+tau10 = 16;
+tau11 = 20;
+
+tauS = 80;
+
+t = 0:length(bits)-1;
+
+for i = 1:length(bits)/tauS
+   z00(i) = cos(2*pi/tau00*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+   z01(i) = cos(2*pi/tau01*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+   z10(i) = cos(2*pi/tau10*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+   z11(i) = cos(2*pi/tau11*t((i-1)*tauS+1:i*tauS))*bits((i-1)*tauS+1:i*tauS)';
+end
+% 
+
+r00 = repdecode(z00,2);
+r01 = repdecode(z01,2);
+r10 = repdecode(z10,2);
+r11 = repdecode(z11,2);
+
+z00 = round(r00/(tauS/2));
+z01 = round(r01/(tauS/2));
+z10 = round(r10/(tauS/2));
+z11 = round(r11/(tauS/2));
+
+for i=1:length(z00)
+    z(2*i-1:2*i) = z00(i).*[0,0]+z01(i).*[0,1]+z10(i).*[1,0]+z11(i).*[1,1];
+end
+
+subplot(4,1,1)
+plot(z00,'ro')
+subplot(4,1,2)
+plot(z01,'ro')
+subplot(4,1,3)
+plot(z10,'ro')
+subplot(4,1,4)
+plot(z11,'ro')
+
+
+
+
+
+
+
+end

fm4modulate.asv

@@ -0,0 +1,32 @@
+function y= fm4modulate(data)
+
+tau00 = 4;
+tau01 = 8;
+tau10 = 16;
+tau11 = 20;
+
+tauS = 80;
+
+% repencode
+
+t = linspace(0,tauS*length(data),length(data));
+
+y00 = findbits(data,[0,0]);
+y01 = findbits(data,[0,1]);
+y10 = findbits(data,[1,0]);
+y11 = findbits(data,[1,1]);
+
+z00 = repencode(y00,tauS).*cos(2*pi/tau00*t);
+z01 = repencode(y01,tauS).*cos(2*pi/tau01*t);
+z10 = repencode(y10,tauS).*cos(2*pi/tau10*t);
+z11 = repencode(y11,tauS).*cos(2*pi/tau11*t);
+
+z = z00 + z01 + z10 + z11;
+
+plot(t,z)
+
+
+
+
+
+end

fm4modulate.m

@@ -0,0 +1,38 @@
+function z = fm4modulate(data)
+
+tau00 = 4;
+tau01 = 8;
+tau10 = 16;
+tau11 = 20;
+
+tauS = 80;
+
+% repencode
+
+t = linspace(0,tauS*length(data)-1,tauS*length(data));
+
+y00 = findbits(data,[0,0]);
+y01 = findbits(data,[0,1]);
+y10 = findbits(data,[1,0]);
+y11 = findbits(data,[1,1]);
+
+z00 = repencode(y00,2*tauS).*cos(2*pi/tau00*t);
+z01 = repencode(y01,2*tauS).*cos(2*pi/tau01*t);
+z10 = repencode(y10,2*tauS).*cos(2*pi/tau10*t);
+z11 = repencode(y11,2*tauS).*cos(2*pi/tau11*t);
+
+% subplot(5,1,1)
+% plot(z00)
+% subplot(5,1,2)
+% plot(z01)
+% subplot(5,1,3)
+% plot(z10)
+% subplot(5,1,4)
+% plot(z11)
+% subplot(5,1,5)
+% plot(data,'ro')
+z = z00 + z01 + z10 + z11;
+% 
+% plot(t,z)
+
+end

fmdemodulate.m

@@ -0,0 +1,46 @@
+function bhat = demodulate(r)
+%UNTITLED2 Summary of this function goes here
+%   Detailed explanation goes here
+
+tau0 = 20;
+tau1 = 8;
+tauS = 80;
+
+t = 0:length(r)-1;
+
+% y0 = cos(2*pi/tau0*t).*r;
+% y1 = cos(2*pi/tau1*t).*r;
+
+for i = 1:length(r)/tauS
+   z0(i) = cos(2*pi/tau0*t((i-1)*tauS+1:i*tauS))*r((i-1)*tauS+1:i*tauS)'
+   z1(i) = cos(2*pi/tau1*t((i-1)*tauS+1:i*tauS))*r((i-1)*tauS+1:i*tauS)'
+end
+
+
+bhat = (z1-z0+tauS/2)/tauS;
+% 
+% subplot(7,1,1);
+% plot(t,r);
+% axis([0 length(r) -1 1]);
+% subplot(7,1,2);
+% plot(t,y0);
+% axis([0 length(r) -1 1]);
+% subplot(7,1,3);
+% plot([1:length(r)/tauS],z0,'o');
+% axis([1 length(r)/tauS 0 40]);
+% subplot(7,1,4);
+% plot(t,y1);
+% axis([0 length(r) -1 1]);
+% subplot(7,1,5);
+% plot([1:length(r)/tauS],z1,'o');
+% axis([1 length(r)/tauS 0 40]);
+% subplot(7,1,6);
+% plot([1:length(r)/tauS],z1-z0,'o');
+% axis([1 length(r)/tauS -40 40]);
+% subplot(7,1,7);
+% plot([1:length(r)/tauS],bhat,'o');
+% axis([1 length(r)/tauS 0 1]);
+
+end
+
+

fmmodulate.m

@@ -0,0 +1,27 @@
+function y = modulate(b)
+%UNTITLED Summary of this function goes here
+%   Detailed explanation goes here
+tau0 = 20;
+tau1 = 8;
+tauS = 80;
+
+b_rep = repencode(b, tauS);
+t = linspace(0,length(b_rep)-1,length(b_rep));
+
+y1 = b_rep.*cos(2*pi/tau1*t);
+y2 = (1-b_rep).*cos(2*pi/tau0*t);
+% subplot(5,1,1)
+% x1 = 1:length(b);
+% plot(x1,b,'o');
+% subplot(5,1,2)
+% plot(t, b_rep)
+% subplot(5,1,3)
+% plot(t, y1)
+% subplot(5,1,4)
+% plot(t,y2)
+% subplot(5,1,5)
+% plot(t,y1+y2);
+
+y = y1+y2;
+end
+

repdecode.m

@@ -0,0 +1,16 @@
+function outputsequence = repdecode(bit_sequence, n)
+% function output = repdecode(seq, n)
+%
+% Zweck: Wandelt eine repetierte Bitsequenz wieder in die normale
+%        Bitsequenz zurück in dem es den Mittelwert aus den n-mal repetitiven Bits
+%        nimmt.
+%
+% Parameters:   bit_sequence    ... Bitsequenz
+%               n               ... Repetitionszahl jedes Bits
+%
+% Returns:      outputsequence  ... ausgegebene Bitsequenz
+%
+outputsequence = round(sum(reshape(bit_sequence, n, length(bit_sequence)/n))/n);
+
+end
+

repencode.m

@@ -0,0 +1,20 @@
+function outputsequence = repencode(bit_sequence, n)
+% function outputsequence = repencode(seq, n)
+%
+% Zweck: Repetiercode der jedes Bit in einer Bitsequenz n-mal repetiert
+%        z.B. 101 -> 111000111
+%
+% Parameters:   bit_sequence    ... Bitsequenz
+%               n               ... Repetitionszahl jedes Bits
+%
+% Returns:      outputsequence  ... ausgegebene Bitsequenz
+%
+
+% Macht aus einer aus einem [1xm]-Vektor eine [nxm]-Matrix
+SEQ = repmat(bit_sequence, n,1);
+
+% Wandelt eine [nxm]-Matrix in ein [1x(n*m)]-Vektor um
+outputsequence = reshape(SEQ, 1, length(bit_sequence)*n);
+
+end
+