Научная статья на тему 'Метод аппроксимации синусоидального равноамплитудного полинома'

Метод аппроксимации синусоидального равноамплитудного полинома Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
661
58
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Фролов С. С., Шевеленко В. Д., Гусаров А. А.

В статье предложен один из методов аппроксимации равноамплитудных тригонометрических полиномов, обладающих равномерным неизменным дискретным спектром, метод аппроксимации синусоидального полинома. Получены аналитические выражения для дискретного гармонического спектра аппроксимирующей функции, просчитаны его определяющие характеристики -коэффициент гармоник паразитного спектра и неравномерность амплитуд гармоник в рабочей области спектра.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Метод аппроксимации синусоидального равноамплитудного полинома»

Фролов C.C., Шевеленко В.Д., Гусаров А.А.

Оренбургский государственный университет

МЕТОД АППРОКСИМАЦИИ СИНУСОИДАЛЬНОГО РАВНОАМПЛИТУДНОГО ПОЛИНОМА

В статье предложен один из методов аппроксимации равноамплитудных тригонометрических полиномов, обладающих равномерным неизменным дискретным спектром, - метод аппроксимации синусоидального полинома. Получены аналитические выражения для дискретного гармонического спектра аппроксимирующей функции, просчитаны его определяющие характеристики -коэффициент гармоник паразитного спектра и неравномерность амплитуд гармоник в рабочей области спектра.

Статья посвящена одному из методов аппроксимации равноамплитудных полиномов, частично освещенных в работе [1, с. 166175], - методу аппроксимации синусоидального полинома вида

sin2 (х)

DsN(х) = —(2),N = 3,4,... . (1)

sinlf)

Результаты исследований, приведенные в статье, предполагается использовать для разработки методов формирования испытательных сигналов (ИС) в устройствах контроля и измерения частотных характеристик (ЧХ) узлов радиоэлектронной аппаратуры и средств телекоммуникаций в области инфра-низких частот (ОИНЧ).

Измерение, контроль и идентификация ЧХ четырехполюсников в указанном диапазоне современными способами малопроизводительны, трудоемки и не отличаются высокой точностью. Так, например, высокоточные ручные методы анализа, в которых ЧХ анализируется путем последовательного воздействия на исследуемый четырехполюсник гармоническими сигналами разных частот,

- требуют сложных систем стабилизации параметров ИС и значительных трудозатрат -до нескольких десятков человеко-часов.

Известны также методы автоматизированного анализа и контроля ЧХ, также основанные на последовательном воздействии гармоническими сигналами генератора с автоматической перестройкой частоты [2, 3, 4, 5]. Но в ОИНЧ они не позволяют повысить производительность измерений ЧХ при той же высокой точности предыдущих методов.

Рабочий диапазон частот современных автоматизированных измерителей амплитудно-частотных характеристик (АЧХ) с приме-

нением генераторов качающейся частоты (ГКЧ) ограничен снизу частотой 20 Гц [6]. Ограниченность вызвана значительной величиной неравномерности спектра испытательного сигнала на выходе ГКЧ - частотно-мо-дулированного импульса в полосе качания. Известно, что степень неравномерности спектральной плотности одиночного линейно частотно-модулированного (ЛЧМ) импульса напрямую связана с величиной так называемой базы [7]

В = А/ -ти, (2)

где А/ - девиация частоты ЛЧМ-импульса

или ширина полосы измерения АЧХ;

Ти - длительность ЛЧМ-импульса.

С увеличением величины базы В спектральная плотность ЛЧМ-импульса становится более равномерной в полосе качания [7]. В одном из лучших измерителе АЧХ Х1-41 минимально возможная величина ширины полосы качания - 100 Гц, нижняя граничная частота полосы - / = 20 Гц, максимальное значение длительности рабочего хода развертки (длительности ЛЧМ-импульса) - 10 с [6]. Соответственно величина базы не превысит значения В = 1000, при котором обеспечивается паспортная величина неравномерности уровня выходного напряжения ГКЧ 5и = 0,5 Дб, или 5,92%. При измерениях ЧХ в ОИНЧ частотный диапазон наблюдения и соответственно полоса качания А/ сузятся. Для сохранения той же равномерности напряжения на выходе ГКЧ или ее улучшения необходимо обеспечить величину базы В > 1000, то есть в соответствии с (2) - увеличить длительность ЛЧМ-импульса Ти. Это потребует на частотах /<20 Гц величин Ти порядка сотен - нескольких тысяч секунд. То есть для измерения потребуется время, соизмеримое

