Pronto habrá tiempo de profesionalizar un poco más el blog...por ahora... van los test de mecánica...acaato!
buen provecho con la mecánica.
FELICITACIONES A STAVROS Y LEANDRO POR SER LOS PRIMEROS GALARDONADOS CON EL PREMIO ALUMNO Y ANTIALUMNO DEL MES ...RESPECTIVAMENTE.
5 comentarios:
Matlab reducido a 300mb
http://rapidshare.com/users/KW32Q1/1?byfilename=1
TEST 1 TEMA 1
CHOLESKY
function parte1(n)
% a)Parte i genere la matriz para n
P1=n*eye(n-1,n-1);
P2=ones(1,n-1);
P3=ones(n-1,1);
A=[P1 P3; P2 n];
clear P1 P2 P3
b=ones(n,1);
save Work
clc
A
input('Presiona Enter para Seguir')
% a)Parte ii calcule la matriz R
R=chol(A);
save Work
clc
R
input('Presiona Enter para Seguir')
% a)Parte iii resuelva mediante Cholesky
V_p=eig(R); %Valores Propios de R
Check=true; %Chequeo
clc
%Aqui se decide si se aplica el metodo o no
for i=1:(n-1)
if V_p(i)<0
disp('error')
Check=false
break
else
continue
end
end
if Check==false
clc
disp('no se puede aplicar el metodo de cholesky')
disp('matriz no es definida positiva')
clear
else
clc
disp('se puede aplicar el metodo de cholesky')
disp('ya que la matriz es definida positiva')
input('Presiona Enter para Seguir')
clear V_p
load Work
%Cholesky ==> A=R'*R => R'*(R*x)=b
% R*x=y ; R'*y=b
y=R'\b;
x=R\y;
NORMA_INF=norm(x,inf);
clear y
save Work
clc
x
NORMA_INF
end
if Check~=false
clc
disp('soluciones pedidas')
% b) Soluciones
r11=R(1,1)
rnn=R(n,n)
% c) Norma Infinito
NORMA_INF
save Work
end
TEST 1 TEMA 1
GRAFICAR LAS FUNCIONES
function parte2
%Definir x,f(x)=y,g(x)=g
x=-2:0.2:2;
D=size(x); %Tamaño de D para generar f(x) y g(x) correctos
D=D(2);
y=zeros(1,D); %generar f(x)
g=zeros(1,D); %generar g(x)
% y=f(x)
for i=1:D
if (x(i))>=(-2) & (x(i))<(0)
y(i)=(-1*((x(i)).^2));
elseif (x(i))>=(0) & (x(i))<=(2)
y(i)=(1*((x(i)).^2));
end
end
clc
y
disp('resultado de f(x)')
input('presiona Enter para continuar')
%ahora g=g(x)
g=x.^3
disp('resultado de g(x)')
input('presiona Enter para continuar y mostrar el grafico')
plot(x,y,'r',x,g,'b')
save Work2
TEST 1 TEMA 2
Matriz y Normas
function parte1(n)
A=ones(n,n);
b=ones(n,1);
A=A+diag(n*ones(n,1));
x=A\b;
N_A=norm(A,2);
C_A=cond(A,2);
N_A_inv=C_A/N_A;
N_x=norm(x,2);
save Work
clc
TEST 1 TEMA 2
GRAFICAR
function parte2
%Definir x,f(x)=y,g(x)=g
x=-2:0.2:2;
D=size(x); %Tamaño de D para generar f(x) y g(x) correctos
D=D(2);
y=zeros(1,D); %generar f(x)
g=zeros(1,D); %generar g(x)
% y=f(x)
for i=1:D
if (x(i))>=(-2) & (x(i))<(0)
y(i)=1;
elseif (x(i))>=(0) & (x(i))<=(2)
y(i)=(0.5);
end
end
clc
y
disp('resultado de f(x)')
input('presiona Enter para continuar')
%ahora g=g(x)
for i=1:D
g(i)=abs(1-abs(x(i)-1));
end
g
disp('resultado de g(x)')
input('presiona Enter para continuar y mostrar el grafico')
plot(x,y,'r',x,g,'b')
save Work2
Publicar un comentario