Найдено 8 результатов

iliaW
27 дек 2017, 15:38
Форум: Цифровые фильтры
Тема: Нужно переделать под БИХ фильтр
Ответы: 0
Просмотры: 2224

Нужно переделать под БИХ фильтр

Привет, задача такая: нужно переделать цикл в этом коде под БИХ фильтр, сейчас он сделан под КИХ, и почему то вообще не хочет пропускать синусы через заданное разностное уравнение. Помогите пожалуйста Fs = 44100; t = (0:2047)/Fs; x = sin(2*pi*1000*t)+sin(2*pi*2000*t)+sin(2*pi*3000*t); b =[10.23 6 7 ...
iliaW
26 дек 2017, 16:27
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Re: Проверка фильтра

Нашёл на другом форуме вот такой код: clear all; close all; clc; Fs = 44100; t = (0:2047)/Fs; x = 5*sin(2*pi*5000*t)+5*sin(2*pi*10000*t)+5*sin(2*pi*15000*t); sp = abs(fft(x)); plot(sp(1:end/2)); b=[0.8 0.2]; figure; freqz(b); % b = [ 1 3 -2 4 -3 -2 1 ]; n = 1; y = zeros(1, length(x) ); % свертка for...
iliaW
25 дек 2017, 12:47
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Re: Проверка фильтра

А как сделать, что бы вот это работало: Fs = 8000; dt = 1/Fs; StopTime = 0.25; t = (0:dt:StopTime-dt)'; Fc = 4000; x = sin(2*pi*Fc*t); y(n)= x(n-4)+4x(n-3)+7x(n-2)+6x(n-1)+10,23x-0,02964y(n-4)+0,11856y(n-3)-1,19019y(n-2)+2,14326y(n-1)-1,04199y; z=conv(x,y) figure; plot(t,z); Как бы коэффициенты я др...
iliaW
21 дек 2017, 04:38
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Re: Проверка фильтра

Т.е надо заново пересчитывать коэффициенты?
iliaW
20 дек 2017, 15:55
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Re: Проверка фильтра

b =[10.23 6 7 4 1]; a =[1.04199 2.13426 -1.19019 0.11856 -0.02964]; y=filter(b, a, x) Здесь x вектор входного сигнала y - выход фильтра Чёт у меня не понятная штука получилась, какая-то, когда я в матлаб это дело забил: Fs = 8000; dt = 1/Fs; StopTime = 0.25; t = (0:dt:StopTime-dt)'; Fc = 4000; x = ...
iliaW
19 дек 2017, 04:17
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Re: Проверка фильтра

Если через корреляцию делать, вот код есть, можно ли в него записать это уравнение и сгенерировать синусоидальный сигнал? x=5; y=3; N=x+y-1; a=zeros(1,N); b=zeros(1,N); cor=zeros(1,N); a(1:x)=1; b(1:y)=1; for i=0:1:N-1; c=a.*b; tmp=sum(c)/N; j=i+1; cor(j)=tmp; b=circshift(b,[0,-1]); end; cor=circshi...
iliaW
18 дек 2017, 14:14
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Re: Проверка фильтра

Мне нужно именно код написать, что бы можно было потом параметры входного синусоидального сигнала менять а по поводу -1,04199y - не знаю...Вот если что передаточная
iliaW
18 дек 2017, 09:38
Форум: Цифровые фильтры
Тема: Проверка фильтра
Ответы: 10
Просмотры: 5456

Проверка фильтра

Добрый день, ситуация такая: Сделал расчёт коэффициентов БИХ фильтра, получился фильтр Батрворта 4го порядка, далее получил разностное уравнение, теперь нужно пропустить через это разностное уравнение синусоидальный сигнал, и проверить что он выдаст. Всё это дело нужно провернуть в матлабе Вот получ...