совет с фильтром

@апывпыпыап · Регистрация: 10.05.2012

Студворк — интернет-сервис помощи студентам

1)расчет фильтра прототипа(см.вариант 11)

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 Rp=0.039; 
Rs=42; 
 
Wp=[ 7.08e3 8.67e3]; 
Ws=[7.4e3 8.4e3];
 
 
 
[n,Wn]=cheb2ord(Wp, Ws, Rp, Rs,  's')
% Далее рассчитываем коэффициенты передаточной характеристики 
% с учетом порядка фильтра, частот среза и уровня ослабления сигнала: 
 
w0=  [7.3696]
 
[b,a]=cheby2(n, Rs, w0,'high',  's');
% построение графиков (рис.6): 
freqs(b,a) 
title('cheb2ord')
figure;
nyquist(tf([b],[a])), grid;

2)бих -фильтр

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
function bih
f0=[300 7.4e3];
n=6;
 [b,a]=cheby2(6,80,[.2,.5]);
fs=44000;
[bz,az]=impinvar(b,a,fs);%импульсная инвариантная характеристика
[hz,wz]=freqz(bz,az,[],fs,'whole');
h=freqs(b,a,2*pi*wz);
plot(wz,abs(h),':')
xlim([-200 1400])
xlabel('Частота (Гц)')
ylabel('Амплитуда')
title('Аналоговый фильтр-прототип')
grid
figure
subplot(3,1,1);
plot(wz,abs(hz))
xlim([-200 1400])
xlabel('Частота (Гц)')
ylabel('Амплитуда')
title('АЧХ БИХ-фильтра')
grid
subplot(3,1,2);
plot(wz,angle(hz))
xlim([-200 1400])
xlabel('Частота (Гц)')
ylabel('Фаза')
title('ФЧХ')
grid
subplot(3,1,3);
[H,nT]=impz(bz,az,[],fs);
subplot(3,1,3)
stem(nT,H),grid
xlabel('Отсчеты')
ylabel('Амплитуда')
title('ИХ');
end

3)ких фильтр

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
function kih
n=9
winda=blackman(n+1);
fs=11025;
fn=fs/2;
Wn=[300/fn 700/fn];
kih=fir1(n,Wn,'bandpass',winda);
[h, f] = freqz(kih, 1, [], fs);
subplot(3,1,1)
plot(f, abs(h))
title('АЧХ');
grid
xlim([0 1000])
subplot(3,1,2)
fi=angle(h);
fi=unwrap(fi);
plot(f,fi)
title('ФЧХ')
grid
xlim([0 1000])
[H,nT]=impz(kih,1,[],fs);
subplot(3,1,3)
stem(nT,H),grid
title('ИХ');
end

вопрос в том как их соединить в один файл(я имею в виду начальные значения) для расчета (а то вылазат ошибки)?

@R2D2 · 07.12.2013, 15:11

А Вы не хотите создать для КИХ и БИХ фильтра два скрипта? Или, если уже городить функцию, то задайте ей входные параметры.

@апывпыпыап · 07.12.2013, 18:14 **[ТС]**

у меня вопрос заключается в следующем :
насколько я понимаю сначала рассчитываться фильтр прототип, а потом эти параметры используются для КИХ и БИХ или они рассчитываться/задаются отдельно( я про порядок фильтра, частоты)?

Добавлено через 26 минут
грубо говоря какие внести изменения для КИХ и БИХ для их реализации ?

Добавлено через 1 час 46 минут
так кто-то может мне объяснит про связь коэффициентов ?

@апывпыпыап · 10.12.2013, 13:53 **[ТС]**

кто то может проверить правильность программы

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
function plottanh
 close all, clear,clc
 
 Rp=0.039; 
Rs=42; 
fs=44000;
Wp=[ 7.08e3]; 
Ws=[ 8.67e3];
 
 
[n,Wn]=cheb2ord(Wp, Ws, Rp, Rs,  's')
w0= Wn;
 
