??#include?“stdlib.h“
??#include?“math.h“
??#include?“stdio.h“
??void?brmul（abmnkc）
??int?mnk;
??double?a[]b[]c[];
??{?int?ijlu;
????for?（i=0;?i<=m-1;?i++）
????for?（j=0;?j<=k-1;?j++）
??????{?u=i*k+j;?c[u]=0.0;
????????for?（l=0;?l<=n-1;?l++）
??????????c[u]=c[u]+a[i*n+l]*b[l*k+j];
??????}
??}
??int?brinv（double?a[]int?n）
??{?int?*is*jsijkluv;
????double?dp;
????is=malloc（n*sizeof（int））;
????js=malloc（n*sizeof（int））;
????for?（k=0;?k<=n-1;?k++）
??????{?d=0.0;
????????for?（i=k;?i<=n-1;?i++）
????????for?（j=k;?j<=n-1;?j++）
??????????{?l=i*n+j;?p=fabs（a[l]）;
????????????if?（p>d）?{?d=p;?is[k]=i;?js[k]=j;}
??????????}
????????if?（d+1.0==1.0）
??????????{?free（is）;?free（js）;?printf（“err**not?inv\n“）;
????????????return（0）;
??????????}
????????if?（is[k]!=k）
??????????for?（j=0;?j<=n-1;?j++）
????????????{?u=k*n+j;?v=is[k]*n+j;
??????????????p=a[u];?a[u]=a[v];?a[v]=p;
????????????}
????????if?（js[k]!=k）
??????????for?（i=0;?i<=n-1;?i++）
????????????{?u=i*n+k;?v=i*n+js[k];
??????????????p=a[u];?a[u]=a[v];?a[v]=p;
????????????}
????????l=k*n+k;
????????a[l]=1.0/a[l];
????????for?（j=0;?j<=n-1;?j++）
??????????if?（j!=k）
????????????{?u=k*n+j;?a[u]=a[u]*a[l];}
????????for?（i=0;?i<=n-1;?i++）
??????????if?（i!=k）
????????????for?（j=0;?j<=n-1;?j++）
??????????????if?（j!=k）
????????????????{?u=i*n+j;
??????????????????a[u]=a[u]-a[i*n+k]*a[k*n+j];
????????????????}
????????for?（i=0;?i<=n-1;?i++）
??????????if?（i!=k）
????????????{?u=i*n+k;?a[u]=-a[u]*a[l];}
??????}
????for?（k=n-1;?k>=0;?k--）
??????{?if?（js[k]!=k）
??????????for?（j=0;?j<=n-1;?j++）
????????????{?u=k*n+j;?v=js[k]*n+j;
??????????????p=a[u];?a[u]=a[v];?a[v]=p;
????????????}
????????if?（is[k]!=k）
??????????for?（i=0;?i<=n-1;?i++）
????????????{?u=i*n+k;?v=i*n+is[k];
??????????????p=a[u];?a[u]=a[v];?a[v]=p;
????????????}
??????}
????free（is）;?free（js）;
????return（1）;
??}

