integrales

y=ax^b
v=a(dp)^b
ln(v)=ln(a)+bln(dp)
y=a°+bx
clear all, clc
dp =[30 35.5 50.5 75 92 105 115 130 153.5 180 199.5]';
v = [3.83 4.17 4.97 6.06 6.71 7.17 7.51 7.98 8.67 9.39 9.89]';
y = log(v);x= log(dp)
plot(x,y,'o')
A=[ones(size(x))  x];
coef = A\y;
a= exp(coef(1))
b=coef(2)
Dp=(3.2:0.5:5.7)';
V=[ones(size(Dp))  Dp]*coef;
plot(x,y,'o', Dp,V,'-'); grid on
En ventana de comandos:a = 0.6983 b = 0.5005
La ecuación nos queda: y=0.6983x^0.5005
25 DE OCTUBRE
y=ae^bx
ln(y)=ln(a)+bx
y=a°+bx
clear all, clc
x=[1 3 5 7 9]';
y=[115 105 95 85 80]';
Y=log(y);A=[ones(size(x)) x];
coef=A\Y;
a=exp(coef(1))
b=coef(2);
X=(1:0.1:9)';
YY=[ones(size(X)) X]*coef;
plot(x,Y,'o',X,YY,'-')
y=120.2736e-0.0469x
x= [1 3 5 7 9];
y=120.2736.*exp(-0.0469*x)
caso 3:
y=ab^xLn(y)=ln(a)+ln(b)x
y=a°+mx
LA PENDIENTE ES EL ln(b)
clear all,clc
x=[1 2 3 4 5]';
y=[132 150 175 195 220]';
Y=log(y);
A=[ones(size(x)) x];
coef=A\Y
a=exp(coef(1))
b=exp(coef(2))X=(1:0.1:5)';
YY=[ones(size(X)) X]*coef;
plot(x,Y,'o',X,YY,'-')
coef =    4.7598
    0.1284
a =   116.7185
b =    1.1370
>> x=[1 2 3 4 5];
>> y=116.7185*(1.1370).^x
y = 132.7089 150.8901 171.5620195.0660 221.7900
clear all,clc
x=[1 2 3 4 5]';
y=[132 150 175 195 220]';
coef=polyfit(x,y,1)
X=(1:0.1:5)';
Y=polyval(coef,X);
plot(x,y,'o',X,Y,'-')
coef =
   22.1000  108.1000
y=   22.1000 x+108.1000
Caso 4:
Elabore un programa que grafique tanto los datos dados como la ecuación ajustada.
xy^a=b
ln(x)+aln(y)=ln(b)
ln(y)=ln(b)/a – (1/a)ln(x)
y=Ao + mx
ln(b)/a=coef(1)ln(b)=a*coef(1)
b=exp(a*coef(1))
-(1/a)=coef 2
a=coef(2)/-1

clear all, clc
x=[200 150 100 60 40 10]';
y=[1 1.5 1.8 2.4 4.1 6.5]';
X=log(x);
Y=log(y);
A=[ones(size(X)) X];
coef=A\Y;
a= -1/(coef(2))b= exp(a*coef(1))
XX=(5.2:-0.2:2.3)';
YY=[ones(size(XX)) XX]*coef;
plot(X,Y,'d',XX,YY,'-')

en comandos:
a = 1.6251
b = 260.3178
la ecuación queda xy^1.6251=260.3178
clear all, clc...