Matlab - Ayuda matlab

hola amigos me dejaron un proyecto en matlab para realizar en muy poco tiempo y de verdad estoy perdido no se si me puedan ayudar.

Se conoce que la intensidad, I, de cierta fuente extraterrestre varía de forma sinusoidal alrededor de cierta constante C0:, con período igual a 1000 años. La información que se tiene es la siguiente:

Fecha Intensidad(I)
1/1/1900 167
1/1/1910 152
1/1/1920 136
1/1/1930 121
1/1/1940 105
1/1/1950 89
1/1/1960 74
1/1/1970 58
1/1/1980 42
1/1/1990 27
1/1/2000 11

A partir de estos datos:

a. Si el tiempo (t) es medido en años, ¿ cuál es la frecuencia fundamental (ϖ0) de I en rad/año? .

b. Calcule los coeficientes c0, c1 y c2 tales que la función:

fhat(t)= c0 + c1 sinϖ0t + c2 cosϖ0t

Sea el ajuste por mínimos cuadrados de los datos dados.

c. Grafique la función de ajuste calculada, desde el año 1900 hasta el año 2999. En el mismo gráfico represente los datos dados en la tabla.

Solo se que se resuelve por series de Fourier.
cualquier ayuda la vale oro gracias.

Hola Rooger;

Tu consulta es interesante.

Con el fin de dar una solución, debería realizar algo parecido a cuando encuentras los coeficientes de fourier de una señal cualquier aproximandolo a una sumatoria de seno o cosenos.

Espero que sea de alguna ayuda.


Saludos
Dave Correa
Servicios de Programación
[email protected]
http://fismatlab.blogspot.com/

hola amigo me sirvió muchísimo tu ayuda. mira aquí te dejo el programa ya terminado.

---------------------------------------------------------------------------------------------
f = [167 152 136 121 105 89 74 58 42 27 11 ]
n = [1900:10:2000]


w0 = (2*pi*(1/1000));

G1 = n.^0;
G2 = sin (w0*n);
G3 = cos (w0*n);

G = [G1' G2' G3'];
c = (G'*G)\(f*G)'

fhat =(G*c)
figure(1)

plot (n,f,'gO',n,fhat,'b')


n = [1900:10:3000]
xhat=c(1)+c(2)*(sin(w0*n))+c(3)*(cos(w0*n));
figure(2)
plot (n,xhat,'b')
grid on


% es para que veas los periodos que cada 1000 años.

n = [1000:10:3000]
xhat=c(1)+c(2)*(sin(w0*n))+c(3)*(cos(w0*n));
figure(3)
plot (n,xhat,'b')
grid on

----------------------------------------------------------------------------------------------------
funcional ya esta solo me falta graficar en 1 misma grafica todas las variables.
si bien en la figura 1 solo se ve el ajuste entre los años 1900 a 2000
y en la figura 2 se ve la forma que toma la funcion entre los años 1000 a 3000

eh aqui mi ultimo problema como graficarlas en una sola me marca error porque me dice que no son del mismo largo.

Ayuda.

De nuevo por aquí le cuento que ya acabe mi proyecto hasta hice una interfaz y me quedo de pelos.

un concejo para los amigos que buscan ayuda, debemos buscar información, trabajar en ejemplos, leer libros y lo mas importante entender , no busquemos que nos pongan en bandeja de plata nuestro trabajo desarrollenlo ustedes mismos.