с временными затратами в ручных методах, а в структуре ИАЧХ потребуются дополнительные устройства стабилизации внутренних параметров измерителя для уменьшения влияния внешних дестабилизирующих факторов - колебаний температуры, влажности, напряжения внешней питающей сети и так далее. Кроме того, ЛЧМ-импульсы, другие частотно модулированные сигналы имеют ненулевой нелинейный фазовый спектр, что делает их малопригодными для измерения ФЧХ.

Кроме ЛЧМ-сигналов в современных измерителях ЧХ в качестве зондирующих сигналов применяются импульсные испытательные сигналы (ИИС) с короткой длительностью прямоугольной или произвольной формы [8], а также многочастотные сигналы [9, 10].

Использование ИИС короткой длительности основано на отсутствии зависимости спектральной характеристики дельта-функции от частоты [11]. Однако неидеальное соответствие реальных ИИС короткой длительности дельта-импульсам приводит к искажению формы их спектральной функции. Например, спектральная характеристика прямоугольного импульса короткой длительности равномерна с точностью до 5% лишь в небольшой части ее главного лепестка [12] - в диапазоне частот

А/ = 0,...,0.17 • —. (3)

Ти

При зондировании исследуемых четырехполюсников кратковременными ИИС или ИИС произвольной формы для измерения ЧХ возможно использование более широкого участка спектральной характеристики. Но в таком случае необходима дополнительная операция нормирования спектра [8] сигнала, принятого с выхода исследуемого объекта, а в конструкции измерителя ЧХ необходимо предусмотреть дополнительное запоминающее устройство с информацией о спектре ИС.

Многочастотный сигнал представляет из себя суммарный сигнал с выходов нескольких синхронизированных генераторов гармонических сигналов разных. Соответственно, количество измеряемых точек ЧХ ограничено количеством используемых генераторов. Кроме того, для изготовления одного

синусоидального генератора инфранизкой частоты требуются реактивные элементы (в основном конденсаторы) с большими значениями параметров, а соответственно и габаритными размерами, а также дополнительное устройство стабилизации частоты и амплитуды. Таким образом, структура многочастотного генератора в ОИНЧ получится очень громоздкой и ограниченой в плане универсальности применения.

В работе [1, с. 166-175] в целях усовершенствования измерителей ЧХ в ОИНЧ в плане снижения трудоемкости, повышения производительности измерений с сохранением приемлемой точности предлагается использовать в качестве генераторов ИС генераторы равноамплитудных полиномов (РАП)

- функций, разлагающихся в тригонометрический ряд Фурье с одинаковыми амплитудами и нулевыми начальными фазами [13].

В связи с невозможностью реализации генераторов РАП на аналоговых схемах из-за присутствия в полиномах точек с неопределенностью вида «ноль на ноль» в работе [1, с. 166-175] предлагается вариант решения задачи воспроизведения РАП - разработка методов аппроксимации, а также предложены общие принципы упрощенной и уточненной аппроксимации равноамплитудного косинусоидального полинома вида

Ом (х) = ~±(4, N = 3,4, 8іп1і)

(4)

Упрощенная методика характеризуется относительной простотой математического выражения аппроксимирующей функции по сравнению со второй, в которой для повышения точности необходимо увеличивать количество операций умножения. Однако и качественные характеристики спектра аппроксимирующей функции в первой методике невысоки: коэффициент гармоник кг м паразитного спектра приближается к 5% по сравнению с 0,5% в «уточненном» методе, равномерность спектра в рабочей полосе приближается к 10% против 2,2% в «уточненном».

Поэтому в настоящей статье продолжается поиск методов аппроксимации РАП, не уступающих по точности первому методу, но более простых, чем второй, и предлагается

упомянутый вначале метод аппроксимации функции (1), которая представляется суммой синусоидальных функций с одинаковыми амплитудами и нулевыми начальными фазами

N , >

(х) = £ зт((2и - 1) )N = 3 , 4 ,

(5)

1. Суть метода заключается в следующем.

