first check user inputs
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