In　order　to　understand　the　evolution　of　chemical　structure　in　glasses,　Raman　features　were　fitted　into　different　peak-fitting　functions,　including　Gaussian　and　Lorentzian　types.　It　was　found　that,　while　single　function　was　used　for　peak-fitting,　Gaussian　function　showed　better　fitted　curves　in　overlapping　bands　than　Lorentzian　function,　and　Lorentzian　function　was　more　suitable　for　single　broad　band.　On　the　other　hand,　while　multiple　functions　were　utilized　to　fit　a　vibrational　peak,　the　results　showed　better　matching　between　fitted　curve　and　raw　data.　Therefore　the　Raman　spectra　were　decomposed　into　different　structural　units.　It　was　found　that　the　evolution　of　different　structural　units　was　correlated　with　that　of　the　physical　parameters　reported　in　literatures　in　a　series　of　GeAsSe　glasses.　Therefore,　it　was　concluded　that　the　proposed　method　could　be　employed　to　decompose　broad　Raman　spectra　in　amorphous　materials　so　that　the　chemical　structure　of　the　material　could　be　understood　better.