Código (c) [Seleccionar]
/*Declaración de funciones*/
void conv(int dim1,int x[],int dim2,int h[],int y[]);
/*Funcion main*/
int main()
{
int i;
printf("Introduzca la dimension del vector de entrada: ");
scanf("%d",&dim1);
xn = malloc(dim1 * sizeof (int));
printf("Introduce los valores del primer vector:\n");
for (i=0; i<dim1; i++) {
printf("\nElemento %d: ", i+1);
scanf("%d", &aux);
*(xn+i)= aux;
}
printf("\n");
printf("Introduzca la dimension del segundo vector: ");
scanf("%d",&dim2);
hn = malloc(dim2 * sizeof (int));
printf("Introduce los valores del segundo vector:\n");
for (i=0; i<dim2; i++) {
printf("\nElemento %d: ", i+1);
scanf("%d", &aux);
*(hn+i)=aux;
}
yn = malloc((dim1+dim2-1) * sizeof (int));
if(dim1<=dim2)
{
int dimcirc=dim2;
conv(dim1,xn,dim2,hn,yn);
printf("Dimension de la Convolucion Circular=%d \n" , dimcirc);
zn = malloc(dimcirc * sizeof (int));
zi=0;
wn = malloc(dimcirc * sizeof (int));
k=dim2-dim1-2;
q=dim1+dim2-2;
for(j=dimcirc-1;j>=0;j--)
{
*(zn+j)=0;
*(wn+j)=0;
*(wn+j)=*((yn+(q%dimcirc))+j);
*(zn+j)=*((yn+(k%dimcirc))-j);
}
for(zi=0;zi<dimcirc;zi++)
{
printf("zn:[%d]=%d \n",zi,(*(zn-zi)+*(wn-zi)));
}
}
else
{
int dimcirc=dim1;
paux = xn;
xn = hn;
hn = paux;
conv(dim2,xn,dim1,hn,yn);
printf("Dimension de la Convolucion Circular=%d \n" , dimcirc);
zn = malloc(dimcirc * sizeof (int));
zi=0;
wn = malloc(dimcirc * sizeof (int));
k=dim1-dim2-2;
q=dim1+dim2-2;
for(j=dimcirc-1;j>=0;j--)
{
*(zn+j)=0;
*(wn+j)=0;
*(wn+j)=*((yn+(q%dimcirc))+j);
*(zn+j)=*((yn+(k%dimcirc))-j);
}
for(zi=0;zi<dimcirc;zi++)
{
printf("zn:[%d]=%d \n",zi,(*(zn-zi)+*(wn-zi)));
}
}
}
void conv(int dim1,int x[],int dim2,int h[],int y[])
{
/*Declaro e Inicializo las variables contador de los vectores*/
int xi=0,yi=0,hi=0;
int dimres=dim1+dim2-1;
printf("Dimensión de la Convolución Lineal=%d \n",dimres);
for(xi=0;xi<dim1;xi++)
{
for(hi=0;hi<dim2;hi++)
{
y[xi+hi]=(y[xi+hi]+x[xi]*h[hi]);
}
}
for(yi=0;yi<dimres;yi++)
{ printf("yn:[%d]=%d \n",yi,y[yi]);
}
}