TP02 prog

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

TP02/exercicio_01_05_05.m
clc
clear
p = 1000;
r = 0;
theta = 0:pi/20:2*pi;
e = 0.5;
r = p./(1- e* cos (theta));
polar(theta,r,'b')
title('\bfÓrbita do Satélite');
hold on
e = 0.25;
r = p./(1- e* cos (theta));
polar(theta,r,'--g')
hold on
e = 0;
r = p./(1- e* cos (theta));
polar(theta,r,'-r')
hold off
legend('\epsilon = 0.5','\epsilon = 0.25','\epsilon = 0');
TP02/exercicio_01_12_05.m
x=0:0.1:2;
y=exp(-x).*sin(x);
figure;
%gráfico de haste
subplot(2,2,1)
stem(x,y);
title('\bfExemplo de um Stem plot');
xlabel('\bf\itx');
ylabel('\bf\ity');
axis([0 2 0 0.5]); 
%gráfico de escada
subplot(2,2,2)
stairs(x,y);
title('\bfExemplo de um Stem plot');
xlabel('\bf\itx');
ylabel('\bf\ity');
axis([0 2 0 0.5]); 
%gráfico de barras
subplot(2,2,3)
bar(x,y);
title('\bfExemplo de um Stem plot');
xlabel('\bf\itx');
ylabel('\bf\ity');
axis([0 2 0 0.5]); 
%gráfico compasso
subplot(2,2,4)
compass(x,y);
title('\bfExemplo de um Stem plot');
xlabel('\bf\itx');
ylabel('\bf\ity');
axis([0 2 0 0.5]);
TP02/exercicio_01_28_04.m
clc
clear
conv = pi/180;
g = -9,81;
vo=20;
distancia=zeros(1,91); 
for ii=1:91
 theta=ii-1;
 vxo=vo*cos(theta*conv);
 vyo=vo*sin(theta*conv);
 max_tempo=-2*vyo/g;
 distancia(ii)=vxo*max_tempo;
end 
fprintf('Distancia versus angulo theta :\n');
for ii=1:91
 theta=ii-1;
 fprintf('%2d %8.4f\n',theta, distancia(ii));
end 
[maxdist indice]=max(distancia);
maxangulo=indice - 1;
fprintf('\nDistancia maxima = %8.4f a %2d graus. \n',maxdist, maxangulo);
for ii=5:10:85
 theta=ii;
 vxo=vo*cos(theta*conv);
 vyo=vo*sin(theta*conv);
 max_tempo=-2*vyo/g; 
 x=zeros(1,21);
 y=zeros(1,21);
 for jj=1:21
 tempo=(jj-1)*max_tempo/20;
 x(jj)=vxo*tempo;
 y(jj)=vyo*tempo+0.5*g*tempo^2;
 end 
 plot(x,y,'b');
 if ii==5
 hold on;
 end
end 
title('Trajetoria da bola vs angulo inicial theta');
xlabel('x(metros');
ylabel('y(metros');
axis([0 45 0 25]);
grid on; 
vxo=vo*cos(maxangulo*conv);
vyo=vo*sin(maxangulo*conv);
max_tempo=-2*vyo/g; 
x=zeros(1,21);
y=zeros(1,21);
for jj=1:21
 tempo=(jj-1)*max_tempo/20;
 x(jj)=vxo*tempo;
 y(jj)=vyo*tempo+0.5*g*tempo^2;
end
plot(x,y, 'r','LineWidth',3.0);
hold off 
TP02/exercicio_02_12_05.m
G=6.672*(10^(-11));%N m2 / Kg2
M=5.98*(10^24); %Kg
R=6371000; %m
h=0:500000:40000000; %m
g=-G.*M./((R+h).^2);
semilogx (h,g,'r--','LineWidth',3.0);% (variaveis, cor e tipo de linha, espessura, valor da espessura)
xlabel('\bf\ith');
ylabel('\bf\itg');
fprintf('\n ALTURA = %8.4f aceleração %2d . \n', h, g)
TP02/exercicio_02_28_04.m
n=input('Entre com a quantidade de numeros :');
cont=0;
soma=0;
while cont<n
 x=input('Entre com um valor:');
 soma=soma+x^2;
 cont=cont+1;
end
media=sqrt(soma);
fprintf('A media da raiz quadrada (rms) =%f',media); 
TP02/exercicio_03_05_05.m
n=input('Entre com a quantidade de numeros :');
cont=0;
soma=0;
soma1=0;
while cont<n
 x=input('Entre com um valor:');
 soma=soma+x;
 soma1 = soma1 + x^2;
 cont=cont+1;
end
media= (1/n)*soma;
desvio = sqrt(n*soma1 - soma^2/(n * (n-1)));
fprintf('A media aritmética é =%f\n O desvio padrão é =%f',media,desvio); 
TP02/exercicio_03_28_04.m
clear all; close all; clc;
% Início do programa
n = input('Digite a quantidade de elementos da sua amostra: \n');
vlr = zeros(n,1);
q = 0;
for i = 1:n
 vlr(n) = input('Digite o valor: \n');
end
% Aritmética
for i = 1:n
 q = vlr(n)+q;
end
mdart = q/n
% RMS
for i = 1:n
 q = (vlr(n))^2+q;
end
mdrms = q/n
% Geométrica
for i = 1:n
 q = vlr(n)+q;
