数字信号处理DSP毕业论文 第27页

k=1;

for(j=1;j<=nz;j++)

{l=iext[j];

dtemp=ad[j]*des[l];

dnum+=dtemp;

dtemp=((double)k)*ad[j]/wt[l];

dden+=dtemp;

k=0-k;}

*pdev=dnum/dden;

nu=1;

if(*pdev>0.0) nu=-1;

*pdev*=(-(double)nu);

k=nu;

for(j=1;j<nz;j++)

{l=iext[j];

dtemp=((double)k)*(*pdev)/wt[l];

y[j]=des[l]+dtemp;

k=-k;}

if(*pdev<=dev1)

{ouch(niter);

goto 19;}

dev1=*pdev;

jchnge=0;

k1=iext[1];

kn2=iext[nz];

klow=0;nut=-nu;

j=1;

while(1)

{if(j==nzz)  yna=comp;

if(j>=nzz)   break;

kup=iext[j+1];

l=iext[j]+1;

nut=-nut;

if(j==2) y1=comp;

comp=*pdev;

if(l>=kup) goto 13;

err=gee(l,nz,ad,x,y,grid);

err=(err-des[l])*wt[l];

dtemp=((double)nut)*err-comp;

if(dtemp<=0.0) goto 13;

comp=(double)nut*err;}}}}

 << 上一页  [21] [22] [23] [24] [25] [26] [27] [28] 下一页