command of alignment for general sequences
2 ビュー (過去 30 日間)
古いコメントを表示
hi, is there a command of sequence alignment for any sequence generally ? I did not mean swalign or nwalign, these just for biological sequences. if not are there anyone wrote a code of sequence alignment for non biological sequences in matlab?
thanks
4 件のコメント
回答 (1 件)
Hin Kwan Wong
2011 年 11 月 28 日
What sequence? If you mean any two numerical sequences which are correlated. Their alignment can be done by using cross correlation and finding where the peak is in terms of the delay, then shifting one of the sequence by the delay.
xcorr = ifft(fft(signalA).*conj(fft(signalB));
delay = find((xcorr==max(xcorr))==1);
AlignedA = circshift(signalA,1-delay);
0 件のコメント
参考
カテゴリ
Help Center および File 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!