Which toolbox use for Phase Analysis

>> close all

%Reading input sequence

[Header Sequence]=fastaread('D:\KERALA PROJECT\DCB-training-kerala\KERALA PROJECT\HIV -1 Sequence\hiv-1 B (K03455).fasta');

Num=length(Sequence);

S=Sequence;

% nt 2 int

gene=nt2int(Sequence);

% Complex Mapping

u(find(S=='A'))=1+j;

u(find(S=='G'))=-1+j;

u(find(S=='C'))=-1-j;

u(find(S=='T'))=1-j;

u(find(S=='W'))=1;

u(find(S=='Y'))=-j;

u(find(S=='S'))=-1;

u(find(S=='R'))=j;

u(find(S=='K'=='M'=='N'))=0;

u(find(S=='D'))=1/3+j/3;

u(find(S=='H'))=1/3-j/3;

u(find(S=='B'))=-1/3-j/3;

u(find(S=='V'))=-1/3+j/3;

% Phase Analysis

z=phase(u);

P=angle(u);

plot(P);

figure(1);

hold all;

% Cumulated Phase

B=cumsum(P);

plot(B);

figure(1);

hold all;

% Unwrapped Phase Analysis

Q =unwrap(P);

plot(Q);

figure(1);

hold all;

title('Phase Analysis');

xlabel('Bases');

ylabel('Radian');

Unrecognized function or variable 'phase'.

[Anmol Dhiman]

Can you share sample file and what is exact purpose of phase function as you are not using it in your code anywhere else.

[AMITABH CHATURVEDI]

Sir, Actual I will convert the Nucleotide Sequence (Biological Sequence) to digital sequence. I have used complex mapping for it. Please help me on urgent basis.

[Anmol Dhiman]

HI ,

There is no phase function defined in MATLAB, that you are trying to use. What exactly do you want to do with the phase function .