Question

卷积积分的C语言程序

#include"stdio.h"
#include"math.h"
#define H(x) (x)/3
int kn;
static float ss[301],ai,c,a[3],g;
//打印子程序
void gw11(g)
float g[301];
{
char ch[61];
float t,gmax,gmin,d;
int i,lp,nx;
ch[1]='|';
for(i=2;i<=60;i++)
ch[i]='';
t=a[1]+c;
gmin=gmax=g[1];
for(i=1;i<=300;i++)
{ if(g[i]>gmax)
gmax=g[i];
if(g[i]<gmin)
gmin=g[i];
}
d=gmax-gmin;
printf("\n");
for(i=1;i<=61;i++)
printf("=");
printf("T G(t)\n");
for(lp=1;lp<=kn;lp++)
{
nx=(g[lp]-gmin)/d*59+1;
ch[nx]='*';
for(i=1;i<=60;i++)
printf("%c",ch[i]);
printf("%6.2f%8.3f\n",t,g[lp]);
ch[nx]=' ';
ch[1]='|';
t+=ai;
}
}
/*主函数*/
main()
{
static float t[301],s[301],b[3],hm,bi,d,p,q,z,fm;
int nf,i,kk,l,ll,km,k,m,kb;
printf("\nF(t)的段数，卷积步长，积分步长 \n");
scanf("%d%f%f,&nf,&ai,&bi");
for(i=1;i<=nf;i++)
{
printf("\ nf(t)的第%d段的下.上界："，i);
scanf("%f%f",&a[i],&b[i]);
}
printf(" \ nH(t)的上下界：")；
scanf("%f%f",&c,&d);
kk=1;
for(l=1;l<=nf;l++)
{
km=((b[l]+d)-(a[l]-c))/ai;
//确定每一步积分的上.下限
for(k=1;k<=km;k++)
{
t[k]=k*ai+a[l];
if(t[k]>b[l])
{
if((t[k]-d+c)>a[l])
{p=t[k]-d+c;
q=b[l];
}
else
{p=a[l];
q=b[l];
}
}
else
{
if((t[k]-d+c)>a[l])
{
p=t[k]-d+c;
q=t[k];
}
else
{
p=a[l];
q=t[k];
}
}
/*积分计算*/
ll=(q-p)/bi;
s[k]=0;
for(m=1;m<=ll;m++)
{
z=p+m*bi;
fm=1.0;
hm=H(t[k]-z);
sk+=fm*hm*bi;
}
/*卷积结果累加*/
ss[kk]+=s[k];
kk++;
}
kn=kk;
kb=((b[l]+d)-(a[l+1]+c))/ai;
kk-=kb;
}
printf(" \ n打印卷积结果及其曲线 \ n");
gw11(ss);
}
哪里出错了，谢谢，小妹谢了

Answer 1

代码并不长,建议单步运行去看看,看走到哪一步时结果与你想象的不一致,再分析那一段代码...