1) Определяются локальные точки экстремума X 0,м и Х1 м РАП (1) соответственно на интервалах (рисунок 1а) х е (о,-2^) и х е ( ■4п), а также точка перегиба Хы на интервале хе (х0м,Х 1м). Подробнее об особенностях определения указанных точек будет сказано ниже. Ввиду нечетности РАП (1) экстремумами для интервала х е (- (п _ ж) и х е (- Хп ,о) являются точки х = -X 1м , х = — X0,м, а точкой перегиба на интервале х е (- Х1м, -Х0м) точка х = -Хм.

2) Одиночную реализацию аппроксимирующей функции на промежутке хе(- к, к\ формируем следующим образом:

а) Дляхе (-Х0 ,м,Х0 ,X) функцию (1) аппроксимируем отрезком синусоиды

D.s00N(х) = ^^0,N 81п(®0,мх), (6)

б) Для множества х е (- м,-Х0Л \и

и хе [Х0,м ,-2п) полином (1) аппроксимируется соответственно четверть-волнами синусоид

Аш^ зтХ^х + Аф1,м)- )01,м, при х е (- Хм,-Х0,м ]

Аш1,м 81п(ю1,мх -Аф1,м)+ )01,м, при х е [х0,м, Хм )

Аш2,м 81п(о)х,Xх + АФ2,М )- )01,М , пРи х е(- м ,-хх] (7)

Аш2,Х ^*1п(ю2,Xх - Аф2,N )+ )01,Х , пРи х е [хх ,тК)

в) На множестве хе (-N,-Х^и-Хп,-4п) применяются отрезки синусоид

Ds10 N (х) =

'2, N>

'2, N>

АШ3 N зтХз^х + Афз,N)- )02

при х е (- Х1^,- тЛ]

АШ3,N зтХз^х - Афз,N)+ )02 при х е , X, ы )

Б820N (х) Н N 1 \ (8)

Аш3,N 81пХ4,NX + Аф4,N )- )02,N ,

пРи х е (- N,-Х 1^ ]

Аш3, N я1п(ю (х -Аф 4, N )+ ) 0х, N ,

п^и.х е [х 1^ ,тп).

3) Параметры Аш^, ..., Аш^, 70^, У02„, ю0,№ ., ю4,„, определяются из условий

Ds00N (х0,N ) = (х0,N ) = -Ds00N (- Хо,N ) =

= -Ds10N (- Х0, N ) = DsN (Х 0,N )

-О^1^ (ХN ) = -Ds10N (- XN )= DsN (ХN )

Ds20N (X1,N ) = -Ds20N (- X1,N ) = DsN (Х 1^ )

Ds10N (2п)= Ds20N (хп)= Ds10N (- 2§)=

= Ds20N (- $= Ds20 N ((*)= Ds20N (- ^)= 0

ю,

l0,N 'Х0,N = Ю1,NX0,N АфЩ =

= Ю3,NX1,N Аф3,N = Ю4,NX1,N Аф4,N =

Ю1,NXN - АфЩ = Ю2,NXN - Аф2,N =п

,(9)

ю

2п

N

2,^ N ' 2п

^ N'

4п

'4^ N '

Аф 4,N =-:Х

Рисунок 1. Графики аппроксимирующей функции

из которых вытекает следующее решение:

П=1

Am0lN = DsN (X0iN )

Am1lN = Dsn (xo,n )— Dsn (xn )=

Am2lN = Y01іN = DsN (XN )= Dsv (xil ж ).

Am3l N = Y 02lN =

®0lN = '

2X

2

“і, N =

0і N

2 •(Xn — X o,n ) ’

2•(If — Xn Ї

“3і N ='

W

П

лл _ 2k '

X 1і N n

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

4п — X

N X 1і N

(10)

Аф2,N = Ю2,NXN -Аф3,N = Ю3,NX1,N -'|;

Аф4,N = Ю4^Х 1^ - "2 .

4) Остальные полуволны функции (1) на промежутке

х еХ-п,- ^\и[ ,п\ (11)

аппроксимируются с помощью выражения

Ds30 n (x) =

(—1)

sln

2к1

< < 2n(i + l) при------< X <----------l

N V

(— 1)+

sin2 f)

" +05) \l

2n(l +1) . . 2П

при-------- --- < x <-------.

V N '

o, при x й (—п, п]

(12)

где i =

ж '~2

V-І

2

- І при чётныхVі при нечётныхV.

Суммируя выражения (6)-(8) и (12), получим функцию, аппроксимирующую полином (1) на промежутке х е (-пп)

Ds00 (х) =

= Ds000 (х)+ Ds100 (х)+ Ds200 (х)+ Ds300 (х). (13)

Зеркально отобразив отрезок кривой функции Ds0 N (х) на промежутке хеХ-пД! относительно оси х=-п и отрезок той же функции на промежутке хе (п,0\ относительно оси х=П , достроим одиночную реализа-

цию аппроксимирующей функции на промежутке хє(-2п,2п) (рисунок 1б):

- (х)= ООМ (- (х + 2п)),

+(х)= О^ОМ (2п- х), (14)

ММ (х) = -(х)+ ОяО N (х) + +(х). (15)

5) Периодически повторяя функцию (15) с периодом Т=4п на всем множестве значений х, получим общую аппроксимирующую функцию

(х)= Ё°51М(х - 4П). (16)

2. Методика получения значений локальных экстремумов X0М, X 1,м и точка перегиба Хм.

1) Для функции (1) изменен масштаб аргумента в N раз

2 (5)

Ds

ж

/ \

x

, V

V У

sln

(17)

Длительность каждой полуволны Ах теперь не зависит от N - Ах=2п. От N будет зависеть теперь период повторения.

2) Определены выражения первой и второй производной функции (16)

