viernes, 6 de febrero de 2009

jueves, 5 de febrero de 2009

Seguimos??


Despues de este periplo estudiantil, en breve iré añadiendo todos los trabajos matlaberos que tengo en la retaguardia...........

viernes, 28 de noviembre de 2008

integrales mediante rectángulos, cambiando el numero de ptos tomados

clear all;
close all;
clc;
%------------
NUM=[256, 512, 1024, 2048, 4096]; %diferentes numeros de ptos tomados
numPUN=length(NUM);
tic; % indica el tiempo que tarda en realizar el ejercicio
for m=1:numPUN
ndata=NUM(m);%número de puntos de xx
xx=linspace(0,10,ndata);
fun=xx.^3 + 0.25*xx +3;
integral=0;
paso=xx(2)-xx(1);
for k=1:ndata-1
integral=integral + fun(k)*paso;
end%fin del bucle con contador 'k'
resul(m)=integral;
end%fin del bucle con contador 'm'
timeE=toc;

clacula la integral mediante rectángulos

%integral
clear all;clc; close all;
num=1024;
xx=linspace(0,10,num);
fun=xx.^3+0.25*xx+3;
plot(xx,fun);
deltax=xx(2)-xx(1); %almacena el incremento
integral=0;
for k=1:(num-1);
integral=integral+fun(k)*deltax; % realiza la integral por la aproximación mediante barras
end;

Función y derivada numérica (método Euler adelante)

% AUTOR:
% FECHA:
ndata=4096;
x=linspace(-10,10,ndata);
fun=x.^3+x+2;
fun1=fun+rand(1,ndata);
hold on;
plot(x,fun,'k-');
plot(x,fun1,'r+');
%graf2=figure(2);
%graf2=plot(x,fun-fun1,'r+'); % gráfica de la derivada de la función
%afectada por el random
%DERIVADA NUMERICA
deltax= x(2)-x(1);
for k=1:(ndata-1);
deri(k)= ((fun(k+1)-fun(k))/deltax);
end;
graf3=figure(3);
graf3=plot(x(1:ndata-1),deri,'b-','LineWidth',1.5);
pos=[30 100 550 350]; %determino la posición de la gráfica
set(figure(3),'Position',pos); % constituye la posición elegida

miércoles, 19 de noviembre de 2008

segundo bucle, paso a paso con funcion pause

num=4
resul=1
for k=1:num
[k resul]
pause; % espera a que toquemos cualquier tecla
resul=resul*k
[k resul]
pause;
end

utilización de For, creamos primer bucle muy sencillo.

%autor:
%18/11/08

% realizar bucle determinar factorial de un numero (menor de 15)
num=8; %numero del cual voy a realizar el factoria
lresul=1 %variable que almacena el resultado

for k=1:num
resul=resul*k
end % end del for