Arquivos para OCTAVE
Site: | Instituto Federal de Educação, Ciência e Tecnologia Sul-rio-grandense - Campus Passo Fundo |
Curso: | VIBRAÇÕES |
Livro: | 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