Borrar filtros
Borrar filtros

how to have character replacement

2 visualizaciones (últimos 30 días)
Sharen H
Sharen H el 27 de Feb. de 2013
every digit is converted to binary and are grouped with two bits each and i have to replace it with characters.
example A=11000100
i have to replace as
if
11 then 'A'
00 'C'
01 'G'
10 'T'
final result should be a=ACGT
I have written a code it takes long time to get executed
bi = dec2bin(secret,8) - '0';
out = reshape(mat2cell(bi,ones(size(bi,1),1),2*ones(size(bi,2)/2,1))',1,[]);
[r,c]=size(out);
for a = 1:c
if( out{a} == [0 0])
out{a} = 'C';
else
if( out{a} == [0 1])
out{a} = 'T';
else
if( out{a} == [1 0])
out{a} = 'A';
else
if( out{a} == [1 1])
out{a} = 'G';
end
end
end
end
end
  1 comentario
Jan
Jan el 27 de Feb. de 2013
Editada: Jan el 27 de Feb. de 2013
Is A the double 11000100? If so, why do you run DEC2BIN? Or is the "A" in the text equal to "bi" in the code?

Iniciar sesión para comentar.

Respuesta aceptada

Azzi Abdelmalek
Azzi Abdelmalek el 27 de Feb. de 2013
A='11000100'
s='ACGT'
f=['11';'00';'01';'10']
b=reshape(A,2,[])'
for k=1:numel(s)
idx=find(ismember(b,f(k,:),'rows'))
c(idx)=s(k)
end
out=c(:)'
  3 comentarios
Azzi Abdelmalek
Azzi Abdelmalek el 27 de Feb. de 2013
clear
s='GCTA'
f=['11';'00';'01';'10'];
secret=[70 77];
bi = dec2bin(secret,8);
for p=1:size(bi,1)
b=reshape(bi(p,:),2,[])';
for k=1:numel(s)
idx=find(ismember(b,f(k,:),'rows'));
c(idx)=s(k);
end
out{p}=c(:)';
end
out
Sharen H
Sharen H el 27 de Feb. de 2013
Editada: Sharen H el 27 de Feb. de 2013
Sir if it is just reverse how should i do ie the input i have is
s=GCTA
i should replace it as
s = 11000110
Please help

Iniciar sesión para comentar.

Más respuestas (1)

Honglei Chen
Honglei Chen el 27 de Feb. de 2013
Editada: Honglei Chen el 27 de Feb. de 2013
Sounds like a sequencing problem. You can do it like this
x = '11000100';
idx = bin2dec(reshape(x,2,[]).')+1;
matchset = 'CGTA';
matchset(idx)
  2 comentarios
Azzi Abdelmalek
Azzi Abdelmalek el 27 de Feb. de 2013
Editada: Azzi Abdelmalek el 27 de Feb. de 2013
Typo,
matchset(idx)
instead of
matcheset(idx)
Honglei Chen
Honglei Chen el 27 de Feb. de 2013
corrected, thanks.

Iniciar sesión para comentar.

Categorías

Más información sobre Numeric Types en Help Center y File Exchange.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by