viernes, 6 de febrero de 2009
jueves, 5 de febrero de 2009
Seguimos??
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;
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;
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
% 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
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
%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
Suscribirse a:
Entradas (Atom)