Commit dd7370d1 by Russell Luke

### Matlab test code

parent 14a1c51a
 ... @@ -5,30 +5,30 @@ M=max(size(x)); ... @@ -5,30 +5,30 @@ M=max(size(x)); fmag=zeros(M, M/2); fmag=zeros(M, M/2); Fmag=zeros(M,M); Fmag=zeros(M,M); for(i=1:M/2) for(i=1:M/2) j=1; j=1; datapow=x(i,j)+datapow; datapow=x(i,j)+datapow; fmag(i,j)=sqrt(x(i,j))/M; fmag(i,j)=sqrt(x(i,j))/M; Fmag(i,j) = fmag(i,j); Fmag(i,j) = fmag(i,j); for(j=2:M/2) for(j=2:M/2) datapow=2*x(i,j)+datapow; datapow=2*x(i,j)+datapow; fmag(i,j)=sqrt(x(i,j))/M; fmag(i,j)=sqrt(x(i,j))/M; Fmag(i,j) = fmag(i,j); Fmag(i,j) = fmag(i,j); if(i>1) if(i>1) Fmag(M-(j-2),M-(i-2)) = fmag(i,j); Fmag(M-(j-2),M-(i-2)) = fmag(i,j); end end end end end end for(i=M/2+1:M) for(i=M/2+1:M) j=1; j=1; datapow=x(i,j)+datapow; datapow=x(i,j)+datapow; fmag(i,j)=sqrt(x(i,j))/M; fmag(i,j)=sqrt(x(i,j))/M; Fmag(i,j) = fmag(i,j); Fmag(i,j) = fmag(i,j); for(j=2:M/2) for(j=2:M/2) datapow=2*x(i,j)+datapow; datapow=2*x(i,j)+datapow; fmag(i,j)=sqrt(x(i,j))/M; fmag(i,j)=sqrt(x(i,j))/M; Fmag(i,j) = fmag(i,j); Fmag(i,j) = fmag(i,j); Fmag(j,i) = fmag(i,j); Fmag(j,i) = fmag(i,j); end end end end Fmag(1,M/2:M)=fmag(M/2:M,1); Fmag(1,M/2:M)=fmag(M/2:M,1); ... ...
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!