void?main（）
{

	double?z[607]zz[607][1]w[4][1]ww[1][4]p[4][4]I[4][4]k[4][1]kk[4][1]ss[4][1]pp[4][4]X[4][1]XX[1][4]x1[807];
	double?s[4][1]a[4][1]b[1]c[1][1]f[4][4]g[4][4]x2[1][1];//ee[1][1]
????double?u1v1dq;
	int?ijN;
	double?u[607]e[607]v[607];

????FILE?*fp1*fp2*fp3*fp;

	if（（fp1=fopen（“m.txt““r“））==NULL）
????{
??????printf（“ERROR“）;exit（1）;
????}
????if（（fp2=fopen（“wnoise.txt““r“））==NULL）
????{
?????printf（“ERROR“）;exit（1）;
????}
????if（（fp3=fopen（“IV.txt““w“））==NULL）
????{
?????printf（“ERROR“）;exit（1）;
????}
	if（（fp=fopen（“result.txt““w“））==NULL）
????{
?????printf（“ERROR“）;exit（1）;
????}
????for（i=0;i<607;i++）
????{
?????fscanf（fp1“%lf?“&u1）;
?????u[i]=u1;
?????fscanf（fp2“%lf?“&v1）;
?????e[i]=v1;
????}

????e[-1]=e[-2]=0.0;
	for（i=0;i<607;i++）
		v[i]=e[i]-1.0*e[i-1]+0.2*e[i-2];??//產(chǎn)生有色噪聲

	z[0]=v[0];
	z[1]=1.5*z[0]+u[0]+v[1];?
	for（i=2;i<607;i++）
	z[i]=1.5*z[i-1]-0.7*z[i-2]+u[i-1]+0.5*u[i-2]+v[i];


	for（i=0;i<607;i++）
	{
		zz[i][0]=z[i];
	
	}
	for（i=0;i<4;i++）
	{
		s[i

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

????......R?????????0??2003-11-11?08:50??1.白噪聲的產(chǎn)生\Debug\WNoise.txt

????......R????241756??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\White.exe

????......R????290560??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\White.ilk

????......R??????9951??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\White.obj

????......R????281632??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\White.pch

????......R????451584??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\White.pdb

????......R?????41984??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\vc60.idb

????......R?????61440??2003-06-28?10:45??1.白噪聲的產(chǎn)生\Debug\vc60.pdb

????...D..R?????????0??2003-11-11?08:58??1.白噪聲的產(chǎn)生\Debug

????......R???????704??2003-06-23?16:01??1.白噪聲的產(chǎn)生\WHITE.CPP

????......R??????5089??2003-06-28?10:45??1.白噪聲的產(chǎn)生\WNoise.txt

????......R??????3389??2003-06-28?10:45??1.白噪聲的產(chǎn)生\White.dsp

????......R???????535??2003-06-28?10:45??1.白噪聲的產(chǎn)生\White.dsw

????......R?????33792??2003-06-28?10:45??1.白噪聲的產(chǎn)生\White.ncb

????......R?????48640??2003-06-28?10:45??1.白噪聲的產(chǎn)生\White.opt

????......R???????707??2003-06-28?10:45??1.白噪聲的產(chǎn)生\White.plg

????...D..R?????????0??2003-11-11?08:58??1.白噪聲的產(chǎn)生

????......R?????????0??2003-11-06?15:23??10.遞推廣義最小二乘法\Debug\Gauss.txt

????......R????245857??2003-06-28?10:50??10.遞推廣義最小二乘法\Debug\Gls.exe

????......R????291656??2003-06-28?10:50??10.遞推廣義最小二乘法\Debug\Gls.ilk

????......R?????16741??2003-06-28?10:50??10.遞推廣義最小二乘法\Debug\Gls.obj

????......R????297040??2003-06-28?10:50??10.遞推廣義最小二乘法\Debug\Gls.pch

????......R????451584??2003-06-28?10:50??10.遞推廣義最小二乘法\Debug\Gls.pdb

????......R??????7157??2003-03-25?17:33??10.遞推廣義最小二乘法\Debug\PRBS.txt

????......R?????34550??2003-11-07?10:35??10.遞推廣義最小二乘法\Debug\mdata.txt

????......R?????27350??2003-11-07?10:35??10.遞推廣義最小二乘法\Debug\mdata1.txt

????......R?????21348??2003-11-07?10:35??10.遞推廣義最小二乘法\Debug\mdata2.txt

????......R?????21348??2003-11-07?10:35??10.遞推廣義最小二乘法\Debug\mdata3.txt

????......R?????21348??2003-11-07?10:35??10.遞推廣義最小二乘法\Debug\mdata4.txt

????......R?????50176??2003-06-28?10:50??10.遞推廣義最小二乘法\Debug\vc60.idb

............此處省略553個(gè)文件信息