= 1/k;
>> tau = 0;
>> b0 = 100;
>> u = zeros (1, k);
>> up = zeros (1, k);
>> ut = zeros (1, k);
>> utp = zeros (1, k);
>> for i = 2: k;
>> y (i +1) = y (i) + yp (i) * h;
>> yp (i +1) = (y (i +1)-y (i-1))/2 * h;
>> alfa0 = y (i +1)/944;
>> alfa1 = yp (i +1)/944;
>> A = z (i) - alfa0;
>> B = (zp (i) + alfa * A)/beta;
>> z (i +1) = exp (-alfa * h) * (A * cos (beta * h) + B * sin (beta * h)) + alfa0 + alfa1 * h;
>> zp (i +1) = exp (-alfa * h) * (cos (beta * h) * (-A * alfa + B * beta) - sin (beta * h) * (A * beta + B * alfa)) + alfa1;
>> zpp (i +1) = (-1) * exp (-alfa * h) * (cos (beta * h) * (-A * alfa * alfa + 2 * alfa * beta * B + A * beta * beta) +>>>> sin (beta * h) * (-B * alfa * alfa -2 * A * alfa * beta + B * beta * beta)) ;
>> x (i +1) = 8.56 * zpp (i +1) + 240 * zp (i +1) + 16000 * z (i +1);
>> upt (i) = up (i) + b0 * x (i) * tau;
>> ut (i) = u (i) + up (i) * tau + (b0 * x (i) * tau * tau)/2;
>> up (i +1) = upt (i) + b0 * x (i +1) * (h-tau);
>> u (i +1) = ut (i) + upt (i) * (h-tau) + (b0 * x (i +1) * (h-tau) * (h-tau))/2;
>> y (i +1) = 1 - u (i +1);
>> end;
Висновки по розділу:
У даному розділі був побудований алгоритм і програма реалізації цифрового фільтра полуаналітіческім методом з одного похідної.
2.2.1. Побудова перехідних процесів цифрового фільтра, побудованого полуаналітіческім методом з одного похідної, при різних кроках дискретизації
Побудуємо графіки перехідних процесів цифрового фільтра для кроків дискретизації 0.01, 0.005 і 0.002 відповідно.
В
Рис. 2.2.2 Графіки перехідного процесу фільтра. p align="justify"> Перехідний процес цифрового фільтра при h = 0,01 с
Стале значення одно 1,776. p align="justify"> Максимальне значення 6,8.
Перехідний процес цифрового фільтра при h = 0,005 с
Стале значення одно 1,776. p align="justify"> Максимальне значення 7,92.
Перехідний процес цифрового фільтра при h = 0,002 с
Стале значення одно 1,776. p align="justify"> Максимальне значення 8,08.
Порівняємо отримані значення з аналоговим випадком:
Крок діскретізаціі0, 010,0050,002 Аналоговий фільтрМаксімальное знач.6, 87,928,088.56 Стале знач.1, 7761,7761,7761,776