系统辨识与建模：系统辨识与建模

系统辨识与建模：系统辨识与建模 系统辨识与建模：系统辨识与建模 系统辨识与建模：系统辨识与建模

  #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];  //产生有色噪声

	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

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

    ......R         0  2003-11-11 08:50  1.白噪声的产生\Debug\WNoise.txt

    ......R    241756  2003-06-28 10:45  1.白噪声的产生\Debug\White.exe

    ......R    290560  2003-06-28 10:45  1.白噪声的产生\Debug\White.ilk

    ......R      9951  2003-06-28 10:45  1.白噪声的产生\Debug\White.obj

    ......R    281632  2003-06-28 10:45  1.白噪声的产生\Debug\White.pch

    ......R    451584  2003-06-28 10:45  1.白噪声的产生\Debug\White.pdb

    ......R     41984  2003-06-28 10:45  1.白噪声的产生\Debug\vc60.idb

    ......R     61440  2003-06-28 10:45  1.白噪声的产生\Debug\vc60.pdb

    ...D..R         0  2003-11-11 08:58  1.白噪声的产生\Debug

    ......R       704  2003-06-23 16:01  1.白噪声的产生\WHITE.CPP

    ......R      5089  2003-06-28 10:45  1.白噪声的产生\WNoise.txt

    ......R      3389  2003-06-28 10:45  1.白噪声的产生\White.dsp

    ......R       535  2003-06-28 10:45  1.白噪声的产生\White.dsw

    ......R     33792  2003-06-28 10:45  1.白噪声的产生\White.ncb

    ......R     48640  2003-06-28 10:45  1.白噪声的产生\White.opt

    ......R       707  2003-06-28 10:45  1.白噪声的产生\White.plg

    ...D..R         0  2003-11-11 08:58  1.白噪声的产生

    ......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个文件信息