dx

,(ів)

.sin! )^{(( - 0-5) sin((1+V )|)-((+0-5) sin((1-V )|)}

' ^srf ((v)

d2Dv ()

3x2

2V2 coSx)s1r2(-2v)—Vs1r(x)s1r(w )+s1n2fX)• (со^^ 4V2s1n^

(19)

3) В программной среде MathCad числен ным решением уравнения

D (f)

dx

= 0

(20)

с помощью функции root получены корни Xm0N и XmXN на промежутках хе(0,2п) и хе (0,4п) для N = 5...100. Примененная точность решения такова, что

dD,

і Xm0lN ^

Ч V0 І

dx

< 10"

dD,

і Xm1lN \

IZJ

dx

< 2.2 • 10"

2

2

В процессе решения было замечено, что для N > 42 искомые экстремумы отличаются только в 4-м значащем разряде.

4) Также в программной среде MathCad численным решением уравнения

dx2

(21)

с помощью той же функции root найдены точки перегиба XmN на промежутке x е {Xm0N, Xmj, N) для N = 5...100. Примененная точность решения позволила добиться результата

3 2 т-\ |X1,N

д Dn\—

dx2

< іо-

при этом в ходе того расчетного эксперимента замечено, что для N > 31 искомые точки перегиба отличаются только в 4-м значащем разряде.

5) Для конкретного значения N для не-масштабированной функции (1) соответствующие точки экстремумов и точки перегиба пересчитываются при помощи выражений

Хт ы Хтх ^

X о, N ="

'о, N

N

, X1, N =

N

, XN =

XmN

N

-. (22)

- Jit

Av4i L-..

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

А' !<><)

"1 |

і L ——<

- л

и

б

Рисунок 2.

ІГ

В ходе расчетного эксперимента в среде MatchCad с помощью выражения

АN (x)=-

DsN (x) - Dsn

Am

•1оо%

(23)

о, N

просчитано относительное, приведенное к глобальному максимуму отклонение аппроксимирующей функции и равноамплитудного полинома (1) для различных N. Проанализировав графики выражения (30), сделали вывод, что с ростом N максимальное ее значение фактически не увеличивается и стремится к значению 1,3% (рисунок 2).

3. Нас в первую очередь интересуют следующие отличительные особенности дискретного спектра модели (16):

1) степень неравномерности амплитудного дискретного спектра модели (16) в рабочей области - на частотах, для которых определен спектр полинома (1) или (5);

2) отличие от нуля значений фазового спектра для тех же частот;

3) количественные характеристики спектра «паразитных» гармоник.

Так как полином (1) на разных промежутках аппроксимируется функциями разного вида, то определение спектра целесообразно выполнять, придерживаясь следующей последовательности:

- определить выражения спектральных плотностей одиночных аппроксимирующих функций каждого отрезка в выражениях (6)-(8) и (12) для хе(-2п,2п) - то есть для отрезков кривых одиночной реализации функции (16);

- определить выражения для комплексных амплитуд гармоник их дискретных спектров как выборки из соответствующих выражений спектральной плотности на частотах, кратных частоте первой гармоники [11]

2п 1

®i =— = -; 1 T 2 ’

(24)

- определить выражения для комплексных амплитуд синусоидальных составляющих общего дискретного спектра, воспользовавшись принципом суперпозиции.

