Zaloguj się
Blog na Matlablog
Forum polskich użytkowników
 
UŻYTKOWNICY GRUPY PROFIL Zaloguj się, by sprawdzić wiadomości FAQ
 



Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

W jaki sposób dokonać filtracji sygnału?
Forum MATLAB Strona Główna-> Przetwarzanie sygnałów
Post Wysłany: 4 Czerwca 2007, Pon 7:01 pm Temat postu: W jaki sposób dokonać filtracji sygnału? Odpowiedz z cytatem
 
AUTOR:
springer



Dołączył: 04 Cze 2007
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość
Mam problem z pewnym zadaniem i proszę Was o pomoc. Nie miałem dotad styczności z matlabem, a mój angielski jest na tyle słaby że przydatność dokumentacji matlaba nie pomaga mi w rozwiązaniu zadania.

Moje zadanie jest następujące:

1. Zaprojektuj filtr o nieskonczonej odpowiedzi impulsowej metodą optymalną
2. Dokonaj implementacji filtru na wygenerowanej przez siebie próbce sygnału.
3. Przeprowadź analizę widmową testowego sygnału przed i po filtracji

Część pierwszą zrobiłem jednak proszę o pomoc z punktem 2 i 3 naprawdę nie mam pojęcia jak to zrobić.

**************************
Wygenerowany m-pliku mojego filtru:

function Hd = filtr_sk
%FILTR_SK Returns a discrete-time filter object

%
% M-File generated by MATLAB(R) 7.0 and the Signal Processing Toolbox 6.2.
%
% Generated on: 27-May-2007 21:35:03
%

% Butterworth Bandstop filter designed using FDESIGN.BANDSTOP.

% All frequency values are in Hz.
Fs = 60000; % Sampling Frequency

Fpass1 = 15000; % First Passband Frequency
Fstop1 = 16000; % First Stopband Frequency
Fstop2 = 22000; % Second Stopband Frequency
Fpass2 = 25000; % Second Passband Frequency
Apass1 = 2; % First Passband Ripple (dB)
Astop = 60; % Stopband Attenuation (dB)
Apass2 = 2; % Second Passband Ripple (dB)
match = 'stopband'; % Band to match exactly

% Construct an FDESIGN object and call its BUTTER method.
h = fdesign.bandstop(Fpass1, Fstop1, Fstop2, Fpass2, Apass1, Astop, ...
Apass2, Fs);
Hd = butter(h, 'MatchExactly', match);

% [EOF]


 

Post Wysłany: 5 Czerwca 2007, Wto 7:04 am Temat postu: Odpowiedz z cytatem
 
AUTOR:
Vieniava
Może pisać książki


Dołączył: 13 Cze 2006
Posty: 445
Skąd: Warszawa


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email Odwiedź stronę autora
filtroawanie sygnalu
spof=filter(Hd, signal);

...gdzie signal jest sygnalem ktory chcesz filtrowac.

analiza widmowaa:
poszukaj dokładnie na tym forum Wink
 

Post Wysłany: 5 Czerwca 2007, Wto 6:21 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
springer



Dołączył: 04 Cze 2007
Posty: 2


Ogląda profil użytkownika Wyślij prywatną wiadomość
niestety nieudało mi się dużo znaleźć na forum, znalazłem jedynie coś takiego:

Nfft=128;
width=50;
w=[ones(1,width) zeros(1, Nfft-width-1)];
widmoamp=abs( fft(w,Nfft) );
plot(widmoamp);
pause;

i dopisałem do tego funkcję filtracji

%filtracja
spof=filter(Hd, widmoamp);
plot(spof);

co prawda wyszły dwa różne wykresy ale nie jestem pewien co do ich poprawności, oczywiście wszystko to dopisałem do m-pliku wygenerowanego z moim filtrem.

Nie mam pojęcia w jaki sposób wygenerować próbkę sygnału, by później dokonać na niej analizy widmowej?


 

Post Wysłany: 8 Czerwca 2007, Pią 5:24 pm Temat postu: Odpowiedz z cytatem
 
AUTOR:
Vieniava
Może pisać książki


Dołączył: 13 Cze 2006
Posty: 445
Skąd: Warszawa


Ogląda profil użytkownika Wyślij prywatną wiadomość Wyślij email Odwiedź stronę autora
pojedyncza sinusoida:

Kod:
Fs=8000;       % czestotliwosc probkowania 8kHz
t
=0:1/Fs:2.5; % 2 i pol sekundy
f
=1000;          % 1000 Hz
sygnal
=2*sin(2*pi*f*t);


bardziej skomplikowane sygnaly? ... prosze:
Kod:
doc chirp

Kod:
doc pulstran

Kod:
doc sawtooth
 

Forum MATLAB Strona Główna-> Przetwarzanie sygnałów
Wyświetl posty z ostatnich:   

Napisz nowy temat     Odpowiedz do tematu Zobacz poprzedni temat :: Zobacz następny temat

Wszystkie czasy w strefie CET (Europa)

Skocz do:  

Statystyki forum:



Od dnia 08.06.2006 forum odwiedzano 46815418
Najwięcej użytkowników 266 było obecnych 19 Lutego 2015, Czw 7:03 pm

Aktualnie online:




Najnowsze posty na forum:
Uprowadzona (Taken) {S02E16} online napisy PL 2x16 CDA  (19 Czerwca 2018, Wto 11:17 am)
Quantico {S03E08} online napisy/lektor PL 3x08 CDA/Zalukaj  (19 Czerwca 2018, Wto 11:16 am)
12 malp (12 monkeys) {S04E09} ogladaj online 4x09 napisy PL  (19 Czerwca 2018, Wto 11:14 am)
12 małp (12 monkeys) {S04E08} ogladaj online 4x08 napisy PL  (19 Czerwca 2018, Wto 11:13 am)
12 malp {S04E07} ogladaj online 4x07 napisy PL CDA/Zalukaj  (19 Czerwca 2018, Wto 11:11 am)
Paquita Salas {S02E05} ogladaj online PL 2x05 CDA/Zalukaj  (19 Czerwca 2018, Wto 11:09 am)
Paquita Salas {S02E04} ogladaj online PL 2x04 CDA/Zalukaj  (19 Czerwca 2018, Wto 11:08 am)
Paquita Salas {S02E03} ogladaj online PL 2x03 CDA/Zalukaj  (19 Czerwca 2018, Wto 11:07 am)
Paquita Salas {S02E02} ogladaj online PL 2x02 CDA/Zalukaj  (19 Czerwca 2018, Wto 11:06 am)
Paquita Salas {S02E01} ogladaj online PL 2x01 CDA/Zalukaj  (19 Czerwca 2018, Wto 11:04 am)
Twoje prawa:
Nie możesz pisać nowych tematów
Nie możesz odpowiadać w tematach
Nie możesz zmieniać swoich postów
Nie możesz usuwać swoich postów
Nie możesz głosować w ankietach
Nie możesz załączać plików na tym forum
Nie możesz ściągać plików na tym forum