function?[ba]=testinvfreqz（gwvarargin）


narginchk（49）%%輸入變量個數(shù)范圍，根據(jù)變量個數(shù)設(shè)置函數(shù)的初始化參數(shù)
if?ischar（varargin{1}）
????realStr?=?lower（varargin{1}）;??%
????varargin（1）?=?[];??????????????%
else
???realStr?=?‘real‘;
end
gaussFlag?=?length（varargin）>3;
if?length（varargin）<6
????varargin{6}?=?[];
end
%%檢查H是否位可用數(shù)據(jù)
signal.internal.sigcheckfloattype（g‘‘‘infreqz‘‘H‘）;

%強制的精度規(guī)則
w?=?signal.internal.gigcasttofloat（w‘double‘‘invfreqz‘‘W‘‘a(chǎn)llownumeric‘）;
[nbnawfmaxitertolpf]?=?deal（varargin{:}）;


nb?=?signal.internalsigcasttofloat（nb‘double‘‘invfreaz‘‘NB‘...
????‘a(chǎn)llownumeric‘）;
na?=?signal.internal.signcasttofloat（na‘double‘‘invfreqz‘‘NA‘...
????‘a(chǎn)llownumeric‘）;
wf?=?signal.internal.sigcasttorfloat（wf‘double‘‘invfreqz‘‘Wt‘...
????‘a(chǎn)llownumeric‘）;
maxiter?=?signal.internalsigcasttofloat（maxiter‘double‘‘invfreqz‘...
????‘ITER‘‘a(chǎn)llownumeric‘）;
tol?=?signal.internalsigcasttofloat（tol‘double‘‘invfreqz‘‘TOL‘....
????‘a(chǎn)llownumeric‘）;
????
switch?realStr
????case?‘real‘
????????realFlag?=?1;
????case?‘complex‘
????????realFlag?=?0;
????otherwise
????????warning（message（‘signal:invfreqz:InvalidParam‘realStr））;
????????realFlag?=?0;
end



nk=0;T=1;%取樣間隔位T，分子nk從0開始

nb?=?nb?+?nk?+?1;
if?isempty（pf）
????verb=0;
elseif?（strcmp（pf‘trace‘））
????verb?=?1;
else
????error（message（‘signal:invfreqz:NotSupported‘pf））;
end

if?isempty（wf）wf=ones（length（w）1）;end
wf=sqrt（wf）;

if?length（g）~=length（w）error（message（‘signal:invfreqz:InvalidDimensions‘‘H‘‘W‘））end
if?length（wf）~=length（w）error（message（‘signal:invfreqz:InvalidDimensions‘‘H‘‘W‘））end
if?any（（w>pi））|（w<0））&&realFlag
????warning（message（‘signal:invfreqz:InvalindRegion‘‘W‘‘INVFREQZ‘‘‘‘complex‘‘‘））
end

[rwcw]=size（w）;???if?rw>cw??w=w‘;???end
[rgcg]?=?size（g）;?if?cg[rwfcwf]=size;????if?cwf>rwfwf=wf‘;?end


nm=max（nanb+nk-1）;
OM=exp（-li*（0:nm）‘*w*T）;

%最小二乘平方估計
Dva=（OM（2:na+1:）.‘）.*（g*ones（1:na））;????%對
Dvb=-（OM（nk+1:nk+nb:）.‘）;???????????????%
D=[Dva?Dvb].*（wf*ones（1na+nb））;?????????%

if?realFlag
????R=real（D‘*D）;????????????????????????%取結(jié)果實數(shù)部分
????Vd=real（D‘*（-g.*wf））;
else?
????R=D‘*D;?????????????????????????????%
????Vd=D‘*（-g.*wf）;?????????????????????%
end

th=R\Vd;
a=[1?th（1:na）.‘];b=[zeros（1nk）?th（na+1:na+nb）.‘];

if?~gaussFlagreturnend

%設(shè)置默認值
if?isempty（maxiter）maxiter?=?30;end

if?isempty（tol）
????tol?=?0.01;
end
indb=1:length（b）;indg=1:length（a）;
a=polystab（a）;???????????%穩(wěn)定多項式
b=polystab（b）;

%初步估計
GC=（（b*OM（indb:））./（a*OM（indg:）））.‘;
e=（GC-g）.*wf;
Vcap=e‘*e;?t=[a（2:na+1）?b（nk+1:nk+nb）].‘;
if（verb）
????clcdisp（[‘‘?getString（message（‘signal:invfreqs:INTIALESTIMATE‘））]）;
????disp（[getString（message（‘signal:invfreqs:CurrentFit‘））?‘‘?num2str（Vcap）]）;
????disp（getString（message（‘signal:invfreqs:Parvector‘）））;
????disp（t）
end


%最小化循環(huán)

gndir=2*tol+1;?l=0;st=0;
while?all（[norm（gndir）>tol?1????l=l+1;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4803??2018-04-03?09:55??2\invfreqz20180320.m

?????文件???????8931??2018-04-11?11:22??2\myinvfreqz.m

?????文件??????????0??2018-03-22?11:21??2\t1.m

?????文件???????3381??2018-01-22?14:25??2\test0505_3a_lmsfb1_invsz_0117.m

?????文件???????4404??2018-05-25?18:04??2\test0505_3a_lmsfb1_invsz_Soff_0118.m

?????文件???????3032??2018-03-05?17:47??2\test0505_3_lmsff_invsz_Soff_0118.m

?????文件???????4504??2018-01-22?18:37??2\test0505_3_lmsff_lmsfb1_invsz_0122.m

?????文件???????4969??2018-05-23?14:27??2\test0505_3_lmsff_lmsfb1_invsz_Soff_0122.m

?????文件???????1374??2018-05-07?16:43??2\test20180305.m

?????文件????????851??2018-03-20?15:43??2\test20180314.m

?????文件???40935424??2017-04-02?00:38??2\traffic_noise.WAV

?????文件????????787??2018-05-22?18:11??2\voice.m

?????目錄??????????0??2018-05-23?18:16??2

-----------?---------??----------?-----??----

?????????????40972460????????????????????13