Výpočet konvoluce v C – C / C++ – Fórum – Programujte.com

Potřebuji algoritmus pro výpočet konvoluce. Mám 2 pole: Signál např. o n=4 prvcích a jádro např. o x=3 prvcích.

Signál je pole libovol. kladných čísel a jádro je pole lichého počtu prvků. 

Jádro[3] = 1,2,1

Jádro[5] = 1,2,4,2,1

Nevím jak to správně napsat. Vím, že pomocí cyklů jen stále nemohu přijít jak na to.

 V C sa my to nechce robit, ale v c++ to vyzera asi takto:

#include <iostream>
#include <vector>
using namespace std;

int main(){
	vector<int> signal{1, 7, 4, 0};
	vector<int> jadro{1, 2, 1};
	vector<int> vysledok(signal.size() + jadro.size() - 1, 0);
	
	for(size_t s = 0; s < signal.size(); ++s)
		for(size_t j = 0; j < jadro.size(); ++j)
			vysledok[s+j] += signal[s] * jadro[j];
	
	
	for(auto v : vysledok)cout << v << ", ";
}

pouzil som algoritmus z http://cs.wikipedia.org/wiki/Konvoluce

#1 Grommi

#include <stdio.h>
#define delkaSignal 4
#define delkaJadro 3

int main()
{
    int poleSignal[delkaSignal+delkaJadro]={0,1,7,4,0};
    int poleJadro[delkaJadro]={1,2,1};
    int poleVysl[delkaSignal]={0};

    for (short j=0;j<delkaSignal;j++)
    {
        for (short i=0;i<delkaJadro;i++)
        {
            poleVysl[j]+=poleSignal[i+j]*poleJadro[i];
        }
        printf("%d, ",poleVysl[j]);
    }
    return 0;
}

#3 Jiří G
Ten tvoj kod je daky divny a mas v poleSignal jednu nulu navyse.

Zkoušel jsem to a pořád to nefunguje. Koukal by na to někdo prosím ?! Jen kde může být chyba.

Zdroják: http://www11.zippyshare.com/v/46748158/file.html

Děkuju za pomoc.

	//vector<int> signal{1, 7, 4, 0};
	int signal[4] = {1, 7, 4, 0};
	//vector<int> jadro{1, 2, 1};
	int jadro[3] = {1, 2, 1};
	//vector<int> vysledok(signal.size() + jadro.size() - 1, 0);
	int vysledok[4 + 3 - 1];
	for(size_t i = 0; i < (4 + 3 - 1); ++i)vysledok[i] = 0; //vynulovanie vysledku 
	
	//for(size_t s = 0; s < signal.size(); ++s)
	//	for(size_t j = 0; j < jadro.size(); ++j)
	//		vysledok[s+j] += signal[s] * jadro[j];
	
	for(size_t s = 0; s < 4; ++s)
		for(size_t j = 0; j < 3; ++j)
			vysledok[s+j] += signal[s] * jadro[j];
	
	//vypis vysledku, prvy a posledny prvok vysledku sa ignoruje:
	for(size_t i = 1; i < (4 + 3 - 1) -1; ++i)printf("%d, ", vysledok[i]);
}