Как я могу увеличить / уменьшить (частоту / высоту) сигнала - PullRequest
Новая
       16

Как я могу увеличить / уменьшить (частоту / высоту) сигнала

0 голосов
/ 05 апреля 2011

Как я могу увеличить / уменьшить (частоту / высоту) и фазу, используя fft / ifft? Я думаю, что у меня есть базовый код, но я не уверен, что делать дальше

Мне сказали, что repmat и resample могут помочь, да, я хочу линейно сместить все компоненты. Я хочу скорректировать фазу с течением времени, чтобы она создавала стоячую волну. Один сигнал с изменением фазы в одном направлении и другой сигнал с изменением фазы в противоположном направлении

PS: это сделано в коде Octave / Matlab

Пример. У меня есть сигнал, который повторяется 1 раз в секунду, и я хочу, чтобы он повторялся 3 раза в секунду. 

%Voiceprint raise lower freq phase conjugate signal
tic
clear all, clc,clf,tic
%% Sound /beep calculation complete
filerawbeepStr='calculations_complete.wav';
filerawbeeppathStr='/home/rat/Documents/octave/raw/';
filevoiceprepathStr='/home/rat/Documents/octave/eq_research/main/transform/voice/';
filewavpathStr='/home/rat/Documents/octave/eq_research/main/transform/wav/';
[ybeep, Fsbeep, nbitsbeep] = wavread(strcat(filerawbeeppathStr,filerawbeepStr));
%addpath(”/home/rat/Documents/octave/eq_research/main/transform/”); %add path to location of functions

%1a voice print import
[vp_sig_orig, fs_rate, nbitsraw] = wavread(strcat(filevoiceprepathStr,'voice8000fs.wav')); 

%vp_sig_orig=vp_sig_orig’;
vp_sig_len=length(vp_sig_orig);

%2a create frequency domain
ya_fft = fft(vp_sig_orig);
vp_sig_phase_orig = unwrap(angle(ya_fft));

%get Magnitude
ya_fft_mag = abs(ya_fft);

%3a frequency back to time domain
ya_ifft=real(ifft(ya_fft));

%adjust frequency/phase here? How?
vp_sig_new=real(ifft(ya_fft_mag.*exp(i*vp_sig_phase_orig)));

subplot(3,1,1), plot(vp_sig_orig),title('1 original time domain')
subplot(3,1,2), plot(ya_ifft),title('2 rebuild time domain')
subplot(3,1,3), plot(vp_sig_new),title('3 adjusted time')

1 Ответ

0 голосов
/ 27 апреля 2011

Это один из способов сделать это, но если сигнал большой, вам, возможно, придется увеличить количество точек для сигнала 

clear,clc
fs = 44100;                   % Sampling frequency
t=linspace(0,1,fs);
freq=1;
ya = sin(2*pi*freq*t)'; %+ 1*sin(2*pi*250*t); 


num_per_sec=5
yb=repmat(ya,num_per_sec,1);%replicate matrix 
xxo=linspace(0,1,length(yb))'; %go from 0 to 1 sec can change speed by incr/decr 1
xxi=linspace(0,1,length(ya))'; %go from 0 to 1 sec and get total samplerate from total y value
yi_t=interp1(xxo,yb,xxi,'linear');

plot(yi_t)
...