胡广书第二章第11题第三问
| 分类: VC |
代码:
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;
}
}
#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 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;
}
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
前一篇:胡广书第一章第15题第3问

加载中…