4. Итак, определим спектр функции (16). 1) Определение спектра отрезков функции (16) на интервалах

x є

( (о,N + 4п',Xо,N + 4П ) и (25)

іо

х є (- 2п(і - 2])-2п(і - 2]) + Xо^ _|и и [(1 + 2])- X,2п(1 + 2]) ] є г. (26)

а) Выражение спектральной плотности на промежутке (25) при ] = 0

Х 0^

^ (й) = I (х)е-]х<3х =

- Х 0,N

= - 2]ю- Ашо,N од) . (27)

-й2

б) Выражения для комплексных амплитуд составляющих дискретного спектра на интервале (25) определим, сделав выборки из выражения спектральной плотности (26) для частот, кратных частоте первой гармоники [10]

(28)

й = пй , п = 1,2,

А’ =

Л 0,п _

2]8 о N (пй ) 2п • Ашо

пХ п

т

п

4ю„\, -й2

. (29)

Так как функция (16) нечетная, то результирующая постоянная составляющая аппроксимирующей функции Ао = о .

в) Отрезки главных лепестков на промежутке (26) при ] = 0 симметричны относительно начала координат, поэтому их спектральная плотность определится как удвоенная мнимая часть спектральной плотности одного из отрезков с обратным знаком

2п

(й)=-2 ] | DsN (х)зіп(йх )оХ =

2п-Х о, N

- 2]Ашопййо п яіп2пй- йсоз((2п - Хо,N ))}

22 йо N - й

.(3о)

г) Выражение для комплексной амплитуды синусоидальных составляющих дискретного спектра на промежутке (26):

А±,п =

2 (пй) 2п(- 1)п+1 Аш,

cos

Ґ V \

Х о, Nn

0,N

К

4йо2^ -й2

= -(- 1)п А'

о,п

(31)

д) Общее выражение для составляющих дискретного спектра на промежутках (25) и (26)

Ао,п = А'о,п+Ао,п = ( ( 1)! )А’о,п = ^

2А'оп, прип = 2к-1, о, прип = 2к, к = 1,2, к

(32)

2) Определение спектральной характеристики отрезков функции (16) на множестве

х є (-4П + 4П,-Хо,N + 4П]и[хо,N + 4П',-4П + 4П1(33)

хє (-2п(1 - 2])+ X^ ,-2п(1 - 2])+ -4п]и

и[2п(1 + 2])-N ,2п(1 + 2])-X о, N )]є г. (34) а) Различия в поведении одиночной реализации функции (16) на интервалах (33) и (34) (при ] = 0) аналогичны различию ее поведения на интервалах (25), (26) - отрезок кривой аппроксимирующей функции на промежутке х є (- 2п+Xо^ ;-2п + ^] - зеркальное отображение отрезка на промежутке х є (- N,-X] относительно оси х=-п, а отрезок на промежутке х є [2п - N ;2п - Xо,N) -зеркальное отображение отрезка на промежутке х є [Xо,N ,^п). Поэтому выражения для составляющих дискретного спектра каждого отрезка функции (16) на множестве (33) и (34) (А1п - А4п) определены по аналогичной методике пункта 1) и приведены в таблице 1.

3) Выражения п-й гармоники спектра для і-х полуволн множества

[- [+ 4П,- N + 4п] + 4п,^П^+) + 4п]и(35) и [- 2п(1 -2./)+ N ,-2п(1 - 2])+ ^и и [(1 + 2[)-,2п(1 + 2])-N] ] є г ,

хє

где і = -

12, к, N -1, При чётных N,

12, к, -1, При нечётных N,

а также спектра «хвостов» при нечетных N -то есть на интервалах

хє и[-п-п + 4п',-п + п + 4п)и

и(п-П + 4П,П + П + 4П] (36)

(рисунок 3) также получены по аналогичной методике:

а) Для множества (35)

