how to have character replacement

조회 수: 1 (최근 30일)
Sharen H
Sharen H 2013년 2월 27일
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
Jan
Jan 2013년 2월 27일
편집: Jan 2013년 2월 27일
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?

댓글을 달려면 로그인하십시오.

채택된 답변

Azzi Abdelmalek
Azzi Abdelmalek 2013년 2월 27일
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
Azzi Abdelmalek
Azzi Abdelmalek 2013년 2월 27일
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 2013년 2월 27일
편집: Sharen H 2013년 2월 27일
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

댓글을 달려면 로그인하십시오.

추가 답변 (1개)

Honglei Chen
Honglei Chen 2013년 2월 27일
편집: Honglei Chen 2013년 2월 27일
Sounds like a sequencing problem. You can do it like this
x = '11000100';
idx = bin2dec(reshape(x,2,[]).')+1;
matchset = 'CGTA';
matchset(idx)
  댓글 수: 2
Azzi Abdelmalek
Azzi Abdelmalek 2013년 2월 27일
편집: Azzi Abdelmalek 2013년 2월 27일
Typo,
matchset(idx)
instead of
matcheset(idx)
Honglei Chen
Honglei Chen 2013년 2월 27일
corrected, thanks.

댓글을 달려면 로그인하십시오.

카테고리

Help CenterFile Exchange에서 Graphics Object Programming에 대해 자세히 알아보기

Community Treasure Hunt

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

Start Hunting!

Translated by