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

nkce „StepFunction.m“ vytvoříme jednotkový impuls o šířce a středu s využitím operace posunu v čase a rozdílu jednotkových skoků. Opět naprogramujeme jako funkci pojmenovanou UnitImpuls.m neboť ji budeme v dalším používat. Bude:
% function UnitImpuls(t,a,b)
% t=time axis [sec]
% a= half width of impuls [sec]
% b= centre of impuls [sec]

function [ft]= UnitImpuls(t,a,b);
P=max(t)-min(t);
if abs(b)>P/2% limit of impuls center
b=sign(b)*P/2;
end;
ft=StepFunction(t,b+a)-StepFunction(t,b-a);
return

Volání této funkce naprogramujte opět jako funkci s názvem UnitImpulsPlot.m. Bude
% function UnitImpulsPlot(a,b,A)
% a= half width of impuls [sec]
% b= center of impuls [sec]
% A= amplitude of impuls [V]
function UnitImpulsPlot(a,b,A);
P=100; % time interval [sec]
N=512;% number of samples
t=linspace(-P/2,P/2,N); % discrete time axis

name=('Unit impuls');
ft=A*UnitImpuls(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í UnitImpulsPlot(10,0,1) obdržíte následující průběh. Volejte funkci s různými parametry a sledujte výsledný obrázek.
Poznámka: Tento impuls můžeme také vytvořit pomocí operace otočení časové osy (flipping) a součinu takto:
name='Unit impuls(use flip)'
a=10;
A=1;
ftp=A*StepFunction(t,a);% shifted unit step
ftn=fliplr(ftp);% flipping of ftp
ft=ftp.*ftn;% multiplication element-by-element
plot(t,ft)
axis([-P/2 P/2 -0.2*A 1.2*A])
title(name)
xlabel('t [sec]');
ylabel('f(t) [V]');
Vytvořme pomocí jednotkového skoku impuls o jednotkové amplitudě mající obě polarity. Bude:
name='Unit impuls- both polarities'
a=10;
A=1;
ft=A*StepFunction(t,a)-2*A*StepFunction(t,0)+A*StepFunction(t,-a);
plot(t,ft)
axis([-P/2 P/2 -1.2*A 1.2*A])
title(name)
xlabel('t');
ylabel('f(t)');

Vytvořme lineární signál o délce trvání . Bude:
name='Linear impuls'
a=25;
A=1;
ft=(A/a)*t;
ft=ft.*UnitImpuls(t,a,0);
plot(t,ft);
axis([-P/2 P/2 -1.2*A 1.2*A])
title(name)
xlabel('t [sec]');
ylabel('f(t) [V]');
a je ukázán na následujícím obrázku.
Dalším užitečným signálem, který si vytvoříme jako MATLABovskou funkci bude trojúhelníkový impuls, který nazveme TriangleImpuls.m a naprogramujeme ho opět jako funkci. Bude
% function TriangleImpuls(t,a,b)
% t= time axis [sec]
% a= half width of impuls [sec]
% b= centre of impuls [sec]

function [ft] = TriangleImpuls(t,a,b);
P=max(t)-min(t);
if abs(b)>P/2% limit of impuls center
b=sign(b)*P/2;
end;

ftp=(1/a)*(t+b+a);
ftp=ftp.*(StepFunction(t,a+b)-StepFunction(t,b));
ftn=(1/a)*(-t-b+a);
ftn=ftn.*(StepFunction(t,b)-StepFunction(t,-a+b));
ft=ftp+ftn;
return

Jeho volání opět naprogramujeme jako funkci
% function TriangleImpulsPlot(a,b,A)
% a= half width of impuls [sec]
% b= centre of impuls [sec]
% A= amplitude of i

Stáhnout celý tento materiál

Předchozí

1 2 3

Další

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