о л т 2 (пп ^ ' (і+о.5)лп

д =(-(-іУ) 8^ 8тП)8іп А5-і-п I1 ( V П ( лг2 2) (і+о.5)ги . (37)

п- п • \4М - п ^іп ^ ж

б) Для «хвостов» (36)

А ( ( 1)п V {2 -п2}пп-2^

Аб,п = I1 -(-1) 2---1 2 ,ЛГ2\--------------.(38)

' К-п(п2 -4^ ) 4 ’

4) Согласно принципу суперпозиции выражение для п-й гармоники для всего множества (36) получится следующим:

соз

2

2

т

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Таблица 1.

№ Интервал Выражение составляющих дискретного спектра

1. X е (-X + 4П, х 0,« + 4П) и(- 2п(1 - 2 j )-2п(1 - 2 j )+ X о>№ ]и и[(1 + 2 j)-X o,n ,2п(1 + 2 j)] е Z f ^X о,N ^ cos A o,„ -(-(-1)") n 4®o2,N -®2

2. X е (- X N + 4nl,-X 0,N + 4П?]и U[X o,N + 4П, Xn + 4П')и u(-2п(1 -2 j)+X o,N;-2n(1 -2 j)+Xn ]u u[2n(1 + 2j)-Xn;2n(1 + 2j)-Xo,n) j е Z - "XN " o N A =(1 -(-1)" )2 I2 ^ ^ S1I4T-¥Г ^ - ■ 4 , 4ffl12N - n2 Y1,N ( nXN "Xo n \l -— PS -N -coS Tf"j

3. X е (-тл + 4П,-XN + 4^]u[Xn + 4П,2Nr + 4n)u u(- 2n(1 - 2 j )+ Xn , -2n(1 - 2 j )+ f-]u u[2n(1 + 2j)-N,2n(1 + 2 j)-Xn)j е Z A ( ( A- )Y1,n coS "f- coS N ®2,n Sln "if + "coS jf]

A 2,- V1 V 4 / 1 л <2 2 f n [ 2" 4®22,n - -2 J

4. x е (-X!,n + 4П,- ■f + 4nj]u u[f + 4nj, X 1,n + 4n')u u(- 2n(1 - 2 j )+ f ,-2n(1 - 2 j )+ X ,n ]u u[2n(1 + 2 j)-X^,2n(1 + 2 j)-N) j е Z A ( ( Л- ) Am3f . "(coS "lif1 + coS ^^) cos AiffL-coS

“3,n ^ V ' П [ 4ffl32N -n2 - J

5. X е (- ^N + 4П,-X 1,N + 4nj]u u[X 1, n + 4n\-f + 4П )u u(- 2n(1 - 2 j )+ X ,n , -2n(1 - 2 j )+ f]u u[2n(1 + 2 j)-N,2n(1 + 2 j)-X,n)j е Z. A ( ( 1)" ^2 J[m 3,N "X1 N ?mi 1 "X1 N . 2П" 1 co^^— coS^f " \co^^ + coS^fj

4,П ^ V ' ' n [ " 4®32n -"2 [

ев Ю 40 чо Для чётных N: X е(-2п(1 -2 j)+ N,-f + 4nj) u(( + 4П ,2n(1 + 2 j )-N.), для нечётных N: X е(-2п(1 -2 j)+ f,-n(1 -4 j)-f) u(-n(1 - 4 j )+§■, - in + 4П) u(( + 4П', n(1 + 4 j)-N)u u( + 4 j )+-N ,2n(1 + 2 j )-N) I . ,\ 8N 2 Sin f -1Sln ('+os5f)”’ A 5," = ( - (-1)" ) (4 (2 2W 2 1 ) —iJ+^ П" • ^4N - " ) i=2 Sin v 2f N-1 , , , 0 1 t 2 • n" 1 ' (i+o.5 )яй A ( ( Л" ) 8 N S1n 2f ^ S1n 2N A 5," - I1 V U /_ ( 2 2 ) ) • (i+o.5)ni П" • ^4 N -" ) i-2 Sin V 2f

7. Для нечётных N X е u[-n--N + 4nj,-n+n + 4nj)u и(п--П + 4nj' ,n+-N + 4nj' ] a -(-(-1 jr)f; ": }coiJy;2^ fI;2(J-L(1- ^) n • - 4 f j

Рисунок 3.

N1

А5,п А5,і,п

і=2

і \ 8N2 чіп — N1 0'+0-5К

:(і-(- 1)п) ^ 5іП 2 N у8іп -

' пп '(4N2 - п2) 1

=2 ЯіП

(і+0.5)гсп

Ж

І N — 1 чётных N,

(39)

(40)

где N1 = - _

IN_1 _ 1 ирм нечётных N.

5) Результаты вычислений приведем в таблицу 1.

6) Общее выражение для комплексной амплитуды п-й гармоники получаем, используя принцип суперпозиции:

А п =

А0,п + А1,п + А2,п + А3,п +

+ А4 п + А5 п для чётных 14,

А0,п + А1,п + А2,п + А3,п + (41)

+ А4,п + А5 п + А6 п для не чётных

5. По графикам полученного амплитудного (рисунок 4, а) и фазового (рисунок 4, б, где фп = а^(А п) спектров видим:

1) появление в спектре функции (16) паразитных гармоник;

л .

2) искажение равномерности основного амплитудного спектра;

3) фазовый спектр принимает два значения - 0 и 180 градусов, причем в рабочей области фазовый спектр отсутствует.

Вклад паразитных гармоник определен выражением для коэффициента гармоник

kN' N

0.5 I Л,

п=N

N

0.51 Л,

п=1

kN N 0

I Лп2

п=N

N 2

I лп

п=1

(42)

а степень неравномерности - по максимальному относительному отклонению основных гармоник от гармоник исходной функции (1)

А„

= тах

п=1,.., N

Лп — 1 1

•100.

(43)

Указанные характеристики просчитаны в программной среде MathCad для N = 6,...,100. Как видим из графиков на рисунке 5, коэффициент гармоник паразитного спектра не превышает значения 1,22%, а степень неравномерности достигает величины 5,6%.

В заключение можем сказать, что предлагаемый метод по исследованным характеристикам (67), (68) занимает промежуточное положение между методами, описанными работе [1, с. 166-175], но обладает меньшей функциональной сложностью реализации в отличие от «уточненного» метода, так как не требует выполнения дополнительных операций умножения.

=1

- —

/ ■

і і !

■і ,

Рисунок 5

Список использованной литературы:

1. Фролов С.С. Способы реализации равноамплитудных полиномов // Материалы всероссийской научно-практической конференции «Современные информационные технологии в науке, образовании и практике». - Оренбург, 2004.

2. Ax. 1762269 СССР, кл. 5G01R27/28. Устройство для контроля амплитудно-фазочастотных характеристик / К.П. Чух-риенко, С.Г. Лукаш, БЛ. Кучер. - Опубл. 1992, Бюл. №34.

3. Ax. 1675798 СССР, кл. 5G01R27/28. Цифровой анализатор частотных характеристик четырехполюсников / М.Я. Минц, B.H. Чинков, A^. Bоронкин. - Опубл. 1991, Бюл. №33.

4. Ax. 1679415 СССР, кл. 5G01R27/28. Устройство для контроля амплитудно-частотной характеристики / A.B. Aнош-кин, H.H. Шиенок - Опубл. 1991, Бюл. №35.

5. Ax. 2025899 РФ, кл. 5H04B3/46. Устройство для контроля и настройки амплитудно-фазочастотных характеристик / K.A. Семенов, B.A. Марков, B.^ Шаров и др. - Опубл. 1994, Бюл. №24.

6. Прибор для исследования амплитудно-частотных характеристик X1-41 // Техническое описание и инструкция по эксплуатации. - 1982.

7. Баскаков С.И. Радиотехнические цепи и сигналы: Учебник для вузов по специальности «Радиотехника». - М.: BbK-шая школа, 2000. - 462 с.

8. Ax. 1712898 СССР, кл. 5G01R27/28. Способ определения частотных характеристик измерительных каналов информационно-измерительных систем / A.A. Плавильщиков. - Опубл. 1989, Бюл. №6.

9. Ax. 2054684 РФ, кл. 6G01R23/16. Способ измерения амплитудно-частотных характеристик / Слюсар B^. - Опубл. 1996, Бюл. №25.

10. Ax. 1832360 СССР, кл. 5H03B19/00. Формирователь многочастотного сигнала / М.Я. Минц, B.H. Чинков, ЮА. Шмшилов, A^. Савицкий, A.B. Гуров - Опубл. 1991, Бюл. №29.

11. Гоноровский И.С. Радиотехнические цепи и сигналы. Часть 1. Сигналы. Линейные системы с постоянными и переменными параметрами.-М.: «Советское радио»., 1966. - 439 с.

12. Фролов С.С., Шевеленко B^., Бурькова E.B. Метод аппроксимации равноамплитудных полиномов // Bестник ОГУ. - №5.

13. Двайт Г.Б. Таблицы интегралов и другие математические формулы: Пер. с англ. H.B. Леви / Под ред. K.A. Семендяева.- М.: №ука, 1977. - 244 с.

i Надоели баннеры? Вы всегда можете отключить рекламу.