تولید یک رشته بیت (0و1) تصادفی

[supern star]

میخوام تو متلب ی رشته بیت شامل 0 و 1 به صورت تصادفی تولید کنم . مثلا 99 تا . بعد 3 بیت 3 بیت یک سمبل بشه . یعنی به 33 تا سمبل تبدیل کنم و بعدش دستور پلاتشم لازم دارم ک ببینم رشته بیتا رو . ممنون میشم دوستان ک اطلاعات دارند ، کمک کنند ...

[ehsan_faal]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

سلام.

Bits=randsrc(1,99,[0,1]);
Symbols=reshape(Bits,[33,3]);
stem(Bits,'ok')
ylim([-1,2])
grid on

[supern star]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

سلام اقا احسان
ابتدا ممنون از توجهتون
بعدش اینکه این error داد . اول روی randscr که خودم randi زدم ، اوکی شد ، بعدش الان باز error روی reshape میده :
Error using reshape
To RESHAPE the number of elements must not change.

Error in rand (line 5)
Symbols=reshape(Bits,[33,3]);

حالا من خودم این خط reshape رو برداشتم و با randi زدم ، این (0,1) نشون نمیده ، فقط صفحه رو اومده تقسیم بندی کرده . من میخوام اون 33 سمبل ک هر کدوم معرف 3 بیت تصادفی هستش ، رو پلات کنم و ببینم . مثلا 000 ، 001،010 و همینجور تا اخر ( کد گری) ...

[ehsan_faal]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

سلام. برای من که اون کد جواب میده :
http://s5.picofile.com/file/8166142450/Figure1.png

شما میتونی به جای خط تولید بیت از این دو خط هم استفاده کنی:

 MainBits=[0,1];
Bits=MainBits(randi([1,2],[1,99]));

در مورد کد گری هم توی این لینک اطلاعات مفیدی هستش:
http://www.gaussianwaves.com/2012/10...nd-gray-codes/

[supern star]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

