加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

胡广书第二章第11题第三问

(2022-11-21 19:18:32)
分类: VC
代码:
#include "stdafx.h"
#include "MSP.h"
#include "stdio.h"
#include

#include

#define N 40

float b[2]={0,-0.2};
float a[3]={1,0,0.8};
float y[N];
float x[N];
int n=N;

int lb=1;
int la=2;


void  mfitout(float b[],float a[],int lb,int la,float x[],int n,
float y[])
{
        int k,i,ki;
        float sum;
       for(k=0;k
           {
            sum=0.0;
            for(i=0;i<=lb;i++)
               {
                ki=k-i;
                if(ki<0)break;
                if(fabs(b[i])
                sum=sum+b[i]*x[ki];
                if(sum>1.e10)exit(3);
                }
       if(la!=0)
              for(i=1;i<=la;i++)
                {
                 ki=k-i;
                 if(ki<0)break;
         if(fabs(a[i])
                  sum=sum-a[i]*y[ki];
          if(sum>1.e10)exit(3);
                 }
             y[k]=sum;
            }
        }




int main(int argc, char* argv[])
{
FILE *fp;
char i=0;
if((fp=fopen("file2","w"))==NULL)
{
printf("cannot open this file\n");
exit(0);
}
for(i=0;i
{
x[i]=1;
}
mfitout(b,a,lb,la,x,n,y);
for(i=0;i
{
fprintf(fp,"%f\n",y[i] );
}
fclose(fp);
return 0;
}



胡广书第二章第11题第三问
数据:
0.000000
-0.200000
-0.200000
-0.040000
-0.040000
-0.168000
-0.168000
-0.065600
-0.065600
-0.147520
-0.147520
-0.081984
-0.081984
-0.134413
-0.134413
-0.092470
-0.092470
-0.126024
-0.126024
-0.099181
-0.099181
-0.120655
-0.120655
-0.103476
-0.103476
-0.117220
-0.117220
-0.106224
-0.106224
-0.115020
-0.115020
-0.107984
-0.107984
-0.113613
-0.113613
-0.109110
-0.109110
-0.112712
-0.112712
-0.109830

0

阅读 收藏 喜欢 打印举报/Report
  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有