Как построить АЧХ и ФЧХ (ГВЗ) КИХ фильтра?, по известным коэффициентам Опции

[_Ivana]

Есть рассчитанные коэффициенты КИХ-фильтра. Как по ним определить и построить его графики АЧХ и ФЧХ? Можно ссылку если это где-то разжевано.

ЗЫ матлаб навскидку скачать не удалось, есть GNU Octave (с интерфейсом командной строки), GNU Plot.

[_Ivana]

160 значений окна при определенном параметре. АЧХ после центрального лепестка -180дБ, не знаю хорошо или плохо. Желающие могут проанализировать а матлабе.

CODE

h = [
0.000000000549916,
0.000000000943069,
0.000000001605766,
0.00000000271468,
0.000000004556784,
0.000000007594641,
0.000000012568122,
0.000000020651584,
0.00000003369469,
0.000000054588341,
0.000000087816027,
0.000000140277549,
0.000000222509216,
0.000000350475889,
0.000000548180237,
0.000000851428878,
0.000001313220812,
0.000002011388839,
0.000003059339326,
0.000004621010446,
0.000006931515669,
0.000010325369814,
0.000015274720634,
0.000022440638956,
0.000032741260117,
0.00004744141772,
0.000068269357207,
0.000097567138033,
0.000138482389466,
0.000195210115272,
0.00027329416138,
0.000379998654668,
0.000524760046923,
0.000719730182933,
0.000980419855492,
0.001326450392725,
0.001782417715224,
0.002378868779462,
0.003153384196031,
0.004151752937631,
0.005429215379979,
0.007051739516525,
0.009097282277044,
0.011656973868303,
0.014836148559658,
0.018755131198059,
0.023549676015608,
0.029370944244185,
0.036384901078617,
0.044771012107889,
0.054720125909643,
0.066431444370009,
0.080108506446187,
0.095954145083895,
0.114164420793247,
0.134921588210201,
0.158386212246758,
0.184688615724559,
0.213919907421971,
0.246122904238181,
0.28128331912935,
0.319321632748266,
0.360086096503364,
0.40334732364511,
0.44879490940794,
0.496036478819302,
0.544599490731681,
0.593936029961729,
0.643430699125206,
0.692411582795637,
0.74016410574706,
0.785947452494301,
0.829013066329605,
0.868624612161054,
0.904078677963546,
0.934725412809768,
0.959988261804627,
0.979381964063306,
0.992528030695028,
0.999167014177938,
0.999167014177938,
0.992528030695028,
0.979381964063306,
0.959988261804627,
0.934725412809768,
0.904078677963546,
0.868624612161054,
0.829013066329605,
0.785947452494301,
0.74016410574706,
0.692411582795637,
0.643430699125206,
0.593936029961729,
0.544599490731681,
0.496036478819302,
0.44879490940794,
0.40334732364511,
0.360086096503364,
0.319321632748266,
0.28128331912935,
0.246122904238181,
0.213919907421971,
0.184688615724559,
0.158386212246758,
0.134921588210201,
0.114164420793247,
0.095954145083895,
0.080108506446187,
0.066431444370009,
0.054720125909643,
0.044771012107889,
0.036384901078617,
0.029370944244185,
0.023549676015608,
0.018755131198059,
0.014836148559658,
0.011656973868303,
0.009097282277044,
0.007051739516525,
0.005429215379979,
0.004151752937631,
0.003153384196031,
0.002378868779462,
0.001782417715224,
0.001326450392725,
0.000980419855492,
0.000719730182933,
0.000524760046923,
0.000379998654668,
0.00027329416138,
0.000195210115272,
0.000138482389466,
0.000097567138033,
0.000068269357207,
0.00004744141772,
0.000032741260117,
0.000022440638956,
0.000015274720634,
0.000010325369814,
0.000006931515669,
0.000004621010446,
0.000003059339326,
0.000002011388839,
0.000001313220812,
0.000000851428878,
0.000000548180237,
0.000000350475889,
0.000000222509216,
0.000000140277549,
0.000000087816027,
0.000000054588341,
0.00000003369469,
0.000000020651584,
0.000000012568122,
0.000000007594641,
0.000000004556784,
0.00000000271468,
0.000000001605766,
0.000000000943069,
0.000000000549916];

UPD добавил картинку из Октавы. (офф - быстро не смог, какая-то авария с компом, антивирус пытается справиться, жрет 100% проца, остальное работает в час по чайной ложке)

Насчет "определиться что хочу от окна", сначала имхо стОит его проанализировать по всякому, посмотреть возможные плюсы и минусы, а потом уже думать как применять и чего хотеть.

Сообщение отредактировал _Ivana - Jul 20 2012, 18:52

Эскизы прикрепленных изображений

 

[ToR_TDA]

160 значений окна при определенном параметре. АЧХ после центрального лепестка -180дБ, не знаю хорошо или плохо. Желающие могут проанализировать а матлабе.

Может сократите количество тапов в фильтре, а то тяжело со 160 значениями оперировать, да и не нужно это для сравнения. Можно попробовать сравнить например ваше окно и Гаусса. К вопросу о -180Дб: подавление боковых лепестков будет неизбежно приводить к расширению основного. Можно хоть 360Дб окном Гаусса получить, но основной лепесток расползется. Тут нужно искать компромисс.
P.S. Что-то я не понял, по вашим графикам полоса пропускания вашего фильтра равна 0. Может все же сделать 0.5. Фильтр с окном Гаусса ниже.

Эскизы прикрепленных изображений