Arquivos para OCTAVE

Impresso por: Usuário visitante
Data: sexta, 18 abr 2025, 16:14

Descrição

Arquivos dos exemplos para OCTAVE.

1. Exemplo 1.15 - Representação gráfica de batimentos

# ex1_15.m

# Representação gráfica de batimentos

A = 1;

w = 20;

delta = 1;

for i = 1:1001

            t(i)=15*(i/1000);

             x(i)=2*A*cos(delta*t(i)/2)*cos((w*delta/2)*t(i));

end


plot (t,x);

xlabel('t');

ylabel('x(t)');

title(' Fenomeno do Batimento');

2. Exercício 1.80

# Exerc_1.80

clear

for i=1:201

                t(i) = (i-1)*30/200;

                x1(i) = 3 * sin(30*t(i));

                x2(i) = 3 * sin(29*t(i));

                x(i) = x1(i)+x2(i);

endfor

plot(t,x);

xlabel('t');

ylabel('x');


3. Exemplo Angulo de Fase

# Exerc_1.80

clear

for i=1:201

          t(i) = (i-1)*30/200;

           x1(i) = 3 * sin(t(i));

           x2(i) = 3 * sin(t(i)+pi/2);

endfor

plot(t,x1);

hold on

plot(t,x2);

4. Fourier

function[azero,a,b,xsin,xcos, w,theta]=fourier(N,m,time,x,t)

pi = 3.1416;

sumz = 0.0;

w = (2*pi)/0.12;

for i=1:N

  sumz=sumz+x(i);

endfor

azero = (2.0/N)*sumz;

for ii=1:m

  sums = 0.0;

  sumc = 0.0;

  for i=1:N

    theta(ii) = (2.0*pi*t(i)*ii)/time;

    xcos(i)=x(i)*cos(theta(ii));

    xsin(i)=x(i)*sin(theta(ii));

    sums=sums+xsin(i);

    sumc=sumc+xcos(i);

endfor

a(ii)=2.0*sumc/N;

b(ii)=2.0*sums/N;

endfor