end
mdgmt = q^(1/n)
% Harmônica
for i = 1:n
 q = (1/vlr(n))+q;
end
mdhrm = n/q
TP02/exercicio_04_05_05.m
%Desenvolver um programa em MATLAB que desenhe as curvas de seno
%cosseno no mesmo gráfico. Coloque titulo, rótulos, legenda e 
%tipos de cores e traçados para diferenciar as duas curvas.
%Variação do angulo -> -2pi a 2pi com passo de pi/20.
clc
clear
theta = -2*pi:pi/20:2*pi;
y = cos(theta);
plot(theta,y,'--b')
hold on
theta = -2*pi:pi/20:2*pi;
y = sin(theta);
plot(theta,y,'-.r')
hold off
title('\bfSeno e Cosseno');
legend('\bf Cos','\bf Sen');
TP02/exercici_02_05_05.m
clc
clear
x = 0.1:0.1:3;
y = x.^2 - 3*x + 2;
plot(x,y)
xlabel('x');
ylabel('y');
TP02/ex_01_14_04.m
clc
clear
close all
x = 0:pi/20:2*pi;
s = sin(x);
c = cos(2*x);
plot(s,c,'or','LineWidth',2,'MarkerFaceColor','b',...
 'MarkerSize',6);
TP02/ex_01_19_05.m
%% Recomendável
clc
clear
close all
%% Inicialização de variáveis
t = 1;
l = 1.6;
%% Plotando o gráfico
k = 0:1:5; %.^k pq o k é um vetor
P = exp(-l*t)*((l*t).^k)./factorial(k);
plot(k,P,'-ro','LineWidth',3);
%% Acabamento do gráfico
grid on
title('\bfDistribuição de Poisson');
xlabel('\bf\itEventos');
ylabel('\bf\itProbabilidade');
 
TP02/ex_01_26_05.m
clc
clear all
close all
a=0;
b=0;
c=0;
M = [1,3,5; 7,11,-13;17,19,-23];
X = [a,b,c];
B = [22;-10;-14];
X = M^(-1)*B;
disp(X)
TP02/ex_02_14_04.fig
TP02/ex_02_19_05.m
clc
clear
close all
r = input('Escreve o ângulo em radianos');
g = (180*r)/pi;
b = g - fix(g);
a = b*60;
min = fix(a);
c = min - fix(min);
d = c*60;
segundos = fix(d);
fprintf('%d - %d - %f',g,min,segundos);
TP02/ex_02_26_05.m
clc
clear
close all
x=[1.5 1.8 2.2 2.7 3 3.8 4 4.2 4.3 4.4];
y=[1.7 2.56 3.1 4.3 4.7 5 5.2 5.5 6 6.2];
n=10;
soma_x=0;
soma_y=0;
soma_xy=0;
soma_x2=0;
for i=1:n
 soma_x=soma_x+x(i);
 soma_y=soma_y+y(i);
 soma_x2=soma_x2+x(i)^2;
 soma_xy=soma_xy+x(i)*y(i);
end
my=soma_y/n;
mx=soma_x/n;
m=(n*soma_xy-soma_x*soma_y)/(n*soma_x2-(soma_x)^2);
b=my - m*mx;
disp('Regressão linear por minimos quadrados:');
fprintf('Inclinação (m) =%8.3f\n',m);
fprintf('Intersecão da reta com eixo y (b) = %8.3f\n',b);
fprintf('Numero de pontos = %8d\n',n);
plot(x,y,'ro','MarkerFaceColor','r');
hold on
xmin=min(x);
xmax=max(x);
ymin=m*xmin+b;
ymax=m*xmax+b;
plot([xmin xmax],[ymin ymax],'k-','LineWidth',2);
hold off
title('\bfAjuste de Minimos quadrados');
xlabel('\bf\itx');
ylabel('\bf\ity');
legend('Medidas originais','Reta ajustada');
grid on
TP02/ex_03_14_04.m
clc
clear
close all
t = 0:10;
w = 188.5*(1 - exp(-0.2*t));
T = 10*exp(-0.2*t);
P = w.*t;
plot(t,T,'r');
hold on
plot(t,w,'g');
hold on
plot(t,P,'b');
hold off
legend('Torque induzido no eixo de um motor','Velocidade rotacional do eixo de um motor','Potência de saída de um motor')
TP02/ex_03_19_05.m
G=6.672*(10^(-11));%N m2 / Kg2
M=5.98*(10^24); %Kg
R=6371000; %m
h=0:500000:40000000; %m
g=-G.*M./((R+h).^2);
semilogx (h,g,'r--','LineWidth',3.0);% (variaveis, cor e tipo de linha, espessura, valor da espessura)
xlabel('\bf\ith');
ylabel('\bf\itg');
fprintf('\n ALTURA = %8.4f aceleração %2d . \n', h, g)
TP02/ex_03_26_05.m
clc
clear
a = input('Digite 8 valores espaçados entre colchetes\n a=');
b = input('Digite 5 valores espaçados entre colchetes\n b=');
c = intersecao(a,b);
disp(c)
TP02/intersecao.m
function c = intersecao(a,b)
na = length(a);
nb = length(b);
k = 1;
for i = 1:na
 for j = 1:nb
 if a(i) == b(j)
 c(k)=a(i);
 k = k+1;
 end
 end
end