Home > K25104 > PDEs > amp_fac.m

amp_fac

PURPOSE ^

first check user inputs

SYNOPSIS ^

function [ H ] = amp_fac( alpha, beta,r,s )

DESCRIPTION ^

 first check user inputs

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

0001 function [ H ] = amp_fac( alpha, beta,r,s )
0002 
0003 % first check user inputs
0004 [a1,a2]=size(alpha);
0005 [b1,b2]=size(beta);
0006 n=r+s+1;
0007 if a2~=1 || b2~=1;
0008     error('input must be column vectors');
0009 elseif a1~=n || b1~=n;
0010     error('alpha and beta must be of length (r+s+1)');
0011 end
0012 
0013 syms phi;
0014 num=0; denom=0;   % initialise numerator and denominator
0015 for k=1:n;
0016     num=num+beta(k)*exp(1i*(k-r-1)*phi);
0017     denom=denom+alpha(k)*exp(1i*(k-r-1)*phi);
0018 end
0019 H=num./denom;
0020 end

Generated on Mon 18-Jan-2016 10:25:49 by m2html © 2005