Спектр сигнала в Матлабе

23.05.2012, 14:30. **Показов** 22943. **Ответов** 4

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

Здравствуйте! Помогите пожалуйста студенту) Задача в следующем: нужно взять сложный полигармонический сигнал, разложить его в фурье, выкинуть половину гармоник, затем сделать обратное преобразование фурье и получить сигнал с уже обрезанными гармониками. Затем опять его разложить, выкинуть половину гармоник и заново сложить, до тех пор, пока одна гармоника только не останется...
Если взять spectrum=fft2(signal,1024), а затем обратно signal=ifft2(spectrum), то все отлично..а если например в spectrum(512)=0, то после обратного преобразования получается что-то непонятное...похоже чего-то я не понимаю((...
Помогите,пожалуйста, решить данную задачу...

@vital792 · 23.05.2012, 15:45

Сообщение от AlexxxRF

spectrum=fft2(signal,1024)

сигнал двумерный что ли?
Если все таки нет, то так:

Matlab M

clear, clc, close;
 
% сигнал
t = (0:1023)/100;
s = sin(2*pi*t*10)+...
    sin(2*pi*t*20)+...
    sin(2*pi*t*30)+...
    sin(2*pi*t*40);
 
sp = fft(s); % спектр
plot(abs(sp)) % рисуем
 
sp(400:512+512-400) = 0; % обнуляем полосу в которую попадает последняя гармоника
% (с учетом зеркального отображения во второй половине)
 
s1 = real(ifft(sp)); % обратно во временное представление
sp1 = fft(s1);
 
figure
plot(abs(sp1));

Для остальных гармоник также. Доработай сам, это только идея.

raxp · 23.05.2012, 15:49

...а с нуля частот пройтись по спектру и найти первую или нулевую гармонику?

23.05.2012, 21:26

Аааа.....все понятно..все получилось,ура!!!)))Оказывается перед ifft надо real было ставить...
Спасибо огроменное!!!!))

@oleggaaa · 21.04.2015, 22:05

спектр получается размытым ,а не дискретным

AlexxxRF
		1
	Спектр сигнала в Матлабе 23.05.2012, 14:30. Показов 22943. Ответов 4 Метки нет (Все метки) Здравствуйте! Помогите пожалуйста студенту) Задача в следующем: нужно взять сложный полигармонический сигнал, разложить его в фурье, выкинуть половину гармоник, затем сделать обратное преобразование фурье и получить сигнал с уже обрезанными гармониками. Затем опять его разложить, выкинуть половину гармоник и заново сложить, до тех пор, пока одна гармоника только не останется... Если взять spectrum=fft2(signal,1024), а затем обратно signal=ifft2(spectrum), то все отлично..а если например в spectrum(512)=0, то после обратного преобразования получается что-то непонятное...похоже чего-то я не понимаю((... Помогите,пожалуйста, решить данную задачу...

raxp 10231 / 6609 / 498 Регистрация: 28.12.2010 Сообщений: 21,156 Записей в блоге: 1
	23.05.2012, 15:49	3
	...а с нуля частот пройтись по спектру и найти первую или нулевую гармонику? 1

AlexxxRF
	23.05.2012, 21:26	4
	Аааа.....все понятно..все получилось,ура!!!)))Оказывается перед ifft надо real было ставить... Спасибо огроменное!!!!))

@oleggaaa 0 / 0 / 0 Регистрация: 28.02.2015 Сообщений: 30
	21.04.2015, 22:05	5
	спектр получается размытым ,а не дискретным 0

Опции темы