navigation
hledat
Rozdíly
Zde můžete vidět rozdíly mezi vybranou verzí a aktuální verzí dané stránky.
|
courses:a4b33opt:cviceni6 [2009/10/28 14:55] funic |
courses:a4b33opt:cviceni6 [2025/01/03 18:28] (aktuální) |
||
|---|---|---|---|
| Řádek 2: | Řádek 2: | ||
| [[https://cw.felk.cvut.cz/lib/exe/fetch.php/courses/a4b33opt/cviceni/ls.pdf?id=courses%3Aa4b33opt%3Acviceni%3Astart&cache=cache|Zadání]] | [[https://cw.felk.cvut.cz/lib/exe/fetch.php/courses/a4b33opt/cviceni/ls.pdf?id=courses%3Aa4b33opt%3Acviceni%3Astart&cache=cache|Zadání]] | ||
| - | <code> | + | <code matlab> |
| % uloha 1: | % uloha 1: | ||
| Data = load('mzdy.txt','-ascii'); | Data = load('mzdy.txt','-ascii'); | ||
| Řádek 29: | Řádek 29: | ||
| % | % | ||
| + | </code> | ||
| + | |||
| + | |||
| + | ===== druhý příklad ===== | ||
| + | |||
| + | <code matlab> | ||
| + | clear | ||
| + | close all | ||
| + | |||
| + | %nacteme | ||
| + | Data = load('teplota.txt','-ascii'); | ||
| + | days = Data(:,1); | ||
| + | temp = Data(:,2); | ||
| + | COUNT = size(days, 1); | ||
| + | PI = 3.14159265; | ||
| + | |||
| + | %matice podle PDFka | ||
| + | A = [ones(COUNT, 1) days sin(2*PI*days/365) cos(2*PI*days/365) ]; | ||
| + | %brutalne slozity vypocet | ||
| + | x = A\temp | ||
| + | |||
| + | hold on | ||
| + | plot(days, temp, 'ro'); | ||
| + | viz = 1:1:1200; %pouzijeme misto dnu, pro jemnejsi kresbu grafu | ||
| + | plot(viz, x(1,1) + x(2,1)*viz + x(3,1)*sin(2*PI*viz/365) + x(4,1)*cos(2*PI*viz/365), 'b-'); %vykresleni interpolovane funkce | ||
| + | hold off | ||
| + | |||
| + | |||
| + | T = A*x - temp; | ||
| + | T'*T %kvadrat kriterialni funkce | ||
| + | |||
| </code> | </code> | ||
