CONVERT DNA SEQUENCE INTO COMPLEX NUMBER

1 次查看(过去 30 天)
Now I have to convert the DNA sequence into complex number.
what i did was :
s=randseq(20) % to have random sequence of DNA consisting of ATCG s=
ATCGCGCGATATATCGCGCG
next, I am stuck. the thing is, i need to convert each letter to complex. for example, A=-1+1i T=1-1i C=-1=1i G=1+1i
i need to convert and map the sequence into complex to use it after that in daubechies wavelet.
Thank you

回答(1 个)

David Young
David Young 2015-2-21
values('ATCG') = [-1+1i 1-1i -1-1i 1+1i]; % mapping from characters to complex
s = 'ATCGCGCGATATATCGCGCG'; % any specific sequence
scomplex = values(s); % the sequence represented as complex numbers

类别

Help CenterFile Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by