Studijní materiály

Zpět

Hromadně přidat materiály

Cvičení Matlab 1

BASS - Analýza signálů a soustav

Hodnocení materiálu:

Zjednodušená ukázka:

Stáhnout celý tento materiál

mpuls [V]
function TriangleImpulsPlot(a,b,A)
P=100; % time interval [sec]
N=512;% number of samples [-]
t=linspace(-P/2,P/2,N); % discrete time axis

name=('Triangle impuls');
ft=A*TriangleImpuls(t,a,b);
plot(t,ft)
axis([-P/2 P/2 -0.2*A 1.2*A])
title([name,' a=',num2str(a),' b=',num2str(b),' A=',num2str(A)]);
xlabel('t [sec]');
ylabel('f(t) [V]');
return

a po zavolání TriangleImpulsPlot(10,0,1) je výsledek na následujícím obrázku. Volejte funkci s různými parametry a sledujte výsledný obrázek.
Ještě jeden signál, a to tlumený kosinusový signál si vytvoříme jako funkci DampedCosineImpuls.m. Bude
% function DampedCosineImpuls(t,a,b,w0)
% t= time axis [sec]
% a= damping parameter [1/sec]
% b= centre of impuls [sec]
% w0= cosine frequency [rad/sec]
function [ft] = DampedCosineImpuls(t,a,b,w0);
P=max(t)-min(t);
if abs(b)>P/2% limit of impuls centre
b=sign(b)*P/2;
end;
ft=cos(w0*(t+b)).*exp(-a*abs(t+b));
return;
Volání této funkce opět naprogramujeme jako funkci DampedCosineImpulsPlot.m

% function DampedCosineImpulsPlot(a,b,w0,A);
% a= damping parameter [1/sec]
% b= centre of impuls [sec]
% w0= frequency of cosine [rad/sec]
% A= amplitude of impuls [V]

function DampedCosineImpulsPlot(a,b,w0,A);

P=100; % time interval [sec]
N=512;% number of samples [-]
t=linspace(-P/2,P/2,N); % discrete time axis

name=('Damped cosine impuls');
ft=A*DampedCosineImpuls(t,a,b,w0);
plot(t,ft)
axis([-P/2 P/2 -1.2*A 1.2*A])
title([name,' a=',num2str(a),' b=',num2str(b),' w0=',num2str(w0),' A=',num2str(A)]);
xlabel('t [sec]');
ylabel('f(t) [V]');
return;

Volání této funkce DampedCosineImpulsPlot(0.05,0,1,1) generuje následující obrázek.

Experimentujte s parametry volání a sledujte výsledek. Prověřte např., že pro je generována netlumená kosinusovka nebo pro je generován exponenciální signál, nebo pro je generován narůstající kosinový signál.

Cvičení Matlab Analýza periodických signálů
Naprogramujeme nyní výpočet koeficientů Fourierovy řady. Koeficienty jsou dány integrálem
kde je perioda a základní frekvence je . Jak časová osa tak i průběh signálu jsou V Matlabu dány hodnotami v diskrétních časových okamžicích tj. vektory kde . Tedy i funkce bude vektorem o prvcích a integrál bude nahrazen sumou a diferenciál nahradíme přírůstkem . Bude
Jelikož prvky vektorů v Matlabu jsou číslovány od do provedeme v sumě ještě substituci a výraz pro koeficient pak bude
.
Tento výraz naprogramujeme jako funkci FourierCoefficients.m
% FourierCoefficients(t,ft,M)
% t= time axis
% ft= signal
% M= half number of Fourier coefficient

function [cm] = FourierCoefficients(t,ft,M);
P=max(t)-min(t);% fundamental period P
N=length(t);% number of samples
w0=2*pi/P;% fundamental frequency
for m=-M:1:M;
cm(m+M+1)=ft*exp(-j*m*w0*t)'/(N-1);% scalar product
if abs(real(cm(m+M+1)))

Stáhnout celý tento materiál

Předchozí

1 2 3

Vloženo: 28.05.2009

Velikost: 156,18 kB

Stáhnout celý tento materiál

Komentáře

Tento materiál neobsahuje žádné komentáře.

Mohlo by tě zajímat:

Skupina předmětu BASS - Analýza signálů a soustav
Reference vyučujících předmětu BASS - Analýza signálů a soustav

Podobné materiály