[b,a]=cheby2(n, Rs, w0,'high',  's');
freqs(b,a) 
title('cheb2ord')
figure;
 zplane (b,a),grid 
 title(' карта нулей и полюсов');
  figure
impulse(b,a),grid;
 title(' imp');
figure
 
  f=0.01:0.01:100;
  tau_group=grpdelay(a,b,f,1000)/100;
 
plot(f,tau_group,'r')
grid on
title('Групповая  задержка'),
xlabel('Частота (доли от Fs)'),
ylabel('задержка (N)'),
legend('Групповая задержка', 2)
figure
 
 
%ких
winda=blackman(n+1);
fn=fs/2;
Wq=[w0/fn]
kih=fir1(n,Wq,'high',winda);
[h, f] = freqz(kih, 1, [], fs);
subplot(3,1,1)
plot(f, abs(h))
title('АЧХ');
grid
xlim([0 1000])
subplot(3,1,2)
fi=angle(h);
fi=unwrap(fi);
plot(f,fi)
title('ФЧХ')
grid
xlim([0 1000])
[H,nT]=impz(kih,1,[],fs);
subplot(3,1,3)
stem(nT,H),grid
title('КИХ');
figure
  zplane (b,a),grid 
 title(' карта нулей и полюсов'),
 
  f=0.01:0.01:100;%?????
  tau_group=grpdelay(a,b,f,1000)/100;
figure
plot(f,tau_group,'r')
grid on
title('Групповая  задержка'),
xlabel('Частота (доли от Fs)'),
ylabel('задержка (N)'),
legend('Групповая задержка', 2)
figure
                                       
                                    
 
 
%bih
[bz,az]=impinvar(b,a,fs) %импульсная инвариантная характеристика
[hz,wz]=freqz(bz,az,[],fs,'whole');
 
h=freqs(b,a,2*pi*wz);
plot(wz,abs(h),':')
 
xlabel('Частота (Гц)')
ylabel('Амплитуда')
title('Аналоговый фильтр-прототип')
grid
figure
subplot(3,1,1);
plot(wz,abs(hz))
 
xlabel('Частота (Гц)')
ylabel('Амплитуда')
title('АЧХ БИХ-фильтра')
grid
subplot(3,1,2);
plot(wz,angle(hz))
xlim([-200 1400])
xlabel('Частота (Гц)')
ylabel('Фаза')
title('ФЧХ')
grid
subplot(3,1,3);
[H,nT]=impz(bz,az,[],fs);
subplot(3,1,3)
stem(nT,H),grid
xlabel('Отсчеты')
ylabel('Амплитуда')
title('БИХ');
figure
 zplane (b,a),grid 
 title(' карта нулей и полюсов'),
  figure
  
  f=0.01:0.01:100;%?????
  tau_group=grpdelay(a,b,f,1000)/100;
 
plot(f,tau_group,'r')
grid on
title('Групповая  задержка'),
xlabel('Частота (доли от Fs)'),
ylabel('задержка (N)'),
legend('Групповая задержка', 2)
 
 
 
 end

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Логарифм записывается как: (x-2)log(x^2+2) - означает логарифм (x^2+2) по основанию (x-2). Унарный минус обозначается как ! в-строка - входное арифметическое выражение в инфиксной(обычной). . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .

@R2D2 939 / 839 / 121 Регистрация: 23.11.2012 Сообщений: 2,488
	07.12.2013, 15:11
	А Вы не хотите создать для КИХ и БИХ фильтра два скрипта? Или, если уже городить функцию, то задайте ей входные параметры. 0

@апывпыпыап 0 / 0 / 0 Регистрация: 10.05.2012 Сообщений: 69
	07.12.2013, 18:14 [ТС]
	у меня вопрос заключается в следующем : насколько я понимаю сначала рассчитываться фильтр прототип, а потом эти параметры используются для КИХ и БИХ или они рассчитываться/задаются отдельно( я про порядок фильтра, частоты)? Добавлено через 26 минут грубо говоря какие внести изменения для КИХ и БИХ для их реализации ? Добавлено через 1 час 46 минут так кто-то может мне объяснит про связь коэффициентов ? 0