28/08/2017
A=[0 0 0 0 1 3 9; 0 0 0 0 0 1 6; 1 -1 -1 -1 0 0 0; 2 0 -1 -2 0 0 0; 0 1 2 4 -1 -2 -4; 0 0 1 4 0 -1 -4; 0 1 1 1 0 0 0 ]
b=[0 0 0 0 0 0 .5]'
z=A\b;
p1=@(x) z(1)x^2;
p1=@(x) z(1)*x^2;p2=@(x) z(2)+ z(3)*x+ z(4)*x^2;p3=@(x) z(5)+ z(6)*x+ z(7)*x^2;t=[0 1 2 3]hold onfplot(f,0,1)fplot(p1,0,1)fplot(p1,[0,1])
p1=@(x) z(1)*x.^2;
p2=@(x) z(2)+ z(3)*x+ z(4)*x.^2;
p3=@(x) z(5)+ z(6)*x+ z(7)*x.^2;opción 1
hold on
fplot(p1,[0,1])
fplot(p2,[1,2])
fplot(p1,[2,3])
hold offopción 2
x1=p1(linspace(0,1,200));
x2=p2(linspace(1,2,200));
x3=p3(linspace(2,3,200));
plot(linspace(0,1,200),x1, linspace(1,2,200),x2, linspace(2,3,200),x3)
Ajuste de Funciones periódicas
https://www.epa.gov/outdoor-air-quality-data/download-daily-data
y=[282.8815789;
281.6315789;
281.7131579;
282.5368421;
281.7894737;
283.0921053;
285.5105263;
285.4236842;
283.1;
280.3526316;
277.9236842;
278.0447368;
281.25;
275.8894737;
276.9921053;
278.7631579;
277.5157895;
276.4421053;
283.4921053;
280.25;
275.9131579;
279.7973684;
278.8394737;
279.5105263;
276.3368421;
271.0289474;
273.6052632;
272.8684211;
267.7631579;
266.4815789;
267.1368421;
270.4868421;
270.6;
266.6052632;
271.8736842;
273.4657895;
267.5210526;
268.7368421;
269.1473684;
267.6315789;
268.6289474;
270.3815789;
270.6263158;
271.45;
271.3368421;
266.1631579;
265.1789474;
265.2763158;
268.2236842;
268.8684211;
269.3078947;
269.3315789;
266.7736842;
272.2815789;
268.8894737;
267.2552632;
265.6815789;
267.3263158;
271.5763158;
266.8552632;
270.2447368;
272.3631579;
266.0868421;
261.8289474;
263.4026316;
265.5921053;
262.2315789;
261.8342105;
261.7710526;
259.3921053;
258.6;
262.4026316;
264.0763158;
263.0947368;
266.9552632;
270.3605263;
272.1842105;
266.1763158;
267.3684211;
265.9763158;
264.6868421;
262.6078947;
260.4552632;
263.7552632;
262.3868421;
258.1157895;
260.6736842;
265.95;
266.6605263;
265.2605263;
266.7578947;
267.2;
270.4473684;
272.8894737;
270.3657895;
271.4842105;
269.2447368;
274.1315789;
272.8210526;
270.3657895;
269.9868421;
267.9368421;
270.5105263;
267.6815789;
262.2710526;
264.8578947;
266.8868421;
268.0026316;
269.3263158;
266.2868421;
272.0736842;
268.2763158;
271.2552632;
274.3157895;
270.8684211;
270.0289474;
271.6315789;
275.9921053;
275.9052632;
278.0526316;
278.5842105;
278.0921053;
277.1447368;
275.4763158;
277.0473684;
281.2657895;
282.5921053;
284.7526316;
281.0710526;
284.9842105;
280.4815789;
283.8684211;
284.1184211;
287.8;
293.2736842;
294.1052632;
297.6210526;
286.4947368;
283.4421053;
288.2868421;
293.2947368;
294.1026316;
295.95;
291.9578947;
296.5763158;
298.9868421;
299.0763158;
300.9131579;
303.0315789;
303.2789474;
302.9947368;
304.3552632;
305.9868421;
299.6921053;
300.5342105;
298.0815789;
299.7394737;
299.2263158;
295.4789474;
295.1368421;
300.3605263;
300.0315789;
304.5105263;
298.8578947;
302.6868421;
308.9815789;
306.45;
303.5289474;
305.55;
314.4052632;
314.4684211;
314.8447368;
315.7921053;
318.4;
314.2842105;
316.7026316;
315.55;
321.8578947;
324.0315789;
324.3131579;
324.5973684;
325.5263158;
321.2;
318.5763158;
322.3552632;
319.6184211;
319.8605263;
327.1157895;
324.8947368;
319.0473684;
323.9631579;
327.3078947;
328.2552632;
328.2;
324.0447368;
326.9236842;
329.3631579;
324.3578947;
328.6210526;
332.2342105;
334.4894737;
329.4789474;
328.5473684;
333.0657895;
327.3973684;
331.5631579;
333.9421053;
335.0394737;
325.9526316;
329.0526316;
330.6789474;
330.9105263;
332.3763158;
337.2947368;
336.9131579;
336.1078947;
343.25;
346.4894737;
348.2447368;
348.7605263;
339.9657895;
336.2947368;
343.6052632;
346.1842105;
339.2263158;
348.8684211;
345.7736842;
348.5394737;
343.2789474;
349.1815789;
352.3921053;
347.3710526;
348.1105263;
352.6078947;
348.3921053;
346.2131579;
343.2131579;
346.6578947;
354.2947368;
348.2578947;
347.9368421;
341.7684211;
344.8552632;
347.2105263;
342.6052632;
344.2315789;
351.1289474;
350.2368421;
351.5578947;
349.7605263;
351.6973684;
351.2368421;
346.2473684;
351.2868421;
353.4789474;
352.5789474;
357.3289474;
357.2394737;
354.3421053;
356.0868421;
352.1315789;
358.9;
358.9657895;
358.0263158;
356.2526316;
352.8157895;
351.9605263;
352.8105263;
350.0210526;
355.2394737;
357.9315789;
353.3736842;
356.4447368;
354.1842105;
358.7552632;
359.6842105;
363.1789474;
357.7842105;
359.6552632;
361.5947368;
363.8973684;
363.8184211;
358.9736842;
356.4342105;
360.8368421;
362.6921053;
356.7342105;
358.0236842;
357.1447368;
356.5473684;
354.7289474;
353.2631579;
349.9710526;
352.2815789;
355.5684211;
350.4921053;
349.3421053;
343.5921053;
344.6552632;
341.1421053;
333.6947368;
339.3973684;
340.3315789;
339.4052632;
336.5789474;
339.5605263;
340.7921053;
337.3868421;
339.7184211;
333.5657895;
334.1368421;
334.1842105;
330.4710526;
326.8552632;
329.7157895;
335.6789474;
328.6815789;
329.3184211;
324.9342105;
325.9289474;
326.6263158;
328.4078947;
327.4921053;
324.4157895;
323.1657895;
322.1447368;
320.5;
320.5552632;
315.7236842;
317.6184211;
315.0210526;
313.8289474;
312.6552632;
310.9710526;
308.8026316;
304.5921053;
308.0421053;
304.5421053;
302.5236842;
301.7394737;
302.1;
304.9631579;
299.7078947;
297.4052632;
295.3078947;
297.9894737;
297.0394737;
294.6210526;
292.2052632;
294.2315789;
291.1026316;
292.1710526;
293.7605263;
293.6184211;
295.5263158;
294.9368421;
292.5315789;
292.9236842;
291.3815789;
290.3657895;
288.3947368;
287.9157895;
286.8578947;
285.4684211;
287.3394737];
t=(1:365)';
w=1/365;
M=[ones(size(t)) sin(2*pi*w*t) cos(2*pi*w*t)];
p=M\y
tt=(linspace(1,365,3000))';
yy=[ones(size(tt)) sin(2*pi*w*tt) cos(2*pi*w*tt)]*p;
plot(t,y,'o',tt,yy)