Forum ZOTKS

Navadne diferencialne enačbe bomo reševali z ukazom Bolj preproste (pot, hitrost, pospešek) smo že 11. uro, bolj zahtevne
(recimo, raketa v gravitacijskem polju ali Lorenzov čudni atraktor)
pridejo na vrsto 12. uro.

lpa

# Octave 4.2.2, Mon Feb 13 18:06:01 2023 CET <ales@piko> clc uu=sigsum([1,2,3],[30,50,60],[0,0,0]) fs=1000; ts=2; tt=0:1/fs:(ts-1/fs); uu=sigsum(tt,[1,2,3],[30,50,60],[0,0,0]); uu uu=sigsum(tt,[1,2,3],[30,50,60],[0,0,0]); [tt,uu] [tt;uu] [tt;uu]' [tt;uu]' uu=sigsum(tt,[1,2,3],[30,50,60],[0,0,0]); [...

Naredili bomo najbolj osnovno obdelavo signala. Kot primer bomo vzeli signal, ki je vsota sinusnih signalov treh frekvenc. Ena je visoka, relativno majhne amplitude, zato njen prispevek izgleda kot motnja. Postopek je naslednji: signal (s FFT) pretvorimo v spekter, spekter popravimo in transformiram...

% od signala do spektra in nazaj: FFT in IFFT % vsota ve"c signalov poljubnih amplitud % parametri na izbiro (glej ffto1b.m) function fft03a(ffin,vv,tmax,fs) tspause=2.0; tt=0:1/fs:(tmax-1/fs); ss=akordvol(ffin,vv,tmax,fs); plot(tt,ss,"linewidth",3);grid; pause(tspause); fftss=fft(ss...

% produkt ve"cih signalov poljubnih amplitud (akord) % parametri na izbiro function fft02b(ffin,vv,tmax,fs) tt=0:1/fs:(tmax-1/fs); ss=sigprod(ffin,vv,tmax,fs); plot(tt,ss,"linewidth",3);grid; pause(2); fftss=fft(ss); ff=0:1/tmax:(fs-1/tmax); nmax=round(tmax*fs); plot(ff(1:nmax),abs(f...

% produkt dveh signalov poljubnih amplitud (akord) % parametri na izbiro function fft02a(f1,f2,v1,v2,tmax,fs) tt=0:1/fs:(tmax-1/fs); ss=sigprod([f1,f2],[v1,v2],tmax,fs); plot(tt,ss,"linewidth",3);grid; pause(2); fftss=fft(ss); ff=0:1/tmax:(fs-1/tmax); nmax=round(tmax*fs); plot(ff(1:nmax),...

Posebni za ton, akord in produkt: function [yy]=tonvol(f,v,ts,fs) tt=0:1/fs:(ts-1/fs); yy=v.*sin(2*pi*f*tt); end function [aa]=akordvol(ff,vv,ts,fs) n=length(ff); bb=[]; for i=1:n bb=[bb;tonvol(ff(i),vv(i),ts,fs)]; end aa=mean(bb); end function [aa]=sigprod(ff,vv,ts,fs) n=length(ff); aa=tonvol(ff(1)...

% vsota ve"c signalov poljubnih amplitud % parametri na izbiro (glej ffto1b.m) function fft01c(ffin,vv,tmax,fs) tspause=2.0; tt=0:1/fs:(tmax-1/fs); ss=akordvol(ffin,vv,tmax,fs); plot(tt,ss,"linewidth",3);grid; pause(tspause); fftss=fft(ss); ff=0:1/tmax:(fs-1/tmax); nmax=round(tmax*fs...

% vsota dveh signalov poljubnih amplitud (akord) % parametri na izbiro % predlog: ft01b(10,15,1,0.8,1,fs), kjer fs od 100 od 10000 % (probaj tudi 40, kjer je hud aliasing) function fft01b(f1,f2,v1,v2,tmax,fs) tt=0:1/fs:(tmax-1/fs); ss=akordvol([f1,f2],[v1,v2],tmax,fs); plot(tt,ss,"linewidth&qu...

% vsota dveh signalov istih amplitud (akord) % vnaprej pripravljeno function fft01a() f1=10; f2=15; tmax=1.0; fs=1000; tt=0:1/fs:(tmax-1/fs); ss=akordvol([f1,f2],[1,0.5],tmax,fs); plot(tt,ss,"linewidth",3);grid; pause(1); fftss=fft(ss); ff=0:1/tmax:(fs-1/tmax); nmax=round(tmax*fs); plot(f...