convolve.txt

实现信号通过滤波器的功能

#include<graphics.h>
#include<math.h>
#include<stdio.h>
#include<conio.h>
float x[3],h[4],y[7];
pro1()
{
int n,i;
float c,k;
printf("please input 3 words for x(n) end with Enter.\n");
for(i=0;i<=2;i++)
{scanf("%f",&c);
x[i]=c;}
printf("please input 4 words for h(n) end with Enter.\n");
for(i=0;i<=3;i++)
{scanf("%f",&k);h[i]=k;}for(n=0;n<=2;n++)
    {
    printf("%f \n",x[n]);
    }getch();
for(n=0;n<=3;n++)
    {
     printf("%f \n",h[n]);
    }getch();}
void show()
{
int i=0;
int driver=DETECT,MODE;
initgraph(&driver,&MODE,"");
setbkcolor(0);
cleardevice();
for(i=0;i<=2;i++)
{
    setfillstyle(SOLID_FILL,RED);
    bar((i+1)*10+20,(int)(400-x[i]*10),(i+1)*10+22,400);
}
printf("x(n)");getch();
cleardevice();
for(i=0;i<=3;i++)
{
    setfillstyle(SOLID_FILL,GREEN);
    bar((i+1)*10+20,(int)(400-h[i]*10),(i+1)*10+22,400);
}
printf("h(n)");getch();
cleardevice();
for(i=0;i<=5;i++)
{
    setfillstyle(SOLID_FILL,CYAN);
    bar((i+1)*10+20,(int)(400-y[i]),(i+1)*10+22,400);
}
printf("y(n)=x(n)*h(n);");getch();
cleardevice();
printf("y(n)000=x(n)*h(n)\n");
for(i=0;i<=5;i++)
printf("%f    ",y[i]);
getch();}
convolve()
{
    int i,j;
    for(i=0;i<=5;i++)
    {

       float t=0.0;
       for(j=0;j<=2;j++)
         {
         if(i-j>=0&&i-j<=3)
        t+=x[j]*h[i-j];
      }
      y[i]=t;
     }
}
main()
{
pro1();
convolve();
show();
}