سلام اقا احسان
خیلی ممنون ، اوکی شد
حالا میخوام همین سمبل ها که همون طور ک گفتم ، معادله سه بیت هستش ، رو ب ی شکل خاصی در بیارم . نمیدونم شما از مدولاسیون 8QAM در مبحث مخابرات اطلاعاتی دارین یا نه ، اما میخوام انو شبیه سازیش کنم . شکلی ک استادمون داده ، ی شکل استاندارد نیس ( برای اینکه کارمون سخت کنه و تقلب نکنیم :biggrin
http://upload.wikimedia.org/wikipedia/commons/f/f4/Circular_8QAM.png
چینش این سمبل ها باید مثل عکس بالایی باشه . فقط تفاوت کار من اینه ک نقاط دایره کوچیکه روی محوره و دایره بزرگه بین محور ها و تو ناحیه ها قرار دارند!
هر کدوم از اون نقطه ها در واقع معادل یک سمبل مثلا 010 هستش ک باید کد گری چینش بشه و اینم ینی نقاط نزدیک فقط باس تو ی بیت باهم اختلاف داشته باشند . ینی مثلا اگه ی نقطه بود 010 ، باید نزدیک ترین نقطش مثلا باشه 011 یا 110 ! ب همین منظور شعاع دو دایره رو هم داده .
حالا میخوام ببینم این سمبل های تولید شده رو چجوری ب این چنین شکلی تبدیل کنم ک بعدش باید عمل مدولاسیون و دمدولاسیون و بقیه کاراشو انجام بدم ... :smile:

[ehsan_faal]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

سلام.راستش من اصلا توی زمینه هایی مثل مدولاسیون و از این حرفا سر رشته ندارم. شما کدتون رو بنویسین هر جا مشکلی داشتین سعی میکنم واستون رفعش کنم. توی نت هم چیزهای فراوونی راجع به این مباحث هست، ضمن اینکه مثالهای تولباکس comm خیلی عالی این مباحث رو توضیح داده. یه نمونه کد واسه کار شما میتونه این باشه:

%--------------------------------------------------------------------------
% MQAM_Simulation.m -- ETSIT-UPM LTDS 2010-2011
%
% Simulation for a selected QAM modulation.
%
% Authors:
%  Luis Antonio Úbeda Medina (lubeme23@gmail.com)
%  Héctor Veiga Ortiz (hveiga@hawk.iit.edu)
%
% Input:
% longitud - Amount of random bits to be tested. 
% M - Deepness of modulation: 1 = 4-QAM, 2 = 16-QAM, 3 = 64-QAM.
% SNR - SNR - Signal Noise Ratio for the AWGN Channel in dB.
%
% Example: MQAM_Simulation(1000,2,5)
%
%
% Copyright 2010 Héctor Veiga Ortiz and Luis Antonio Úbeda Medina
% 
%
% This file is part of MQAM Simulator.
% 
% MQAM Simulator is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
% 
% MQAM Simulator is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
% 
% You should have received a copy of the GNU General Public License
% along with MQAM Simulator. If not, see <http://www.gnu.org/licenses/>.
%--------------------------------------------------------------------------

function [BER] = MQAM_Simulation(longitud,M,SNR)
close all;
% Setting up parameters.
% Adding zeros to amount of bits to do it 2*M multiple.
resto = rem(longitud,2*M);
% Creation of random bits
bin = randint(longitud+(2*M-resto),1,2);
L=length(bin); 

% Definition of bits with enough time to be observed. Setting up variables.
k=100;
f=10;
bit1=ones(1,k);
bit0=0*bit1;
Nivel=M;
symbol=ones(1,2*Nivel*k);
mbit=[];mx=[];my=[];
bin = fliplr(bin);

% Modulation in 2M bits groups
x=0;
y=0;

% We create groups of 2M bits to modulate
for n=0:2*Nivel:L-2*Nivel
  bit=[];
  xi=0;
  yi=0;
  % We assign an amplitude both 'x' and 'y' to the group
  for m= 1:2:2*Nivel
    if bin(n+m)==0   
      xi=xi+(2^((m-1)/2));
      bit=[bit bit0];
    else
      xi=xi-(2^((m-1)/2));
      bit=[bit bit1];
    end
    if bin(n+m+1)==0
      yi=yi+(2^((m-1)/2));
      bit=[bit bit0];
    else
      yi=yi-(2^((m-1)/2));
      bit=[bit bit1];
    end
  end
  
  x=xi*symbol;
  y=yi*symbol;
  % We store the generated symbol with the calculated amplitude inside mx and
  % my variables. We update the mbit string with the last 2M modulated bits.
  mx=[mx x];
  my=[my y];
  
  mbit=[mbit bit];
end

% Modulating and showing the constellation inside a scatterplot.
v=0:2*pi/k:2*pi*L-2*pi/k;

msync = mx + my*j;
qam=real(msync).*cos(f*v)-imag(msync).*sin(f*v); 
scatterplot(msync),grid,xlabel('I'),ylabel('Q'),title('Constellation before sending');
pause;

% Addition of AWGB to simulate the channel
Vn=awgn(qam,SNR,'measured');

% Demodulation of received signal
Vnx=Vn.*cos(f*v);
Vny=Vn.*-1.*sin(f*v);

% Low pass filtering with a Butterworth filter
[b,a]=butter(2,0.04);
Hx=2.*filter(b,a,Vnx);
Hy=2.*filter(b,a,Vny);
ML=length(Hx);

% Showing the constellation received
msync2=[];
for m=k:2*Nivel*k:ML
  Haux = Hx(m) + Hy(m)*j;
  msync2 = [msync2 Haux];
end
scatterplot(msync2),grid,xlabel('I'),ylabel('Q'),title('Constellation received');
pause;


mdeb=[];
for m=k:2*Nivel*k:ML
  sym=[];
  thx=0;thy=0;
  fx=0;fy=0;
  for n=1:Nivel
    if Hy(m) > thy
      sym = [sym bit0];
 
      fy=1;
    else
      sym = [sym bit1];
      fy=-1;
    end

    if Hx(m) > thx
      sym = [sym bit0];
      fx=1;
    else
      sym = [sym bit1];
      fx=-1;
    end
    thy = thy + fy*(2^(Nivel-n));
    thx = thx + fx*(2^(Nivel-n));
  end
  
  mdeb=[mdeb fliplr(sym)];
end

% Calculation of wrong bits
err = 0;
bin2 = [];

for n=1:length(bin)
  if(bin(n) == 1) 
    bin2 = [bin2 bit1];
  elseif bin(n) == 0
    bin2 = [bin2 bit0];
  end
end

for n=1:length(bin2)
  if(bin2(n) ~= mdeb(n)) 
    err = err + 1;
  end
end

error_bits_total = round(err/k);
% Calculation of BER to return the result
ber = error_bits_total/length(bin);
BER = ber;

% Showing final results
disp(['Total wrong bits = ' num2str(error_bits_total)]);
disp(['BER = ' num2str(ber)]);

figure(1);
subplot(4,1,1),plot(mbit,'r','linewidth',2),axis([0 k*L -0.5 1.5]),grid on,legend('Input data');
subplot(4,1,2),plot(qam,'m','linewidth',1.5),axis([0 k*L min(qam)*1.1 max(qam)*1.1]),grid on,legend('QAM Modulation');
subplot(4,1,3),plot(Vn,'g','linewidth',1.5),axis([0 k*L min(Vn)*1.1 max(Vn)*1.1]),grid on,legend('QAM Modulation with AWGN');
subplot(4,1,4),plot(mdeb,'k','linewidth',1.5),axis([0 k*L -0.5 1.5]),grid on,legend('Output data');

[supern star]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

اقا این ی سری فرمولا داره مدولاسیونش ک فک میکنم با اون فرمولا خودش فک کنم شکلش دربیاد .
بعدش اینکه این توضیحاتی ک فرمودین ، دقیقا توی کجای متلب هستش ؟ چون من زیاد با متلب اشنا نیستمو و فقط شاید ی ماه باشه ک برای پروژه هام دارم باهاش کار میکنم ...

[ehsan_faal]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

راهنمای تولباکس Comm:

http://s5.picofile.com/file/8166215792/comm.pdf.html

[supern star]

پاسخ : تولید یک رشته بیت (0و1) تصادفی

دستتون درد نکنه اقا احسان
دمتون گرم
باز کاری داشتم ، مزاحمتون خواهم شد
لطف کردین ... :smile: