Search All of the Math Forum:
Views expressed in these public forums are not endorsed by
NCTM or The Math Forum.


aniket
Posts:
23
Registered:
1/15/13


How to implement a multiple sets of oscillator in matlab ?
Posted:
Feb 10, 2013 7:49 AM


I have adaptive oscillator function and giving external input to this oscillator. I want to create multiple oscillator function sets like more than 1 function and the external input will be same for all oscillator funtion, because of multiple oscillator i can estimate frequencies in the frequency spectrum of external input signal and their output will be summed and plot against time . How should i do this ?
i am doing simulation for real world recorded signal using one oscillator function and it is taking too much time. i want to simulate using multiple oscillator .
the following code is having only one oscillator like this i want to simulate using 2 or more oscillators. anyone is having idea how to do simulate with multiple oscillator ?
///////////////////////////////////////////
adaptive hopf oscillator function
///////////////////////////////////////////
function dz = myeqd(t,y,ti,xd) dz = zeros(3,1); mu=0.2; r= sqrt(y(1)^2 + y(2)^2); K=100; F=interp1(ti,xd,t); dz(1)= (mu  r^2)*y(1)  y(3)*y(2) +K*F; dz(2) = (mu  r^2)*y(2) + y(3)*y(1); dz(3) = (K*F) * (y(2)/sqrt(y(1)^2 + y(2)^2)); ///////////////////////////////////////////////////
main code
///////////////////////////////////////////////////
fin='Zuege_037_Engine_44.1k.wav'; R=50; %decimation by a factor of R to reduce unnecessary information N=512; %length of the filter used for decimation [x,Fs,nbit]=wavread(fin); xd=decimate(x,R,N,'FIR'); Fsd=Fs/R; td=(1:length(xd))/Fsd; ti=td; [T,Y]=ode45(@(t,y) myeqd(t,y,ti,xd),td,[10;10;1000]); plot (T,Y)



